Search

Delivering large-scale routine pathogen genomics surveillance for public health is of considerable interest, although translational research models that promote national-level implementation are not well defined. We describe the development and deployment of the Australian Pathogen Genomics Program (AusPathoGen), a comprehensive national partnership between academia, public health laboratories, and public health agencies that commenced in January, 2021. Successfully establishing and delivering a national programme requires inclusive and transparent collaboration between stakeholders, defined and clear focus on public health priorities, and support for strengthening national genomics capacity. Major enablers for delivering such a programme include technical solutions for data integration and analysis, such as the genomics surveillance platform AusTrakka, standard bioinformatic analysis methods, and national ethics and data sharing agreements that promote nationally integrated surveillance systems. Training of public health officials to interpret and act on genomic data is crucial, and evaluation and cost-effectiveness programmes will provide a benchmark and evidence for sustainable investment in genomics nationally and globally., Competing Interests: Declaration of interests We declare no competing interests., (Copyright © 2024 The Author(s). Published by Elsevier Ltd.. All rights reserved.)

Pangenome graphs can represent all variation between multiple reference genomes, but current approaches to build them exclude complex sequences or are based upon a single reference. In response, we developed the PanGenome Graph Builder (PGGB), a pipeline for constructing pangenome graphs without bias or exclusion. PGGB uses all-to-all alignments to build a variation graph in which we can identify variation, measure conservation, detect recombination events, and infer phylogenetic relationships.

Human respiratory syncytial virus (RSV) is a major cause of acute respiratory infection. In 2020, RSV was effectively eliminated from the community in New Zealand due to non-pharmaceutical interventions (NPI) used to control the spread of COVID-19. However, in April 2021, following a brief quarantine-free travel agreement with Australia, there was a large-scale nationwide outbreak of RSV that led to reported cases more than five times higher, and hospitalisations more than three times higher, than the typical seasonal pattern. In this study, we generated 1,471 viral genomes of both RSV-A and RSV-B sampled between 2015 and 2022 from across New Zealand. Using a phylodynamics approach, we used these data to better understand RSV transmission patterns in New Zealand prior to 2020, and how RSV became re-established in the community following the relaxation of COVID-19 restrictions. We found that in 2021, there was a large epidemic of RSV in New Zealand that affected a broader age group range compared to the usual pattern of RSV infections. This epidemic was due to an increase in RSV importations, leading to several large genomic clusters of both RSV-A ON1 and RSV-B BA9 genotypes in New Zealand. However, while a number of importations were detected, there was also a major reduction in RSV genetic diversity compared to pre-pandemic seasonal outbreaks. These genomic clusters were temporally associated with the increase of migration in 2021 due to quarantine-free travel from Australia at the time. The closest genetic relatives to the New Zealand RSV genomes, when sampled, were viral genomes sampled in Australia during a large, off-season summer outbreak several months prior, rather than cryptic lineages that were sustained but not detected in New Zealand. These data reveal the impact of NPI used during the COVID-19 pandemic on other respiratory infections and highlight the important insights that can be gained from viral genomes.

Background: New Zealand's (NZ) complete absence of community transmission of influenza and respiratory syncytial virus (RSV) after May 2020, likely due to COVID-19 elimination measures, provided a rare opportunity to assess the impact of border restrictions on common respiratory viral infections over the ensuing 2 years., Methods: We collected the data from multiple surveillance systems, including hospital-based severe acute respiratory infection surveillance, SHIVERS-II, -III and -IV community cohorts for acute respiratory infection (ARI) surveillance, HealthStat sentinel general practice (GP) based influenza-like illness surveillance and SHIVERS-V sentinel GP-based ARI surveillance, SHIVERS-V traveller ARI surveillance and laboratory-based surveillance. We described the data on influenza, RSV and other respiratory viral infections in NZ before, during and after various stages of the COVID related border restrictions., Results: We observed that border closure to most people, and mandatory government-managed isolation and quarantine on arrival for those allowed to enter, appeared to be effective in keeping influenza and RSV infections out of the NZ community. Border restrictions did not affect community transmission of other respiratory viruses such as rhinovirus and parainfluenza virus type-1. Partial border relaxations through quarantine-free travel with Australia and other countries were quickly followed by importation of RSV in 2021 and influenza in 2022., Conclusion: Our findings inform future pandemic preparedness and strategies to model and manage the impact of influenza and other respiratory viral threats., (© 2024 The Authors. Influenza and Other Respiratory Viruses published by John Wiley & Sons Ltd.)

Background: The Global Typhoid Genomics Consortium was established to bring together the typhoid research community to aggregate and analyse Salmonella enterica serovar Typhi (Typhi) genomic data to inform public health action. This analysis, which marks 22 years since the publication of the first Typhi genome, represents the largest Typhi genome sequence collection to date (n=13,000)., Methods: This is a meta-analysis of global genotype and antimicrobial resistance (AMR) determinants extracted from previously sequenced genome data and analysed using consistent methods implemented in open analysis platforms GenoTyphi and Pathogenwatch., Results: Compared with previous global snapshots, the data highlight that genotype 4.3.1 (H58) has not spread beyond Asia and Eastern/Southern Africa; in other regions, distinct genotypes dominate and have independently evolved AMR. Data gaps remain in many parts of the world, and we show the potential of travel-associated sequences to provide informal 'sentinel' surveillance for such locations. The data indicate that ciprofloxacin non-susceptibility (>1 resistance determinant) is widespread across geographies and genotypes, with high-level ciprofloxacin resistance (≥3 determinants) reaching 20% prevalence in South Asia. Extensively drug-resistant (XDR) typhoid has become dominant in Pakistan (70% in 2020) but has not yet become established elsewhere. Ceftriaxone resistance has emerged in eight non-XDR genotypes, including a ciprofloxacin-resistant lineage (4.3.1.2.1) in India. Azithromycin resistance mutations were detected at low prevalence in South Asia, including in two common ciprofloxacin-resistant genotypes., Conclusions: The consortium's aim is to encourage continued data sharing and collaboration to monitor the emergence and global spread of AMR Typhi, and to inform decision-making around the introduction of typhoid conjugate vaccines (TCVs) and other prevention and control strategies., Funding: No specific funding was awarded for this meta-analysis. Coordinators were supported by fellowships from the European Union (ZAD received funding from the European Union's Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie grant agreement No 845681), the Wellcome Trust (SB, Wellcome Trust Senior Fellowship), and the National Health and Medical Research Council (DJI is supported by an NHMRC Investigator Grant [GNT1195210])., Competing Interests: MC, ZD, DI, AA, MA, MC, KC, JC, BH, KK, MM, CP, SV, HW, AA, AA, SA, JA, PA, BB, AB, JC, Kd, AD, Jd, PD, CD, SD, ME, LF, DG, GG, MG, AG, CG, MG, RH, RH, YH, JH, OI, JI, JJ, CJ, DJ, JJ, GK, AK, AK, AK, SK, RK, RK, AL, ML, RL, SL, GM, TM, CM, AM, GN, SN, SN, TN, SN, EN, IO, SP, AP, FQ, FQ, SR, SR, DR, PR, RR, TR, JR, SS, SS, KS, MS, JS, JS, VS, JS, RS, SS, MS, AS, AS, KT, DT, AT, MT, MT, RT, NT, ST, KV, MV, BV, FW, JW, GD, SA, JK, DA, SB, KH No competing interests declared, NF NAF chairs the Wellcome Surveillance and Epidemiology of Drug Resistant Infections (SEDRIC) group, which has a focus on antimicrobial resistance. This could be perceived as relevant although not a direct conflict, IB IB has consulted to BlueDot and the NHL Players' Association, AP AJP is chair of the UK Department of Health and Social Care's (DHSC) Joint Committee on Vaccination and Immunisation (JCVI) but does not take part in the JCVI COVID-19 committee. He was a member of WHO SAGE until 2022. AJPs employer, Oxford University has entered into a partnership with AstraZeneca for development of a COVID-19 vaccine. AJP has provided advice to Shionogi & Co., Ltd on development of a COVID19 vaccine, (© 2023, Carey et al.)

Whole genome sequencing has revolutionized infectious disease surveillance for tracking and monitoring the spread and evolution of pathogens. However, using a linear reference genome for genomic analyses may introduce biases, especially when studies are conducted on highly variable bacterial genomes of the same species. Pangenome graphs provide an efficient model for representing and analyzing multiple genomes and their variants as a graph structure that includes all types of variations. In this study, we present a practical bioinformatics pipeline that employs the PanGenome Graph Builder and the Variation Graph toolkit to build pangenomes from assembled genomes, align whole genome sequencing data and call variants against a graph reference. The pangenome graph enables the identification of structural variants, rearrangements, and small variants (e.g., single nucleotide polymorphisms and insertions/deletions) simultaneously. We demonstrate that using a pangenome graph, instead of a single linear reference genome, improves mapping rates and variant calling for both simulated and real datasets of the pathogen Neisseria meningitidis . Overall, pangenome graphs offer a promising approach for comparative genomics and comprehensive genetic variation analysis in infectious disease. Moreover, this innovative pipeline, leveraging pangenome graphs, can bridge variant analysis, genome assembly, population genetics, and evolutionary biology, expanding the reach of genomic understanding and applications., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Yang, Guarracino, Biggs, Black, Ismail, Wold, Merriman, Prins, Garrison and de Ligt.)

Background: Genomic technologies have become routine in the surveillance and monitoring of the coronavirus disease 2019 (COVID-19) pandemic, as evidenced by the millions of SARS-CoV-2 sequences uploaded to international databases. Yet the ways in which these technologies have been applied to manage the pandemic are varied., Main Text: Aotearoa New Zealand was one of a small number of countries to adopt an elimination strategy for COVID-19, establishing a managed isolation and quarantine system for all international arrivals. To aid our response, we rapidly set up and scaled our use of genomic technologies to help identify community cases of COVID-19, to understand how they had arisen, and to determine the appropriate action to maintain elimination. Once New Zealand pivoted from elimination to suppression in late 2021, our genomic response changed to focusing on identifying new variants arriving at the border, tracking their incidence around the country, and examining any links between specific variants and increased disease severity. Wastewater detection, quantitation and variant detection were also phased into the response. Here, we explore New Zealand's genomic journey through the pandemic and provide a high-level overview of the lessons learned and potential future capabilities to better prepare for future pandemics., Conclusions: Our commentary is aimed at health professionals and decision-makers who might not be familiar with genetic technologies, how they can be used, and why this is an area with great potential to assist in disease detection and tracking now and in the future., (© 2023. The Author(s).)

In the second quarter of 2022, there was a global surge of emergent SARS-CoV-2 lineages that had a distinct growth advantage over then-dominant Omicron BA.1 and BA.2 lineages. By generating 10,403 Omicron genomes, we show that Aotearoa New Zealand observed an influx of these immune-evasive variants (BA.2.12.1, BA.4, and BA.5) through the border. This is explained by the return to significant levels of international travel following the border's reopening in March 2022. We estimate one Omicron transmission event from the border to the community for every ~5,000 passenger arrivals at the current levels of travel and restriction. Although most of these introductions did not instigate any detected onward transmission, a small minority triggered large outbreaks. Genomic surveillance at the border provides a lens on the rate at which new variants might gain a foothold and trigger new waves of infection., (© 2022. The Author(s).)

New Zealand's COVID-19 elimination strategy heavily relied on the use of genomics to inform contact tracing, linking cases to the border and to clusters during community outbreaks. In August 2021, New Zealand entered its second nationwide lockdown after the detection of a single community case with no immediately apparent epidemiological link to the border. This incursion resulted in the largest outbreak seen in New Zealand caused by the Delta Variant of Concern. Here we generated 3806 high quality SARS-CoV-2 genomes from cases reported in New Zealand between 17 August and 1 December 2021, representing 43% of reported cases. We detected wide geographical spread coupled with undetected community transmission, characterised by the apparent extinction and reappearance of genomically linked clusters. We also identified the emergence, and near replacement, of genomes possessing a 10-nucleotide frameshift deletion that caused the likely truncation of accessory protein ORF7a. By early October, New Zealand moved from an elimination strategy to a suppression strategy and the role of genomics changed markedly from being used to track and trace, towards population-level surveillance., (© 2022. The Author(s).)

In New Zealand, international arrivals are quarantined and undergo severe acute respiratory syndrome coronavirus 2 screening; those who test positive are transferred to a managed isolation facility (MIF). Solo traveler A and person E from a 5-person travel group (BCDEF) tested positive. After transfer to the MIF, person A and group BCDEF occupied rooms >2 meters apart across a corridor. Persons B, C, and D subsequently tested positive; viral sequences matched A and were distinct from E. The MIF was the only shared location of persons A and B, C, and D, and they had no direct contact. Security camera footage revealed 4 brief episodes of simultaneous door opening during person A's infectious period. This public health investigation demonstrates transmission from A to B, C, and D while in the MIF, with airborne transmission the most plausible explanation. These findings are of global importance for coronavirus disease public health interventions and infection control practices.

Calcium nitrate has been reported to benefit reproductive outcomes in sows and their offspring when administered via the feed (15 to 19 mg/kg-body weight [bw]/day) during the periparturient period. Traditionally, dietary nitrate had been considered a methemoglobinemia (MetHb) risk in swine. Similar hazard concerns have existed in humans, but a recent benefit/risk analysis established that nitrate levels associated with well-recognized health benefits outweigh potential risks. A similar benefit/risk perspective in swine was lacking and challenged by sparse published hazard data, often referenced within larger reviews related to all livestock. The objective of this review was to better characterize the potential for adverse health and performance effects reported in the literature for swine consuming nitrate and to provide metrics for evaluating the reliability of the studies reviewed. Supplemental exposure via feed or drinking water was considered for any life stage, dose, and exposure duration. More than 30 relevant studies, including case reports and reviews, examined calcium, potassium, sodium, or unspecified nitrate salts at doses up to 1,800 mg nitrate/kg-bw/day for exposures ranging from 1 to 105 d. The studies primarily evaluated weight gain, blood methemoglobin levels, or vitamin A homeostasis in sows or growing swine. An extensive review of the literature showed reports of adverse effects at low nitrate doses to be of low reliability. Conversely, reliable studies corroborate nitrate intake from feed or drinking water at levels equal to or greater than the European Food Safety Authority's no-observed-adverse-effect level (NOAEL) for swine of 410 mg nitrate/kg-bw/day, with no MetHb or other adverse effects on reproduction, growth, or vitamin A levels. Using a weight-of-evidence evaluation, we have moderate-to-high confidence that the NOAEL for nitrate supplementation in swine is likely between 600 and 800 mg/kg-bw/day. These levels are several-fold higher than dietary nitrate concentrations (19 mg/kg-bw/day) that are known to benefit birth outcomes in sows. This review elucidates the quality and reliability of the information sources historically used to characterize nitrate in swine feed as a contaminant of concern. Results from this evaluation can assist risk managers (e.g., regulatory officials and veterinarians) in consideration of proposed benefits as well as reassuring swine producers that low-level nitrate supplementation is not anticipated to be a concern., (© The Author(s) 2021. Published by Oxford University Press on behalf of the American Society of Animal Science.)

Background: COVID-19 elimination measures, including border closures have been applied in New Zealand. We have modelled the potential effect of vaccination programmes for opening borders. Methods: We used a deterministic age-stratified Susceptible, Exposed, Infectious, Recovered (SEIR) model. We minimised spread by varying the age-stratified vaccine allocation to find the minimum herd immunity requirements (the effective reproduction number R eff <1 with closed borders) under various vaccine effectiveness (VE) scenarios and R 0 values. We ran two-year open-border simulations for two vaccine strategies: minimising R eff and targeting high-risk groups. Findings: Targeting of high-risk groups will result in lower hospitalisations and deaths in most scenarios. Reaching the herd immunity threshold (HIT) with a vaccine of 90% VE against disease and 80% VE against infection requires at least 86•5% total population uptake for R 0 =4•5 (with high vaccination coverage for 30-49-year-olds) and 98•1% uptake for R 0 =6. In a two-year open-border scenario with 10 overseas cases daily and 90% total population vaccine uptake (including 0-15 year olds) with the same vaccine, the strategy of targeting high-risk groups is close to achieving HIT, with an estimated 11,400 total hospitalisations (peak 324 active and 36 new daily cases in hospitals), and 1,030 total deaths. Interpretation: Targeting high-risk groups for vaccination will result in fewer hospitalisations and deaths with open borders compared to targeting reduced transmission. With a highly effective vaccine and a high total uptake, opening borders will result in increasing cases, hospitalisations, and deaths. Other public health and social measures will still be required as part of an effective pandemic response. Funding: This project was funded by the Health Research Council [20/1018]. Research in context ., Competing Interests: DM is a member of COVID-19 Vaccine Strategy Taskforce, NZ Government; COVID-19 Vaccine Strategy Scientific and Technical Advisory Group, NZ Government; Advisory Group, Vaccine Alliance Aotearoa New Zealand (VAANZ); COVID-19 Expert Advisory Network, NZ Ministry of Health; and an independent member of the Clinical Trials Steering Committee, University of Oxford COVID-19 Vaccine trials. CRS received COVID-19 related grant funding from the NZ Health Research Council, NZ Ministry of Business, Innovation and Employment, Chief Scientist Office Scotland, UK National Institute for Health Research and UK Medical Research Council. JdL, NF and TN received COVID-19 related grant funding from the NZ Health Research Council, and NZ Ministry of Business, Innovation and Employment. BN and DW received COVID-19 related grant funding from the NZ Health Research Council. DW received grant funding from the Ministry of Business, Innovation and Employment as part of the Te Pūnaha Matatini COVID-19 Modelling Programme. SJ received COVID-19 related grant funding from the NZ Ministry of Business, Innovation and Employment. SJ and JdL's institute (ESR) received funding from the New Zealand Ministry of Health to undertake national infectious disease surveillance.The authors declare no conflict of interest., (© 2021 The Author(s). Published by Elsevier Ltd.)

Late in 2020, two genetically-distinct clusters of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) with mutations of biological concern were reported, one in the United Kingdom and one in South Africa. Using a combination of data from routine surveillance, genomic sequencing and international travel we track the international dispersal of lineages B.1.1.7 and B.1.351 (variant 501Y-V2). We account for potential biases in genomic surveillance efforts by including passenger volumes from location of where the lineage was first reported, London and South Africa respectively. Using the software tool grinch (global report investigating novel coronavirus haplotypes), we track the international spread of lineages of concern with automated daily reports, Further, we have built a custom tracking website (cov-lineages.org/global&#95;report.html) which hosts this daily report and will continue to include novel SARS-CoV-2 lineages of concern as they are detected., Competing Interests: No competing interests were disclosed., (Copyright: © 2021 O'Toole Á et al.)

Since severe acute respiratory syndrome coronavirus 2 was first eliminated in New Zealand in May 2020, a total of 13 known coronavirus disease (COVID-19) community outbreaks have occurred, 2 of which led health officials to issue stay-at-home orders. These outbreaks originated at the border via isolating returnees, airline workers, and cargo vessels. Because a public health system was informed by real-time viral genomic sequencing and complete genomes typically were available within 12 hours of community-based positive COVID-19 test results, every outbreak was well-contained. A total of 225 community cases resulted in 3 deaths. Real-time genomics were essential for establishing links between cases when epidemiologic data could not do so and for identifying when concurrent outbreaks had different origins.

New Zealand, Australia, Iceland, and Taiwan all saw success in controlling their first waves of Coronavirus Disease 2019 (COVID-19). As islands, they make excellent case studies for exploring the effects of international travel and human movement on the spread of COVID-19. We employed a range of robust phylodynamic methods and genome subsampling strategies to infer the epidemiological history of Severe acute respiratory syndrome coronavirus 2 in these four countries. We compared these results to transmission clusters identified by the New Zealand Ministry of Health by contact tracing strategies. We estimated the effective reproduction number of COVID-19 as 1-1.4 during early stages of the pandemic and show that it declined below 1 as human movement was restricted. We also showed that this disease was introduced many times into each country and that introductions slowed down markedly following the reduction of international travel in mid-March 2020. Finally, we confirmed that New Zealand transmission clusters identified via standard health surveillance strategies largely agree with those defined by genomic data. We have demonstrated how the use of genomic data and computational biology methods can assist health officials in characterising the epidemiology of viral epidemics and for contact tracing., (© The Author(s) 2021. Published by Oxford University Press.)

Canonical Wnt signaling plays a key role during organ development, homeostasis and regeneration and these processes are conserved between invertebrates and vertebrates. Mutations in Wnt pathway components are commonly found in various types of cancer. Upon activation of canonical Wnt signaling, β-catenin binds in the nucleus to members of the TCF-LEF family and activates the transcription of target genes. Multiple Wnt target genes, including Lgr5/LGR5 and Axin2/AXIN2, have been identified in mouse models and human cancer cell lines. Here we set out to identify the transcriptional targets of Wnt signaling in five human tissues using organoid technology. Organoids are derived from adult stem cells and recapitulate the functionality as well as the structure of the original tissue. Since the Wnt pathway is critical to maintain the organoids from the human intestine, colon, liver, pancreas and stomach, organoid technology allows us to assess Wnt target gene expression in a human wildtype situation. We performed bulk mRNA sequencing of organoids immediately after inhibition of Wnt pathway and identified 41 genes as commonly regulated genes in these tissues. We also identified large numbers of target genes specific to each tissue. One of the shared target genes is TEAD4, a transcription factor driving expression of YAP/TAZ signaling target genes. In addition to TEAD4, we identified a variety of genes which encode for proteins that are involved in Wnt-independent pathways, implicating the possibility of direct crosstalk between Wnt signaling and other pathways. Collectively, this study identified tissue-specific and common Wnt target gene signatures and provides evidence for a conserved role for these Wnt targets in different tissues., Competing Interests: Declaration of competing interest The authors declare that they have no conflict of interest., (Copyright © 2021 The Authors. Published by Elsevier Inc. All rights reserved.)

Real-time genomic sequencing has played a major role in tracking the global spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), contributing greatly to disease mitigation strategies. In August 2020, after having eliminated the virus, New Zealand experienced a second outbreak. During that outbreak, New Zealand used genomic sequencing in a primary role, leading to a second elimination of the virus. We generated genomes from 78% of the laboratory-confirmed samples of SARS-CoV-2 from the second outbreak and compared them with the available global genomic data. Genomic sequencing rapidly identified that virus causing the second outbreak in New Zealand belonged to a single cluster, thus resulting from a single introduction. However, successful identification of the origin of this outbreak was impeded by substantial biases and gaps in global sequencing data. Access to a broader and more heterogenous sample of global genomic data would strengthen efforts to locate the source of any new outbreaks.

The strategy in New Zealand (Aotearoa) to eliminate coronavirus disease requires that international arrivals undergo managed isolation and quarantine and mandatory testing for severe acute respiratory syndrome coronavirus 2. Combining genomic and epidemiologic data, we investigated the origin of an acute case of coronavirus disease identified in the community after the patient had spent 14 days in managed isolation and quarantine and had 2 negative test results. By combining genomic sequence analysis and epidemiologic investigations, we identified a multibranched chain of transmission of this virus, including on international and domestic flights, as well as a probable case of aerosol transmission without direct person-to-person contact. These findings show the power of integrating genomic and epidemiologic data to inform outbreak investigations.

Since the first wave of coronavirus disease in March 2020, citizens and permanent residents returning to New Zealand have been required to undergo managed isolation and quarantine (MIQ) for 14 days and mandatory testing for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). As of October 20, 2020, of 62,698 arrivals, testing of persons in MIQ had identified 215 cases of SARS-CoV-2 infection. Among 86 passengers on a flight from Dubai, United Arab Emirates, that arrived in New Zealand on September 29, test results were positive for 7 persons in MIQ. These passengers originated from 5 different countries before a layover in Dubai; 5 had negative predeparture SARS-CoV-2 test results. To assess possible points of infection, we analyzed information about their journeys, disease progression, and virus genomic data. All 7 SARS-CoV-2 genomes were genetically identical, except for a single mutation in 1 sample. Despite predeparture testing, multiple instances of in-flight SARS-CoV-2 transmission are likely.

New Zealand, a geographically remote Pacific island with easily sealable borders, implemented a nationwide 'lockdown' of all non-essential services to curb the spread of COVID-19. Here, we generate 649 SARS-CoV-2 genome sequences from infected patients in New Zealand with samples collected during the 'first wave', representing 56% of all confirmed cases in this time period. Despite its remoteness, the viruses imported into New Zealand represented nearly all of the genomic diversity sequenced from the global virus population. These data helped to quantify the effectiveness of public health interventions. For example, the effective reproductive number, R e of New Zealand's largest cluster decreased from 7 to 0.2 within the first week of lockdown. Similarly, only 19% of virus introductions into New Zealand resulted in ongoing transmission of more than one additional case. Overall, these results demonstrate the utility of genomic pathogen surveillance to inform public health and disease mitigation.

The SARS-CoV-2 epidemic has rapidly spread outside China with major outbreaks occurring in Italy, South Korea, and Iran. Phylogenetic analyses of whole-genome sequencing data identified a distinct SARS-CoV-2 clade linked to travellers returning from Iran to Australia and New Zealand. This study highlights potential viral diversity driving the epidemic in Iran, and underscores the power of rapid genome sequencing and public data sharing to improve the detection and management of emerging infectious diseases., (© The Author(s) 2020. Published by Oxford University Press.)

N-(n-butyl) thiophosphoric triamide (NBPT) (Figure 1) is an active ingredient in nitrogen stabilizer (urease inhibitor), which temporarily inhibits the action of the urease enzyme to improve the efficiency of urea-containing fertilizers. Given the potential for NBPT residues to be present in milk and tissues of dairy cattle, due diligence is needed to demonstrate the safety of NBPT in urea-based fertilizers used on forages and crops intended for consumption by Holstein dairy cows. This study used controlled dosing of NBPT in capsule form to dairy cattle for 28 d, followed by a 14-d depuration phase to assess the potential for residues to exist in milk and tissues of dairy cattle at exaggerated use levels. Fourteen lactating cows were selected for the dosing and depuration phases of the study, based on health, body weight (BW), and milk production. There were four treatment groups: 0 mg NBPT/kg BW (Control) ( n = 2 cows), 1 mg NBPT/kg BW (1×) ( n = 3 cows), 3 mg NBPT/kg BW (3×) ( n = 3 cows), and 10 mg NBPT/kg/BW (10×) ( n = 6 cows); levels were based on maximum tolerable amount of urea that a cow can ingest on a daily basis (1×) and the maximum concentration of NBPT commercially used when treating urea (0.1 wt% NBPT in urea). At the end of the 28-d dosing phase, cows were randomly selected for the 14-d depuration phase of the study (one control and three 10× cows). The results showed no NBPT residue is detectable at all dose levels, except that a residue level was above the lower limit of quantitation in a single milk and subcutaneous fat sample in the highest (10×) treatment group, which represents the level of NBPT that would be theoretically present in 10× the lethal dosing of daily consumable urea to a cow. Overall, the study demonstrated that it is unlikely for NBPT residues to be present in cattle milk or edible tissues or to cause negative effects on animal health under good agricultural practice., (© The Author(s) 2019. Published by Oxford University Press on behalf of the American Society of Animal Science.)

Male breast cancer (MBC) is extremely rare and accounts for less than 1% of all breast malignancies. Therefore, clinical management of MBC is currently guided by research on the disease in females. In this study, DNA obtained from 45 formalin-fixed paraffin-embedded (FFPE) MBCs with and 90 MBCs (52 FFPE and 38 fresh-frozen) without matched normal tissues was subjected to massively parallel sequencing targeting all exons of 1943 cancer-related genes. The landscape of mutations and copy number alterations was compared to that of publicly available estrogen receptor (ER)-positive female breast cancers (smFBCs) and correlated to prognosis. From the 135 MBCs, 90% showed ductal histology, 96% were ER-positive, 66% were progesterone receptor (PR)-positive, and 2% HER2-positive, resulting in 50, 46 and 4% luminal A-like, luminal B-like and basal-like cases, respectively. Five patients had Klinefelter syndrome (4%) and 11% of patients harbored pathogenic BRCA2 germline mutations. The genomic landscape of MBC to some extent recapitulated that of smFBC, with recurrent PIK3CA (36%) and GATA3 (15%) somatic mutations, and with 40% of the most frequently amplified genes overlapping between both sexes. TP53 (3%) somatic mutations were significantly less frequent in MBC compared to smFBC, whereas somatic mutations in genes regulating chromatin function and homologous recombination deficiency-related signatures were more prevalent. MDM2 amplifications were frequent (13%), correlated with protein overexpression (P = 0.001) and predicted poor outcome (P = 0.007). In conclusion, despite similarities in the genomic landscape between MBC and smFBC, MBC is a molecularly unique and heterogeneous disease requiring its own clinical trials and treatment guidelines.

Nucleotide excision repair (NER) is one of the main DNA repair pathways that protect cells against genomic damage. Disruption of this pathway can contribute to the development of cancer and accelerate aging. Mutational characteristics of NER-deficiency may reveal important diagnostic opportunities, as tumors deficient in NER are more sensitive to certain treatments. Here, we analyzed the genome-wide somatic mutational profiles of adult stem cells (ASCs) from NER-deficient Ercc1 -/Δ mice. Our results indicate that NER-deficiency increases the base substitution load twofold in liver but not in small intestinal ASCs, which coincides with the tissue-specific aging pathology observed in these mice. Moreover, NER-deficient ASCs of both tissues show an increased contribution of Signature 8 mutations, which is a mutational pattern with unknown etiology that is recurrently observed in various cancer types. The scattered genomic distribution of the base substitutions indicates that deficiency of global-genome NER (GG-NER) underlies the observed mutational consequences. In line with this, we observe increased Signature 8 mutations in a GG-NER-deficient human organoid culture, in which XPC was deleted using CRISPR-Cas9 gene-editing. Furthermore, genomes of NER-deficient breast tumors show an increased contribution of Signature 8 mutations compared with NER-proficient tumors. Elevated levels of Signature 8 mutations could therefore contribute to a predictor of NER-deficiency based on a patient's mutational profile., (© 2019 Jager et al.; Published by Cold Spring Harbor Laboratory Press.)

Organoids are self-organizing 3D structures grown from stem cells that recapitulate essential aspects of organ structure and function. Here, we describe a method to establish long-term-expanding human airway organoids from broncho-alveolar resections or lavage material. The pseudostratified airway organoids consist of basal cells, functional multi-ciliated cells, mucus-producing secretory cells, and CC10-secreting club cells. Airway organoids derived from cystic fibrosis (CF) patients allow assessment of CFTR function in an organoid swelling assay. Organoids established from lung cancer resections and metastasis biopsies retain tumor histopathology as well as cancer gene mutations and are amenable to drug screening. Respiratory syncytial virus (RSV) infection recapitulates central disease features, dramatically increases organoid cell motility via the non-structural viral NS2 protein, and preferentially recruits neutrophils upon co-culturing. We conclude that human airway organoids represent versatile models for the in vitro study of hereditary, malignant, and infectious pulmonary disease., (© 2019 The Authors. Published under the terms of the CC BY NC ND 4.0 license.)

Objective: The nuclear receptor PPARγ is the master regulator of adipocyte differentiation, distribution, and function. In addition, PPARγ induces terminal differentiation of several epithelial cell lineages, including colon epithelia. Loss-of-function mutations in PPARG result in familial partial lipodystrophy subtype 3 (FPDL3), a rare condition characterized by aberrant adipose tissue distribution and severe metabolic complications, including diabetes. Mutations in PPARG have also been reported in sporadic colorectal cancers, but the significance of these mutations is unclear. Studying these natural PPARG mutations provides valuable insights into structure-function relationships in the PPARγ protein. We functionally characterized a novel FPLD3-associated PPARγ L451P mutation in helix 9 of the ligand binding domain (LBD). Interestingly, substitution of the adjacent amino acid K450 was previously reported in a human colon carcinoma cell line., Methods: We performed a detailed side-by-side functional comparison of these two PPARγ mutants., Results: PPARγ L451P shows multiple intermolecular defects, including impaired cofactor binding and reduced RXRα heterodimerisation and subsequent DNA binding, but not in DBD-LBD interdomain communication. The K450Q mutant displays none of these functional defects. Other colon cancer-associated PPARγ mutants displayed diverse phenotypes, ranging from complete loss of activity to wildtype activity., Conclusions: Amino acid changes in helix 9 can differently affect LBD integrity and function. In addition, FPLD3-associated PPARγ mutations consistently cause intra- and/or intermolecular defects; colon cancer-associated PPARγ mutations on the other hand may play a role in colon cancer onset and progression, but this is not due to their effects on the most well-studied functional characteristics of PPARγ., (Copyright © 2018 The Authors. Published by Elsevier GmbH.. All rights reserved.)

Biological, clinical, and pharmacological research now often involves analyses of genomes, transcriptomes, proteomes, and interactomes, within and between individuals and across species. Due to large volumes, the analysis and integration of data generated by such high-throughput technologies have become computationally intensive, and analysis can no longer happen on a typical desktop computer.In this chapter we show how to describe and execute the same analysis using a number of workflow systems and how these follow different approaches to tackle execution and reproducibility issues. We show how any researcher can create a reusable and reproducible bioinformatics pipeline that can be deployed and run anywhere. We show how to create a scalable, reusable, and shareable workflow using four different workflow engines: the Common Workflow Language (CWL), Guix Workflow Language (GWL), Snakemake, and Nextflow. Each of which can be run in parallel.We show how to bundle a number of tools used in evolutionary biology by using Debian, GNU Guix, and Bioconda software distributions, along with the use of container systems, such as Docker, GNU Guix, and Singularity. Together these distributions represent the overall majority of software packages relevant for biology, including PAML, Muscle, MAFFT, MrBayes, and BLAST. By bundling software in lightweight containers, they can be deployed on a desktop, in the cloud, and, increasingly, on compute clusters.By bundling software through these public software distributions, and by creating reproducible and shareable pipelines using these workflow engines, not only do bioinformaticians have to spend less time reinventing the wheel but also do we get closer to the ideal of making science reproducible. The examples in this chapter allow a quick comparison of different solutions.

Intrinsic and/or acquired resistance represents one of the great challenges in targeted cancer therapy. A deeper understanding of the molecular biology of cancer has resulted in more efficient strategies, where one or multiple drugs are adopted in novel therapies to tackle resistance. This beneficial effect of using combination treatments has also been observed in colorectal cancer patients harboring the BRAF(V600E) mutation, whereby dual inhibition of BRAF(V600E) and EGFR increases antitumor activity. Notwithstanding this success, it is not clear whether this combination treatment is the only or most effective treatment to block intrinsic resistance to BRAF inhibitors. Here, we investigate molecular responses upon single and multi-target treatments, over time, using BRAF(V600E) mutant colorectal cancer cells as a model system. Through integration of transcriptomic, proteomic and phosphoproteomics data we obtain a comprehensive overview, revealing both known and novel responses. We primarily observe widespread up-regulation of receptor tyrosine kinases and metabolic pathways upon BRAF inhibition. These findings point to mechanisms by which the drug-treated cells switch energy sources and enter a quiescent-like state as a defensive response, while additionally compensating for the MAPK pathway inhibition., (© 2018 by The American Society for Biochemistry and Molecular Biology, Inc.)