Somatic variant analysis of linked-reads sequencing data with Lancet.

SUMMARY: We present a new version of the popular somatic variant caller, Lancet, that supports the analysis of linked-reads sequencing data. By seamlessly integrating barcodes and haplotype read assignments within the colored De Bruijn graph local-assembly framework, Lancet computes a barcode-aware coverage and identifies variants that disagree with the local haplotype structure. AVAILABILITY AND IMPLEMENTATION: Lancet is implemented in C++ and available for academic and non-commercial research purposes as an open-source package at https://github.com/nygenome/lancet. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Citizen Science Surveillance of Triazole-Resistant Aspergillus fumigatus in United Kingdom Residential Garden Soils.

Compost is an ecological niche for Aspergillus fumigatus due to its role as a decomposer of organic matter and its ability to survive the high temperatures associated with the composting process. Subsequently, composting facilities are associated with high levels of A. fumigatus spores that are aerosolized from compost and cause respiratory illness in workers. In the UK, gardening is an activity enjoyed by individuals of all ages, and it is likely that they are being exposed to A. fumigatus spores when handling commercial compost or compost they have produced themselves. In the present study, 246 citizen scientists collected 509 soil samples from locations in their gardens in the UK, from which were cultured 5,174 A. fumigatus isolates. Of these isolates, 736 (14%) were resistant to tebuconazole: the third most-sprayed triazole fungicide in the UK, which confers cross-resistance to the medical triazoles used to treat A. fumigatus lung infections in humans. These isolates were found to contain the common resistance mechanisms in the A. fumigatus cyp51A gene TR34/L98H or TR46/Y121F/T289A, as well as the less common resistance mechanisms TR34, TR53, TR46/Y121F/T289A/S363P/I364V/G448S, and (TR46)2/Y121F/M172I/T289A/G448S. Regression analyses found that soil samples containing compost were significantly more likely to grow tebuconazole-susceptible and tebuconazole-resistant A. fumigatus strains than those that did not and that compost samples grew significantly higher numbers of A. fumigatus than other samples. IMPORTANCE The findings presented here highlight compost as a potential health hazard to individuals with predisposing factors to A. fumigatus lung infections and as a potential health hazard to immunocompetent individuals who could be exposed to sufficiently high numbers of spores to develop infection. Furthermore, we found that 14% of A. fumigatus isolates in garden soils were resistant to an agricultural triazole, which confers cross-resistance to medical triazoles used to treat A. fumigatus lung infections. This raises the question of whether compost bags should carry additional health warnings regarding inhalation of A. fumigatus spores, whether individuals should be advised to wear facemasks while handling compost, or whether commercial producers should be responsible for sterilizing compost before shipping. The findings support increasing public awareness of the hazard posed by compost and investigating measures that can be taken to reduce the exposure risk.

Post-epizootic microbiome associations across communities of neotropical amphibians.

Microbiome-pathogen interactions are increasingly recognized as an important element of host immunity. While these host-level interactions will have consequences for community disease dynamics, the factors which influence host microbiomes at larger scales are poorly understood. We here describe landscape-scale pathogen-microbiome associations within the context of post-epizootic amphibian chytridiomycosis, a disease caused by the panzootic chytrid fungus Batrachochytrium dendrobatidis. We undertook a survey of Neotropical amphibians across altitudinal gradients in Ecuador ~30 years following the observed amphibian declines and collected skin swab-samples which were metabarcoded using both fungal (ITS-2) and bacterial (r16S) amplicons. The data revealed marked variation in patterns of both B. dendrobatidis infection and microbiome structure that are associated with host life history. Stream breeding amphibians were most likely to be infected with B. dendrobatidis. This increased probability of infection was further associated with increased abundance and diversity of non-Batrachochytrium chytrid fungi in the skin and environmental microbiome. We also show that increased alpha diversity and the relative abundance of fungi are lower in the skin microbiome of adult stream amphibians compared to adult pond-breeding amphibians, an association not seen for bacteria. Finally, stream tadpoles exhibit lower proportions of predicted protective microbial taxa than pond tadpoles, suggesting reduced biotic resistance. Our analyses show that host breeding ecology strongly shapes pathogen-microbiome associations at a landscape scale, a trait that may influence resilience in the face of emerging infectious diseases.

Enhanced genome assembly and a new official gene set for Tribolium castaneum.

BACKGROUND: The red flour beetle Tribolium castaneum has emerged as an important model organism for the study of gene function in development and physiology, for ecological and evolutionary genomics, for pest control and a plethora of other topics. RNA interference (RNAi), transgenesis and genome editing are well established and the resources for genome-wide RNAi screening have become available in this model. All these techniques depend on a high quality genome assembly and precise gene models. However, the first version of the genome assembly was generated by Sanger sequencing, and with a small set of RNA sequence data limiting annotation quality. RESULTS: Here, we present an improved genome assembly (Tcas5.2) and an enhanced genome annotation resulting in a new official gene set (OGS3) for Tribolium castaneum, which significantly increase the quality of the genomic resources. By adding large-distance jumping library DNA sequencing to join scaffolds and fill small gaps, the gaps in the genome assembly were reduced and the N50 increased to 4753kbp. The precision of the gene models was enhanced by the use of a large body of RNA-Seq reads of different life history stages and tissue types, leading to the discovery of 1452 novel gene sequences. We also added new features such as alternative splicing, well defined UTRs and microRNA target predictions. For quality control, 399 gene models were evaluated by manual inspection. The current gene set was submitted to Genbank and accepted as a RefSeq genome by NCBI. CONCLUSIONS: The new genome assembly (Tcas5.2) and the official gene set (OGS3) provide enhanced genomic resources for genetic work in Tribolium castaneum. The much improved information on transcription start sites supports transgenic and gene editing approaches. Further, novel types of information such as splice variants and microRNA target genes open additional possibilities for analysis.

Adaptive Radiation Genomics of Two Ecologically Divergent Hawai'ian Honeycreepers: The 'akiapola'au and the Hawai'i 'amakihi.

The Hawai'ian honeycreepers (drepanids) are a classic example of adaptive radiation: they adapted to a variety of novel dietary niches, evolving a wide range of bill morphologies. Here we investigated genomic diversity, demographic history, and genes involved in bill morphology phenotypes in 2 honeycreepers: the 'akiapola'au (Hemignathus wilsoni) and the Hawai'i 'amakihi (Chlorodrepanis virens). The 'akiapola'au is an endangered island endemic, filling the "woodpecker" niche by using a unique bill morphology, while the Hawai'i 'amakihi is a dietary generalist common on the islands of Hawai'i and Maui. We de novo sequenced the 'akiapola'au genome and compared it to the previously sequenced 'amakihi genome. The 'akiapola'au is far less heterozygous and has a smaller effective population size than the 'amakihi, which matches expectations due to its smaller census population and restricted ecological niche. Our investigation revealed genomic islands of divergence, which may be involved in the honeycreeper radiation. Within these islands of divergence, we identified candidate genes (including DLK1, FOXB1, KIF6, MAML3, PHF20, RBP1, and TIMM17A) that may play a role in honeycreeper adaptations. The gene DLK1, previously shown to influence Darwin's finch bill size, may be related to honeycreeper bill morphology evolution, while the functions of the other candidates remain unknown.

Elevated Prevalence of Azole-Resistant Aspergillus fumigatus in Urban versus Rural Environments in the United Kingdom.

Azole resistance in the opportunistic pathogen Aspergillus fumigatus is increasing, dominated primarily by the following two environmentally associated resistance alleles: TR34/L98H and TR46/Y121F/T289A. By sampling soils across the South of England, we assess the prevalence of azole-resistant A. fumigatus (ARAf) in samples collected in both urban and rural locations. We characterize the susceptibility profiles of the resistant isolates to three medical azoles, identify the underlying genetic basis of resistance, and investigate their genetic relationships. ARAf was detected in 6.7% of the soil samples, with a higher prevalence in urban (13.8%) than rural (1.1%) locations. Twenty isolates were confirmed to exhibit clinical breakpoints for resistance to at least one of three medical azoles, with 18 isolates exhibiting resistance to itraconazole, 6 to voriconazole, and 2 showing elevated minimum inhibitory concentrations to posaconazole. Thirteen of the resistant isolates harbored the TR34/L98H resistance allele, and six isolates carried the TR46/Y121F/T289A allele. The 20 azole-resistant isolates were spread across five csp1 genetic subtypes, t01, t02, t04B, t09, and t18 with t02 being the predominant subtype. Our study demonstrates that ARAf can be easily isolated in the South of England, especially in urban city centers, which appear to play an important role in the epidemiology of environmentally linked drug-resistant A. fumigatus.

MARDy: Mycology Antifungal Resistance Database.

Summary: The increase of antifungal drug resistance is a major global human health concern and threatens agriculture and food security; in order to tackle these concerns, it is important to understand the mechanisms that cause antifungal resistance. The curated Mycology Antifungal Resistance Database (MARDy) is a web-service of antifungal drug resistance mechanisms, including amino acid substitutions, tandem repeat sequences and genome ploidy. MARDy is implemented on a Linux, Apache, MySQL and PHP web development platform and includes a local installation of BLASTn of the database of curated genes. Availability and implementation: MARDy can be accessed at http://www.mardy.net and is free to use. The complete database can be retrieved, ordered by organism, gene and drug. Missing or new mycological antifungal resistance data can be relayed to the development team through a contribute entry form. Updates and news will be publicized via a dedicated Twitter feed: @MARDYfungi.

Membrane topology and identification of key residues of EaDAcT, a plant MBOAT with unusual substrate specificity.

Euonymus alatus diacylglycerol acetyltransferase (EaDAcT) catalyzes the transfer of an acetyl group from acetyl-CoA to the sn-3 position of diacylglycerol to form 3-acetyl-1,2-diacyl-sn-glycerol (acetyl-TAG). EaDAcT belongs to a small, plant-specific subfamily of the membrane bound O-acyltransferases (MBOAT) that acylate different lipid substrates. Sucrose gradient density centrifugation revealed that EaDAcT colocalizes to the same fractions as an endoplasmic reticulum (ER)-specific marker. By mapping the membrane topology of EaDAcT, we obtained an experimentally determined topology model for a plant MBOAT. The EaDAcT model contains four transmembrane domains (TMDs), with both the N- and C-termini orientated toward the lumen of the ER. In addition, there is a large cytoplasmic loop between the first and second TMDs, with the MBOAT signature region of the protein embedded in the third TMD close to the interface between the membrane and the cytoplasm. During topology mapping, we discovered two cysteine residues (C187 and C293) located on opposite sides of the membrane that are important for enzyme activity. In order to identify additional amino acid residues important for acetyltransferase activity, we isolated and characterized acetyltransferases from other acetyl-TAG-producing plants. Among them, the acetyltransferase from Euonymus fortunei possessed the highest activity in vivo and in vitro. Mutagenesis of conserved amino acids revealed that S253, H257, D258 and V263 are essential for EaDAcT activity. Alteration of residues unique to the acetyltransferases did not alter the unique acyl donor specificity of EaDAcT, suggesting that multiple amino acids are important for substrate recognition.

Hybrid de novo genome assembly and centromere characterization of the gray mouse lemur (Microcebus murinus).

BACKGROUND: The de novo assembly of repeat-rich mammalian genomes using only high-throughput short read sequencing data typically results in highly fragmented genome assemblies that limit downstream applications. Here, we present an iterative approach to hybrid de novo genome assembly that incorporates datasets stemming from multiple genomic technologies and methods. We used this approach to improve the gray mouse lemur (Microcebus murinus) genome from early draft status to a near chromosome-scale assembly. METHODS: We used a combination of advanced genomic technologies to iteratively resolve conflicts and super-scaffold the M. murinus genome. RESULTS: We improved the M. murinus genome assembly to a scaffold N50 of 93.32 Mb. Whole genome alignments between our primary super-scaffolds and 23 human chromosomes revealed patterns that are congruent with historical comparative cytogenetic data, thus demonstrating the accuracy of our de novo scaffolding approach and allowing assignment of scaffolds to M. murinus chromosomes. Moreover, we utilized our independent datasets to discover and characterize sequences associated with centromeres across the mouse lemur genome. Quality assessment of the final assembly found 96% of mouse lemur canonical transcripts nearly complete, comparable to other published high-quality reference genome assemblies. CONCLUSIONS: We describe a new assembly of the gray mouse lemur (Microcebus murinus) genome with chromosome-scale scaffolds produced using a hybrid bioinformatic and sequencing approach. The approach is cost effective and produces superior results based on metrics of contiguity and completeness. Our results show that emerging genomic technologies can be used in combination to characterize centromeres of non-model species and to produce accurate de novo chromosome-scale genome assemblies of complex mammalian genomes.

Inferential considerations for low-count RNA-seq transcripts: a case study on the dominant prairie grass Andropogon gerardii.

BACKGROUND: Differential expression (DE) analysis of RNA-seq data still poses inferential challenges, such as handling of transcripts characterized by low expression levels. In this study, we use a plasmode-based approach to assess the relative performance of alternative inferential strategies on RNA-seq transcripts, with special emphasis on transcripts characterized by a small number of read counts, so-called low-count transcripts, as motivated by an ecological application in prairie grasses. Big bluestem (Andropogon gerardii) is a wide-ranging dominant prairie grass of ecological and agricultural importance to the US Midwest while edaphic subspecies sand bluestem (A. gerardii ssp. Hallii) grows exclusively on sand dunes. Relative to big bluestem, sand bluestem exhibits qualitative phenotypic divergence consistent with enhanced drought tolerance, plausibly associated with transcripts of low expression levels. Our dataset consists of RNA-seq read counts for 25,582 transcripts (60% of which are classified as low-count) collected from leaf tissue of individual plants of big bluestem (n = 4) and sand bluestem (n = 4). Focused on low-count transcripts, we compare alternative ad-hoc data filtering techniques commonly used in RNA-seq pipelines and assess the inferential performance of recently developed statistical methods for DE analysis, namely DESeq2 and edgeR robust. These methods attempt to overcome the inherently noisy behavior of low-count transcripts by either shrinkage or differential weighting of observations, respectively. RESULTS: Both DE methods seemed to properly control family-wise type 1 error on low-count transcripts, whereas edgeR robust showed greater power and DESeq2 showed greater precision and accuracy. However, specification of the degree of freedom parameter under edgeR robust had a non-trivial impact on inference and should be handled carefully. When properly specified, both DE methods showed overall promising inferential performance on low-count transcripts, suggesting that ad-hoc data filtering steps at arbitrary expression thresholds may be unnecessary. A note of caution is in order regarding the approximate nature of DE tests under both methods. CONCLUSIONS: Practical recommendations for DE inference are provided when low-count RNA-seq transcripts are of interest, as is the case in the comparison of subspecies of bluestem grasses. Insights from this study may also be relevant to other applications focused on transcripts of low expression levels.

Tools and pipelines for BioNano data: molecule assembly pipeline and FASTA super scaffolding tool.

BACKGROUND: Genome assembly remains an unsolved problem. Assembly projects face a range of hurdles that confound assembly. Thus a variety of tools and approaches are needed to improve draft genomes. RESULTS: We used a custom assembly workflow to optimize consensus genome map assembly, resulting in an assembly equal to the estimated length of the Tribolium castaneum genome and with an N50 of more than 1 Mb. We used this map for super scaffolding the T. castaneum sequence assembly, more than tripling its N50 with the program Stitch. CONCLUSIONS: In this article we present software that leverages consensus genome maps assembled from extremely long single molecule maps to increase the contiguity of sequence assemblies. We report the results of applying these tools to validate and improve a 7x Sanger draft of the T. castaneum genome.

Genetic determinants of anti-malarial acquired immunity in a large multi-centre study.

BACKGROUND: Many studies report associations between human genetic factors and immunity to malaria but few have been reliably replicated. These studies are usually country-specific, use small sample sizes and are not directly comparable due to differences in methodologies. This study brings together samples and data collected from multiple sites across Africa and Asia to use standardized methods to look for consistent genetic effects on anti-malarial antibody levels. METHODS: Sera, DNA samples and clinical data were collected from 13,299 individuals from ten sites in Senegal, Mali, Burkina Faso, Sudan, Kenya, Tanzania, and Sri Lanka using standardized methods. DNA was extracted and typed for 202 Single Nucleotide Polymorphisms with known associations to malaria or antibody production, and antibody levels to four clinical grade malarial antigens [AMA1, MSP1, MSP2, and (NANP)4] plus total IgE were measured by ELISA techniques. Regression models were used to investigate the associations of clinical and genetic factors with antibody levels. RESULTS: Malaria infection increased levels of antibodies to malaria antigens and, as expected, stable predictors of anti-malarial antibody levels included age, seasonality, location, and ethnicity. Correlations between antibodies to blood-stage antigens AMA1, MSP1 and MSP2 were higher between themselves than with antibodies to the (NANP)4 epitope of the pre-erythrocytic circumsporozoite protein, while there was little or no correlation with total IgE levels. Individuals with sickle cell trait had significantly lower antibody levels to all blood-stage antigens, and recessive homozygotes for CD36 (rs321198) had significantly lower anti-malarial antibody levels to MSP2. CONCLUSION: Although the most significant finding with a consistent effect across sites was for sickle cell trait, its effect is likely to be via reducing a microscopically positive parasitaemia rather than directly on antibody levels. However, this study does demonstrate a framework for the feasibility of combining data from sites with heterogeneous malaria transmission levels across Africa and Asia with which to explore genetic effects on anti-malarial immunity.

Citizen science reveals landscape-scale exposures to multiazole-resistant Aspergillus fumigatus bioaerosols.

Using a citizen science approach, we identify a country-wide exposure to aerosolized spores of a human fungal pathogen, Aspergillus fumigatus, that has acquired resistance to the agricultural fungicide tebuconazole and first-line azole clinical antifungal drugs. Genomic analysis shows no distinction between resistant genotypes found in the environment and in patients, indicating that at least 40% of azole-resistant A. fumigatus infections are acquired from environmental exposures. Hotspots and coldspots of aerosolized azole-resistant spores were not stable between seasonal sampling periods. This suggests a high degree of atmospheric mixing resulting in an estimated per capita cumulative annual exposure of 21 days (±2.6). Because of the ubiquity of this measured exposure, it is imperative that we determine sources of azole-resistant A. fumigatus to reduce treatment failure in patients with aspergillosis.

Environment is associated with chytrid infection and skin microbiome richness on an amphibian rich island (Taiwan).

Growing evidence suggests that the origins of the panzootic amphibian pathogens Batrachochytrium dendrobatidis (Bd) and Batrachochytrium salamandrivorans (Bsal) are in Asia. In Taiwan, an island hotspot of high amphibian diversity, no amphibian mass mortality events linked to Bd or Bsal have been reported. We conducted a multi-year study across this subtropical island, sampling 2517 individuals from 30 species at 34 field sites, between 2010 and 2017, and including 171 museum samples collected between 1981 and 2009. We analyzed the skin microbiome of 153 samples (6 species) from 2017 in order to assess any association between the amphibian skin microbiome and the probability of infection amongst different host species. We did not detect Bsal in our samples, but found widespread infection by Bd across central and northern Taiwan, both taxonomically and spatially. Museum samples show that Bd has been present in Taiwan since at least 1990. Host species, geography (elevation), climatic conditions and microbial richness were all associated with the prevalence of infection. Host life-history traits, skin microbiome composition and phylogeny were associated with lower prevalence of infection for high altitude species. Overall, we observed low prevalence and burden of infection in host populations, suggesting that Bd is enzootic in Taiwan where it causes subclinical infections. While amphibian species in Taiwan are currently threatened by habitat loss, our study indicates that Bd is in an endemic equilibrium with the populations and species we investigated. However, ongoing surveillance of the infection is warranted, as changing environmental conditions may disturb the currently stable equilibrium.

Microbiome function predicts amphibian chytridiomycosis disease dynamics.

BACKGROUND: The fungal pathogen Batrachochytrium dendrobatidis (Bd) threatens amphibian biodiversity and ecosystem stability worldwide. Amphibian skin microbial community structure has been linked to the clinical outcome of Bd infections, yet its overall functional importance is poorly understood. METHODS: Microbiome taxonomic and functional profiles were assessed using high-throughput bacterial 16S rRNA and fungal ITS2 gene sequencing, bacterial shotgun metagenomics and skin mucosal metabolomics. We sampled 56 wild midwife toads (Alytes obstetricans) from montane populations exhibiting Bd epizootic or enzootic disease dynamics. In addition, to assess whether disease-specific microbiome profiles were linked to microbe-mediated protection or Bd-induced perturbation, we performed a laboratory Bd challenge experiment whereby 40 young adult A. obstetricans were exposed to Bd or a control sham infection. We measured temporal changes in the microbiome as well as functional profiles of Bd-exposed and control animals at peak infection. RESULTS: Microbiome community structure and function differed in wild populations based on infection history and in experimental control versus Bd-exposed animals. Bd exposure in the laboratory resulted in dynamic changes in microbiome community structure and functional differences, with infection clearance in all but one infected animal. Sphingobacterium, Stenotrophomonas and an unclassified Commamonadaceae were associated with wild epizootic dynamics and also had reduced abundance in laboratory Bd-exposed animals that cleared infection, indicating a negative association with Bd resistance. This was further supported by microbe-metabolite integration which identified functionally relevant taxa driving disease outcome, of which Sphingobacterium and Bd were most influential in wild epizootic dynamics. The strong correlation between microbial taxonomic community composition and skin metabolome in the laboratory and field is inconsistent with microbial functional redundancy, indicating that differences in microbial taxonomy drive functional variation. Shotgun metagenomic analyses support these findings, with similar disease-associated patterns in beta diversity. Analysis of differentially abundant bacterial genes and pathways indicated that bacterial environmental sensing and Bd resource competition are likely to be important in driving infection outcomes. CONCLUSIONS: Bd infection drives altered microbiome taxonomic and functional profiles across laboratory and field environments. Our application of multi-omics analyses in experimental and field settings robustly predicts Bd disease dynamics and identifies novel candidate biomarkers of infection. Video Abstract.

Somatic whole genome dynamics of precancer in Barrett's esophagus reveals features associated with disease progression.

While the genomes of normal tissues undergo dynamic changes over time, little is understood about the temporal-spatial dynamics of genomes in premalignant tissues that progress to cancer compared to those that remain cancer-free. Here we use whole genome sequencing to contrast genomic alterations in 427 longitudinal samples from 40 patients with stable Barrett's esophagus compared to 40 Barrett's patients who progressed to esophageal adenocarcinoma (ESAD). We show the same somatic mutational processes are active in Barrett's tissue regardless of outcome, with high levels of mutation, ESAD gene and focal chromosomal alterations, and similar mutational signatures. The critical distinction between stable Barrett's versus those who progress to cancer is acquisition and expansion of TP53-/- cell populations having complex structural variants and high-level amplifications, which are detectable up to six years prior to a cancer diagnosis. These findings reveal the timing of common somatic genome dynamics in stable Barrett's esophagus and define key genomic features specific to progression to esophageal adenocarcinoma, both of which are critical for cancer prevention and early detection strategies.

A Low-Cost Tebuconazole-Based Screening Test for Azole-Resistant Aspergillus fumigatus.

The global emergence of azole resistance in Aspergillus fumigatus is resulting in health and food security concerns. Rapid diagnostics and environmental surveillance methods are key to understanding the distribution and prevalence of azole resistance. However, such methods are often associated with high costs and are not always applicable to laboratories based in the least-developed countries. Here, we present and validate a low-cost screening protocol that can be used to differentiate between azole-susceptible "wild-type" and azole-resistant A. fumigatus isolates. © 2020 The Authors. Basic Protocol 1: Preparation of Tebucheck multi-well plates Basic Protocol 2: Inoculation of Tebucheck multi-well plates.

Cross-Disciplinary Genomics Approaches to Studying Emerging Fungal Infections.

Emerging fungal pathogens pose a serious, global and growing threat to food supply systems, wild ecosystems, and human health. However, historic chronic underinvestment in their research has resulted in a limited understanding of their epidemiology relative to bacterial and viral pathogens. Therefore, the untargeted nature of genomics and, more widely, -omics approaches is particularly attractive in addressing the threats posed by and illuminating the biology of these pathogens. Typically, research into plant, human and wildlife mycoses have been largely separated, with limited dialogue between disciplines. However, many serious mycoses facing the world today have common traits irrespective of host species, such as plastic genomes; wide host ranges; large population sizes and an ability to persist outside the host. These commonalities mean that -omics approaches that have been productively applied in one sphere and may also provide important insights in others, where these approaches may have historically been underutilised. In this review, we consider the advances made with genomics approaches in the fields of plant pathology, human medicine and wildlife health and the progress made in linking genomes to other -omics datatypes and sets; we identify the current barriers to linking -omics approaches and how these are being underutilised in each field; and we consider how and which -omics methodologies it is most crucial to build capacity for in the near future.

Deep whole-genome sequencing of 3 cancer cell lines on 2 sequencing platforms.

To test the performance of a new sequencing platform, develop an updated somatic calling pipeline and establish a reference for future benchmarking experiments, we performed whole-genome sequencing of 3 common cancer cell lines (COLO-829, HCC-1143 and HCC-1187) along with their matched normal cell lines to great sequencing depths (up to 278x coverage) on both Illumina HiSeqX and NovaSeq sequencing instruments. Somatic calling was generally consistent between the two platforms despite minor differences at the read level. We designed and implemented a novel pipeline for the analysis of tumor-normal samples, using multiple variant callers. We show that coupled with a high-confidence filtering strategy, the use of combination of tools improves the accuracy of somatic variant calling. We also demonstrate the utility of the dataset by creating an artificial purity ladder to evaluate the somatic pipeline and benchmark methods for estimating purity and ploidy from tumor-normal pairs. The data and results of the pipeline are made accessible to the cancer genomics community.

Captivity and Infection by the Fungal Pathogen Batrachochytrium salamandrivorans Perturb the Amphibian Skin Microbiome.

The emerging fungal pathogen, Batrachochytrium salamandrivorans (Bsal) is responsible for the catastrophic decline of European salamanders and poses a threat to amphibians globally. The amphibian skin microbiome can influence disease outcome for several host-pathogen systems, yet little is known of its role in Bsal infection. In addition, many experimental in-vivo amphibian disease studies to date have relied on specimens that have been kept in captivity for long periods without considering the influence of environment on the microbiome and how this may impact the host response to pathogen exposure. We characterized the impact of captivity and exposure to Bsal on the skin bacterial and fungal communities of two co-occurring European newt species, the smooth newt, Lissotriton vulgaris and the great-crested newt, Triturus cristatus. We show that captivity led to significant losses in bacterial and fungal diversity of amphibian skin, which may be indicative of a decline in microbe-mediated protection. We further demonstrate that in both L. vulgaris and T. cristatus, Bsal infection was associated with changes in the composition of skin bacterial communities with possible negative consequences to host health. Our findings advance current understanding of the role of host-associated microbiota in Bsal infection and highlight important considerations for ex-situ amphibian conservation programmes.