Search

Application of whole genome sequencing (WGS) has allowed monitoring of the emergence of variants of concern (VOC) of severe acute respiratory syndrome-related coronavirus 2 (SARS-CoV-2) globally. Genomic investigation of emerging variants and surveillance of clinical progress has reduced the public health impact of infection during the COVID-19 pandemic. These steps required developing and implementing a proficiency testing program (PTP), as WGS has been incorporated into routine reference laboratory practice. In this study, we describe how the PTP evaluated the capacity and capability of one New Zealand and 14 Australian public health laboratories to perform WGS of SARS-CoV-2 in 2022. The participants' performances in characterising a specimen panel of known SARS-CoV-2 isolates in the PTP were assessed based on: (1) genome coverage, (2) Pango lineage, and (3) sequence quality, with the choice of assessment metrics refined based on a previously reported assessment conducted in 2021. The participants' performances in 2021 and 2022 were also compared after reassessing the 2021 results using the more stringent metrics adopted in 2022. We found that more participants would have failed the 2021 assessment for all survey samples and a significantly higher fail rate per sample in 2021 compared to 2022. This study highlights the importance of choosing appropriate performance metrics to reflect better the laboratories' capacity to perform SARS-CoV-2 WGS, as was done in the 2022 PTP. It also displays the need for a PTP for WGS of SARS-CoV-2 to be available to public health laboratories ongoing, with continuous refinements in the design and provision of the PTP to account for the dynamic nature of the COVID-19 pandemic as SARS-CoV-2 continues to evolve., (Copyright © 2024 Royal College of Pathologists of Australasia. Published by Elsevier B.V. All rights reserved.)

Human cytomegalovirus (CMV) infection is the leading non-genetic aetiology of congenital malformation in developed countries, causing significant fetal neurological injury. This study investigated potential CMV pathogenetic mechanisms of fetal neural malformation using in vitro human cerebral organoids. Cerebral organoids were permissive to CMV replication, and infection dysregulated cellular pluripotency and differentiation pathways. Aberrant expression of dual-specificity tyrosine phosphorylation-regulated kinases (DYRK), sonic hedgehog (SHH), pluripotency, neurodegeneration, axon guidance, hippo signalling and dopaminergic synapse pathways were observed in CMV-infected organoids using immunofluorescence and RNA-sequencing. Infection with CMV resulted in dysregulation of 236 Autism Spectrum Disorder (ASD)-related genes (p = 1.57E-05) and pathways. This notable observation suggests potential links between congenital CMV infection and ASD. Using DisGeNET databases, 103 diseases related to neural malformation or mental disorders were enriched in CMV-infected organoids. Cytomegalovirus infection-related dysregulation of key cerebral cellular pathways potentially provides important, modifiable pathogenetic mechanisms for congenital CMV-induced neural malformation and ASD., (© 2024. The Author(s).)

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., (© 2024 The Authors. Insect Science published by John Wiley & Sons Australia, Ltd on behalf of Institute of Zoology, Chinese Academy of Sciences.)

Objectives: The 2022 seasonal respiratory syncytial virus (RSV) epidemic in Sydney, Australia saw an unprecedented number of RSV detections. We aimed to characterize genomic and immunologic factors associated with the surge in RSV cases., Methods: Whole genome sequences of RSV were generated from 264 RSV-infected infants and linked to case-matched clinical data from the 2022 southern hemisphere RSV season. We then performed an immunologic analysis of baseline RSV-specific humoral immunity in women of childbearing age before and throughout the coronavirus disease 2019 pandemic., Results: Clinical analysis revealed a high burden of disease across patients of all health backgrounds. More than one-half of RSV-related health care visits by infants resulted in hospitalization, and one-quarter required high-flow respiratory support or a higher level of care. Viral phylogenetic analyses revealed that 2022 Sydney RSV sequences were closely related to viruses that had been circulating globally since 2017, including those detected in recent US outbreaks. Nonsynonymous mutations within the palivizumab and nirsevimab binding sites were detected at low frequencies. There was no difference in baseline RSV-neutralizing antibody titers between 2020 and 2022., Conclusions: Collectively, these findings suggest that neither the emergence of a novel RSV genotype nor hypothesized immune debt was associated with the surge of RSV cases and hospitalizations in 2022. Continued genomic and immunologic surveillance is required to further understand the factors driving outbreaks of RSV globally, and to inform guidelines for the rollout and ongoing use of recently developed immunotherapeutics and vaccines.

The COVID-19 pandemic has necessitated the rapid development and implementation of whole-genome sequencing (WGS) and bioinformatic methods for managing the pandemic. However, variability in methods and capabilities between laboratories has posed challenges in ensuring data accuracy. A national working group comprising 18 laboratory scientists and bioinformaticians from Australia and New Zealand was formed to improve data concordance across public health laboratories (PHLs). One effort, presented in this study, sought to understand the impact of the methodology on consensus genome concordance and interpretation. SARS-CoV-2 WGS proficiency testing programme (PTP) data were retrospectively obtained from the 2021 Royal College of Pathologists of Australasia Quality Assurance Programmes (RCPAQAP), which included 11 participating Australian laboratories. The submitted consensus genomes and reads from eight contrived specimens were investigated, focusing on discordant sequence data and findings were presented to the working group to inform best practices. Despite using a variety of laboratory and bioinformatic methods for SARS-CoV-2 WGS, participants largely produced concordant genomes. Two participants returned five discordant sites in a high-Cτ replicate, which could be resolved with reasonable bioinformatic quality thresholds. We noted ten discrepancies in genome assessment that arose from nucleotide heterogeneity at three different sites in three cell-culture-derived control specimens. While these sites were ultimately accurate after considering the participants' bioinformatic parameters, it presented an interesting challenge for developing standards to account for intrahost single nucleotide variation (iSNV). Observed differences had little to no impact on key surveillance metrics, lineage assignment and phylogenetic clustering, while genome coverage <90 % affected both. We recommend PHLs bioinformatically generate two consensus genomes with and without ambiguity thresholds for quality control and downstream analysis, respectively, and adhere to a minimum 90 % genome coverage threshold for inclusion in surveillance interpretations. We also suggest additional PTP assessment criteria, including primer efficiency, detection of iSNVs and minimum genome coverage of 90 %. This study underscores the importance of multidisciplinary national working groups in informing guidelines in real time for bioinformatic quality acceptance criteria. It demonstrates the potential for enhancing public health responses through improved data concordance and quality control in SARS-CoV-2 genomic analysis during pandemic surveillance.

Background: The Omicron era of the COVID-19 pandemic commenced at the beginning of 2022 and whilst it started with primarily BA.1, it was latter dominated by BA.2 and the related sub-lineage BA.5. Following resolution of the global BA.5 wave, a diverse grouping of Omicron sub-lineages emerged derived from BA.2, BA.5 and recombinants thereof. Whilst emerging from distinct lineages, all shared similar changes in the Spike glycoprotein affording them an outgrowth advantage through evasion of neutralising antibodies., Methods: Over the course of 2022, we monitored the potency and breadth of antibody neutralization responses to many emerging variants in the Australian community at three levels: (i) we tracked over 420,000 U.S. plasma donors over time through various vaccine booster roll outs and Omicron waves using sequentially collected IgG pools; (ii) we mapped the antibody response in individuals using blood from stringently curated vaccine and convalescent cohorts. (iii) finally we determine the in vitro efficacy of clinically approved therapies Evusheld and Sotrovimab., Findings: In pooled IgG samples, we observed the maturation of neutralization breadth to Omicron variants over time through continuing vaccine and infection waves. Importantly, in many cases, we observed increased antibody breadth to variants that were yet to be in circulation. Determination of viral neutralization at the cohort level supported equivalent coverage across prior and emerging variants with isolates BQ.1.1, XBB.1, BR.2.1 and XBF the most evasive. Further, these emerging variants were resistant to Evusheld, whilst increasing neutralization resistance to Sotrovimab was restricted to BQ.1.1 and XBF. We conclude at this current point in time that dominant variants can evade antibodies at levels equivalent to their most evasive lineage counterparts but sustain an entry phenotype that continues to promote an additional outgrowth advantage. In Australia, BR.2.1 and XBF share this phenotype and, in contrast to global variants, are uniquely dominant in this region in the later months of 2022., Interpretation: Whilst the appearance of a diverse range of omicron lineages has led to primary or partial resistance to clinically approved monoclonal antibodies, the maturation of the antibody response across both cohorts and a large donor pools importantly observes increasing breadth in the antibody neutralisation responses over time with a trajectory that covers both current and known emerging variants., Funding: This work was primarily supported by Australian Medical Foundation research grants MRF2005760 (SGT, GM & WDR), Medical Research Future Fund Antiviral Development Call grant (WDR), the New South Wales Health COVID-19 Research Grants Round 2 (SGT & FB) and the NSW Vaccine Infection and Immunology Collaborative (VIIM) (ALC). Variant modeling was supported by funding from SciLifeLab's Pandemic Laboratory Preparedness program to B.M. (VC-2022-0028) and by the European Union's Horizon 2020 research and innovation programme under grant agreement no. 101003653 (CoroNAb) to B.M., Competing Interests: Declaration of interests Funding bodies did not contribute to study design, data collection, data analysis or writing of the manuscript. Study design, data collection, data analysis and data interpretation was performed by the listed authors. Writing and review was performed by the listed authors. A.L.C. Participated on the AstraZeneca COVID-19 Advisory Board (2021), Seqirus COVID Advisory Board (2020) & CSL Seqirus APAC Advisory Council (2022). N.R., T.B., S.A.C., S.M., and T.H. are employees and shareholders of CSL LTD, a manufacturer of the polyclonal immunoglobulin preparations used as a reagent within the study., (Copyright © 2023 The Author(s). Published by Elsevier B.V. All rights reserved.)

Australia experienced widespread COVID-19 outbreaks from infection with the SARS-CoV-2 Delta variant between June 2021 and February 2022. A 17-nucleotide frameshift-inducing deletion in ORF7a rapidly became represented at the consensus level (Delta-ORF7a Δ17del ) in most Australian outbreak cases. Studies from early in the COVID-19 pandemic suggest that frameshift-inducing deletions in ORF7a do not persist for long in the population; therefore, Delta-ORF7a Δ17del genomes should have disappeared early in the Australian outbreak. In this study, we conducted a retrospective analysis of global Delta genomes to characterise the dynamics of Delta-ORF7a Δ17del over time, determined the frequency of all ORF7a deletions worldwide, and compared global trends with those of the Australian Delta outbreak. We downloaded all GISAID clade GK Delta genomes and scanned them for deletions in ORF7a. For each deletion we identified, we characterised its frequency, the number of countries it was found in, and how long it persisted. Of the 4,018,216 Delta genomes identified globally, 134,751 (~3.35%) possessed an ORF7a deletion, and ORF7a Δ17del was the most common. ORF7a Δ17del was the sole deletion in 28,014 genomes, of which 27,912 (~99.6%) originated from the Australian outbreak. During the outbreak, ~87% of genomes were Delta-ORF7a Δ17del , and genomes with this deletion were sampled until the outbreak's end. These data demonstrate that, contrary to suggestions early in the COVID-19 pandemic, genomes with frameshifting deletions in ORF7a can persist over long time periods. We suggest that the proliferation of Delta-ORF7a Δ17del genomes was likely a chance founder effect. Nonetheless, the frequency of ORF7a deletions in SARS-CoV-2 genomes worldwide suggests they might have some benefit for virus transmission.

Background: Genetically distinct viral variants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) have been recorded since January 2020. The introduction of global vaccine programs has contributed to lower COVID-19 hospitalisation and mortality rates, particularly in developed countries. In late 2021, Omicron BA.1 emerged, with substantially altered genetic differences and clinical effects from other variants of concern. Shortly after dominating global spread in early 2022, BA.1 was supplanted by the genetically distinct Omicron lineage BA.2. A sub-lineage of BA.2, designated BA.5, presently has an outgrowth advantage over BA.2 and other BA.2 sub-lineages. Here we study the neutralisation of Omicron BA.1, BA.2 and BA.5 and pre-Omicron variants using a range of vaccine and convalescent sera and therapeutic monoclonal antibodies using a live virus neutralisation assay. Using primary nasopharyngeal swabs, we also tested the relative fitness of BA.5 compared to pre-Omicron and Omicron viral lineages in their ability to use the ACE2-TMPRSS2 pathway., Methods: Using low passage clinical isolates of Clade A.2.2, Beta, Delta, BA.1, BA.2 and BA.5, we determined humoral neutralisation in vitro in vaccinated and convalescent cohorts, using concentrated human IgG pooled from thousands of plasma donors, and licensed monoclonal antibody therapies. We then determined infectivity to particle ratios in primary nasopharyngeal samples and expanded low passage isolates in a genetically engineered ACE2/TMPRSS2 cell line in the presence and absence of the TMPRSS2 inhibitor Nafamostat., Findings: Peak responses to 3 doses of BNT162b2 vaccine were associated with a 9-fold reduction in neutralisation for Omicron lineages BA.1, BA.2 and BA.5. Concentrated pooled human IgG from convalescent and vaccinated donors and BNT162b2 vaccination with BA.1 breakthrough infections were associated with greater breadth of neutralisation, although the potency was still reduced 7-fold across all Omicron lineages. Testing of clinical grade antibodies revealed a 14.3-fold reduction using Evusheld and 16.8-fold reduction using Sotrovimab for the BA.5. Whilst the infectivity of BA.1 and BA.2 was attenuated in ACE2/TMPRSS2 entry, BA.5 was observed to be equivalent to that of an early 2020 circulating clade and had greater sensitivity to the TMPRSS2 inhibitor Nafamostat., Interpretation: Observations support all Omicron variants to significantly escape neutralising antibodies across a range of vaccination and/or convalescent responses. Potency of therapeutic monoclonal antibodies is also reduced and differs across Omicron lineages. The key difference of BA.5 from other Omicron sub-variants is the reversion in tropism back to using the well-known ACE2-TMPRSS2 pathway, utilised efficiently by pre-Omicron lineages. Monitoring if these changes influence transmission and/or disease severity will be key for ongoing tracking and management of Omicron waves globally., Funding: This work was primarily supported by Australian Medical Foundation research grants MRF2005760 (ST, GM & WDR), MRF2001684 (ADK and ST) and Medical Research Future Fund Antiviral Development Call grant (WDR), Medical Research Future Fund COVID-19 grant (MRFF2001684, ADK & SGT) and the New South Wales Health COVID-19 Research Grants Round 2 (SGT)., Competing Interests: Declaration of interests A.O.S contributed to this manuscript in his capacity as adjunct associate lecturer (from January 2022) at the University of New South Wales. Since January 2022, A.O.S is an employee of GlaxoSmithKline Australia (medical science liaison for COVID-19)., (Copyright © 2022 The Author(s). Published by Elsevier B.V. All rights reserved.)

Subterranean habitats are generally very stable environments, and as such evolutionary transitions of organisms from surface to subterranean lifestyles may cause considerable shifts in physiology, particularly with respect to thermal tolerance. In this study we compared responses to heat shock at the molecular level in a geographically widespread, surface-dwelling water beetle to a congeneric subterranean species restricted to a single aquifer (Dytiscidae: Hydroporinae). The obligate subterranean beetle Paroster macrosturtensis is known to have a lower thermal tolerance compared to surface lineages (CT max 38 °C cf. 42-46 °C), but the genetic basis of this physiological difference has not been characterized. We experimentally manipulated the thermal environment of 24 individuals to demonstrate that both species can mount a heat shock response at high temperatures (35 °C), as determined by comparative transcriptomics. However, genes involved in these responses differ between species and a far greater number were differentially expressed in the surface taxon, suggesting it can mount a more robust heat shock response; these data may underpin its higher thermal tolerance compared to subterranean relatives. In contrast, the subterranean species examined not only differentially expressed fewer genes in response to increasing temperatures, but also in the presence of the experimental setup employed here alone. Our results suggest P. macrosturtensis may be comparatively poorly equipped to respond to both thermally induced stress and environmental disturbances more broadly. The molecular findings presented here have conservation implications for P. macrosturtensis and contribute to a growing narrative concerning weakened thermal tolerances in obligate subterranean organisms at the molecular level., (© 2022. The Author(s).)

Cytomegalovirus infection during antiviral prophylaxis occurs in transplant recipients despite individualized regimens based on renal function. Fifty kidney transplant recipients were assessed between 2016 and 2019 for valganciclovir dosing, ganciclovir exposure, cytomegalovirus infection, and genotypic resistance markers during the first year posttransplant. Ganciclovir plasma concentrations were measured using mass spectrometry. Population pharmacokinetics was used to determine individual ganciclovir exposure and to evaluate the ability of manufacturer dosing guidelines to meet therapeutic target daily area under the curve (AUC 24 ) of 40 to 50 μg·h/mL. Full-length UL54 and UL97 were assessed using high-throughput sequencing in cytomegalovirus DNA-positive patient specimens. Valganciclovir doses administered to recipients with creatinine clearance of <40 mL/min were higher than specified by guidelines, and they were lower for recipients with creatinine clearance of ≥40 mL/min. The mean ganciclovir AUC 24 was 33 ± 13 μg·h/mL, and 82% of subjects did not attain the therapeutic target. Pharmacokinetic simulations showed that the guidelines similarly could not attain the therapeutic target in 79% of individuals. Cytomegalovirus breakthrough occurred in 6% (3/50) of recipients, while 12% (6/50) developed late-onset infection. The mean AUC 24 s of recipients with ( n  = 3) and without ( n  = 47) infection were not significantly different ( P =  0.528). However, one recipient with an AUC 24 of 20 μg·h/mL acquired two UL97 ganciclovir resistance mutations. Current prophylaxis guidelines resulted in subtherapeutic ganciclovir exposure in several study recipients, including the emergence of resistance genotypes. IMPORTANCE This study examined the pharmacokinetics and viral genomic data from a prospective cohort of kidney transplant recipients undergoing valganciclovir prophylaxis for cytomegalovirus (CMV) prevention. We showed for the first time using high-throughput sequencing the detection of ganciclovir resistance mutations in breakthrough CMV infection during subtherapeutic plasma ganciclovir as indicated by the pharmacokinetic parameter daily area under the curve (AUC 24 ). In addition, we found that current valganciclovir dosing guidelines for CMV prophylaxis are predicted to attain therapeutic targets in only 21% of recipients, which is consistent with previous pharmacokinetic studies. The novel findings of resistance mutations during subtherapeutic ganciclovir exposure presented here can inform future studies investigating the dynamics of drug selection pressure and the emergence of resistance mutations in vivo .

Genetically distinct variants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) have emerged since the start of the COVID-19 pandemic. Over this period, we developed a rapid platform (R-20) for viral isolation and characterization using primary remnant diagnostic swabs. This, combined with quarantine testing and genomics surveillance, enabled the rapid isolation and characterization of all major SARS-CoV-2 variants circulating in Australia in 2021. Our platform facilitated viral variant isolation, rapid resolution of variant fitness using nasopharyngeal swabs and ranking of evasion of neutralizing antibodies. In late 2021, variant of concern Omicron (B1.1.529) emerged. Using our platform, we detected and characterized SARS-CoV-2 VOC Omicron. We show that Omicron effectively evades neutralization antibodies and has a different entry route that is TMPRSS2-independent. Our low-cost platform is available to all and can detect all variants of SARS-CoV-2 studied so far, with the main limitation being that our platform still requires appropriate biocontainment., (© 2022. The Author(s).)

The repeated evolution of the same traits in distantly related groups (convergent evolution) raises a key question in evolutionary biology: do the same genes underpin convergent phenotypes? Here, we explore one such trait, viviparity (live birth), which, qualitative studies suggest, may indeed have evolved via genetic convergence. There are >150 independent origins of live birth in vertebrates, providing a uniquely powerful system to test the mechanisms underpinning convergence in morphology, physiology, and/or gene recruitment during pregnancy. We compared transcriptomic data from eight vertebrates (lizards, mammals, sharks) that gestate embryos within the uterus. Since many previous studies detected qualitative similarities in gene use during independent origins of pregnancy, we expected to find significant overlap in gene use in viviparous taxa. However, we found no more overlap in uterine gene expression associated with viviparity than we would expect by chance alone. Each viviparous lineage exhibits the same core set of uterine physiological functions. Yet, contrary to prevailing assumptions about this trait, we find that none of the same genes are differentially expressed in all viviparous lineages, or even in all viviparous amniote lineages. Therefore, across distantly related vertebrates, different genes have been recruited to support the morphological and physiological changes required for successful pregnancy. We conclude that redundancies in gene function have enabled the repeated evolution of viviparity through recruitment of different genes from genomic "toolboxes", which are uniquely constrained by the ancestries of each lineage., (© The Author(s) 2022. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

Whole-genome sequencing of viral isolates is critical for informing transmission patterns and for the ongoing evolution of pathogens, especially during a pandemic. However, when genomes have low variability in the early stages of a pandemic, the impact of technical and/or sequencing errors increases. We quantitatively assessed inter-laboratory differences in consensus genome assemblies of 72 matched SARS-CoV-2-positive specimens sequenced at different laboratories in Sydney, Australia. Raw sequence data were assembled using two different bioinformatics pipelines in parallel, and resulting consensus genomes were compared to detect laboratory-specific differences. Matched genome sequences were predominantly concordant, with a median pairwise identity of 99.997%. Identified differences were predominantly driven by ambiguous site content. Ignoring these produced differences in only 2.3% (5/216) of pairwise comparisons, each differing by a single nucleotide. Matched samples were assigned the same Pango lineage in 98.2% (212/216) of pairwise comparisons, and were mostly assigned to the same phylogenetic clade. However, epidemiological inference based only on single nucleotide variant distances may lead to significant differences in the number of defined clusters if variant allele frequency thresholds for consensus genome generation differ between laboratories. These results underscore the need for a unified, best-practices approach to bioinformatics between laboratories working on a common outbreak problem.

We designed and tested a protocol for measuring the performance of individuals in small-sided soccer games. We tested our protocol on three different groups of youth players from elite Brazilian football academies. Players in each group played a series of 3v3 games, in which individuals were randomly assigned into new teams and against new opponents for each game. We calculated each individual's average individual goals scored, goals scored by teammates, goals conceded, and net team goals per game. Our protocol was consistent across days and repeatable across groups, with intraclass correlation coefficients (ICCs) of 0.57-0.69 for average net goals per game across testing days. Players could achieve high success by scoring goals or ensuring their team concede few goals. We also calculated the first and second dimension of a principal component analysis based on each player's number of goals scored, goals scored by teammates, and number of goals conceded per game. Players that were overall high performers had higher PC 1 scores, while PC 2 scores represented the type of contribution made by a player to overall performance. Positive PC 2 values were indicative of high number of individual goals while negative values were associated with more goals from teammates and fewer conceded goals. Our design allows coaches and scouts to easily collect a robust metric of individual performance using randomly designed, small-sided games. We also provide simulations that allow one to apply our methodology for individual talent identification to other small-sided games in any team sport., (© 2021 John Wiley & Sons A/S. Published by John Wiley & Sons Ltd.)

Our understanding of the evolution of complex biological traits is greatly advanced by examining taxa with intermediate phenotypes. The transition from oviparity (egg-laying) to viviparity (live-bearing) has occurred independently in many animal lineages, but there are few phenotypic intermediates. The lizard Saiphos equalis exhibits bimodal reproduction, with some viviparous populations, and other oviparous populations with long egg-retention, a rare trait where most of embryonic development occurs inside the mother prior to late ovipositioning. We posit that oviparous S. equalis represent an intermediate form between "true" oviparity and viviparity. We used transcriptomics to compare uterine gene expression in these two phenotypes, and provide a molecular model for the genetic control and evolution of reproductive mode. Many genes are differentially expressed throughout the reproductive cycle of both phenotypes, which have clearly different gene expression profiles overall. The differentially expressed genes within oviparous and viviparous individuals have broadly similar biological functions putatively important for sustaining embryos, including uterine remodelling, respiratory gas and water exchange, and immune regulation. These functional similarities indicate either that long egg-retention is an exaptation for viviparity, or might reflect parallel evolution of similar gravidity-related changes in gene expression in long egg-retention oviparity. In contrast, gene expression changes across the reproductive cycle of long egg-retaining oviparous S. equalis are dramatically different from those of "true" oviparous skinks (such as Lampropholis guichenoti), supporting our assertion that oviparous S. equalis exhibit an intermediate phenotype between "true" oviparity and viviparity., (© 2020 John Wiley & Sons Ltd.)

The Kinship Theory of Genomic Imprinting (KTGI) posits that, in species where females mate with multiple males, there is selection for a male to enhance the reproductive success of his offspring at the expense of other males and his mating partner. Reciprocal crosses between honey bee subspecies show parent-of-origin effects for reproductive traits, suggesting that males modify the expression of genes related to female function in their female offspring. This effect is likely to be greater in the Cape honey bee (Apis mellifera capensis), because a male's daughters have the unique ability to produce female offspring that can develop into reproductive workers or the next queen without mating. We generated reciprocal crosses between Capensis and another subspecies and used RNA-seq to identify transcripts that are over- or underexpressed in the embryos, depending on the parental origin of the gene. As predicted, 21 genes showed expression bias towards the Capensis father's allele in colonies with a Capensis father, with no such bias in the reciprocal cross. A further six genes showed a consistent bias towards expression of the father's allele across all eight colonies examined, regardless of the direction of the cross. Consistent with predictions of the KTGI, six of the 21 genes are associated with female reproduction. No gene consistently showed overexpression of the maternal allele., (© 2020 John Wiley & Sons Ltd.)

Australia is an excellent setting to explore relationships between climate change and diversification dynamics. Aridification since the Eocene has resulted in spectacular radiations within one or more Australian biomes. Acacia is the largest plant genus on the Australian continent, with around 1000 species, and is present in all biomes. We investigated the macroevolutionary dynamics of Acacia within climate space. We analysed phylogenetic and climatic data for 503 Acacia species to estimate a time-calibrated phylogeny and central climatic tendencies for BioClim layers from 132 000 herbarium specimens. Diversification rate heterogeneity and rates of climate space exploration were tested. We inferred two diversification rate increases, both associated with significantly higher rates of climate space exploration. Observed spikes in climate disparity within the Pleistocene correspond with onset of Pleistocene glacial-interglacial cycling. Positive time dependency in environmental disparity applies in the basal grade of Acacia, though climate space exploration rates were lower. Incongruence between rates of climate space exploration and disparity suggests different Acacia lineages have experienced different macroevolutionary processes. The second diversification rate increase is associated with a south-east Australian mesic lineage, suggesting adaptations to progressively aridifying environments and ability to transition into mesic environments contributed to Acacia's dominance across Australia., (© 2019 National Herbarium of New South Wales New Phytologist © 2019 New Phytologist Trust.)

Evolution leaves heterogeneous patterns of nucleotide variation across the genome, with different loci subject to varying degrees of mutation, selection, and drift. In phylogenetics, the potential impacts of partitioning sequence data for the assignment of substitution models are well appreciated. In contrast, the treatment of branch lengths has received far less attention. In this study, we examined the effects of linking and unlinking branch-length parameters across loci or subsets of loci. By analyzing a range of empirical data sets, we find consistent support for a model in which branch lengths are proportionate between subsets of loci: gene trees share the same pattern of branch lengths, but form subsets that vary in their overall tree lengths. These models had substantially better statistical support than models that assume identical branch lengths across gene trees, or those in which genes form subsets with distinct branch-length patterns. We show using simulations and empirical data that the complexity of the branch-length model with the highest support depends on the length of the sequence alignment and on the numbers of taxa and loci in the data set. Our findings suggest that models in which branch lengths are proportionate between subsets have the highest statistical support under the conditions that are most commonly seen in practice. The results of our study have implications for model selection, computational efficiency, and experimental design in phylogenomics., (© The Author(s) 2019. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Angiosperms (flowering plants) are the most diverse and species-rich group of plants. The vast majority (∼99.95%) of angiosperms form a clade called Mesangiospermae, which is subdivided into five major groups: eudicots, monocots, magnoliids, Chloranthales, and Ceratophyllales. The relationships among these Mesangiospermae groups have been the subject of long debate. In this study, we assembled a phylogenomic dataset of 1594 genes from 151 angiosperm taxa, including representatives of all five lineages, to investigate the phylogeny of major angiosperm lineages under both coalescent- and concatenation-based methods. We dissected the phylogenetic signal and found that more than half of the genes lack phylogenetic information for the backbone of angiosperm phylogeny. We further removed the genes with weak phylogenetic signal and showed that eudicots, Ceratophyllales, and Chloranthales form a clade, with magnoliids and monocots being the next successive sister lineages. Similar frequencies of gene tree conflict are suggestive of incomplete lineage sorting along the backbone of the angiosperm phylogeny. Our analyses suggest that a fully bifurcating species tree may not be the best way to represent the early radiation of angiosperms. Meanwhile, we inferred that the crown-group angiosperms originated approximately between 255.1 and 222.2 million years ago, and Mesangiospermae diversified into the five extant groups in a short time span (∼27 million years) at the Early to Late Jurassic., (© 2020 The Author(s).)

Establishing an accurate evolutionary timescale for green plants (Viridiplantae) is essential to understanding their interaction and coevolution with the Earth's climate and the many organisms that rely on green plants. Despite being the focus of numerous studies, the timing of the origin of green plants and the divergence of major clades within this group remain highly controversial. Here, we infer the evolutionary timescale of green plants by analyzing 81 protein-coding genes from 99 chloroplast genomes, using a core set of 21 fossil calibrations. We test the sensitivity of our divergence-time estimates to various components of Bayesian molecular dating, including the tree topology, clock models, clock-partitioning schemes, rate priors, and fossil calibrations. We find that the choice of clock model affects date estimation and that the independent-rates model provides a better fit to the data than the autocorrelated-rates model. Varying the rate prior and tree topology had little impact on age estimates, with far greater differences observed among calibration choices and clock-partitioning schemes. Our analyses yield date estimates ranging from the Paleoproterozoic to Mesoproterozoic for crown-group green plants, and from the Ediacaran to Middle Ordovician for crown-group land plants. We present divergence-time estimates of the major groups of green plants that take into account various sources of uncertainty. Our proposed timeline lays the foundation for further investigations into how green plants shaped the global climate and ecosystems, and how embryophytes became dominant in terrestrial environments., (© The Author(s) 2019. Published by Oxford University Press, on behalf of the Society of Systematic Biologists. All rights reserved. For permissions, please email: journals.permissions@oup.com.)

Data sets comprising small numbers of genetic markers are not always able to resolve phylogenetic relationships. This has frequently been the case in molecular systematic studies of plants, with many analyses being based on sequence data from only two or three chloroplast genes. An example of this comes from the riceflowers Pimelea Banks & Sol. ex Gaertn. (Thymelaeaceae), a large genus of flowering plants predominantly distributed in Australia. Despite the considerable morphological variation in the genus, low sequence divergence in chloroplast markers has led to the phylogeny of Pimelea remaining largely uncertain. In this study, we resolve the backbone of the phylogeny of Pimelea in comprehensive Bayesian and maximum-likelihood analyses of plastome sequences from 41 taxa. However, some relationships received only moderate to poor support, and the Pimelea clade contained extremely short internal branches. By using topology-clustering analyses, we demonstrate that conflicting phylogenetic signals can be found across the trees estimated from individual chloroplast protein-coding genes. A relaxed-clock dating analysis reveals that Pimelea arose in the mid-Miocene, with most divergences within the genus occurring during a subsequent rapid diversification. Our new phylogenetic estimate offers better resolution and is more strongly supported than previous estimates, providing a platform for future taxonomic revisions of both Pimelea and the broader subfamily. Our study has demonstrated the substantial improvements in phylogenetic resolution that can be achieved using plastome-scale data sets in plant molecular systematics., (Copyright © 2018 Elsevier Inc. All rights reserved.)

Evolutionary timescales can be inferred from molecular sequence data using a Bayesian phylogenetic approach. In these methods, the molecular clock is often calibrated using fossil data. The uncertainty in these fossil calibrations is important because it determines the limiting posterior distribution for divergence-time estimates as the sequence length tends to infinity. Here, we investigate how the accuracy and precision of Bayesian divergence-time estimates improve with the increased clock-partitioning of genome-scale data into clock-subsets. We focus on a data set comprising plastome-scale sequences of 52 angiosperm taxa. There was little difference among the Bayesian date estimates whether we chose clock-subsets based on patterns of among-lineage rate heterogeneity or relative rates across genes, or by random assignment. Increasing the degree of clock-partitioning usually led to an improvement in the precision of divergence-time estimates, but this increase was asymptotic to a limit presumably imposed by fossil calibrations. Our clock-partitioning approaches yielded highly precise age estimates for several key nodes in the angiosperm phylogeny. For example, when partitioning the data into 20 clock-subsets based on patterns of among-lineage rate heterogeneity, we inferred crown angiosperms to have arisen 198-178 Ma. This demonstrates that judicious clock-partitioning can improve the precision of molecular dating based on phylogenomic data, but the meaning of this increased precision should be considered critically., (© The Author 2017. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Recent advances in molecular phylogenetics and a series of important palaeobotanical discoveries have revolutionized our understanding of angiosperm diversification. Yet, the origin and early evolution of their most characteristic feature, the flower, remains poorly understood. In particular, the structure of the ancestral flower of all living angiosperms is still uncertain. Here we report model-based reconstructions for ancestral flowers at the deepest nodes in the phylogeny of angiosperms, using the largest data set of floral traits ever assembled. We reconstruct the ancestral angiosperm flower as bisexual and radially symmetric, with more than two whorls of three separate perianth organs each (undifferentiated tepals), more than two whorls of three separate stamens each, and more than five spirally arranged separate carpels. Although uncertainty remains for some of the characters, our reconstruction allows us to propose a new plausible scenario for the early diversification of flowers, leading to new testable hypotheses for future research on angiosperms.

The evolutionary timescale of angiosperms has long been a key question in biology. Molecular estimates of this timescale have shown considerable variation, being influenced by differences in taxon sampling, gene sampling, fossil calibrations, evolutionary models, and choices of priors. Here, we analyze a data set comprising 76 protein-coding genes from the chloroplast genomes of 195 taxa spanning 86 families, including novel genome sequences for 11 taxa, to evaluate the impact of models, priors, and gene sampling on Bayesian estimates of the angiosperm evolutionary timescale. Using a Bayesian relaxed molecular-clock method, with a core set of 35 minimum and two maximum fossil constraints, we estimated that crown angiosperms arose 221 (251-192) Ma during the Triassic. Based on a range of additional sensitivity and subsampling analyses, we found that our date estimates were generally robust to large changes in the parameters of the birth-death tree prior and of the model of rate variation across branches. We found an exception to this when we implemented fossil calibrations in the form of highly informative gamma priors rather than as uniform priors on node ages. Under all other calibration schemes, including trials of seven maximum age constraints, we consistently found that the earliest divergences of angiosperm clades substantially predate the oldest fossils that can be assigned unequivocally to their crown group. Overall, our results and experiments with genome-scale data suggest that reliable estimates of the angiosperm crown age will require increased taxon sampling, significant methodological changes, and new information from the fossil record. [Angiospermae, chloroplast, genome, molecular dating, Triassic.]., (© The Author(s) 2016. Published by Oxford University Press, on behalf of the Society of Systematic Biologists. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)