Mobster: Accurate detection of mobile element insertions in next generation sequencing data

Mobile elements are major drivers in changing genomic architecture and can cause disease. The detection of mobile elements is hindered due to the low mappability of their highly repetitive sequences. We have developed an algorithm, called Mobster, to detect non-reference mobile element insertions in next generation sequencing data from both whole genome and whole exome studies. Mobster uses discordant read pairs and clipped reads in combination with consensus sequences of known active mobile elements. Mobster has a low false discovery rate and high recall rate for both L1 and Alu elements. Mobster is available at http://sourceforge.net/projects/mobster. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1186/s13059-014-0488-x) contains supplementary material, which is available to authorized users

Identification of novel human Wnt target genes using adult endodermal tissue-derived organoids

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

Use of CRISPR-modified human stem cell organoids to study the origin of mutational signatures in cancer.

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

COVID-19 vaccine strategies for Aotearoa New Zealand:a mathematical modelling study

Summary: 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 Reff<1 with closed borders) under various vaccine effectiveness (VE) scenarios and R0 values. We ran two-year open-border simulations for two vaccine strategies: minimising Reff 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 R0=4•5 (with high vaccination coverage for 30–49-year-olds) and 98•1% uptake for R0=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 contex

Long-term culture of genome-stable bipotent stem cells from adult human liver.

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.This work was supported by grants to MH (EU/236954) and to HC (The United European Gastroenterology Federation (UEGF) Research Prize 2010, EU/232814-StemCellMark and NWO/116002008). MH is supported by The Wellcome Trust Sir Henry Dale fellowship. The Rspo cell line was kindly provided by Dr. Calvin Kuo.This is the final published version. It first appeared at http://www.cell.com/abstract/S0092-8674%2814%2901566-9

Global diversity and antimicrobial resistance of typhoid fever pathogens : insights from a meta-analysis of 13,000 Salmonella Typhi genomes

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 becomedominant 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. © Carey et al. **Please note that there are multiple authors for this article therefore only the name of the first 30 including Federation University Australia affiliate “Andrew Greenhill” is provided in this record*

Genomic landscape of rat strain and substrain variation

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

Global disparities in SARS-CoV-2 genomic surveillance

Genomic sequencing is essential to track the evolution and spread of SARS-CoV-2, optimize molecular tests, treatments, vaccines, and guide public health responses. To investigate the global SARS-CoV-2 genomic surveillance, we used sequences shared via GISAID to estimate the impact of sequencing intensity and turnaround times on variant detection in 189 countries. In the first two years of the pandemic, 78% of high-income countries sequenced >0.5% of their COVID-19 cases, while 42% of low- and middle-income countries reached that mark. Around 25% of the genomes from high income countries were submitted within 21 days, a pattern observed in 5% of the genomes from low- and middle-income countries. We found that sequencing around 0.5% of the cases, with a turnaround time <21 days, could provide a benchmark for SARS-CoV-2 genomic surveillance. Socioeconomic inequalities undermine the global pandemic preparedness, and efforts must be made to support low- and middle-income countries improve their local sequencing capacity

Deficiency of nucleotide excision repair is associated with mutational signature observed in cancer

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