FINE-SCALE PHYLOGEOGRAPHIC CONGRUENCE DESPITE DEMOGRAPHIC INCONGRUENCE IN TWO LOW-MOBILITY SAPROXYLIC SPRINGTAILS

Tropical Logistic Regression Model on Space of Phylogenetic Trees

Classification of gene trees is an important task both in the analysis of multi-locus phylogenetic data, and assessment of the convergence of Markov Chain Monte Carlo (MCMC) analyses used in Bayesian phylogenetic tree reconstruction. The logistic regression model is one of the most popular classification models in statistical learning, thanks to its computational speed and interpretability. However, it is not appropriate to directly apply the standard logistic regression model to a set of phylogenetic trees, as the space of phylogenetic trees is non-Euclidean and thus contradicts the standard assumptions on covariates. It is well-known in tropical geometry and phylogenetics that the space of phylogenetic trees is a tropical linear space in terms of the max-plus algebra. Therefore, in this paper, we propose an analogue approach of the logistic regression model in the setting of tropical geometry. Our proposed method outperforms classical logistic regression in terms of Area under the ROC Curve in numerical examples, including with data generated by the multi-species coalescent model. Theoretical properties such as statistical consistency have been proved and generalization error rates have been derived. Finally, our classification algorithm is proposed as an MCMC convergence criterion for Mr Bayes. Unlike the convergence metric used by Mr Bayes which is only dependent on tree topologies, our method is sensitive to branch lengths and therefore provides a more robust metric for convergence. In a test case, it is illustrated that the tropical logistic regression can differentiate between two independently run MCMC chains, even when the standard metric cannot.

Shared phylogeographic patterns and environmental responses of co-distributed soybean pests: Insights from comparative phylogeographic studies of Riptortus pedestris and Riptortus linearis in the subtropics of East Asia

Comparative phylogeographic studies of closely related species sharing co-distribution areas can elucidate the role of shared historical factors and environmental changes in shaping their phylogeographic pattern. The bean bugs, Riptortus pedestris and Riptortus linearis, which both inhabit subtropical regions in East Asia, are recognized as highly destructive soybean pests. Many previous studies have investigated the biological characteristics, pheromones, chemicals and control mechanisms of these two pests, but few studies have explored their phylogeographic patterns and underlying factors. In this study, we generated a double-digest restriction site-associated DNA sequencing (ddRAD-seq) dataset to investigate phylogeographic patterns and construct ecological niche models (ENM) for both Riptortus species. Our findings revealed similar niche occupancies and population genetic structures between the two species, with each comprising two phylogeographic lineages (i.e., the mainland China and the Indochina Peninsula clades) that diverged approximately 0.1 and 0.3 million years ago, respectively. This divergence likely resulted from the combined effects of temperatures variation and geographical barriers in the mountainous regions of Southwest China. Further demographic history and ENM analyses suggested that both pests underwent rapid expansion prior to the Last Glacial Maximum (LGM). Furthermore, ENM predicts a northward shift of both pests into new soybean-producing regions due to global warming. Our study indicated that co-distribution soybean pests with overlapping ecological niches and similar life histories in subtropical regions of East Asia exhibit congruent phylogeographic and demographic patterns in response to shared historical biogeographic drivers.

Tooth replacement in the early-diverging neornithischian Jeholosaurus shangyuanensis and implications for dental evolution and herbivorous adaptation in Ornithischia

BACKGROUND

Tooth replacement patterns of early-diverging ornithischians, which are important for understanding the evolution of the highly specialized dental systems in hadrosaurid and ceratopsid dinosaurs, are poorly known. The early-diverging neornithischian Jeholosaurus, a small, bipedal herbivorous dinosaur from the Early Cretaceous Jehol Biota, is an important taxon for understanding ornithischian dental evolution, but its dental morphology was only briefly described previously and its tooth replacement is poorly known.

RESULTS

CT scanning of six specimens representing different ontogenetic stages of Jeholosaurus reveals significant new information regarding the dental system of Jeholosaurus, including one or two replacement teeth in nearly all alveoli, relatively complete tooth resorption, and an increase in the numbers of alveoli and replacement teeth during ontogeny. Reconstructions of Zahnreihen indicate that the replacement pattern of the maxillary dentition is similar to that of the dentary dentition but with a cyclical difference. The maxillary tooth replacement rate in Jeholosaurus is probably 46 days, which is faster than that of most other early-diverging ornithischians. During the ontogeny of Jeholosaurus, the premaxillary tooth replacement rate slows from 25 days to 33 days with similar daily dentine formation.

CONCLUSIONS

The tooth replacement rate exhibits a decreasing trend with ontogeny, as in Alligator. In a phylogenetic context, fast tooth replacement and multi-generation replacement teeth have evolved at least twice independently in Ornithopoda, and our analyses suggest that the early-diverging members of the major ornithischian clades exhibit different tooth replacement patterns as an adaption to herbivory.

Common origin of sterol biosynthesis points to a feeding strategy shift in Neoproterozoic animals

Steranes preserved in sedimentary rocks serve as molecular fossils, which are thought to record the expansion of eukaryote life through the Neoproterozoic Era ( ~ 1000-541 Ma). Scientists hypothesize that ancient C27 steranes originated from cholesterol, the major sterol produced by living red algae and animals. Similarly, C28 and C29 steranes are thought to be derived from the sterols of prehistoric fungi, green algae, and other microbial eukaryotes. However, recent work on annelid worms-an advanced group of eumetazoan animals-shows that they are also capable of producing C28 and C29 sterols. In this paper, we explore the evolutionary history of the 24-C sterol methyltransferase (smt) gene in animals, which is required to make C28+ sterols. We find evidence that the smt gene was vertically inherited through animals, suggesting early eumetazoans were capable of C28+ sterol synthesis. Our molecular clock of the animal smt gene demonstrates that its diversification coincides with the rise of C28 and C29 steranes in the Neoproterozoic. This study supports the hypothesis that early eumetazoans were capable of making C28+ sterols and that many animal lineages independently abandoned its biosynthesis around the end-Neoproterozoic, coinciding with the rise of abundant eukaryotic prey.

Comparative mitogenomics and phylogenetic analyses of the genus Menida (Hemiptera, Heteroptera, Pentatomidae)

In order to explore the genetic diversity and phylogenetic relationship of the genus Menida Motschulsky, 1861 and reveal the molecular evolution of the family Pentatomidae, subfamily Pentatominae, complete mitochondrial genomes of three species of Menida were sequenced, and the phylogenetic relationships of tribes within the subfamily Pentatominae were studied based on these results. The mitochondrial genomes of Menidamusiva (Jakovlev, 1876), M.lata Yang, 1934, and M.metallica Hsiao & Cheng, 1977 were 16,663 bp, 16,463 bp, and 16,418 bp, respectively, encoding 37 genes and including 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes, and a control region. The mitochondrial genome characteristics of Menida were compared and analyzed, and the phylogenetic tree of the Pentatominae was constructed based on the mitochondrial genome datasets using Bayesian inference (BI) and maximum likelihood (MI) methods. The results showed that gene arrangements, nucleotide composition, codon preference, gene overlaps, and RNA secondary structures were highly conserved within the Menida and had more similar characteristics in Pentatominae. The phylogenetic analysis shows a highly consistent topological structure based on BI and ML methods, which supported that the genus Menida belongs to the Pentatominae and is closely related to Hoplistoderini. The examined East Asian species of Menida form a monophyletic group with the internal relationships: (M.musiva + (M.lata + (M.violacea + M.metallica))). In addition, these results support the monophyly of Eysarcorini and Strachiini. Placosternum and Cappaeini are stable sister groups in the evolutionary branch of Pentatominae. The results of this study enrich the mitochondrial genome databases of Pentatominae and have significance for further elucidation of the phylogenetic relationships within the Pentatominae.

Diversity of Fusarium spp. Associated with Wheat Node and Grain in Representative Sites Across the Western Canadian Prairies

Fusarium head blight (FHB) and Fusarium crown and root rot (FCRR) are major wheat diseases. Populations of FHB and FCRR pathogens are highly dynamic, and shifts in these populations in different regions is reported. Analyzing fungal populations associated with wheat node and grain tissues collected from different regions can provide useful information and predict diseases that might affect subsequent crops and effective disease management practices. In this study, wheat node and grain samples were collected from four representative sites across the western Canadian prairies in the 2018 growing season to characterize the major Fusarium spp. and other mycobiota associated with wheat in these regions. In total, 994 fungal isolates were recovered, and based on culture and molecular diagnostic methods, three genera constituted over 90% of all fungal isolates, namely Alternaria (39.6%), Fusarium (27.8%), and Parastagonospora (23.9%). A quantitative PCR (qPCR) diagnostic toolkit was developed to quantify the most frequently isolated Fusarium spp. in infected wheat tissues: Fusarium avenaceum, F. culmorum, F. graminearum, and F. poae. This qPCR specificity was validated in silico, in vitro, and in planta and proved specific to the target species. The qPCR results showed that F. graminearum was not detected frequently from wheat node and grain samples collected from four locations in this study. F. poae was the most abundant Fusarium species in grain samples in all tested locations. However, in node samples, F. culmorum (Beaverlodge and Scott) and F. avenaceum (Lacombe and Lethbridge) were the most abundant species. Trichothecene genotyping showed that the 3ADON is the most dominant trichothecene genotype (68%), followed by type-A trichothecenes (29.5%), whereas the 15ADON trichothecene genotype was least dominant (2.5%) and the NIV genotype was not detected. Moreover, a total of 129 translation elongation factor 1-alpha (TEF1α) sequences from nine Fusarium spp. were compared at the haplotype level to evaluate genetic variability and distribution. F. avenaceum and F. poae exhibited higher diversity as reflected by higher number of haplotypes present in these two species compared with the rest.

Subaqueous foraging among carnivorous dinosaurs

Secondary aquatic adaptations evolved independently more than 30 times from terrestrial vertebrate ancestors^1,2. For decades, non-avian dinosaurs were believed to be an exception to this pattern. Only a few species have been hypothesized to be partly or predominantly aquatic^3-11. However, these hypotheses remain controversial^12,13, largely owing to the difficulty of identifying unambiguous anatomical adaptations for aquatic habits in extinct animals. Here we demonstrate that the relationship between bone density and aquatic ecologies across extant amniotes provides a reliable inference of aquatic habits in extinct species. We use this approach to evaluate the distribution of aquatic adaptations among non-avian dinosaurs. We find strong support for aquatic habits in spinosaurids, associated with a marked increase in bone density, which precedes the evolution of more conspicuous anatomical modifications, a pattern also observed in other aquatic reptiles and mammals^14-16. Spinosaurids are revealed to be aquatic specialists with surprising ecological disparity, including subaqueous foraging behaviour in Spinosaurus and Baryonyx, and non-diving habits in Suchomimus. Adaptation to aquatic environments appeared in spinosaurids during the Early Cretaceous, following their divergence from other tetanuran theropods during the Early Jurassic¹⁷.

DNA barcoding ferns in an unexplored tropical montane cloud forest area of southeast Oaxaca, Mexico

DNA barcoding can be useful for species identification and phylogenetic analysis, but its effectivity has not been verified in most neotropical cloud forest plants. We tested three plastid barcodes, rbcLa, matK, and trnH-psbA, in selected pteridophytes, a well-represented group in these forests, from a little-explored area in Oaxaca, Mexico, applying the CBOL criteria for barcoding. We used BLASTn, genetic distance, and monophyly tree-based analyses employing neighbor-joining (NJ), maximum likelihood (ML), and Bayesian inference methods. Universal primers for rbcLa and trnH-psbA were successfully amplified and bi-directionally sequenced, but matK could not be amplified for most species. rbcLa showed the highest species discrimination in BLASTn (66.67%). trnH-psbA exhibited higher significant interspecific divergence values than rbcL and rbcLa + trnH-psbA (two-sample sign test, P value < 2.2e-16). Using NJ and ML phylogenetic trees, monophyletic species were successfully resolved (100%), differing only in support values and displaying full agreement with the most recent fern classification. ML trees showed the highest mean support value (80.95%). trnH-psbA was the only barcode that could detect the Elaphoglossoideae subfamily. Species discrimination did not increase using rbcLa + trnH-psbA. rbcLa is useful for fern barcoding, trnH-psbA is most helpful for phylogenetic analyses, and matK may not work as a universal barcoding marker.

The Mitochondrial Genomes of 18 New Pleurosticti (Coleoptera: Scarabaeidae) Exhibit a Novel trnQ-NCR-trnI-trnM Gene Rearrangement and Clarify Phylogenetic Relationships of Subfamilies within Scarabaeidae

The availability of next-generation sequencing (NGS) in recent years has facilitated a revolution in the availability of mitochondrial (mt) genome sequences. The mt genome is a powerful tool for comparative studies and resolving the phylogenetic relationships among insect lineages. The mt genomes of phytophagous scarabs of the subfamilies Cetoniinae and Dynastinae were under-represented in GenBank. Previous research found that the subfamily Rutelinae was recovered as a paraphyletic group because the few representatives of the subfamily Dynastinae clustered into Rutelinae, but the subfamily position of Dynastinae was still unclear. In the present study, we sequenced 18 mt genomes from Dynastinae and Cetoniinae using next-generation sequencing (NGS) to re-assess the phylogenetic relationships within Scarabaeidae. All sequenced mt genomes contained 37 sets of genes (13 protein-coding genes, 22 tRNAs, and two ribosomal RNAs), with one long control region, but the gene order was not the same between Cetoniinae and Dynastinae species. All mt genomes of Dynastinae species showed the same gene rearrangement of trnQ-NCR-trnI-trnM, whereas all mt genomes of Cetoniinae species showed the ancestral insect gene order of trnI-trnQ-trnM. Phylogenetic analyses (IQ-tree and MrBayes) were conducted using 13 protein-coding genes based on nucleotide and amino acid datasets. In the ML and BI trees, we recovered the monophyly of Rutelinae, Cetoniinae, Dynastinae, and Sericinae, and the non-monophyly of Melolonthinae. Cetoniinae was shown to be a sister clade to (Dynastinae + Rutelinae).

Museomics and phylogenomics with protein-encoding ultraconserved elements illuminate the evolution of life history and phallic morphology of flesh flies (Diptera: Sarcophagidae)

BACKGROUND

The common name of the Flesh flies (Sarcophagidae) usually relates them with organisms feeding on decomposing organic matter, although the biology of one of the largest radiations among insects also includes predation, coprophagy, and even kleptoparasitism. The question of whether the ancestor of all sarcophagids was a predator or a decomposer, or in association to which host have sarcophagids evolved, has thus always piqued the curiosity of flesh fly specialists. Such curiosity has often been hindered by both the impossibility of having a well-supported phylogeny of Sarcophagidae and its sister group to trace live habits and the scarcity of information on the biology of the group. Using a phylogenomic dataset of protein-encoding ultraconserved elements from representatives of all three subfamilies of Sarcophagidae as ingroup and a large Calyptratae outgroup, a robust phylogenetic framework and timescale are generated to understand flesh fly systematics and the evolution of their life histories.

RESULTS

The evolutionary history for Sarcophagidae reconstructed here differs considerably from previous hypotheses. Within subfamily Sarcophaginae, a group of predatory flies, including genera Lepidodexia and Boettcheria, emerged as sister-group to the rest of Sarcophaginae. The genera Oxysarcodexia, Ravinia, and Tricharaea, long considered archaic and early-branching coprophagous and sarcosaprophagous lineages, were found nested well within the Sarcophaginae as sister-group to the sarcosaprophagous Microcerella. Predation on invertebrates is suggested as the ancestral and dominant strategy throughout the early evolution of flesh flies. Several transitions from predation to sarcosaprophagy and coprophagy occur across the sarcophagid phylogenetic tree, in contrast with almost no transitions from sarcosaprophagy or coprophagy to predatory habits. Regarding the morphological evolution of flesh flies, there might be a concerted evolution of male genitalia traits, such as the phallotrema position and the juxta, or the vesica and the folding of the phallotrema. One diversification rate shift was inferred in the evolution of sarcophagids, which is related to the origin of genus Sarcophaga.

CONCLUSIONS

This study has a significant impact on understanding sarcophagid evolution and highlights the importance of having a robust phylogenetic framework to reconstruct the ancestral character state of biological and morphological characters. I discuss the evolution of life histories of the family in relation to their hosts or substrates and outline how sarcosaprophagy, coprophagy, and kleptoparasitism behavior on various hosts may have evolved from predation on invertebrates. This study provides a phylogenetic framework for further physiological and comparative genomic work between predatory, sarcosaprophagous, coprophagous, and kleptoparasitic lineages, which could also have significant implications for the evolution of diverse life histories in other Diptera.

Inferring Evolutionary Timescales without Independent Timing Information: An Assessment of "Universal" Insect Rates to Calibrate a Collembola (Hexapoda) Molecular Clock

Previous estimates of nucleotide substitution rates are routinely applied as secondary or “universal” molecular clock calibrations for estimating evolutionary timescales in groups that lack independent timing information. A major limitation of this approach is that rates can vary considerably among taxonomic groups, but the assumption of rate constancy is rarely evaluated prior to using secondary rate calibrations. Here I evaluate whether an insect mitochondrial DNA clock is appropriate for estimating timescales in Collembola—a group of insect-like arthropods characterized by high levels of cryptic diversity. Relative rates of substitution in cytochrome oxidase subunit 1 (COI) were inferred via Bayesian analysis across a topologically constrained Hexapod phylogeny using a relaxed molecular clock model. Rates for Collembola did not differ significantly from the average rate or from the rates estimated for most other groups (25 of 30), suggesting that (1) their apparent cryptic diversity cannot be explained by accelerated rates of molecular evolution and (2) clocks calibrated using “universal” insect rates may be appropriate for estimating evolutionary timescales in this group. However, of the 31 groups investigated, 10 had rates that deviated significantly from the average (6 higher, 4 lower), underscoring the need for caution and careful consideration when applying secondary insect rate calibrations. Lastly, this study exemplifies a relatively simple approach for evaluating rate constancy within a taxonomic group to determine whether the use of secondary rates are appropriate for molecular clock calibrations.

Dispersal barriers and opportunities drive multiple levels of phylogeographic concordance in the Southern Alps of New Zealand

Phylogeographic concordance, or the sharing of phylogeographic patterns among codistributed species, suggests similar responses to topography or climatic history. While the orientation and timing of breaks between lineages are routinely compared, spatial dynamics within regions occupied by individual lineages provide a second opportunity for comparing responses to past events. In environments with complex topography and glacial history, such as New Zealand's South Island, geographically nested comparisons can identify the processes leading to phylogeographic concordance between and within regional genomic clusters. Here, we used single nucleotide polymorphisms (obtained via ddRADseq) for two codistributed forest beetle species, Agyrtodes labralis (Leiodidae) and Brachynopus scutellaris (Staphylinidae), to evaluate the role of climate change and topography in shaping phylogeographic concordance at two, nested spatial scales: do species diverge over the same geographic barriers, with similar divergence times? And within regions delimited by these breaks, do species share similar spatial dynamics of directional expansion or isolation-by-distance? We found greater congruence of phylogeographic breaks between regions divided by the strongest dispersal barriers (i.e., the Southern Alps). However, these shared breaks were not indicative of shared spatial dynamics within the regions they delimit, and the most similar spatial dynamics between species occurred within regions with the strongest gradients in historical climatic stability. Our results indicate that lack of concordance as traditionally detected by lineage turnover does not rule out the possibility of shared histories, and variation in the presence and type of concordance may provide insights into the different processes shaping phylogeographic patterns across geologically dynamic regions.

Extreme genetic diversity among springtails (Collembola) in subterranean calcretes of arid Australia

The subterranean islands hypothesis for calcretes of the Yilgarn region in Western Australia applies to many stygobitic (subterranean-aquatic) species that are "trapped" evolutionarily within isolated aquifers due to their aquatic lifestyles. In contrast, little is known about the distribution of terrestrial-subterranean invertebrates associated with the calcretes. We used subterranean Collembola from the Yilgarn calcretes to test the hypothesis that troglobitic species, those inhabiting the subterranean unsaturated (non-aquatic) zone of calcretes, are also restricted in their distribution and represent reciprocally monophyletic and endemic lineages. We used the barcoding fragment of the mtDNA cytochrome c oxidase subunit 1 (COI) gene from 183 individuals to reconstruct the phylogenetic history of the genus Pseudosinella Schäffer (Collembola, Lepidocyrtidae) from 10 calcretes in the Yilgarn. These calcretes represent less than 5% of the total possible calcretes in this region, yet we show that their diversity for subterranean Collembola comprises a minimum of 25 new species. Regionally, multiple levels of diversity exist in Pseudosinella, indicative of a complex evolutionary history for this genus in the Yilgarn. These species have probably been impacted by climatic oscillations, facilitating their dispersal across the landscape. The results represent a small proportion of the undiscovered diversity in Australia's arid zone.

Genomic analysis for virulence determinants in feline herpesvirus type-1 isolates

Feline herpesvirus type 1 (FHV-1) is a widespread cause of respiratory and ocular disease in domestic cats. A spectrum of disease severity is observed in host animals, but there has been limited prior investigation into viral genome factors which could be responsible. Stocks of FHV-1 were established from oropharyngeal swabs obtained from twenty-five cats with signs of infection housed in eight animal shelters around the USA. A standardized numerical host clinical disease severity scoring scheme was used for each cat from which an isolate was obtained. Illumina MiSeq was used to sequence the genome of each isolate. Genomic homogeneity among isolates was relatively high. A general linear model for fixed effects determined that only two synonymous single nucleotide polymorphisms across two genes (UL37/39) in the same isolate (from one host animal with a low disease severity score) were significantly associated (p ≤ 0.05) with assigned host respiratory and total disease severity score. No variants in any isolate were found to be significantly associated with assigned host ocular disease severity score. A concurrent analysis of missense mutations among the viral isolates identified three genes as being primarily involved in the observed genomic variation, but none were significantly associated with host disease severity scores. An ancestral state likelihood reconstruction was performed and determined that there was no evidence of a connection between host disease severity score and viral evolutionary state. We conclude from our results that the spectrum of host disease severity observed with FHV-1 is unlikely to be primarily related to viral genomic variations, and is instead due to host response and/or other factors.

Variation in the Internal Transcribed Spacer Region of Phakopsora pachyrhizi and Implications for Molecular Diagnostic Assays

Phakopsora pachyrhizi, the causal agent of soybean rust (SBR), is a global threat to soybean production. Since the discovery of SBR in the continental United States, quantitative polymerase chain reaction assays based on the internal transcribed spacer (ITS) ribosomal DNA locus were established for its rapid detection. However, insufficient data were initially available to test assays against factors that could give rise to misidentification. This study aimed to reevaluate current assays for (i) the potential for false-positive detection caused by nontarget Phakopsora species and (ii) the potential for false-negative detection caused by intraspecific variation within the ITS locus of P. pachyrhizi. A large amount of intraspecific and intragenomic variation in ITS was detected, including the presence of polymorphic ITS copies within single leaf samples and within single rust sori. The diagnostic assays were not affected by polymorphisms in the ITS region; however, current assays are at risk of false positives when screened against other species of Phakopsora. This study raises caveats to the use of multicopy genes (e.g., ITS) in single-gene detection assays and discusses the pitfalls of inferences concerning the aerobiological pathways of disease spread made in the absence of an evaluation of intragenomic ITS heterogeneity.

First report of anatoxin-a producing cyanobacteria in Australia illustrates need to regularly up-date monitoring strategies in a shifting global distribution

Routine monitoring of toxic cyanobacteria depends on up-to-date epidemiological information about their distribution. In Australia, anatoxin producing cyanobacteria are not regularly tested for and thought to be rare if not absent from the continent. Our study investigated the presence of anatoxin-a (ATX-a) producing cyanobacteria in surface water samples (n = 226 from 67 sampling locations) collected from 2010 to 2017 across the state of Victoria, Australia. We (1) detected the presence and distribution of anaC (anatoxin-a synthetase C) gene sequences previously associated with various cyanobacteria, including Cuspidothrix issatschenkoi, Aphanizomenon sp., D. circinale, Anabaena sp., and Oscillatoria sp., from 31 sampling locations, and (2) determined the concentration of ATX-a in samples tested using ELISA, in two instances detected at >4 µg · L-1. These data present the first confirmation of ATX-a producers in Australia. Our study indicates that ATX-a should be included in regular testing of cyanobacterial blooms in Australia and highlights the importance of regular investigation of the distributions of toxic cyanobacteria worldwide, particularly amid the known expanding distribution of many cyanobacterial taxa in a period of increased eutrophication and rising surface water temperatures.

At the confluence of vicariance and dispersal: Phylogeography of cavernicolous springtails (Collembola: Arrhopalitidae, Tomoceridae) codistributed across a geologically complex karst landscape in Illinois and Missouri

The processes of vicariance and dispersal are central to our understanding of diversification, yet determining the factors that influence these processes remains a significant challenge in evolutionary biology. Caves offer ideal systems for examining the mechanisms underlying isolation, divergence, and speciation. Intrinsic ecological differences among cavernicolous organisms, such as the degree of cave dependence, are thought to be major factors influencing patterns of genetic isolation in caves. Using a comparative phylogeographic approach, we employed mitochondrial and nuclear markers to assess the evolutionary history of two ecologically distinct groups of terrestrial cave-dwelling springtails (Collembola) in the genera Pygmarrhopalites (Arrhopalitidae) and Pogonognathellus (Tomoceridae) that are codistributed in caves throughout the Salem Plateau-a once continuous karst region, now bisected by the Mississippi River Valley in Illinois and Missouri. Contrasting phylogeographic patterns recovered for troglobiotic Pygmarrhopalites sp. and eutroglophilic Pogonognathellus sp. suggests that obligate associations with cave habitats can restrict dispersal across major geographic barriers such as rivers and valleys, but may also facilitate subterranean dispersal between neighboring cave systems. Pygmarrhopalites sp. populations spanning the Mississippi River Valley were estimated to have diverged 2.9-4.8 Ma, which we attribute to vicariance resulting from climatic and geological processes involved in Mississippi River Valley formation beginning during the late Pliocene/early Pleistocene. Lastly, we conclude that the detection of many deeply divergent, morphologically cryptic, and microendemic lineages highlights our poor understanding of microarthropod diversity in caves and exposes potential conservation concerns.

Genotyping by sequencing reveals contrasting patterns of population structure, ecologically mediated divergence, and long-distance dispersal in North American palms

Comparative studies can provide powerful insights into processes that affect population divergence and thereby help to elucidate the mechanisms by which contemporary populations may respond to environmental change. Furthermore, approaches such as genotyping by sequencing (GBS) provide unprecedented power for resolving genetic differences among species and populations. We therefore used GBS to provide a genomewide perspective on the comparative population structure of two palm genera, Washingtonia and Brahea, on the Baja California peninsula, a region of high landscape and ecological complexity. First, we used phylogenetic analysis to address taxonomic uncertainties among five currently recognized species. We resolved three main clades, the first corresponding to W. robusta and W. filifera, the second to B. brandegeei and B. armata, and the third to B. edulis from Guadalupe Island. Focusing on the first two clades, we then delved deeper by investigating the underlying population structure. Striking differences were found, with GBS uncovering four distinct Washingtonia populations and identifying a suite of loci associated with temperature, consistent with ecologically mediated divergence. By contrast, individual mountain ranges could be resolved in Brahea and few loci were associated with environmental variables, implying a more prominent role of neutral divergence. Finally, evidence was found for long-distance dispersal events in Washingtonia but not Brahea, in line with knowledge of the dispersal mechanisms of these palms including the possibility of human-mediated dispersal. Overall, our study demonstrates the power of GBS together with a comparative approach to elucidate markedly different patterns of genomewide divergence mediated by multiple effectors.

Molecular characterization of Hepatozoon canis from farm dogs in Pakistan

Hepatozoon canis is a tick-borne pathogen of canids, which is distributed worldwide. However, very little is known about this protozoan parasite in Pakistan. This study provides the first molecular evidence of H. canis from farm dogs from three agro-ecological zones of Punjab, Pakistan. A conventional PCR targeting the 18S rRNA gene was used to characterize H. canis from farm dogs from three districts, namely Kasur, Rawalpindi, and Muzaffargarh, in Punjab. Of 341 blood samples tested, 155 (45.5%) were positive for H. canis, 73 (61.3%) from Kasur, 46 (42.5%) from Rawalpindi, and 36 (31.5%) from Muzaffargarh. Phylogenetic analyses revealed that 18S rRNA sequences of H. canis from this study clustered in three clades with those of H. canis from previously published studies to the exclusion of all other Hepatozoon spp. included in the analysis. This study provides the first insight into H. canis from farm dogs in Pakistan. Furthermore, it lays a foundation for future studies of the parasite to assess the impact of canine hepatozoonosis in dogs from various agro-ecological zones in Pakistan where pet ownership of dogs is increasing.

Union of phylogeography and landscape genetics

Phylogeography and landscape genetics have arisen within the past 30 y. Phylogeography is said to be the bridge between population genetics and systematics, and landscape genetics the bridge between landscape ecology and population genetics. Both fields can be considered as simply the amalgamation of classic biogeography with genetics and genomics; however, they differ in the temporal, spatial, and organismal scales addressed and the methodology used. I begin by briefly summarizing the history and purview of each field and suggest that, even though landscape genetics is a younger field (coined in 2003) than phylogeography (coined in 1987), early studies by Dobzhansky on the "microgeographic races" of Linanthus parryae in the Mojave Desert of California and Drosophila pseudoobscura across the western United States presaged the fields by over 40 y. Recent advances in theory, models, and methods have allowed researchers to better synthesize ecological and evolutionary processes in their quest to answer some of the most basic questions in biology. I highlight a few of these novel studies and emphasize three major areas ripe for investigation using spatially explicit genomic-scale data: the biogeography of speciation, lineage divergence and species delimitation, and understanding adaptation through time and space. Examples of areas in need of study are highlighted, and I end by advocating a union of phylogeography and landscape genetics under the more general field: biogeography.

Longwing (Heliconius) butterflies combine a restricted set of pigmentary and structural coloration mechanisms

BACKGROUND

Longwing butterflies, Heliconius sp., also called heliconians, are striking examples of diversity and mimicry in butterflies. Heliconians feature strongly colored patterns on their wings, arising from wing scales colored by pigments and/or nanostructures, which serve as an aposematic signal.

RESULTS

Here, we investigate the coloration mechanisms among several species of Heliconius by applying scanning electron microscopy, (micro)spectrophotometry, and imaging scatterometry. We identify seven kinds of colored scales within Heliconius whose coloration is derived from pigments, nanostructures or both. In yellow-, orange- and red-colored wing patches, both cover and ground scales contain wavelength-selective absorbing pigments, 3-OH-kynurenine, xanthommatin and/or dihydroxanthommatin. In blue wing patches, the cover scales are blue either due to interference of light in the thin-film lower lamina (e.g., H. doris) or in the multilayered lamellae in the scale ridges (so-called ridge reflectors, e.g., H. sara and H. erato); the underlying ground scales are black. In the white wing patches, both cover and ground scales are blue due to their thin-film lower lamina, but because they are stacked upon each other and at the wing substrate, a faint bluish to white color results. Lastly, green wing patches (H. doris) have cover scales with blue-reflecting thin films and short-wavelength absorbing 3-OH-kynurenine, together causing a green color.

CONCLUSIONS

The pigmentary and structural traits are discussed in relation to their phylogenetic distribution and the evolution of vision in this highly interesting clade of butterflies.

Genetic insights into family group co-occurrence in Cryptocercus punctulatus, a sub-social woodroach from the southern Appalachian Mountains

The wood-feeding cockroach Cryptocercus punctulatus Scudder (Blattodea: Cryptocercidae) is an important member of the dead wood (saproxylic) community in montane forests of the southeastern United States. However, its population biology remains poorly understood. Here, aspects of family group co-occurrence were characterized to provide basic information that can be extended by studies on the evolution and maintenance of sub-sociality. Broad sampling across the species’ range was coupled with molecular data (mitochondrial DNA (mtDNA) sequences). The primary questions were: (1) what proportion of rotting logs contain two or more different mtDNA haplotypes and how often can this be attributed to multiple families inhabiting the same log, (2) are multi-family logs spatially clustered, and (3) what levels of genetic differentiation among haplotypes exist within a log, and how genetically similar are matrilines of co-occurring family groups? Multi-family logs were identified on the premise that three different mtDNA haplotypes, or two different haplotypes among adult females, is inconsistent with a single family group founded by one male–female pair. Results showed that of the 88 rotting logs from which multiple adult C. punctulatus were sampled, 41 logs (47%) contained two or more mtDNA haplotypes, and at least 19 of these logs (22% overall) were inferred to be inhabited by multiple families. There was no strong evidence for spatial clustering of the latter class of logs. The frequency distribution of nucleotide differences between co-occurring haplotypes was strongly right-skewed, such that most haplotypes were only one or two mutations apart, but more substantial divergences (up to 18 mutations, or 1.6% uncorrected sequence divergence) do occasionally occur within logs. This work represents the first explicit investigation of family group co-occurrence in C. punctulatus, providing a valuable baseline for follow-up studies.

MtDNA metagenomics reveals large-scale invasion of belowground arthropod communities by introduced species

Using a series of standardized sampling plots within forest ecosystems in remote oceanic islands, we reveal fundamental differences between the structuring of aboveground and belowground arthropod biodiversity that are likely due to large-scale species introductions by humans. Species of beetle and spider were sampled almost exclusively from single islands, while soil-dwelling Collembola exhibited more than tenfold higher species sharing among islands. Comparison of Collembola mitochondrial metagenomic data to a database of more than 80 000 Collembola barcode sequences revealed almost 30% of sampled island species are genetically identical, or near identical, to individuals sampled from often very distant geographic regions of the world. Patterns of mtDNA relatedness among Collembola implicate human-mediated species introductions, with minimum estimates for the proportion of introduced species on the sampled islands ranging from 45% to 88%. Our results call for more attention to soil mesofauna to understand the global extent and ecological consequences of species introductions.

Comparative phylogeography in the Atlantic forest and Brazilian savannas: pleistocene fluctuations and dispersal shape spatial patterns in two bumblebees

Background

Bombus morio and B. pauloensis are sympatric widespread bumblebee species that occupy two major Brazilian biomes, the Atlantic forest and the savannas of the Cerrado. Differences in dispersion capacity, which is greater in B. morio, likely influence their phylogeographic patterns. This study asks which processes best explain the patterns of genetic variation observed in B. morio and B. pauloensis, shedding light on the phenomena that shaped the range of local populations and the spatial distribution of intra-specific lineages.

Results

Results suggest that Pleistocene climatic oscillations directly influenced the population structure of both species. Correlative species distribution models predict that the warmer conditions of the Last Interglacial contributed to population contraction, while demographic expansion happened during the Last Glacial Maximum. These results are consistent with physiological data suggesting that bumblebees are well adapted to colder conditions. Intra-specific mitochondrial genealogies are not congruent between the two species, which may be explained by their documented differences in dispersal ability.

Conclusions

While populations of the high-dispersal B. morio are morphologically and genetically homogeneous across the species range, B. pauloensis encompasses multiple (three) mitochondrial lineages, and show clear genetic, geographic, and morphological differences. Because the lineages of B. pauloensis are currently exposed to distinct climatic conditions (and elevations), parapatric diversification may occur within this taxon. The eastern portion of the state of São Paulo, the most urbanized area in Brazil, represents the center of genetic diversity for B. pauloensis.

Electronic supplementary material

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

Pangenome and immuno-proteomics analysis of Acinetobacter baumannii strains revealed the core peptide vaccine targets

BACKGROUND

Acinetobacter baumannii has emerged as a significant nosocomial pathogen during the last few years, exhibiting resistance to almost all major classes of antibiotics. Alternative treatment options such as vaccines tend to be most promising and cost effective approaches against this resistant pathogen. In the current study, we have explored the pan-genome of A. baumannii followed by immune-proteomics and reverse vaccinology approaches to identify potential core vaccine targets.

RESULTS

The pan-genome of all available A. baumannii strains (30 complete genomes) is estimated to contain 7,606 gene families and the core genome consists of 2,445 gene families (~32 % of the pan-genome). Phylogenetic tree, comparative genomic and proteomic analysis revealed both intra- and inter genomic similarities and evolutionary relationships. Among the conserved core genome, thirteen proteins, including P pilus assembly protein, pili assembly chaperone, AdeK, PonA, OmpA, general secretion pathway protein D, FhuE receptor, Type VI secretion system OmpA/MotB, TonB dependent siderophore receptor, general secretion pathway protein D, outer membrane protein, peptidoglycan associated lipoprotein and peptidyl-prolyl cis-trans isomerase are identified as highly antigenic. Epitope mapping of the target proteins revealed the presence of antigenic surface exposed 9-mer T-cell epitopes. Protein-protein interaction and functional annotation have shown their involvement in significant biological and molecular processes. The pipeline is validated by predicting already known immunogenic targets against Gram negative pathogen Helicobacter pylori as a positive control.

CONCLUSION

The study, based upon combinatorial approach of pan-genomics, core genomics, proteomics and reverse vaccinology led us to find out potential vaccine candidates against A. baumannii. The comprehensive analysis of all the completely sequenced genomes revealed thirteen putative antigens which could elicit substantial immune response. The integration of computational vaccinology strategies would facilitate in tackling the rapid dissemination of resistant A.baumannii strains. The scarcity of effective antibiotics and the global expansion of sequencing data making this approach desirable in the development of effective vaccines against A. baumannii and other bacterial pathogens.

Toward a paradigm shift in comparative phylogeography driven by trait-based hypotheses

For three decades, comparative phylogeography has conceptually and methodologically relied on the concordance criterion for providing insights into the historical/biogeographic processes driving population genetic structure and divergence. Here we discuss how this emphasis, and the corresponding lack of methods for extracting information about biotic/intrinsic contributions to patterns of genetic variation, may bias our general understanding of the factors driving genetic structure. Specifically, this emphasis has promoted a tendency to attribute discordant phylogeographic patterns to the idiosyncracies of history, as well as an adherence to generic null expectations of concordance with reduced predictive power. We advocate that it is time for a paradigm shift in comparative phylogeography, especially given the limited utility of the concordance criterion as genomic data provide ever-increasing levels of resolution. Instead of adhering to the concordance-discordance dichotomy, comparative phylogeography needs to emphasize the contribution of taxon-specific traits that will determine whether concordance is a meaningful criterion for evaluating hypotheses or may predict discordant phylogeographic structure. Through reference to some case studies we illustrate how refined hypotheses based on taxon-specific traits can provide improved predictive frameworks to forecast species responses to climatic change or biogeographic barriers while gaining unique insights about the taxa themselves and their interactions with their environment. We outline a potential avenue toward a synthetic comparative phylogeographic paradigm that includes addressing some important conceptual and methodological challenges related to study design and application of model-based approaches for evaluating support of trait-based hypotheses under the proposed paradigm.

Community trees: Identifying codiversification in the Páramo dipteran community

Groups of codistributed species that responded in a concerted manner to environmental events are expected to share patterns of evolutionary diversification. However, the identification of such groups has largely been based on qualitative, post hoc analyses. We develop here two methods (posterior predictive simulation [PPS], Kuhner-Felsenstein [K-F] analysis of variance [ANOVA]) for the analysis of codistributed species that, given a group of species with a shared pattern of diversification, allow empiricists to identify those taxa that do not codiversify (i.e., "outlier" species). The identification of outlier species makes it possible to jointly estimate the evolutionary history of co-diversifying taxa. To evaluate the approaches presented here, we collected data from Páramo dipterans, identified outlier species, and estimated a "community tree" from species that are identified as having codiversified. Our results demonstrate that dipteran communities from different Páramo habitats in the same mountain range are more closely related than communities in other ranges. We also conduct simulation testing to evaluate this approach. Results suggest that our approach provides a useful addition to comparative phylogeographic methods, while identifying aspects of the analysis that require careful interpretation. In particular, both the PPS and K-F ANOVA perform acceptably when there are one or two outlier species, but less so as the number of outliers increases. This is likely a function of the corresponding degradation of the signal of community divergence; without a strong signal from a codiversifying community, there is no dominant pattern from which to detect an outlier species. For this reason, both the magnitude of K-F distance distribution and outside knowledge about the phylogeographic history of each putative member of the community should be considered when interpreting the results.

Phylogeographic concordance factors quantify phylogeographic congruence among co-distributed species in the Sarracenia alata pitcher plant system

Comparative phylogeographic investigations have identified congruent phylogeographic breaks in co-distributed species in nearly every region of the world. The qualitative assessments of phylogeographic patterns traditionally used to identify such breaks, however, are limited because they rely on identifying monophyletic groups across species and do not account for coalescent stochasticity. Only long-standing phylogeographic breaks are likely to be obvious; many species could have had a concerted response to more recent landscape events, yet possess subtle signs of phylogeographic congruence because ancestral polymorphism has not completely sorted. Here, we introduce Phylogeographic Concordance Factors (PCFs), a novel method for quantifying phylogeographic congruence across species. We apply this method to the Sarracenia alata pitcher plant system, a carnivorous plant with a diverse array of commensal organisms. We explore whether a group of ecologically associated arthropods have co-diversified with the host pitcher plant, and identify if there is a positive correlation between ecological interaction and PCFs. Results demonstrate that multiple arthropods share congruent phylogeographic breaks with S. alata, and provide evidence that the level of ecological association can be used to predict the degree of similarity in the phylogeographic pattern. This study outlines an approach for quantifying phylogeographic congruence, a central concept in biogeographic research.

Identification of Eastern United States Reticulitermes Termite Species via PCR-RFLP, Assessed Using Training and Test Data

Reticulitermes termites play key roles in dead wood decomposition and nutrient cycling in forests. They also damage man-made structures, resulting in considerable economic loss. In the eastern United States, five species (R. flavipes, R. virginicus, R. nelsonae, R. hageni and R. malletei) have overlapping ranges and are difficult to distinguish morphologically. Here we present a molecular tool for species identification. It is based on polymerase chain reaction (PCR) amplification of a section of the mitochondrial cytochrome oxidase subunit II gene, followed by a three-enzyme restriction fragment length polymorphism (RFLP) assay, with banding patterns resolved via agarose gel electrophoresis. The assay was designed using a large set of training data obtained from a public DNA sequence database, then evaluated using an independent test panel of Reticulitermes from the Southern Appalachian Mountains, for which species assignments were determined via phylogenetic comparison to reference sequences. After refining the interpretive framework, the PCR-RFLP assay was shown to provide accurate identification of four co-occurring species (the fifth species, R. hageni, was absent from the test panel, so accuracy cannot yet be extended to training data). The assay is cost- and time-efficient, and will help improve knowledge of Reticulitermes species distributions.

Naturally rare versus newly rare: demographic inferences on two timescales inform conservation of Galápagos giant tortoises

Long-term population history can influence the genetic effects of recent bottlenecks. Therefore, for threatened or endangered species, an understanding of the past is relevant when formulating conservation strategies. Levels of variation at neutral markers have been useful for estimating local effective population sizes (N e ) and inferring whether population sizes increased or decreased over time. Furthermore, analyses of genotypic, allelic frequency, and phylogenetic information can potentially be used to separate historical from recent demographic changes. For 15 populations of Galápagos giant tortoises (Chelonoidis sp.), we used 12 microsatellite loci and DNA sequences from the mitochondrial control region and a nuclear intron, to reconstruct demographic history on shallow (past ∽100 generations, ∽2500 years) and deep (pre-Holocene, >10 thousand years ago) timescales. At the deep timescale, three populations showed strong signals of growth, but with different magnitudes and timing, indicating different underlying causes. Furthermore, estimated historical N e of populations across the archipelago showed no correlation with island age or size, underscoring the complexity of predicting demographic history a priori. At the shallow timescale, all populations carried some signature of a genetic bottleneck, and for 12 populations, point estimates of contemporary N e were very small (i.e., < 50). On the basis of the comparison of these genetic estimates with published census size data, N e generally represented ∽0.16 of the census size. However, the variance in this ratio across populations was considerable. Overall, our data suggest that idiosyncratic and geographically localized forces shaped the demographic history of tortoise populations. Furthermore, from a conservation perspective, the separation of demographic events occurring on shallow versus deep timescales permits the identification of naturally rare versus newly rare populations; this distinction should facilitate prioritization of management action.

Environmental complexity and biodiversity: the multi-layered evolutionary history of a log-dwelling velvet worm in Montane Temperate Australia

Phylogeographic studies provide a framework for understanding the importance of intrinsic versus extrinsic factors in shaping patterns of biodiversity through identifying past and present microevolutionary processes that contributed to lineage divergence. Here we investigate population structure and diversity of the Onychophoran (velvet worm) Euperipatoides rowelli in southeastern Australian montane forests that were not subject to Pleistocene glaciations, and thus likely retained more forest cover than systems under glaciation. Over a ~100 km transect of structurally-connected forest, we found marked nuclear and mitochondrial (mt) DNA genetic structuring, with spatially-localised groups. Patterns from mtDNA and nuclear data broadly corresponded with previously defined geographic regions, consistent with repeated isolation in refuges during Pleistocene climatic cycling. Nevertheless, some E. rowelli genetic contact zones were displaced relative to hypothesized influential landscape structures, implying more recent processes overlying impacts of past environmental history. Major impacts at different timescales were seen in the phylogenetic relationships among mtDNA sequences, which matched geographic relationships and nuclear data only at recent timescales, indicating historical gene flow and/or incomplete lineage sorting. Five major E. rowelli phylogeographic groups were identified, showing substantial but incomplete reproductive isolation despite continuous habitat. Regional distinctiveness, in the face of lineages abutting within forest habitat, could indicate pre- and/or postzygotic gene flow limitation. A potentially functional phenotypic character, colour pattern variation, reflected the geographic patterns in the molecular data. Spatial-genetic patterns broadly match those in previously-studied, co-occurring low-mobility organisms, despite a variety of life histories. We suggest that for E. rowelli, the complex topography and history of the region has led to interplay among limited dispersal ability, historical responses to environmental change, local adaptation, and some resistance to free admixture at geographic secondary contact, leading to strong genetic structuring at fine spatial scale.

Comparative phylogeography of two sympatric beeches in subtropical China: Species-specific geographic mosaic of lineages

In subtropical China, large-scale phylogeographic comparisons among multiple sympatric plants with similar ecological preferences are scarce, making generalizations about common response to historical events necessarily tentative. A phylogeographic comparison of two sympatric Chinese beeches (Fagus lucida and F. longipetiolata, 21 and 28 populations, respectively) was conducted to test whether they have responded to historical events in a concerted fashion and to determine whether their phylogeographic structure is exclusively due to Quaternary events or it is also associated with pre-Quaternary events. Twenty-three haplotypes were recovered for F. lucida and F. longipetiolata (14 each one and five shared). Both species exhibited a species-specific mosaic distribution of haplotypes, with many of them being range-restricted and even private to populations. The two beeches had comparable total haplotype diversity but F. lucida had much higher within-population diversity than F. longipetiolata. Molecular dating showed that the time to most recent common ancestor of all haplotypes was 6.36 Ma, with most haplotypes differentiating during the Quaternary. [Correction added on 14 October 2013, after first online publication: the timeunit has been corrected to '6.36'.] Our results support a late Miocene origin and southwards colonization of Chinese beeches when the aridity in Central Asia intensified and the monsoon climate began to dominate the East Asia. During the Quaternary, long-term isolation in subtropical mountains of China coupled with limited gene flow would have lead to the current species-specific mosaic distribution of lineages.

Pleistocene survival on central Alpine nunataks: genetic evidence from the jumping bristletail Machilis pallida

Mechanisms of survival during the Pleistocene glaciation periods have been studied for more than a century. Until now, molecular studies that confirmed animal survival on Alpine nunataks, that is, ice-free summits surrounded by glaciers, were restricted to peripheral areas. Here, we search for molecular signatures of inner-Alpine survival of the narrow endemic and putatively parthenogenetic Alpine jumping bristletail Machilis pallida combining mitochondrial and AFLP data from its three known populations. The mitochondrial data indicate survival on both peripheral and central nunataks, the latter suggesting that refugia in the centre of the Alpine main ridge were more widespread than previously recognized. Incongruences between mitochondrial and AFLP patterns suggest a complex evolutionary history of the species and may be explained via parallel fixation of parthenogenesis of different origins during the last glacial maximum. We suggest that the inferred parthenogenesis may have been essential for central nunatak survival, but may pose a serious threat for M. pallida in consideration of the present climatic changes.

A species assemblage approach to comparative phylogeography of birds in southern Australia

We present a novel approach to investigating the divergence history of biomes and their component species using single-locus data prior to investing in multilocus data. We use coalescent-based hierarchical approximate Bayesian computation (HABC) methods (MsBayes) to estimate the number and timing of discrete divergences across a putative barrier and to assign species to their appropriate period of co-divergence. We then apply a coalescent-based full Bayesian model of divergence (IMa) to suites of species shown to have simultaneously diverged. The full Bayesian model results in reduced credibility intervals around divergence times and allows other parameters associated with divergence to be summarized across species assemblages. We apply this approach to 10 bird species that are wholly or patchily discontinuous in semi-arid habitats between Australia's southwest (SW) and southeast (SE) mesic zones. There was substantial support for up to three discrete periods of divergence. HABC indicates that two species wholly restricted to more mesic habitats diverged earliest, between 594,382 and 3,417,699 years ago, three species from semi-arid habitats diverged between 0 and 1,508,049 years ago, and four diverged more recently, between 0 and 396,843 years ago. Eight species were assigned to three periods of co-divergence with confidence. For full Bayesian analyses, we accounted for uncertainty in the two remaining species by analyzing all possible suites of species. Estimates of divergence times from full Bayesian divergence models ranged between 429,105 and 2,006,355; 67,172 and 663,837; and 24,607 and 171,085 for the earliest, middle, and most recent periods of co-divergence, respectively. This single-locus approach uses the power of multitaxa coalescent analyses as an efficient means of generating a foundation for further, targeted research using multilocus and genomic tools applied to an understudied biome.

Large-scale analysis of microRNA evolution

BACKGROUND

In animals, microRNAs (miRNA) are important genetic regulators. Animal miRNAs appear to have expanded in conjunction with an escalation in complexity during early bilaterian evolution. Their small size and high-degree of similarity makes them challenging for phylogenetic approaches. Furthermore, genomic locations encoding miRNAs are not clearly defined in many species. A number of studies have looked at the evolution of individual miRNA families. However, we currently lack resources for large-scale analysis of miRNA evolution.

RESULTS

We addressed some of these issues in order to analyse the evolution of miRNAs. We perform syntenic and phylogenetic analysis for miRNAs from 80 animal species. We present synteny maps, phylogenies and functional data for miRNAs across these species. These data represent the basis of our analyses and also act as a resource for the community.

CONCLUSIONS

We use these data to explore the distribution of miRNAs across phylogenetic space, characterise their birth and death, and examine functional relationships between miRNAs and other genes. These data confirm a number of previously reported findings on a larger scale and also offer novel insights into the evolution of the miRNA repertoire in animals, and it's genomic organization.

Assessment of phylogenetic structure in genome size – gene content correlations

Gene content and gene-coding percentage can be predicted from genome size in newly sequenced organisms. Here, we investigate whether these predictions are influenced by phylogenetic relationships between the involved species. Combining a highly resolved phylogenetic tree with a large compilation of gene content data, our results reveal the presence of significant phylogenetic structure in the correlations between genome size and gene content in both bacteria and eukaryotes. The variation in log(gene content) explained by log(genome size) in combination with phylogeny was found to be 97% in bacteria and 55% in eukaryotes. Further, in bacteria, gene-coding percentages are only significantly correlated to genome size if phylogenetic information is taken into account in the analyses. These findings support the usage of phylogenetic correlation models for gene content predictions.

Phylogeography of Saproxylic and Forest Floor Invertebrates from Tallaganda, South-eastern Australia

The interaction between physiogeographic landscape context and certain life history characteristics, particularly dispersal ability, can generate predictable outcomes for how species responded to Pleistocene (and earlier) climatic changes. Furthermore, the extent to which impacts of past landscape-level changes ‘scale-up’ to whole communities has begun to be addressed via comparative phylogeographic analyses of co-distributed species. Here we present an overview of a body of research on flightless low-mobility forest invertebrates, focusing on two springtails and two terrestrial flatworms, from Tallaganda on the Great Dividing Range of south-eastern Australia. These species are distantly-related, and represent contrasting trophic levels (i.e., slime-mold-grazers vs. higher-level predators). However, they share an association with the dead wood (saproxylic) habitat. Spatial patterns of intraspecific genetic diversity partly conform to topography-based divisions that circumscribe five ‘microgeographic regions’ at Tallaganda. In synthesizing population processes and past events that generated contemporary spatial patterns of genetic diversity in these forest floor invertebrates, we highlight cases of phylogeographic congruence, pseudo-congruence, and incongruence. Finally, we propose conservation-oriented recommendations for the prioritisation of areas for protection.