Incomplete lineage sorting and phenotypic evolution in marsupials.

Incomplete lineage sorting (ILS) makes ancestral genetic polymorphisms persist during rapid speciation events, inducing incongruences between gene trees and species trees. ILS has complicated phylogenetic inference in many lineages, including hominids. However, we lack empirical evidence that ILS leads to incongruent phenotypic variation. Here, we performed phylogenomic analyses to show that the South American monito del monte is the sister lineage of all Australian marsupials, although over 31% of its genome is closer to the Diprotodontia than to other Australian groups due to ILS during ancient radiation. Pervasive conflicting phylogenetic signals across the whole genome are consistent with some of the morphological variation among extant marsupials. We detected hundreds of genes that experienced stochastic fixation during ILS, encoding the same amino acids in non-sister species. Using functional experiments, we confirm how ILS may have directly contributed to hemiplasy in morphological traits that were established during rapid marsupial speciation ca. 60 mya.

Robust and prototypical immune responses toward COVID-19 vaccine in First Nations peoples are impacted by comorbidities.

High-risk groups, including Indigenous people, are at risk of severe COVID-19. Here we found that Australian First Nations peoples elicit effective immune responses to COVID-19 BNT162b2 vaccination, including neutralizing antibodies, receptor-binding domain (RBD) antibodies, SARS-CoV-2 spike-specific B cells, and CD4+ and CD8+ T cells. In First Nations participants, RBD IgG antibody titers were correlated with body mass index and negatively correlated with age. Reduced RBD antibodies, spike-specific B cells and follicular helper T cells were found in vaccinated participants with chronic conditions (diabetes, renal disease) and were strongly associated with altered glycosylation of IgG and increased interleukin-18 levels in the plasma. These immune perturbations were also found in non-Indigenous people with comorbidities, indicating that they were related to comorbidities rather than ethnicity. However, our study is of a great importance to First Nations peoples who have disproportionate rates of chronic comorbidities and provides evidence of robust immune responses after COVID-19 vaccination in Indigenous people.

A Pungent and Painful Toxin.

In this issue of Cell, King et al. (2019) have discovered a cell penetrating peptide isolated from the venom of the Australian Black Rock scorpion that activates the TRPA1 receptor in a unique way to induce pain. Their findings offer new insights into how animals evolved venoms to target specific ion channel functions.

An Atypical Parvovirus Drives Chronic Tubulointerstitial Nephropathy and Kidney Fibrosis.

The occurrence of a spontaneous nephropathy with intranuclear inclusions in laboratory mice has puzzled pathologists for over 4 decades, because its etiology remains elusive. The condition is more severe in immunodeficient animals, suggesting an infectious cause. Using metagenomics, we identify the causative agent as an atypical virus, termed "mouse kidney parvovirus" (MKPV), belonging to a divergent genus of Parvoviridae. MKPV was identified in animal facilities in Australia and North America, is transmitted via a fecal-oral or urinary-oral route, and is controlled by the adaptive immune system. Detailed analysis of the clinical course and histopathological features demonstrated a stepwise progression of pathology ranging from sporadic tubular inclusions to tubular degeneration and interstitial fibrosis and culminating in renal failure. In summary, we identify a widely distributed pathogen in laboratory mice and establish MKPV-induced nephropathy as a new tool for elucidating mechanisms of tubulointerstitial fibrosis that shares molecular features with chronic kidney disease in humans.

Elevated Southern Hemisphere moisture availability during glacial periods.

Late Pleistocene ice-age climates are routinely characterized as having imposed moisture stress on low- to mid-latitude ecosystems1-5. This idea is largely based on fossil pollen evidence for widespread, low-biomass glacial vegetation, interpreted as indicating climatic dryness6. However, woody plant growth is inhibited under low atmospheric CO2 (refs. 7,8), so understanding glacial environments requires the development of new palaeoclimate indicators that are independent of vegetation9. Here we show that, contrary to expectations, during the past 350 kyr, peaks in southern Australian climatic moisture availability were largely confined to glacial periods, including the Last Glacial Maximum, whereas warm interglacials were relatively dry. By measuring the timing of speleothem growth in the Southern Hemisphere subtropics, which today has a predominantly negative annual moisture balance, we developed a record of climatic moisture availability that is independent of vegetation and extends through multiple glacial-interglacial cycles. Our results demonstrate that a cool-moist response is consistent across the austral subtropics and, in part, may result from reduced evaporation under cool glacial temperatures. Insofar as cold glacial environments in the Southern Hemisphere subtropics have been portrayed as uniformly arid3,10,11, our findings suggest that their characterization as evolutionary or physiological obstacles to movement and expansion of animal, plant and, potentially, human populations10 should be reconsidered.

Pathogen spillover driven by rapid changes in bat ecology.

During recent decades, pathogens that originated in bats have become an increasing public health concern. A major challenge is to identify how those pathogens spill over into human populations to generate a pandemic threat1. Many correlational studies associate spillover with changes in land use or other anthropogenic stressors2,3, although the mechanisms underlying the observed correlations have not been identified4. One limitation is the lack of spatially and temporally explicit data on multiple spillovers, and on the connections among spillovers, reservoir host ecology and behaviour and viral dynamics. We present 25 years of data on land-use change, bat behaviour and spillover of Hendra virus from Pteropodid bats to horses in subtropical Australia. These data show that bats are responding to environmental change by persistently adopting behaviours that were previously transient responses to nutritional stress. Interactions between land-use change and climate now lead to persistent bat residency in agricultural areas, where periodic food shortages drive clusters of spillovers. Pulses of winter flowering of trees in remnant forests appeared to prevent spillover. We developed integrative Bayesian network models based on these phenomena that accurately predicted the presence or absence of clusters of spillovers in each of the 25 years. Our long-term study identifies the mechanistic connections between habitat loss, climate and increased spillover risk. It provides a framework for examining causes of bat virus spillover and for developing ecological countermeasures to prevent pandemics.

Tropical forests are approaching critical temperature thresholds.

The critical temperature beyond which photosynthetic machinery in tropical trees begins to fail averages approximately 46.7 °C (Tcrit)1. However, it remains unclear whether leaf temperatures experienced by tropical vegetation approach this threshold or soon will under climate change. Here we found that pantropical canopy temperatures independently triangulated from individual leaf thermocouples, pyrgeometers and remote sensing (ECOSTRESS) have midday peak temperatures of approximately 34 °C during dry periods, with a long high-temperature tail that can exceed 40 °C. Leaf thermocouple data from multiple sites across the tropics suggest that even within pixels of moderate temperatures, upper canopy leaves exceed Tcrit 0.01% of the time. Furthermore, upper canopy leaf warming experiments (+2, 3 and 4 °C in Brazil, Puerto Rico and Australia, respectively) increased leaf temperatures non-linearly, with peak leaf temperatures exceeding Tcrit 1.3% of the time (11% for more than 43.5 °C, and 0.3% for more than 49.9 °C). Using an empirical model incorporating these dynamics (validated with warming experiment data), we found that tropical forests can withstand up to a 3.9 ± 0.5 °C increase in air temperatures before a potential tipping point in metabolic function, but remaining uncertainty in the plasticity and range of Tcrit in tropical trees and the effect of leaf death on tree death could drastically change this prediction. The 4.0 °C estimate is within the 'worst-case scenario' (representative concentration pathway (RCP) 8.5) of climate change predictions2 for tropical forests and therefore it is still within our power to decide (for example, by not taking the RCP 6.0 or 8.5 route) the fate of these critical realms of carbon, water and biodiversity3,4.

Chlorine activation and enhanced ozone depletion induced by wildfire aerosol.

Remarkable perturbations in the stratospheric abundances of chlorine species and ozone were observed over Southern Hemisphere mid-latitudes following the 2020 Australian wildfires1,2. These changes in atmospheric chemical composition suggest that wildfire aerosols affect stratospheric chlorine and ozone depletion chemistry. Here we propose that wildfire aerosol containing a mixture of oxidized organics and sulfate3-7 increases hydrochloric acid solubility8-11 and associated heterogeneous reaction rates, activating reactive chlorine species and enhancing ozone loss rates at relatively warm stratospheric temperatures. We test our hypothesis by comparing atmospheric observations to model simulations that include the proposed mechanism. Modelled changes in 2020 hydrochloric acid, chlorine nitrate and hypochlorous acid abundances are in good agreement with observations1,2. Our results indicate that wildfire aerosol chemistry, although not accounting for the record duration of the 2020 Antarctic ozone hole, does yield an increase in its area and a 3-5% depletion of southern mid-latitude total column ozone. These findings increase concern2,12,13 that more frequent and intense wildfires could delay ozone recovery in a warming world.

The landscape of genomic structural variation in Indigenous Australians.

Indigenous Australians harbour rich and unique genomic diversity. However, Aboriginal and Torres Strait Islander ancestries are historically under-represented in genomics research and almost completely missing from reference datasets1-3. Addressing this representation gap is critical, both to advance our understanding of global human genomic diversity and as a prerequisite for ensuring equitable outcomes in genomic medicine. Here we apply population-scale whole-genome long-read sequencing4 to profile genomic structural variation across four remote Indigenous communities. We uncover an abundance of large insertion-deletion variants (20-49 bp; n = 136,797), structural variants (50 b-50 kb; n = 159,912) and regions of variable copy number (>50 kb; n = 156). The majority of variants are composed of tandem repeat or interspersed mobile element sequences (up to 90%) and have not been previously annotated (up to 62%). A large fraction of structural variants appear to be exclusive to Indigenous Australians (12% lower-bound estimate) and most of these are found in only a single community, underscoring the need for broad and deep sampling to achieve a comprehensive catalogue of genomic structural variation across the Australian continent. Finally, we explore short tandem repeats throughout the genome to characterize allelic diversity at 50 known disease loci5, uncover hundreds of novel repeat expansion sites within protein-coding genes, and identify unique patterns of diversity and constraint among short tandem repeat sequences. Our study sheds new light on the dimensions and dynamics of genomic structural variation within and beyond Australia.

Indigenous Australian genomes show deep structure and rich novel variation.

The Indigenous peoples of Australia have a rich linguistic and cultural history. How this relates to genetic diversity remains largely unknown because of their limited engagement with genomic studies. Here we analyse the genomes of 159 individuals from four remote Indigenous communities, including people who speak a language (Tiwi) not from the most widespread family (Pama-Nyungan). This large collection of Indigenous Australian genomes was made possible by careful community engagement and consultation. We observe exceptionally strong population structure across Australia, driven by divergence times between communities of 26,000-35,000 years ago and long-term low but stable effective population sizes. This demographic history, including early divergence from Papua New Guinean (47,000 years ago) and Eurasian groups1, has generated the highest proportion of previously undescribed genetic variation seen outside Africa and the most extended homozygosity compared with global samples. A substantial proportion of this variation is not observed in global reference panels or clinical datasets, and variation with predicted functional consequence is more likely to be homozygous than in other populations, with consequent implications for medical genomics2. Our results show that Indigenous Australians are not a single homogeneous genetic group and their genetic relationship with the peoples of New Guinea is not uniform. These patterns imply that the full breadth of Indigenous Australian genetic diversity remains uncharacterized, potentially limiting genomic medicine and equitable healthcare for Indigenous Australians.

Continent-wide declines in shallow reef life over a decade of ocean warming.

Human society is dependent on nature1,2, but whether our ecological foundations are at risk remains unknown in the absence of systematic monitoring of species' populations3. Knowledge of species fluctuations is particularly inadequate in the marine realm4. Here we assess the population trends of 1,057 common shallow reef species from multiple phyla at 1,636 sites around Australia over the past decade. Most populations decreased over this period, including many tropical fishes, temperate invertebrates (particularly echinoderms) and southwestern Australian macroalgae, whereas coral populations remained relatively stable. Population declines typically followed heatwave years, when local water temperatures were more than 0.5 °C above temperatures in 2008. Following heatwaves5,6, species abundances generally tended to decline near warm range edges, and increase near cool range edges. More than 30% of shallow invertebrate species in cool latitudes exhibited high extinction risk, with rapidly declining populations trapped by deep ocean barriers, preventing poleward retreat as temperatures rise. Greater conservation effort is needed to safeguard temperate marine ecosystems, which are disproportionately threatened and include species with deep evolutionary roots. Fundamental among such efforts, and broader societal needs to efficiently adapt to interacting anthropogenic and natural pressures, is greatly expanded monitoring of species' population trends7,8.

Tropical tree mortality has increased with rising atmospheric water stress.

Evidence exists that tree mortality is accelerating in some regions of the tropics1,2, with profound consequences for the future of the tropical carbon sink and the global anthropogenic carbon budget left to limit peak global warming below 2 °C. However, the mechanisms that may be driving such mortality changes and whether particular species are especially vulnerable remain unclear3-8. Here we analyse a 49-year record of tree dynamics from 24 old-growth forest plots encompassing a broad climatic gradient across the Australian moist tropics and find that annual tree mortality risk has, on average, doubled across all plots and species over the last 35 years, indicating a potential halving in life expectancy and carbon residence time. Associated losses in biomass were not offset by gains from growth and recruitment. Plots in less moist local climates presented higher average mortality risk, but local mean climate did not predict the pace of temporal increase in mortality risk. Species varied in the trajectories of their mortality risk, with the highest average risk found nearer to the upper end of the atmospheric vapour pressure deficit niches of species. A long-term increase in vapour pressure deficit was evident across the region, suggesting that thresholds involving atmospheric water stress, driven by global warming, may be a primary cause of increasing tree mortality in moist tropical forests.

Fossil evidence unveils an early Cambrian origin for Bryozoa.

Bryozoans (also known as ectoprocts or moss animals) are aquatic, dominantly sessile, filter-feeding lophophorates that construct an organic or calcareous modular colonial (clonal) exoskeleton1-3. The presence of six major orders of bryozoans with advanced polymorphisms in lower Ordovician rocks strongly suggests a Cambrian origin for the largest and most diverse lophophorate phylum2,4-8. However, a lack of convincing bryozoan fossils from the Cambrian period has hampered resolution of the true origins and character assembly of the earliest members of the group. Here we interpret the millimetric, erect, bilaminate, secondarily phosphatized fossil Protomelission gatehousei9 from the early Cambrian of Australia and South China as a potential stem-group bryozoan. The monomorphic zooid capsules, modular construction, organic composition and simple linear budding growth geometry represent a mixture of organic Gymnolaemata and biomineralized Stenolaemata character traits, with phylogenetic analyses identifying P. gatehousei as a stem-group bryozoan. This aligns the origin of phylum Bryozoa with all other skeletonized phyla in Cambrian Age 3, pushing back its first occurrence by approximately 35 million years. It also reconciles the fossil record with molecular clock estimations of an early Cambrian origination and subsequent Ordovician radiation of Bryozoa following the acquisition of a carbonate skeleton10-13.

Genomic insights into population history and biological adaptation in Oceania.

The Pacific region is of major importance for addressing questions regarding human dispersals, interactions with archaic hominins and natural selection processes1. However, the demographic and adaptive history of Oceanian populations remains largely uncharacterized. Here we report high-coverage genomes of 317 individuals from 20 populations from the Pacific region. We find that the ancestors of Papuan-related ('Near Oceanian') groups underwent a strong bottleneck before the settlement of the region, and separated around 20,000-40,000 years ago. We infer that the East Asian ancestors of Pacific populations may have diverged from Taiwanese Indigenous peoples before the Neolithic expansion, which is thought to have started from Taiwan around 5,000 years ago2-4. Additionally, this dispersal was not followed by an immediate, single admixture event with Near Oceanian populations, but involved recurrent episodes of genetic interactions. Our analyses reveal marked differences in the proportion and nature of Denisovan heritage among Pacific groups, suggesting that independent interbreeding with highly structured archaic populations occurred. Furthermore, whereas introgression of Neanderthal genetic information facilitated the adaptation of modern humans related to multiple phenotypes (for example, metabolism, pigmentation and neuronal development), Denisovan introgression was primarily beneficial for immune-related functions. Finally, we report evidence of selective sweeps and polygenic adaptation associated with pathogen exposure and lipid metabolism in the Pacific region, increasing our understanding of the mechanisms of biological adaptation to island environments.

Vast CO2 release from Australian fires in 2019-2020 constrained by satellite.

Southeast Australia experienced intensive and geographically extensive wildfires during the 2019-2020 summer season1,2. The fires released substantial amounts of carbon dioxide into the atmosphere3. However, existing emission estimates based on fire inventories are uncertain4, and vary by up to a factor of four for this event. Here we constrain emission estimates with the help of satellite observations of carbon monoxide5, an analytical Bayesian inversion6 and observed ratios between emitted carbon dioxide and carbon monoxide7. We estimate emissions of carbon dioxide to be 715 teragrams (range 517-867) from November 2019 to January 2020. This is more than twice the estimate derived by five different fire inventories8-12, and broadly consistent with estimates based on a bottom-up bootstrap analysis of this fire episode13. Although fires occur regularly in the savannas in northern Australia, the recent episodes were extremely large in scale and intensity, burning unusually large areas of eucalyptus forest in the southeast13. The fires were driven partly by climate change14,15, making better-constrained emission estimates particularly important. This is because the build-up of atmospheric carbon dioxide may become increasingly dependent on fire-driven climate-carbon feedbacks, as highlighted by this event16.

Genome of a middle Holocene hunter-gatherer from Wallacea.

Much remains unknown about the population history of early modern humans in southeast Asia, where the archaeological record is sparse and the tropical climate is inimical to the preservation of ancient human DNA1. So far, only two low-coverage pre-Neolithic human genomes have been sequenced from this region. Both are from mainland Hòabìnhian hunter-gatherer sites: Pha Faen in Laos, dated to 7939-7751 calibrated years before present (yr cal BP; present taken as AD 1950), and Gua Cha in Malaysia (4.4-4.2 kyr cal BP)1. Here we report, to our knowledge, the first ancient human genome from Wallacea, the oceanic island zone between the Sunda Shelf (comprising mainland southeast Asia and the continental islands of western Indonesia) and Pleistocene Sahul (Australia-New Guinea). We extracted DNA from the petrous bone of a young female hunter-gatherer buried 7.3-7.2 kyr cal BP at the limestone cave of Leang Panninge2 in South Sulawesi, Indonesia. Genetic analyses show that this pre-Neolithic forager, who is associated with the 'Toalean' technocomplex3,4, shares most genetic drift and morphological similarities with present-day Papuan and Indigenous Australian groups, yet represents a previously unknown divergent human lineage that branched off around the time of the split between these populations approximately 37,000 years ago5. We also describe Denisovan and deep Asian-related ancestries in the Leang Panninge genome, and infer their large-scale displacement from the region today.

Widespread phytoplankton blooms triggered by 2019-2020 Australian wildfires.

Droughts and climate-change-driven warming are leading to more frequent and intense wildfires1-3, arguably contributing to the severe 2019-2020 Australian wildfires4. The environmental and ecological impacts of the fires include loss of habitats and the emission of substantial amounts of atmospheric aerosols5-7. Aerosol emissions from wildfires can lead to the atmospheric transport of macronutrients and bio-essential trace metals such as nitrogen and iron, respectively8-10. It has been suggested that the oceanic deposition of wildfire aerosols can relieve nutrient limitations and, consequently, enhance marine productivity11,12, but direct observations are lacking. Here we use satellite and autonomous biogeochemical Argo float data to evaluate the effect of 2019-2020 Australian wildfire aerosol deposition on phytoplankton productivity. We find anomalously widespread phytoplankton blooms from December 2019 to March 2020 in the Southern Ocean downwind of Australia. Aerosol samples originating from the Australian wildfires contained a high iron content and atmospheric trajectories show that these aerosols were likely to be transported to the bloom regions, suggesting that the blooms resulted from the fertilization of the iron-limited waters of the Southern Ocean. Climate models project more frequent and severe wildfires in many regions1-3. A greater appreciation of the links between wildfires, pyrogenic aerosols13, nutrient cycling and marine photosynthesis could improve our understanding of the contemporary and glacial-interglacial cycling of atmospheric CO2 and the global climate system.

Ancient genomes reveal over two thousand years of dingo population structure.

Dingoes are culturally and ecologically important free-living canids whose ancestors arrived in Australia over 3,000 B.P., likely transported by seafaring people. However, the early history of dingoes in Australia-including the number of founding populations and their routes of introduction-remains uncertain. This uncertainty arises partly from the complex and poorly understood relationship between modern dingoes and New Guinea singing dogs, and suspicions that post-Colonial hybridization has introduced recent domestic dog ancestry into the genomes of many wild dingo populations. In this study, we analyzed genome-wide data from nine ancient dingo specimens ranging in age from 400 to 2,746 y old, predating the introduction of domestic dogs to Australia by European colonists. We uncovered evidence that the continent-wide population structure observed in modern dingo populations had already emerged several thousand years ago. We also detected excess allele sharing between New Guinea singing dogs and ancient dingoes from coastal New South Wales (NSW) compared to ancient dingoes from southern Australia, irrespective of any post-Colonial hybrid ancestry in the genomes of modern individuals. Our results are consistent with several demographic scenarios, including a scenario where the ancestry of dingoes from the east coast of Australia results from at least two waves of migration from source populations with varying affinities to New Guinea singing dogs. We also contribute to the growing body of evidence that modern dingoes derive little genomic ancestry from post-Colonial hybridization with other domestic dog lineages, instead descending primarily from ancient canids introduced to Sahul thousands of years ago.

Boom-bust population dynamics drive rapid genetic change.

Increasing environmental threats and more extreme environmental perturbations place species at risk of population declines, with associated loss of genetic diversity and evolutionary potential. While theory shows that rapid population declines can cause loss of genetic diversity, populations in some environments, like Australia's arid zone, are repeatedly subject to major population fluctuations yet persist and appear able to maintain genetic diversity. Here, we use repeated population sampling over 13 y and genotype-by-sequencing of 1903 individuals to investigate the genetic consequences of repeated population fluctuations in two small mammals in the Australian arid zone. The sandy inland mouse (Pseudomys hermannsburgensis) experiences marked boom-bust population dynamics in response to the highly variable desert environment. We show that heterozygosity levels declined, and population differentiation (FST) increased, during bust periods when populations became small and isolated, but that heterozygosity was rapidly restored during episodic population booms. In contrast, the lesser hairy-footed dunnart (Sminthopsis youngsoni), a desert marsupial that maintains relatively stable population sizes, showed no linear declines in heterozygosity. These results reveal two contrasting ways in which genetic diversity is maintained in highly variable environments. In one species, diversity is conserved through the maintenance of stable population sizes across time. In the other species, diversity is conserved through rapid genetic mixing during population booms that restores heterozygosity lost during population busts.