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.

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 analysed 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 thirteen 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.

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.

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.

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.

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.

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.

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.