Search

Human respiratory syncytial virus (RSV) is a major cause of acute respiratory infection. In 2020, RSV was eliminated from New Zealand due to non-pharmaceutical interventions (NPI) used to control the spread of SARS-CoV-2. However, in 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 than typical seasonal patterns. We generated 1470 viral genomes of both RSV-A and RSV-B sampled between 2015-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 due to an increase in importations, leading to several large genomic clusters of both RSV-A ON1 and RSV-B BA9 genotypes. However, while a number of viral importations were detected, there was also a major reduction in RSV genetic diversity compared to pre-pandemic years. 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., Competing Interests: Competing interests The authors declare no competing interests., (© 2024. The Author(s).)

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, a pipeline for constructing pangenome graphs without bias or exclusion. The PanGenome Graph Builder 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., (© 2024. The Author(s), under exclusive licence to Springer Nature America, Inc.)

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.)

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.)

Genetic lipodystrophies are a group of rare syndromes associated with major metabolic complications - including severe insulin resistance, type 2 diabetes mellitus, and hypertriglyceridemia - which are classified according to the distribution of adipose tissue. Lipodystrophies can be present at birth or develop during life and can range from local to partial and general. With at least 18 different genes implicated so far, definite diagnosis can be challenging due to clinical and genetic heterogeneity. In an adult female patient with clinical and metabolic features of partial lipodystrophy we identified via whole genome sequencing (WGS) a single complex AGPAT2 allele [V67M;V167A], functionally equivalent to heterozygosity. AGPAT2 encodes for an acyltransferase implicated in the biosynthesis of triacylglycerol and glycerophospholipids. So far homozygous and compound heterozygous mutations in AGPAT2 have only been associated with generalized lipodystrophy. A SNP risk score analysis indicated that the index patient is not predisposed to lipodystrophy based on her genetic background. The partial phenotype in our patient is therefore more likely associated to the genetic variants in AGPAT2. To test whether the resulting double-mutant AGPAT2 protein is functional we analyzed its in vitro enzymatic activity via mass spectrometry. The resulting AGPAT2 double mutant is enzymatically inactive. Our data support the view that the current classification of lipodystrophies as strictly local, partial or generalized may have to be re-evaluated and viewed more as a continuum, both in terms of clinical presentation and underlying genetic causes. Better molecular understanding of lipodystrophies may lead to new therapies to treat adipose tissue dysfunction in common and rare diseases.

Recent data showed that cancer cells from different tumor subtypes with distinct metastatic potential influence each other's metastatic behavior by exchanging biomolecules through extracellular vesicles (EVs). However, it is debated how small amounts of cargo can mediate this effect, especially in tumors where all cells are from one subtype, and only subtle molecular differences drive metastatic heterogeneity. To study this, we have characterized the content of EVs shed in vivo by two clones of melanoma (B16) tumors with distinct metastatic potential. Using the Cre-LoxP system and intravital microscopy, we show that cells from these distinct clones phenocopy their migratory behavior through EV exchange. By tandem mass spectrometry and RNA sequencing, we show that EVs shed by these clones into the tumor microenvironment contain thousands of different proteins and RNAs, and many of these biomolecules are from interconnected signaling networks involved in cellular processes such as migration. Thus, EVs contain numerous proteins and RNAs and act on recipient cells by invoking a multi-faceted biological response including cell migration., (© 2018 The Authors. Published under the terms of the CC BY 4.0 license.)

Breast cancer (BC) comprises multiple distinct subtypes that differ genetically, pathologically, and clinically. Here, we describe a robust protocol for long-term culturing of human mammary epithelial organoids. Using this protocol, >100 primary and metastatic BC organoid lines were generated, broadly recapitulating the diversity of the disease. BC organoid morphologies typically matched the histopathology, hormone receptor status, and HER2 status of the original tumor. DNA copy number variations as well as sequence changes were consistent within tumor-organoid pairs and largely retained even after extended passaging. BC organoids furthermore populated all major gene-expression-based classification groups and allowed in vitro drug screens that were consistent with in vivo xeno-transplantations and patient response. This study describes a representative collection of well-characterized BC organoids available for cancer research and drug development, as well as a strategy to assess in vitro drug response in a personalized fashion., (Copyright © 2017 Elsevier Inc. All rights reserved.)

Despite improvements in genomics technology, the detection of structural variants (SVs) from short-read sequencing still poses challenges, particularly for complex variation. Here we analyse the genomes of two patients with congenital abnormalities using the MinION nanopore sequencer and a novel computational pipeline-NanoSV. We demonstrate that nanopore long reads are superior to short reads with regard to detection of de novo chromothripsis rearrangements. The long reads also enable efficient phasing of genetic variations, which we leveraged to determine the parental origin of all de novo chromothripsis breakpoints and to resolve the structure of these complex rearrangements. Additionally, genome-wide surveillance of inherited SVs reveals novel variants, missed in short-read data sets, a large proportion of which are retrotransposon insertions. We provide a first exploration of patient genome sequencing with a nanopore sequencer and demonstrate the value of long-read sequencing in mapping and phasing of SVs for both clinical and research applications.

Recognition of individuals with a genetic predisposition to gastric cancer (GC) enables preventive measures. However, the underlying cause of genetic susceptibility to gastric cancer remains largely unexplained. We performed germline whole-exome sequencing on leukocyte DNA of 54 patients from 53 families with genetically unexplained diffuse-type and intestinal-type GC to identify novel GC-predisposing candidate genes. As young age at diagnosis and familial clustering are hallmarks of genetic tumor susceptibility, we selected patients that were diagnosed below the age of 35, patients from families with two cases of GC at or below age 60 and patients from families with three GC cases at or below age 70. All included individuals were tested negative for germline CDH1 mutations before or during the study. Variants that were possibly deleterious according to in silico predictions were filtered using several independent approaches that were based on gene function and gene mutation burden in controls. Despite a rigorous search, no obvious candidate GC predisposition genes were identified. This negative result stresses the importance of future research studies in large, homogeneous cohorts.

Mutational processes underlie cancer initiation and progression. Signatures of these processes in cancer genomes may explain cancer etiology and could hold diagnostic and prognostic value. We developed a strategy that can be used to explore the origin of cancer-associated mutational signatures. We used CRISPR-Cas9 technology to delete key DNA repair genes in human colon organoids, followed by delayed subcloning and whole-genome sequencing. We found that mutation accumulation in organoids deficient in the mismatch repair gene MLH1 is driven by replication errors and accurately models the mutation profiles observed in mismatch repair-deficient colorectal cancers. Application of this strategy to the cancer predisposition gene NTHL1 , which encodes a base excision repair protein, revealed a mutational footprint (signature 30) previously observed in a breast cancer cohort. We show that signature 30 can arise from germline NTHL1 mutations., (Copyright © 2017 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.)

In the adenoma-carcinoma sequence, it is proposed that intestinal polyps evolve through a set of defined mutations toward metastatic colorectal cancer (CRC). Here, we dissect this adenoma-carcinoma sequence in vivo by using an orthotopic organoid transplantation model of human colon organoids engineered to harbor different CRC mutation combinations. We demonstrate that sequential accumulation of oncogenic mutations in Wnt, EGFR, P53, and TGF-β signaling pathways facilitates efficient tumor growth, migration, and metastatic colonization. We show that reconstitution of specific niche signals can restore metastatic growth potential of tumor cells lacking one of the oncogenic mutations. Our findings imply that the ability to metastasize-i.e., to colonize distant sites-is the direct consequence of the loss of dependency on specific niche signals.

Background: Germline chromothripsis causes complex genomic rearrangements that are likely to affect multiple genes and their regulatory contexts. The contribution of individual rearrangements and affected genes to the phenotypes of patients with complex germline genomic rearrangements is generally unknown., Methods: To dissect the impact of germline chromothripsis in a relevant developmental context, we performed trio-based RNA expression analysis on blood cells, induced pluripotent stem cells (iPSCs), and iPSC-derived neuronal cells from a patient with de novo germline chromothripsis and both healthy parents. In addition, Hi-C and 4C-seq experiments were performed to determine the effects of the genomic rearrangements on transcription regulation of genes in the proximity of the breakpoint junctions., Results: Sixty-seven genes are located within 1 Mb of the complex chromothripsis rearrangements involving 17 breakpoints on four chromosomes. We find that three of these genes (FOXP1, DPYD, and TWIST1) are both associated with developmental disorders and differentially expressed in the patient. Interestingly, the effect on TWIST1 expression was exclusively detectable in the patient's iPSC-derived neuronal cells, stressing the need for studying developmental disorders in the biologically relevant context. Chromosome conformation capture analyses show that TWIST1 lost genomic interactions with several enhancers due to the chromothripsis event, which likely led to deregulation of TWIST1 expression and contributed to the patient's craniosynostosis phenotype., Conclusions: We demonstrate that a combination of patient-derived iPSC differentiation and trio-based molecular profiling is a powerful approach to improve the interpretation of pathogenic complex genomic rearrangements. Here we have applied this approach to identify misexpression of TWIST1, FOXP1, and DPYD as key contributors to the complex congenital phenotype resulting from germline chromothripsis rearrangements.

Purpose: We aimed to identify a novel genetic cause of tooth agenesis (TA) and/or orofacial clefting (OFC) by combining whole-exome sequencing (WES) and targeted resequencing in a large cohort of TA and OFC patients., Methods: WES was performed in two unrelated patients: one with severe TA and OFC and another with severe TA only. After deleterious mutations were identified in a gene encoding low-density lipoprotein receptor-related protein 6 (LRP6), all its exons were resequenced with molecular inversion probes in 67 patients with TA, 1,072 patients with OFC, and 706 controls., Results: We identified a frameshift (c.4594delG, p.Cys1532fs) and a canonical splice-site mutation (c.3398-2A>C, p.?) in LRP6, respectively, in the patient with TA and OFC and in the patient with severe TA only. The targeted resequencing showed significant enrichment of unique LRP6 variants in TA patients but not in nonsyndromic OFC patients. Of the five variants in patients with TA, two affected the canonical splice site and three were missense variants; all variants segregated with the dominant phenotype, and in one case the missense mutation occurred de novo., Conclusion: Mutations in LRP6 cause TA in humans.Genet Med 18 11, 1158-1162.

The gradual accumulation of genetic mutations in human adult stem cells (ASCs) during life is associated with various age-related diseases, including cancer. Extreme variation in cancer risk across tissues was recently proposed to depend on the lifetime number of ASC divisions, owing to unavoidable random mutations that arise during DNA replication. However, the rates and patterns of mutations in normal ASCs remain unknown. Here we determine genome-wide mutation patterns in ASCs of the small intestine, colon and liver of human donors with ages ranging from 3 to 87 years by sequencing clonal organoid cultures derived from primary multipotent cells. Our results show that mutations accumulate steadily over time in all of the assessed tissue types, at a rate of approximately 40 novel mutations per year, despite the large variation in cancer incidence among these tissues. Liver ASCs, however, have different mutation spectra compared to those of the colon and small intestine. Mutational signature analysis reveals that this difference can be attributed to spontaneous deamination of methylated cytosine residues in the colon and small intestine, probably reflecting their high ASC division rate. In liver, a signature with an as-yet-unknown underlying mechanism is predominant. Mutation spectra of driver genes in cancer show high similarity to the tissue-specific ASC mutation spectra, suggesting that intrinsic mutational processes in ASCs can initiate tumorigenesis. Notably, the inter-individual variation in mutation rate and spectra are low, suggesting tissue-specific activity of common mutational processes throughout life.

Structural variation (SV) represents a major source of differences between individual human genomes and has been linked to disease phenotypes. However, the majority of studies provide neither a global view of the full spectrum of these variants nor integrate them into reference panels of genetic variation. Here, we analyse whole genome sequencing data of 769 individuals from 250 Dutch families, and provide a haplotype-resolved map of 1.9 million genome variants across 9 different variant classes, including novel forms of complex indels, and retrotransposition-mediated insertions of mobile elements and processed RNAs. A large proportion are previously under reported variants sized between 21 and 100 bp. We detect 4 megabases of novel sequence, encoding 11 new transcripts. Finally, we show 191 known, trait-associated SNPs to be in strong linkage disequilibrium with SVs and demonstrate that our panel facilitates accurate imputation of SVs in unrelated individuals.

Unmapped next-generation sequencing reads are typically ignored while they contain biologically relevant information. We systematically analyzed unmapped reads from whole genome sequencing of 33 inbred rat strains. High quality reads were selected and enriched for biologically relevant sequences; similarity-based analysis revealed clustering similar to previously reported phylogenetic trees. Our results demonstrate that on average 20% of all unmapped reads harbor sequences that can be used to improve reference genomes and generate hypotheses on potential genotype-phenotype relationships. Analysis pipelines would benefit from incorporating the described methods and reference genomes would benefit from inclusion of the genomic segments obtained through these efforts.

Cerebral visual impairment (CVI) is a major cause of low vision in children due to impairment in projection and/or interpretation of the visual input in the brain. Although acquired causes for CVI are well known, genetic causes underlying CVI are largely unidentified. DNAs of 25 patients with CVI and intellectual disability, but without acquired (eg, perinatal) damage, were investigated by whole-exome sequencing. The data were analyzed for de novo, autosomal-recessive, and X-linked variants, and subsequently classified into known, candidate, or unlikely to be associated with CVI. This classification was based on the Online Mendelian Inheritance in Man database, literature reports, variant characteristics, and functional relevance of the gene. After classification, variants in four genes known to be associated with CVI (AHDC1, NGLY1, NR2F1, PGAP1) in 5 patients (20%) were identified, establishing a conclusive genetic diagnosis for CVI. In addition, in 11 patients (44%) with CVI, variants in one or more candidate genes were identified (ACP6, AMOT, ARHGEF10L, ATP6V1A, DCAF6, DLG4, GABRB2, GRIN1, GRIN2B, KCNQ3, KCTD19, RERE, SLC1A1, SLC25A16, SLC35A2, SOX5, UFSP2, UHMK1, ZFP30). Our findings show that diverse genetic causes underlie CVI, some of which will provide insight into the biology underlying this disease process.

Homozygous variants in PGAP1 (post-GPI attachment to proteins 1) have recently been identified in two families with developmental delay, seizures and/or spasticity. PGAP1 is a member of the glycosylphosphatidylinositol anchor biosynthesis and remodeling pathway and defects in this pathway are a subclass of congenital disorders of glycosylation. Here we performed whole-exome sequencing in an individual with cerebral visual impairment (CVI), intellectual disability (ID), and factor XII deficiency and revealed compound heterozygous variants in PGAP1, c.274&#95;276del (p.(Pro92del)) and c.921&#95;925del (p.(Lys308Asnfs*25)). Subsequently, PGAP1-deficient Chinese hamster ovary (CHO)-cell lines were transfected with either mutant or wild-type constructs and their sensitivity to phosphatidylinositol-specific phospholipase C (PI-PLC) treatment was measured. The mutant constructs could not rescue the PGAP1-deficient CHO cell lines resistance to PI-PLC treatment. In addition, lymphoblastoid cell lines (LCLs) of the affected individual showed no sensitivity to PI-PLC treatment, whereas the LCLs of the heterozygous carrier parents were partially resistant. In conclusion, we report novel PGAP1 variants in a boy with CVI and ID and a proven functional loss of PGAP1 and show, to our knowledge, for the first time this genetic association with CVI.

Linkage analysis combined with whole-exome sequencing in a large family with congenital and stable non-syndromic unilateral and asymmetric hearing loss (NS-UHL/AHL) revealed a heterozygous truncating mutation, c.286&#95;303delinsT (p.Ser96Ter), in KITLG. This mutation co-segregated with NS-UHL/AHL as a dominant trait with reduced penetrance. By screening a panel of probands with NS-UHL/AHL, we found an additional mutation, c.200&#95;202del (p.His67&#95;Cys68delinsArg). In vitro studies revealed that the p.His67&#95;Cys68delinsArg transmembrane isoform of KITLG is not detectable at the cell membrane, supporting pathogenicity. KITLG encodes a ligand for the KIT receptor. Also, KITLG-KIT signaling and MITF are suggested to mutually interact in melanocyte development. Because mutations in MITF are causative of Waardenburg syndrome type 2 (WS2), we screened KITLG in suspected WS2-affected probands. A heterozygous missense mutation, c.310C>G (p.Leu104Val), that segregated with WS2 was identified in a small family. In vitro studies revealed that the p.Leu104Val transmembrane isoform of KITLG is located at the cell membrane, as is wild-type KITLG. However, in culture media of transfected cells, the p.Leu104Val soluble isoform of KITLG was reduced, and no soluble p.His67&#95;Cys68delinsArg and p.Ser96Ter KITLG could be detected. These data suggest that mutations in KITLG associated with NS-UHL/AHL have a loss-of-function effect. We speculate that the mechanism of the mutation underlying WS2 and leading to membrane incorporation and reduced secretion of KITLG occurs via a dominant-negative or gain-of-function effect. Our study unveils different phenotypes associated with KITLG, previously associated with pigmentation abnormalities, and will thereby improve the genetic counseling given to individuals with KITLG variants., (Copyright © 2015 The American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

Background: Since the completion of the rat reference genome in 2003, whole-genome sequencing data from more than 40 rat strains have become available. These data represent the broad range of strains that are used in rat research including commonly used substrains. Currently, this wealth of information cannot be used to its full extent, because the variety of different variant calling algorithms employed by different groups impairs comparison between strains. In addition, all rat whole genome sequencing studies to date used an outdated reference genome for analysis (RGSC3.4 released in 2004)., Results: Here we present a comprehensive, multi-sample and uniformly called set of genetic variants in 40 rat strains, including 19 substrains. We reanalyzed all primary data using a recent version of the rat reference assembly (RGSC5.0 released in 2012) and identified over 12 million genomic variants (SNVs, indels and structural variants) among the 40 strains. 28,318 SNVs are specific to individual substrains, which may be explained by introgression from other unsequenced strains and ongoing evolution by genetic drift. Substrain SNVs may have a larger predicted functional impact compared to older shared SNVs., Conclusions: In summary we present a comprehensive catalog of uniformly analyzed genetic variants among 40 widely used rat inbred strains based on the RGSC5.0 assembly. This represents a valuable resource, which will facilitate rat functional genomic research. In line with previous observations, our genome-wide analyses do not show evidence for contribution of multiple ancestral founder rat subspecies to the currently used rat inbred strains, as is the case for mouse. In addition, we find that the degree of substrain variation is highly variable between strains, which is of importance for the correct interpretation of experimental data from different labs.

Noonan syndrome (NS) is a developmental disorder characterized by short stature, facial dysmorphisms and congenital heart defects. To date, all mutations known to cause NS are dominant, activating mutations in signal transducers of the RAS/mitogen-activated protein kinase (MAPK) pathway. In 25% of cases, however, the genetic cause of NS remains elusive, suggesting that factors other than those involved in the canonical RAS/MAPK pathway may also have a role. Here, we used family-based whole exome sequencing of a case-parent trio and identified a de novo mutation, p.(Arg802His), in A2ML1, which encodes the secreted protease inhibitor α-2-macroglobulin (A2M)-like-1. Subsequent resequencing of A2ML1 in 155 cases with a clinical diagnosis of NS led to the identification of additional mutations in two families, p.(Arg802Leu) and p.(Arg592Leu). Functional characterization of these human A2ML1 mutations in zebrafish showed NS-like developmental defects, including a broad head, blunted face and cardiac malformations. Using the crystal structure of A2M, which is highly homologous to A2ML1, we identified the intramolecular interaction partner of p.Arg802. Mutation of this residue, p.Glu906, induced similar developmental defects in zebrafish, strengthening our conclusion that mutations in A2ML1 cause a disorder clinically related to NS. This is the first report of the involvement of an extracellular factor in a disorder clinically related to RASopathies, providing potential new leads for better understanding of the molecular basis of this family of developmental diseases.

Despite the enormous replication potential of the human liver, there are currently no culture systems available that sustain hepatocyte replication and/or function in vitro. We have shown previously that single mouse Lgr5+ liver stem cells can be expanded as epithelial organoids in vitro and can be differentiated into functional hepatocytes in vitro and in vivo. We now describe conditions allowing long-term expansion of adult bile duct-derived bipotent progenitor cells from human liver. The expanded cells are highly stable at the chromosome and structural level, while single base changes occur at very low rates. The cells can readily be converted into functional hepatocytes in vitro and upon transplantation in vivo. Organoids from α1-antitrypsin deficiency and Alagille syndrome patients mirror the in vivo pathology. Clonal long-term expansion of primary adult liver stem cells opens up experimental avenues for disease modeling, toxicology studies, regenerative medicine, and gene therapy., (Copyright © 2015 The Authors. Published by Elsevier Inc. All rights reserved.)