Extant and extinct bilby genomes combined with Indigenous knowledge improve conservation of a unique Australian marsupial.

Ninu (greater bilby, Macrotis lagotis) are desert-dwelling, culturally and ecologically important marsupials. In collaboration with Indigenous rangers and conservation managers, we generated the Ninu chromosome-level genome assembly (3.66 Gbp) and genome sequences for the extinct Yallara (lesser bilby, Macrotis leucura). We developed and tested a scat single-nucleotide polymorphism panel to inform current and future conservation actions, undertake ecological assessments and improve our understanding of Ninu genetic diversity in managed and wild populations. We also assessed the beneficial impact of translocations in the metapopulation (N = 363 Ninu). Resequenced genomes (temperate Ninu, 6; semi-arid Ninu, 6; and Yallara, 4) revealed two major population crashes during global cooling events for both species and differences in Ninu genes involved in anatomical and metabolic pathways. Despite their 45-year captive history, Ninu have fewer long runs of homozygosity than other larger mammals, which may be attributable to their boom-bust life history. Here we investigated the unique Ninu biology using 12 tissue transcriptomes revealing expression of all 115 conserved eutherian chorioallantoic placentation genes in the uterus, an XY1Y2 sex chromosome system and olfactory receptor gene expansions. Together, we demonstrate the holistic value of genomics in improving key conservation actions, understanding unique biological traits and developing tools for Indigenous rangers to monitor remote wild populations.

Pervasive relaxed selection in termite genomes.

Genetic changes that enabled the evolution of eusociality have long captivated biologists. More recently, attention has focussed on the consequences of eusociality on genome evolution. Studies have reported higher molecular evolutionary rates in eusocial hymenopteran insects compared with their solitary relatives. To investigate the genomic consequences of eusociality in termites, we analysed nine genomes, including newly sequenced genomes from three non-eusocial cockroaches. Using a phylogenomic approach, we found that termite genomes have experienced lower rates of synonymous substitutions than those of cockroaches, possibly as a result of longer generation times. We identified higher rates of non-synonymous substitutions in termite genomes than in cockroach genomes, and identified pervasive relaxed selection in the former (24-31% of the genes analysed) compared with the latter (2-4%). We infer that this is due to reductions in effective population size, rather than gene-specific effects (e.g. indirect selection of caste-biased genes). We found no obvious signature of increased genetic load in termites, and postulate efficient purging of deleterious alleles at the colony level. Additionally, we identified genomic adaptations that may underpin caste differentiation, such as genes involved in post-translational modifications. Our results provide insights into the evolution of termites and the genomic consequences of eusociality more broadly.

Novel typing scheme reveals emergence and genetic diversity of Chlamydia pecorum at the local management scale across two koala populations.

To overcome shortcomings in discriminating Chlamydia pecorum strains infecting the koala (Phascolarctos cinereus) at the local level, we developed a novel genotyping scheme for this pathogen to inform koala management at a fine-scale subpopulation level. We applied this scheme to two geographically distinct koala populations in New South Wales, Australia: the Liverpool Plains and the Southern Highlands to South-west Sydney (SHSWS). Our method provides greater resolution than traditional multi-locus sequence typing, and can be used to monitor strain emergence, movement, and divergence across a range of fragmented habitats. Within the Liverpool Plains population, suspected recent introduction of a novel strain was confirmed by an absence of genetic diversity at the earliest sampling events and limited diversity at recent sampling events. Across the partially fragmented agricultural landscape of the Liverpool Plains, diversity within a widespread sequence type suggests that this degree of fragmentation may hinder but not prevent spread. In the SHSWS population, our results suggest movement of a strain from the south, where diverse strains exist, into a previously Chlamydia-free area in the north, indicating the risk of expansion towards an adjacent Chlamydia-negative koala population in South-west Sydney. In the south of the SHSWS where koala subpopulations appear segregated, we found evidence of divergent strain evolution. Our tool can be used to infer the risks of strain introduction across fragmented habitats in population management, particularly through practices such as wildlife corridor constructions and translocations.

Complexity of avian evolution revealed by family-level genomes.

Despite tremendous efforts in the past decades, relationships among main avian lineages remain heavily debated without a clear resolution. Discrepancies have been attributed to diversity of species sampled, phylogenetic method and the choice of genomic regions1-3. Here we address these issues by analysing the genomes of 363 bird species4 (218 taxonomic families, 92% of total). Using intergenic regions and coalescent methods, we present a well-supported tree but also a marked degree of discordance. The tree confirms that Neoaves experienced rapid radiation at or near the Cretaceous-Palaeogene boundary. Sufficient loci rather than extensive taxon sampling were more effective in resolving difficult nodes. Remaining recalcitrant nodes involve species that are a challenge to model due to either extreme DNA composition, variable substitution rates, incomplete lineage sorting or complex evolutionary events such as ancient hybridization. Assessment of the effects of different genomic partitions showed high heterogeneity across the genome. We discovered sharp increases in effective population size, substitution rates and relative brain size following the Cretaceous-Palaeogene extinction event, supporting the hypothesis that emerging ecological opportunities catalysed the diversification of modern birds. The resulting phylogenetic estimate offers fresh insights into the rapid radiation of modern birds and provides a taxon-rich backbone tree for future comparative studies.

Shrinking in the dark: Parallel endosymbiont genome erosions are associated with repeated host transitions to an underground life.

Microbial symbioses have had profound impacts on the evolution of animals. Conversely, changes in host biology may impact the evolutionary trajectory of symbionts themselves. Blattabacterium cuenoti is present in almost all cockroach species and enables hosts to subsist on a nutrient-poor diet. To investigate if host biology has impacted Blattabacterium at the genomic level, we sequenced and analyzed 25 genomes from Australian soil-burrowing cockroaches (Blaberidae: Panesthiinae), which have undergone at least seven separate subterranean, subsocial transitions from above-ground, wood-feeding ancestors. We find at least three independent instances of genome erosion have occurred in Blattabacterium strains exclusive to Australian soil-burrowing cockroaches. These shrinkages have involved the repeated inactivation of genes involved in amino acid biosynthesis and nitrogen recycling, the core role of Blattabacterium in the host-symbiont relationship. The most drastic of these erosions have occurred in hosts thought to have transitioned underground the earliest relative to other lineages, further suggestive of a link between gene loss in Blattabacterium and the burrowing behavior of hosts. As Blattabacterium is unable to fulfill its core function in certain host lineages, these findings suggest soil-burrowing cockroaches must acquire these nutrients from novel sources. Our study represents one of the first cases, to our knowledge, of parallel host adaptations leading to concomitant parallelism in their mutualistic symbionts, further underscoring the intimate relationship between these two partners.

The spotted parrotfish genome provides insights into the evolution of a coral reef dietary specialist (Teleostei: Labridae: Scarini: Cetoscarus ocellatus).

With over 600 valid species, the wrasses (family Labridae) are among the largest and most successful families of the marine teleosts. They feature prominently on coral reefs where they are known not only for their impressive diversity in colouration and form but also for their functional specialisation and ability to occupy a wide variety of trophic guilds. Among the wrasses, the parrotfishes (tribe Scarini) display some of the most dramatic examples of trophic specialisation. Using abrasion-resistant biomineralized teeth, parrotfishes are able to mechanically extract protein-rich micro-photoautotrophs growing in and among reef carbonate material, a dietary niche that is inaccessible to most other teleost fishes. This ability to exploit an otherwise untapped trophic resource is thought to have played a role in the diversification and evolutionary success of the parrotfishes. In order to better understand the key evolutionary innovations leading to the success of these dietary specialists, we sequenced and analysed the genome of a representative species, the spotted parrotfish (Cetoscarus ocellatus). We find significant expansion, selection and duplications within several detoxification gene families and a novel poly-glutamine expansion in the enamel protein ameloblastin, and we consider their evolutionary implications. Our genome provides a useful resource for comparative genomic studies investigating the evolutionary history of this highly specialised teleostean radiation.

The Neotropical endemic liverwort subfamily Micropterygioideae had circum-Antarctic links to the rest of the Lepidoziaceae during the early Cretaceous.

Lepidoziaceae are the third-largest family of liverworts, with about 860 species distributed on all continents. The evolutionary history of this family has not been satisfactorily resolved, with taxa such as Micropterygioideae yet to be included in phylogenetic analyses. We inferred a dated phylogeny of Lepidoziaceae using a data set consisting of 13 genetic markers, sampled from 147 species. Based on our phylogenetic estimate, we used statistical dispersal-vicariance analysis to reconstruct the biogeographic history of the family. We inferred a crown age of 197 Ma (95% credible interval 157-240 Ma) for the family in the Australian region, with most major lineages also originating in the same region. Micropterygioideae are placed as the sister group to Lembidioideae, with these two lineages diverging from each other about 132 Ma in the South American-Australian region. With South America and Australia being connected through Antarctica at the time, our results suggest a circum-Antarctic link between Micropterygioideae and the rest of the family. Crown Micropterygioideae were inferred to have arisen 45 Ma in South America before the continent separated from Antarctica. Extinction from southern temperate regions might explain the present-day restriction of Micropterygioideae to the Neotropics. Our study reveals the influence of past geological events, such as continental drift, on the evolution and distribution of a widespread and diverse family of liverworts.

ClockstaRX: Testing Molecular Clock Hypotheses With Genomic Data.

Phylogenomic data provide valuable opportunities for studying evolutionary rates and timescales. These analyses require theoretical and statistical tools based on molecular clocks. We present ClockstaRX, a flexible platform for exploring and testing evolutionary rate signals in phylogenomic data. Here, information about evolutionary rates in branches across gene trees is placed in Euclidean space, allowing data transformation, visualization, and hypothesis testing. ClockstaRX implements formal tests for identifying groups of loci and branches that make a large contribution to patterns of rate variation. This information can then be used to test for drivers of genomic evolutionary rates or to inform models for molecular dating. Drawing on the results of a simulation study, we recommend forms of data exploration and filtering that might be useful prior to molecular-clock analyses.

Dating in the Dark: Elevated Substitution Rates in Cave Cockroaches (Blattodea: Nocticolidae) Have Negative Impacts on Molecular Date Estimates.

Rates of nucleotide substitution vary substantially across the Tree of Life, with potentially confounding effects on phylogenetic and evolutionary analyses. A large acceleration in mitochondrial substitution rate occurs in the cockroach family Nocticolidae, which predominantly inhabit subterranean environments. To evaluate the impacts of this among-lineage rate heterogeneity on estimates of phylogenetic relationships and evolutionary timescales, we analyzed nuclear ultraconserved elements (UCEs) and mitochondrial genomes from nocticolids and other cockroaches. Substitution rates were substantially elevated in nocticolid lineages compared with other cockroaches, especially in mitochondrial protein-coding genes. This disparity in evolutionary rates is likely to have led to different evolutionary relationships being supported by phylogenetic analyses of mitochondrial genomes and UCE loci. Furthermore, Bayesian dating analyses using relaxed-clock models inferred much deeper divergence times compared with a flexible local clock. Our phylogenetic analysis of UCEs, which is the first genome-scale study to include all 13 major cockroach families, unites Corydiidae and Nocticolidae and places Anaplectidae as the sister lineage to the rest of Blattoidea. We uncover an extraordinary level of genetic divergence in Nocticolidae, including two highly distinct clades that separated ~115 million years ago despite both containing representatives of the genus Nocticola. The results of our study highlight the potential impacts of high among-lineage rate variation on estimates of phylogenetic relationships and evolutionary timescales.

Evaluating the Accuracy of Methods for Detecting Correlated Rates of Molecular and Morphological Evolution.

Determining the link between genomic and phenotypic change is a fundamental goal in evolutionary biology. Insights into this link can be gained by using a phylogenetic approach to test for correlations between rates of molecular and morphological evolution. However, there has been persistent uncertainty about the relationship between these rates, partly because conflicting results have been obtained using various methods that have not been examined in detail. We carried out a simulation study to evaluate the performance of 5 statistical methods for detecting correlated rates of evolution. Our simulations explored the evolution of molecular sequences and morphological characters under a range of conditions. Of the methods tested, Bayesian relaxed-clock estimation of branch rates was able to detect correlated rates of evolution correctly in the largest number of cases. This was followed by correlations of root-to-tip distances, Bayesian model selection, independent sister-pairs contrasts, and likelihood-based model selection. As expected, the power to detect correlated rates increased with the amount of data, both in terms of tree size and number of morphological characters. Likewise, greater among-lineage rate variation in the data led to improved performance of all 5 methods, particularly for Bayesian relaxed-clock analysis when the rate model was mismatched. We then applied these methods to a data set from flowering plants and did not find evidence of a correlation in evolutionary rates between genomic data and morphological characters. The results of our study have practical implications for phylogenetic analyses of combined molecular and morphological data sets, and highlight the conditions under which the links between genomic and phenotypic rates of evolution can be evaluated quantitatively.

Impacts of Taxon-Sampling Schemes on Bayesian Tip Dating Under the Fossilized Birth-Death Process.

Evolutionary timescales can be inferred by molecular-clock analyses of genetic data and fossil evidence. Bayesian phylogenetic methods such as tip dating provide a powerful framework for inferring evolutionary timescales, but the most widely used priors for tree topologies and node times often assume that present-day taxa have been sampled randomly or exhaustively. In practice, taxon sampling is often carried out so as to include representatives of major lineages, such as orders or families. We examined the impacts of different densities of diversified sampling on Bayesian tip dating on unresolved fossilized birth-death (FBD) trees, in which fossil taxa are topologically constrained but their exact placements are averaged out. We used synthetic data generated by simulations of nucleotide sequence evolution, fossil occurrences, and diversified taxon sampling. Our analyses under the diversified-sampling FBD process show that increasing taxon-sampling density does not necessarily improve divergence-time estimates. However, when informative priors were specified for the root age or when tree topologies were fixed to those used for simulation, the performance of tip dating on unresolved FBD trees maintains its accuracy and precision or improves with taxon-sampling density. By exploring three situations in which models are mismatched, we find that including all relevant fossils, without pruning off those that are incompatible with the diversified-sampling FBD process, can lead to underestimation of divergence times. Our reanalysis of a eutherian mammal data set confirms some of the findings from our simulation study, and reveals the complexity of diversified taxon sampling in phylogenomic data sets. In highlighting the interplay of taxon-sampling density and other factors, the results of our study have practical implications for using Bayesian tip dating to infer evolutionary timescales across the Tree of Life. [Bayesian tip dating; eutherian mammals; fossilized birth-death process; phylogenomics; taxon sampling.].

Genomic analysis reveals strong population structure in the Giant Sydney Crayfish (Euastacus spinifer (Heller, 1865)).

Australia is home to over 140 species of freshwater crayfish (Decapoda: Parastacidae), representing a centre of diversity for this group in the Southern Hemisphere. Species delimitation in freshwater crayfish is difficult because many species show significant variation in colouration and morphology. This is particularly evident in the genus Euastacus, which exhibits large variations in colour and spination throughout its putative range. To understand this variation, we investigated the genetic diversity, population structure, phylogeny, and evolutionary timescale of the Giant Sydney Crayfish (Euastacus spinifer (Heller, 1865)). Our data set is sampled from over 70 individuals from across the ∼600 km range of the species, and includes a combination of two mitochondrial markers and more than 7000 single-nucleotide polymorphisms (SNPs) from the nuclear genome. Data were also obtained for representatives of the close relative, Euastacus vesper McCormack and Ahyong, 2017. Genomic SNP analyses revealed strong population structure, with multiple distinct populations showing little evidence of gene flow or migration. Phylogenetic analyses of mitochondrial data revealed similar structure between populations. Taken together, our analyses suggest that E. spinifer, as currently understood, represents a species complex, of which E. vesper is a member. Molecular clock estimates place the divergences within this group during the Pleistocene. The isolated and highly fragmented populations identified in our analyses probably represent relict populations of a previously widespread ancestral species. Periodic flooding events during the Pleistocene are likely to have facilitated the movement of these otherwise restricted freshwater crayfish within and between drainage basins, including the Murray-Darling and South East Coast Drainages. We present evidence supporting the recognition of populations in the southern parts of the range of E. spinifer as one or two separate species, which would raise the number of species within the E. spinifer complex to at least three. Our results add to the growing body of evidence that many freshwater crayfish exhibit highly fragmented, range-restricted distributions. In combination with the life-history traits of these species, the restricted distributions exacerbate the threats already placed on freshwater crayfish, which are among the five most endangered animal groups globally.

The recombinogenic history of turnip mosaic potyvirus reveals its introduction to Japan in the 19th century.

Characterizing the detailed spatial and temporal dynamics of plant pathogens can provide valuable information for crop protection strategies. However, the epidemiological characteristics and evolutionary trajectories of pathogens can differ markedly from one country to another. The most widespread and important virus of brassica vegetables, turnip mosaic virus (TuMV), causes serious plant diseases in Japan. We collected 317 isolates of TuMV from Raphanus and Brassica plants throughout Japan over nearly five decades. Genomic sequences from these isolates were combined with published sequences. We identified a total of eighty-eight independent recombination events in Japanese TuMV genomes and found eighty-two recombination-type patterns in Japan. We assessed the evolution of TuMV through space and time using whole and partial genome sequences of both nonrecombinants and recombinants. Our results suggest that TuMV was introduced into Japan after the country emerged from its isolationist policy (1639-1854) in the Edo period and then dispersed to other parts of Japan in the 20th century. The results of our analyses reveal the complex structure of the TuMV population in Japan and emphasize the importance of identifying recombination events in the genome. Our study also provides an example of surveying the epidemiology of a virus that is highly recombinogenic.

Early diversifications of angiosperms and their insect pollinators: were they unlinked?

The present-day ubiquity of angiosperm-insect pollination has led to the hypothesis that these two groups coevolved early in their evolutionary history. However, recent fossil discoveries and fossil-calibrated molecular dating analyses challenge the notion that early diversifications of angiosperms and insects were inextricably linked. In this article, we examine (i) the discrepancies between dates of emergence for angiosperms and major clades of insects; (ii) the long history of gymnosperm-insect pollination modes, which likely shaped early angiosperm-insect pollination mutualisms; and (iii) how the K-Pg (Cretaceous-Paleogene) mass extinction event was vital in propelling modern angiosperm-insect mutualisms. We posit that the early diversifications of angiosperms and their insect pollinators were largely decoupled until the end of the Cretaceous.

Redescription of Euastacus clydensis Riek, 1969 (Crustacea: Parastacidae), a valid species of spiny crayfish from southern New South Wales, Australia.

The Giant Sydney Crayfish (Euastacus spinifer (Heller, 1865)) was thought to have a wide range in New South Wales, Australia, spanning some 600 km north-south. A recent extensive molecular phylogenetic and population genomic analysis of E. spinifer across its geographical range revealed strong population structure corresponding to several major geographically correlated clades, the southernmost clade being the most genetically divergent and clearly a separate species. This southern clade corresponds to the junior synonym E. clydensis Riek, 1969 and is sister to the clade comprising the remaining populations of E. spinifer and Euastacus vesper. We formally remove E. clydensis from the synonymy of E. spinifer, increasing the recognised number of species of Euastacus to 54. Euastacus clydensis is redescribed based on type and other material, and is distinguished from E. spinifer by differences in abdominal spination and the form of the antennal scaphocerite. Euastacus clydensis has a restricted southern New South Wales range in the Shoalhaven and Jervis Bay-Clyde River catchments, from Moss Vale south to the vicinity of Clyde Mountain; much of the known range of E. clydensis was burnt in the 2019-2020 eastern Australian megafires.

Evolutionary Rate Variation among Lineages in Gene Trees has a Negative Impact on Species-Tree Inference.

Phylogenetic analyses of genomic data provide a powerful means of reconstructing the evolutionary relationships among organisms, yet such analyses are often hindered by conflicting phylogenetic signals among loci. Identifying the signals that are most influential to species-tree estimation can help to inform the choice of data for phylogenomic analysis. We investigated this in an analysis of 30 phylogenomic data sets. For each data set, we examined the association between several branch-length characteristics of gene trees and the distance between these gene trees and the corresponding species trees. We found that the distance of each gene tree to the species tree inferred from the full data set was positively associated with variation in root-to-tip distances and negatively associated with mean branch support. However, no such associations were found for gene-tree length, a measure of the overall substitution rate at each locus. We further explored the usefulness of the best-performing branch-based characteristics for selecting loci for phylogenomic analyses. We found that loci that yield gene trees with high variation in root-to-tip distances have a disproportionately distant signal of tree topology compared with the complete data sets. These results suggest that rate variation across lineages should be taken into consideration when exploring and even selecting loci for phylogenomic analysis.[Branch support; data filtering; nucleotide substitution model; phylogenomics; substitution rate; summary coalescent methods.].

Excluding Loci With Substitution Saturation Improves Inferences From Phylogenomic Data.

The historical signal in nucleotide sequences becomes eroded over time by substitutions occurring repeatedly at the same sites. This phenomenon, known as substitution saturation, is recognized as one of the primary obstacles to deep-time phylogenetic inference using genome-scale data sets. We present a new test of substitution saturation and demonstrate its performance in simulated and empirical data. For some of the 36 empirical phylogenomic data sets that we examined, we detect substitution saturation in around 50% of loci. We found that saturation tends to be flagged as problematic in loci with highly discordant phylogenetic signals across sites. Within each data set, the loci with smaller numbers of informative sites are more likely to be flagged as containing problematic levels of saturation. The entropy saturation test proposed here is sensitive to high evolutionary rates relative to the evolutionary timeframe, while also being sensitive to several factors known to mislead phylogenetic inference, including short internal branches relative to external branches, short nucleotide sequences, and tree imbalance. Our study demonstrates that excluding loci with substitution saturation can be an effective means of mitigating the negative impact of multiple substitutions on phylogenetic inferences. [Phylogenetic model performance; phylogenomics; substitution model; substitution saturation; test statistics.].

Philippine Ayta possess the highest level of Denisovan ancestry in the world.

Multiple lines of evidence show that modern humans interbred with archaic Denisovans. Here, we report an account of shared demographic history between Australasians and Denisovans distinctively in Island Southeast Asia. Our analyses are based on ∼2.3 million genotypes from 118 ethnic groups of the Philippines, including 25 diverse self-identified Negrito populations, along with high-coverage genomes of Australopapuans and Ayta Magbukon Negritos. We show that Ayta Magbukon possess the highest level of Denisovan ancestry in the world-∼30%-40% greater than that of Australians and Papuans-consistent with an independent admixture event into Negritos from Denisovans. Together with the recently described Homo luzonensis, we suggest that there were multiple archaic species that inhabited the Philippines prior to the arrival of modern humans and that these archaic groups may have been genetically related. Altogether, our findings unveil a complex intertwined history of modern and archaic humans in the Asia-Pacific region, where distinct Islander Denisovan populations differentially admixed with incoming Australasians across multiple locations and at various points in time.