Environmental drivers and cryptic biodiversity hotspots define endophytes in Earth's largest terrestrial biome.

Understanding how symbiotic associations differ across environmental gradients is key to predicting the fate of symbioses as environments change, and it is vital for detecting global reservoirs of symbiont biodiversity in a changing world.1,2,3 However, sampling of symbiotic partners at the full-biome scale is difficult and rare. As Earth's largest terrestrial biome, boreal forests influence carbon dynamics and climate regulation at a planetary scale. Plants and lichens in this biome host the highest known phylogenetic diversity of fungal endophytes, which occur within healthy photosynthetic tissues and can influence hosts' resilience to stress.4,5 We examined how communities of endophytes are structured across the climate gradient of the boreal biome, focusing on the dominant plant and lichen species occurring across the entire south-to-north span of the boreal zone in eastern North America. Although often invoked for understanding the distribution of biodiversity, neither a latitudinal gradient nor mid-domain effect5,6,7 can explain variation in endophyte diversity at this trans-biome scale. Instead, analyses considering shifts in forest characteristics, Picea biomass and age, and nutrients in host tissues from 46° to 58° N reveal strong and distinctive signatures of climate in defining endophyte assemblages in each host lineage. Host breadth of endophytes varies with climate factors, and biodiversity hotspots can be identified at plant-community transitions across the boreal zone at a global scale. Placed against a backdrop of global circumboreal sampling,4 our study reveals the sensitivity of endophytic fungi, their reservoirs of biodiversity, and their important symbiotic associations, to climate.

An inspection-based assessment of obstacles to salmon, trout, eel and lamprey migration and river channel connectivity in Ireland.

Knowledge of the location, physical attributes and impacts of obstacles on river connectivity is a requirement for any mitigating action aimed at restoring the connectivity of a river system. Here, we present a study that recorded the numbers and physical diversity of obstacles in 10 river catchments in Ireland, together with the impact these structures had on overall river connectivity. A total of 372 obstacles were recorded, 3 of these were dams, and the remainder were low-head weirs/sluices, obstacles associated with road or rail crossings of rivers and natural structures. The degree of fragmentation was estimated in each catchment by calculating obstacle density and the Dendritic Connectivity Index (DCI). DCI scores were calculated for 4 native Irish fish species with different life-histories, namely diadromous (Atlantic salmon, sea trout, European eel, sea lamprey) and potamodromous (brown trout). Obstacle density ranged between 1.2 and 0.02 obstacles/km of river. Six of the 10 catchments had at least one obstacle located on the mainstem river at least 5 km from its mouth/confluence. These 6 catchments typically had the lowest connectivity scores for diadromous species and ranged between 0.6 and 44.1 (a fully connected river would receive a maximum score of 100). While there was no significant correlation between obstacle density and the DCI score for diadromous fish, a significant negative correlation was detected between obstacle density and the DCI score for potamodromous brown trout. Here, we highlight the merit of these obstacle assessments and associated challenges for decision-making relating to prioritisation of obstacles for removal or modification.

Amplicon sequencing of 42 nuclear loci supports directional gene flow between South Pacific populations of a hydrothermal vent limpet.

In the past few decades, population genetics and phylogeographic studies have improved our knowledge of connectivity and population demography in marine environments. Studies of deep-sea hydrothermal vent populations have identified barriers to gene flow, hybrid zones, and demographic events, such as historical population expansions and contractions. These deep-sea studies, however, used few loci, which limit the amount of information they provided for coalescent analysis and thus our ability to confidently test complex population dynamics scenarios. In this study, we investigated population structure, demographic history, and gene flow directionality among four Western Pacific hydrothermal vent populations of the vent limpet Lepetodrilus aff. schrolli. These vent sites are located in the Manus and Lau back-arc basins, currently of great interest for deep-sea mineral extraction. A total of 42 loci were sequenced from each individual using high-throughput amplicon sequencing. Amplicon sequences were analyzed using both genetic variant clustering methods and evolutionary coalescent approaches. Like most previously investigated vent species in the South Pacific, L. aff. schrolli showed no genetic structure within basins but significant differentiation between basins. We inferred significant directional gene flow from Manus Basin to Lau Basin, with low to no gene flow in the opposite direction. This study is one of the very few marine population studies using >10 loci for coalescent analysis and serves as a guide for future marine population studies.

The application of CRISPR-Cas for single species identification from environmental DNA.

We report the first application of CRISPR-Cas technology to single species detection from environmental DNA (eDNA). Organisms shed and excrete DNA into their environment such as in skin cells and faeces, referred to as environmental DNA (eDNA). Utilising eDNA allows noninvasive monitoring with increased specificity and sensitivity. Current methods primarily employ PCR-based techniques to detect a given species from eDNA samples, posing a logistical challenge for on-site monitoring and potential adaptation to biosensor devices. We have developed an alternative method; coupling isothermal amplification to a CRISPR-Cas12a detection system. This utilises the collateral cleavage activity of Cas12a, a ribonuclease guided by a highly specific single CRISPR RNA. We used the target species Salmo salar as a proof-of-concept test of the specificity of the assay among closely related species and to show the assay is successful at a single temperature of 37°C with signal detection at 535 nM. The specific assay, detects at attomolar sensitivity with rapid detection rates (<2.5 hr). This approach simplifies the challenge of building a biosensor device for rapid target species detection in the field and can be easily adapted to detect any species from eDNA samples from a variety of sources enhancing the capabilities of eDNA as a tool for monitoring biodiversity.

Phylogenetic placement, species delimitation, and cyanobiont identity of endangered aquatic Peltigera species (lichen-forming Ascomycota, Lecanoromycetes).

• Premise of this study: Aquatic cyanolichens from the genus Peltigera section Hydrothyriae are subject to anthropogenic threats and, therefore, are considered endangered. In this study we addressed the phylogenetic placement of section Hydrothyriae within Peltigera. We delimited species within the section and identified their symbiotic cyanobacteria.• Methods: Species delimitation and population structure were explored using monophyly as a grouping criterion (RAxML) and Structurama based on three protein-coding genes in combination with two nuclear ribosomal loci. The 16S and rbcLX sequences for the cyanobionts were analyzed in the broad phylogenetic context of free-living and symbiotic cyanobacteria.• Key results: We confirm with high confidence the placement of section Hydrothyriae within the monophyletic genus Peltigera; however, its phylogenetic position within the genus remains unsettled. We recovered three distinct monophyletic groups corresponding to three species: P. hydrothyria, P. gowardii s.s., and P. aquatica Miadl. & Lendemer, the latter being formally introduced here. Each species was associated with an exclusive set of Nostoc haplotypes.• Conclusions: The ITS region alone provides sufficient genetic information to distinguish the three morphologically cryptic species within section Hydrothyriae. Section Hydrothyriae seems to be associated with a monophyletic lineage of Nostoc, that has not been found in symbiotic association with other members of Peltigera. Capsosira lowei should be transferred to the genus Nostoc. Potential threats to P. aquatica should be re-examined based on the recognition of two aquatic species in western North America.

Arthropod relationships revealed by phylogenomic analysis of nuclear protein-coding sequences.

The remarkable antiquity, diversity and ecological significance of arthropods have inspired numerous attempts to resolve their deep phylogenetic history, but the results of two decades of intensive molecular phylogenetics have been mixed. The discovery that terrestrial insects (Hexapoda) are more closely related to aquatic Crustacea than to the terrestrial centipedes and millipedes (Myriapoda) was an early, if exceptional, success. More typically, analyses based on limited samples of taxa and genes have generated results that are inconsistent, weakly supported and highly sensitive to analytical conditions. Here we present strongly supported results from likelihood, Bayesian and parsimony analyses of over 41 kilobases of aligned DNA sequence from 62 single-copy nuclear protein-coding genes from 75 arthropod species. These species represent every major arthropod lineage, plus five species of tardigrades and onychophorans as outgroups. Our results strongly support Pancrustacea (Hexapoda plus Crustacea) but also strongly favour the traditional morphology-based Mandibulata (Myriapoda plus Pancrustacea) over the molecule-based Paradoxopoda (Myriapoda plus Chelicerata). In addition to Hexapoda, Pancrustacea includes three major extant lineages of 'crustaceans', each spanning a significant range of morphological disparity. These are Oligostraca (ostracods, mystacocarids, branchiurans and pentastomids), Vericrustacea (malacostracans, thecostracans, copepods and branchiopods) and Xenocarida (cephalocarids and remipedes). Finally, within Pancrustacea we identify Xenocarida as the long-sought sister group to the Hexapoda, a result confirming that 'crustaceans' are not monophyletic. These results provide a statistically well-supported phylogenetic framework for the largest animal phylum and represent a step towards ending the often-heated, century-long debate on arthropod relationships.

Resolving arthropod phylogeny: exploring phylogenetic signal within 41 kb of protein-coding nuclear gene sequence.

This study attempts to resolve relationships among and within the four basal arthropod lineages (Pancrustacea, Myriapoda, Euchelicerata, Pycnogonida) and to assess the widespread expectation that remaining phylogenetic problems will yield to increasing amounts of sequence data. Sixty-eight regions of 62 protein-coding nuclear genes (approximately 41 kilobases (kb)/taxon) were sequenced for 12 taxonomically diverse arthropod taxa and a tardigrade outgroup. Parsimony, likelihood, and Bayesian analyses of total nucleotide data generally strongly supported the monophyly of each of the basal lineages represented by more than one species. Other relationships within the Arthropoda were also supported, with support levels depending on method of analysis and inclusion/exclusion of synonymous changes. Removing third codon positions, where the assumption of base compositional homogeneity was rejected, altered the results. Removing the final class of synonymous mutations--first codon positions encoding leucine and arginine, which were also compositionally heterogeneous--yielded a data set that was consistent with a hypothesis of base compositional homogeneity. Furthermore, under such a data-exclusion regime, all 68 gene regions individually were consistent with base compositional homogeneity. Restricting likelihood analyses to nonsynonymous change recovered trees with strong support for the basal lineages but not for other groups that were variably supported with more inclusive data sets. In a further effort to increase phylogenetic signal, three types of data exploration were undertaken. (1) Individual genes were ranked by their average rate of nonsynonymous change, and three rate categories were assigned--fast, intermediate, and slow. Then, bootstrap analysis of each gene was performed separately to see which taxonomic groups received strong support. Five taxonomic groups were strongly supported independently by two or more genes, and these genes mostly belonged to the slow or intermediate categories, whereas groups supported only by a single gene region tended to be from genes of the fast category, arguing that fast genes provide a less consistent signal. (2) A sensitivity analysis was performed in which increasing numbers of genes were excluded, beginning with the fastest. The number of strongly supported nodes increased up to a point and then decreased slightly. Recovery of Hexapoda required removal of fast genes. Support for Mandibulata (Pancrustacea + Myriapoda) also increased, at times to "strong" levels, with removal of the fastest genes. (3) Concordance selection was evaluated by clustering genes according to their ability to recover Pancrustacea, Euchelicerata, or Myriapoda and analyzing the three clusters separately. All clusters of genes recovered the three concordance clades but were at times inconsistent in the relationships recovered among and within these clades, a result that indicates that the a priori concordance criteria may bias phylogenetic signal in unexpected ways. In a further attempt to increase support of taxonomic relationships, sequence data from 49 additional taxa for three slow genes (i.e., EF-1 alpha, EF-2, and Pol II) were combined with the various 13-taxon data sets. The 62-taxon analyses supported the results of the 13-taxon analyses and provided increased support for additional pancrustacean clades found in an earlier analysis including only EF-1 alpha, EF-2, and Pol II.

Diversity in the weapons of sexual selection: horn evolution in the beetle genus Onthophagus (Coleoptera: Scarabaeidae).

Both ornaments and weapons of sexual selection frequently exhibit prolific interspecific diversity of form. Yet, most studies of this diversity have focused on ornaments involved with female mate choice, rather than on the weapons of male competition. With few exceptions, the mechanisms of divergence in weapon morphology remain largely unexplored. Here, we characterize the evolutionary radiation of one type of weapon: beetle horns. We use partial sequences from four nuclear and three mitochondrial genes to develop a phylogenetic hypothesis for a worldwide sample of 48 species from the dung beetle genus Onthophagus (Coleoptera: Scarabaeidae). We then use these data to test for multiple evolutionary origins of horns and to characterize the evolutionary radiation of horns. Although our limited sampling of one of the world's most species-rich genera almost certainly underestimates the number of evolutionary events, our phylogeny reveals prolific evolutionary lability of these exaggerated sexually selected weapons (more than 25 separate gains and losses of five different horn types). We discuss these results in the context of the natural history of these beetles and explore ways that sexual selection and ecology may have interacted to generate this extraordinary diversity of weapon morphology.