Towards complete and error-free genome assemblies of all vertebrate species.

High-quality and complete reference genome assemblies are fundamental for the application of genomics to biology, disease, and biodiversity conservation. However, such assemblies are available for only a few non-microbial species1-4. To address this issue, the international Genome 10K (G10K) consortium5,6 has worked over a five-year period to evaluate and develop cost-effective methods for assembling highly accurate and nearly complete reference genomes. Here we present lessons learned from generating assemblies for 16 species that represent six major vertebrate lineages. We confirm that long-read sequencing technologies are essential for maximizing genome quality, and that unresolved complex repeats and haplotype heterozygosity are major sources of assembly error when not handled correctly. Our assemblies correct substantial errors, add missing sequence in some of the best historical reference genomes, and reveal biological discoveries. These include the identification of many false gene duplications, increases in gene sizes, chromosome rearrangements that are specific to lineages, a repeated independent chromosome breakpoint in bat genomes, and a canonical GC-rich pattern in protein-coding genes and their regulatory regions. Adopting these lessons, we have embarked on the Vertebrate Genomes Project (VGP), an international effort to generate high-quality, complete reference genomes for all of the roughly 70,000 extant vertebrate species and to help to enable a new era of discovery across the life sciences.

Dermonecrosis caused by a spitting cobra snakebite results from toxin potentiation and is prevented by the repurposed drug varespladib.

Snakebite envenoming is a neglected tropical disease that causes substantial mortality and morbidity globally. The venom of African spitting cobras often causes permanent injury via tissue-destructive dermonecrosis at the bite site, which is ineffectively treated by current antivenoms. To address this therapeutic gap, we identified the etiological venom toxins in Naja nigricollis venom responsible for causing local dermonecrosis. While cytotoxic three-finger toxins were primarily responsible for causing spitting cobra cytotoxicity in cultured keratinocytes, their potentiation by phospholipases A2 toxins was essential to cause dermonecrosis in vivo. This evidence of probable toxin synergism suggests that a single toxin-family inhibiting drug could prevent local envenoming. We show that local injection with the repurposed phospholipase A2-inhibiting drug varespladib significantly prevents local tissue damage caused by several spitting cobra venoms in murine models of envenoming. Our findings therefore provide a therapeutic strategy that may effectively prevent life-changing morbidity caused by snakebite in rural Africa.

FAIR Header Reference genome: a TRUSTworthy standard.

The lack of interoperable data standards among reference genome data-sharing platforms inhibits cross-platform analysis while increasing the risk of data provenance loss. Here, we describe the FAIR bioHeaders Reference genome (FHR), a metadata standard guided by the principles of Findability, Accessibility, Interoperability and Reuse (FAIR) in addition to the principles of Transparency, Responsibility, User focus, Sustainability and Technology. The objective of FHR is to provide an extensive set of data serialisation methods and minimum data field requirements while still maintaining extensibility, flexibility and expressivity in an increasingly decentralised genomic data ecosystem. The effort needed to implement FHR is low; FHR's design philosophy ensures easy implementation while retaining the benefits gained from recording both machine and human-readable provenance.

Hydrogen energy futures - foraging or farming?

Exploration for commercially viable natural hydrogen accumulations within the Earth's crust, here compared to 'foraging' for wild food, holds promise. However, a potentially more effective strategy lies in the in situ artificial generation of hydrogen in natural underground reservoirs, akin to 'farming'. Both biotic and abiotic processes can be employed, converting introduced or indigenous components, gases, and nutrients into hydrogen. Through studying natural hydrogen-generating reactions, we can discern pathways for optimized engineering. Some reactions may be inherently slow, allowing for a 'seed and leave' methodology, where sites are infused with gases, nutrients, and specific bacterial strains, then left to gradually produce hydrogen. However, other reactions could offer quicker outcomes to harvest hydrogen. A crucial element of this strategy is our innovative concept of 'X' components-ranging from trace minerals to bioengineered microbes. These designed components enhance biotic and/or abiotic reactions and prove vital in accelerating hydrogen production. Drawing parallels with our ancestors' transition from hunter-gathering to agriculture, we propose a similar paradigm shift in the pursuit of hydrogen energy. As we transition towards a hydrogen-centric energy landscape, the amalgamation of geochemistry, advanced biology, and engineering emerges as a beacon, signalling a pathway towards a sustainable and transformative energy future.

Caecilian Genomes Reveal the Molecular Basis of Adaptation and Convergent Evolution of Limblessness in Snakes and Caecilians.

We present genome sequences for the caecilians Geotrypetes seraphini (3.8 Gb) and Microcaecilia unicolor (4.7 Gb), representatives of a limbless, mostly soil-dwelling amphibian clade with reduced eyes, and unique putatively chemosensory tentacles. More than 69% of both genomes are composed of repeats, with retrotransposons being the most abundant. We identify 1,150 orthogroups that are unique to caecilians and enriched for functions in olfaction and detection of chemical signals. There are 379 orthogroups with signatures of positive selection on caecilian lineages with roles in organ development and morphogenesis, sensory perception, and immunity amongst others. We discover that caecilian genomes are missing the zone of polarizing activity regulatorysequence (ZRS) enhancer of Sonic Hedgehog which is also mutated in snakes. In vivo deletions have shown ZRS is required for limb development in mice, thus, revealing a shared molecular target implicated in the independent evolution of limblessness in snakes and caecilians.

Brassica rapa CURLY LEAF is a major H3K27 methyltransferase regulating flowering time.

MAIN CONCLUSION: In Brassica rapa, the epigenetic modifier BraA.CLF orchestrates flowering by modulating H3K27me3 levels at the floral integrator genes FT, SOC1, and SEP3, thereby influencing their expression. CURLY LEAF (CLF) is the catalytic subunit of the plant Polycomb Repressive Complex 2 that mediates the trimethylation of histone H3 lysine 27 (H3K27me3), an epigenetic modification that leads to gene silencing. While the function of CURLY LEAF (CLF) has been extensively studied in Arabidopsis thaliana, its role in Brassica crops is barely known. In this study, we focused on the Brassica rapa homolog of CLF and found that the loss-of-function mutant braA.clf-1 exhibits an accelerated flowering together with pleiotropic phenotypic alterations compared to wild-type plants. In addition, we carried out transcriptomic and H3K27me3 genome-wide analyses to identify the genes regulated by BraA.CLF. Interestingly, we observed that several floral regulatory genes, including the B. rapa homologs of FT, SOC1 and SEP3, show reduced H3K27me3 levels and increased transcript levels compared to wild-type plants, suggesting that they are direct targets of BraA.CLF and key players in regulating flowering time in this crop. In addition, the results obtained will enhance our understanding of the epigenetic mechanisms regulating key developmental traits and will aid to increase crop yield by engineering new Brassica varieties with different flowering time requirements.

A Retrospective Longitudinal Study of 460 Patients with ABCA4-Associated Retinal Disease.

PURPOSE: To investigate the distribution of genotypes and natural history of ABCA4-associated retinal disease in a large cohort of patients seen at a single institution. DESIGN: Retrospective, single-institution cohort review. PARTICIPANTS: Patients seen at the University of Iowa between November 1986 and August 2022 clinically suspected to have disease caused by sequence variations in ABCA4. METHODS: DNA samples from participants were subjected to a tiered testing strategy progressing from allele-specific screening to whole genome sequencing. Charts were reviewed, and clinical data were tabulated. The pathogenic severity of the most common alleles was estimated by studying groups of patients who shared 1 allele. Groups of patients with shared genotypes were reviewed for evidence of modifying factor effects. MAIN OUTCOME MEASURES: Age at first uncorrectable vision loss, best-corrected visual acuity, and the area of the I2e isopter of the Goldmann visual field. RESULTS: A total of 460 patients from 390 families demonstrated convincing clinical features of ABCA4-associated retinal disease. Complete genotypes were identified in 399 patients, and partial genotypes were identified in 61. The median age at first vision loss was 16 years (range, 4-76 years). Two hundred sixty-five families (68%) harbored a unique genotype, and no more than 10 patients shared any single genotype. Review of the patients with shared genotypes revealed evidence of modifying factors that in several cases resulted in a > 15-year difference in age at first vision loss. Two hundred forty-one different alleles were identified among the members of this cohort, and 161 of these (67%) were found in only a single individual. CONCLUSIONS: ABCA4-associated retinal disease ranges from a very severe photoreceptor disease with an onset before 5 years of age to a late-onset retinal pigment epithelium-based condition resembling pattern dystrophy. Modifying factors frequently impact the ABCA4 disease phenotype to a degree that is similar in magnitude to the detectable ABCA4 alleles themselves. It is likely that most patients in any cohort will harbor a unique genotype. The latter observations taken together suggest that patients' clinical findings in most cases will be more useful for predicting their clinical course than their genotype. FINANCIAL DISCLOSURE(S): Proprietary or commercial disclosure may be found in the Footnotes and Disclosures at the end of this article.

Genomic SNPs resolve the phylogeny of an ancient amphibian island radiation from the Seychelles.

Unusually for oceanic islands, the granitic Seychelles host multiple lineages of endemic amphibians. This includes an ancient (likely ca. 60 million years) radiation of eight caecilian species, most of which occur on multiple islands.These caecilians have a complicated taxonomic history and their phylogenetic inter-species relationships have been difficult to resolve. Double-digest RAD sequencing (ddRADseq) has been applied extensively to phylogeography and increasingly to phylogenetics but its utility for resolving ancient divergences is less well established. To address this, we applied ddRADseq to generate a genome-wide SNP panel for phylogenomic analyses of the Seychelles caecilians, whose phylogeny has so far not been satisfactorily resolved with traditional DNA markers. Based on 129,154 SNPs, we resolved deep and shallow splits, with strong support. Our findings demonstrate the capability of genome-wide SNPs for evolutionary inference at multiple taxonomic levels and support the recently proposed synonymy of Grandisonia Taylor, 1968 with Hypogeophis Peters, 1879. We revealed three clades of Hypogeophis (large-, medium- and short-bodied) and identify a single origin of the diminutive, stocky-bodied and pointy-snouted phenotype.

Comparative morphology of oral glands in snakes of the family Homalopsidae reveals substantial variation and additional independent origins of salt glands within Serpentes.

Using diffusible iodine-based contrast-enhanced computed tomography (diceCT), we examined the morphology of the oral glands of 12 species of the family Homalopsidae. Snakes of this family exhibit substantial interspecific morphological variation in their oral glands. Particular variables are the venom glands, ranging from large (e.g., Subsessor bocourti) to small (e.g., Erpeton tentaculatum). The supra- and infralabial glands are more uniform in morphology, being the second most developed in almost all the sampled species. Premaxillary glands distinct from the supralabial glands were observed in five species (Myron richardsonii, Bitia hydroides, Cantoria violacea, Fordonia leucobalia, and Gerarda prevostiana), in addition to Cerberus rynchops, the only species in which this condition was previously documented associated with the excretion of salt. In the three species of the saltwater group of homalopsids (C. violacea, F. leucobalia, and G. prevostiana), the premaxillary glands also extend posteriorly, occupying a large area above the supralabial gland, a condition not observed in any other species of snake studied thus far. Character evolution analyses indicate that premaxillary glands differentiated from the supralabial gland and evolved independently three or four times in the family, always in lineages that invaded marine habitats. Our results suggest that the differentiated premaxillary glands are likely salt glands, as is the case in C. rynchops. If corroborated, this increases to six or seven the number of independent evolutionary origins of salt glands in snakes that have undergone an evolutionary transition to marine life.

Pan-cancer genome and transcriptome analyses of 1,699 paediatric leukaemias and solid tumours.

Analysis of molecular aberrations across multiple cancer types, known as pan-cancer analysis, identifies commonalities and differences in key biological processes that are dysregulated in cancer cells from diverse lineages. Pan-cancer analyses have been performed for adult but not paediatric cancers, which commonly occur in developing mesodermic rather than adult epithelial tissues. Here we present a pan-cancer study of somatic alterations, including single nucleotide variants, small insertions or deletions, structural variations, copy number alterations, gene fusions and internal tandem duplications in 1,699 paediatric leukaemias and solid tumours across six histotypes, with whole-genome, whole-exome and transcriptome sequencing data processed under a uniform analytical framework. We report 142 driver genes in paediatric cancers, of which only 45% match those found in adult pan-cancer studies; copy number alterations and structural variants constituted the majority (62%) of events. Eleven genome-wide mutational signatures were identified, including one attributed to ultraviolet-light exposure in eight aneuploid leukaemias. Transcription of the mutant allele was detectable for 34% of protein-coding mutations, and 20% exhibited allele-specific expression. These data provide a comprehensive genomic architecture for paediatric cancers and emphasize the need for paediatric cancer-specific development of precision therapies.

The Mechano-Ubiquitinome of Articular Cartilage: Differential Ubiquitination and Activation of a Group of ER-Associated DUBs and ER Stress Regulators.

Understanding how connective tissue cells respond to mechanical stimulation is important to human health and disease processes in musculoskeletal diseases. Injury to articular cartilage is a key risk factor in predisposition to tissue damage and degenerative osteoarthritis. Recently, we have discovered that mechanical injury to connective tissues including murine and porcine articular cartilage causes a significant increase in lysine-63 polyubiquitination. Here, we identified the ubiquitin signature that is unique to injured articular cartilage tissue upon mechanical injury (the "mechano-ubiquitinome"). A total of 463 ubiquitinated peptides were identified, with an enrichment of ubiquitinated peptides of proteins involved in protein processing in the endoplasmic reticulum (ER), also known as the ER-associated degradation response, including YOD1, BRCC3, ATXN3, and USP5 as well as the ER stress regulators, RAD23B, VCP/p97, and Ubiquilin 1. Enrichment of these proteins suggested an injury-induced ER stress response and, for instance, ER stress markers DDIT3/CHOP and BIP/GRP78 were upregulated following cartilage injury on the protein and gene expression levels. Similar ER stress induction was also observed in response to tail fin injury in zebrafish larvae, suggesting a generic response to tissue injury. Furthermore, a rapid increase in global DUB activity following injury and significant activity in human osteoarthritic cartilage was observed using DUB-specific activity probes. Combined, these results implicate the involvement of ubiquitination events and activation of a set of DUBs and ER stress regulators in cellular responses to cartilage tissue injury and in osteoarthritic cartilage tissues. This link through the ER-associated degradation pathway makes this protein set attractive for further investigation in in vivo models of tissue injury and for targeting in osteoarthritis and related musculoskeletal diseases.

Differential diagnosis of contact dermatitis: A practical-approach review by the EADV Task Force on contact dermatitis.

The diagnosis of eczema ('dermatitis') is mostly clinical and depends on the clinical history and exploratory objective findings (primary lesions, patterns). Contact dermatitis remains as an important condition in the group of eczematous disorders, with important socioeconomic and occupational relevance. Although irritant and allergic contact dermatitis have a different pathogenesis, both are characterized by a rather typical morphology, are triggered by external factors and tend to occur primarily in the area of contact with the exogenous agent. In addition, allergic and irritant dermatitis may also co-exist. The importance of diagnosing contact dermatitis, especially when allergic in nature, is both due to the possibility of avoiding the trigger, and due to its role in aggravating other skin conditions. Nevertheless, the heterogeneity of clinical presentations in daily practice may pose an important challenge for the suspicion and correct diagnosis of contact dermatitis. Furthermore, other conditions, with different pathogenesis and treatment, may clinically simulate contact dermatitis. The Task Force aims to conduct a review of the unifying clinical features of contact dermatitis and characterize its main clinical phenotypes, and its simulators, in order to contribute to an early suspicion or recognition of contact dermatitis and enable a correct differential diagnosis.

Severely compromised supply of patch test allergens in Europe hampers adequate diagnosis of occupational and non-occupational contact allergy. A European Society of Contact Dermatitis (ESCD), European Academy of Allergy and Clinical Immunology (EAACI), European Academy of Dermatology and Venereology (EADV) task forces 'Contact Dermatitis' and 'Occupational Skin Disease' position paper.

Patch testing is the only clinically applicable diagnostic method for Type IV allergy. The availability of Type IV patch test (PT) allergens in Europe, however, is currently scarce. This severely compromises adequate diagnostics of contact allergy, leading to serious consequences for the affected patients. Against this background, the European Society of Contact Dermatitis (ESCD) has created a task force (TF) (i) to explore the current availability of PT substances in different member states, (ii) to highlight some of the unique characteristics of Type IV vs. other allergens and (iii) to suggest ways forward to promote and ensure availability of high-quality patch testing substances for the diagnosis of Type IV allergies throughout Europe. The suggestions of the TF on how to improve the availability of PT allergens are supported by the ESCD, the European Academy of Allergy and Clinical Immunology, and the European Academy of Dermatology and Venereology and intend to provide potential means to resolve the present medical crisis.

The anatomy of the head muscles in caecilians (Amphibia: Gymnophiona): Variation in relation to phylogeny and ecology?

In limbless fossorial vertebrates such as caecilians (Gymnophiona), head-first burrowing imposes severe constraints on the morphology and overall size of the head. As such, caecilians developed a unique jaw-closing system involving the large and well-developed m. interhyoideus posterior, which is positioned in such a way that it does not significantly increase head diameter. Caecilians also possess unique muscles among amphibians. Understanding the diversity in the architecture and size of the cranial muscles may provide insights into how a typical amphibian system was adapted for a head-first burrowing lifestyle. In this study, we use dissection and non-destructive contrast-enhanced micro-computed tomography (µCT) scanning to describe and compare the cranial musculature of 13 species of caecilians. Our results show that the general organization of the head musculature is rather constant across extant caecilians. However, the early-diverging Rhinatrema bivittatum mainly relies on the 'ancestral' amphibian jaw-closing mechanism dominated by the m. adductores mandibulae, whereas other caecilians switched to the use of the derived dual jaw-closing mechanism involving the additional recruitment of the m. interhyoideus posterior. Additionally, the aquatic Typhlonectes show a greater investment in hyoid musculature than terrestrial caecilians, which is likely related to greater demands for ventilating their large lungs, and perhaps also an increased use of suction feeding. In addition to three-dimensional interactive models, our study provides the required quantitative data to permit the generation of accurate biomechanical models allowing the testing of further functional hypotheses.

Recommendation to update the British Society for Cutaneous Allergy corticosteroid series.

BACKGROUND: Patch testing is an important investigation when dermatitis is unresponsive to, or worsened by, topical corticosteroid treatment. There is a balance to be struck between testing too many allergens, which is expensive, time consuming and risks causing sensitization, and testing too few, which risks missing the diagnosis. The current British Society for Cutaneous Allergy (BSCA) corticosteroid series comprises eight allergens and was last updated in February 2007. AIM: To review and update the BSCA corticosteroid series. METHODS: We retrospectively analysed data from 16 patch test centres in the UK and Ireland for all patients who were patch tested to a corticosteroid series between August 2017 and July 2019. We recorded the allergens tested, the number and percentage tested to a corticosteroid series and the number of positive results for each allergen. We identified the allergens that test positive in ≥ 0.1% of selectively tested patients. RESULTS: Overall, 3531 patients were tested to a corticosteroid series in the 16 centres. The number of allergens tested ranged from 7 to 18 (mean 10). The proportion of patch test patients who were tested to a corticosteroid series ranged from 1% to 99%. Six allergens in the 2017 BSCA series tested positive in ≥ 0.1% of patients. Nine allergens not in the BSCA corticosteroid series tested positive in ≥ 0.1% of patients. CONCLUSION: This audit demonstrates the importance of regular review of recommended series and the significant variations in practice. The new BSCA corticosteroid series that we recommend contains 13 haptens, with the addition of the patient's own steroid creams as appropriate.

Occupational contact allergy: The European perspective-Analysis of patch test data from ESSCA between 2011 and 2020.

BACKGROUND: Occupational skin diseases have led the occupational disease statistics in Europe for many years. Especially occupational allergic contact dermatitis is associated with a poor prognosis and low healing rates leading to an enormous burden for the affected individual and for society. OBJECTIVES: To present the sensitization frequencies to the most relevant allergens of the European baseline series in patients with occupational contact dermatitis (OCD) and to compare sensitization profiles of different occupations. METHODS: The data of 16 022 patients considered having OCD after patch testing within the European Surveillance System on Contact Allergies (ESSCA) network between January 2011 and December 2020 were evaluated. Patients (n = 46 652) in whom an occupational causation was refuted served as comparison group. RESULTS: The highest percentages of OCD were found among patients working in agriculture, fishery and related workers, metal industry, chemical industry, followed by building and construction industry, health care, food and service industry. Sensitizations to rubber chemicals (thiurams, carbamates, benzothiazoles) and epoxy resins were associated with at least a doubled risk of OCD. After a decline from 2014 onwards, the risks to acquire an occupation-related sensitization to methyl(chloro)isothiazolinone (MCI/MI) and especially to methylisothiazolinone (MI) seem to increase again. Sensitization rates to formaldehyde were stable, and to methyldibromo glutaronitrile (MDBGN) slightly decreasing over time. CONCLUSIONS: Among allergens in the European Baseline Series, occupational relevance is most frequently attributed to rubber accelerators, epoxy resins and preservatives.

Physiological constraints dictate toxin spatial heterogeneity in snake venom glands.

BACKGROUND: Venoms are ecological innovations that have evolved numerous times, on each occasion accompanied by the co-evolution of specialised morphological and behavioural characters for venom production and delivery. The close evolutionary interdependence between these characters is exemplified by animals that control the composition of their secreted venom. This ability depends in part on the production of different toxins in different locations of the venom gland, which was recently documented in venomous snakes. Here, we test the hypothesis that the distinct spatial distributions of toxins in snake venom glands are an adaptation that enables the secretion of venoms with distinct ecological functions. RESULTS: We show that the main defensive and predatory peptide toxins are produced in distinct regions of the venom glands of the black-necked spitting cobra (Naja nigricollis), but these distributions likely reflect developmental effects. Indeed, we detected no significant differences in venom collected via defensive 'spitting' or predatory 'biting' events from the same specimens representing multiple lineages of spitting cobra. We also found the same spatial distribution of toxins in a non-spitting cobra and show that heterogeneous toxin distribution is a feature shared with a viper with primarily predatory venom. CONCLUSIONS: Our findings suggest that heterogeneous distributions of toxins are not an adaptation to controlling venom composition in snakes. Instead, it likely reflects physiological constraints on toxin production by the venom glands, opening avenues for future research on the mechanisms of functional differentiation of populations of protein-secreting cells within adaptive contexts.

Resistance Is Not Futile: Widespread Convergent Evolution of Resistance to Alpha-Neurotoxic Snake Venoms in Caecilians (Amphibia: Gymnophiona).

Predatory innovations impose reciprocal selection pressures upon prey. The evolution of snake venom alpha-neurotoxins has triggered the corresponding evolution of resistance in the post-synaptic nicotinic acetylcholine receptors of prey in a complex chemical arms race. All other things being equal, animals like caecilians (an Order of legless amphibians) are quite vulnerable to predation by fossorial elapid snakes and their powerful alpha-neurotoxic venoms; thus, they are under strong selective pressure. Here, we sequenced the nicotinic acetylcholine receptor alpha-1 subunit of 37 caecilian species, representing all currently known families of caecilians from across the Americas, Africa, and Asia, including species endemic to the Seychelles. Three types of resistance were identified: (1) steric hindrance from N-glycosylated asparagines; (2) secondary structural changes due to the replacement of proline by another amino acid; and (3) electrostatic charge repulsion of the positively charged neurotoxins, through the introduction of a positively charged amino acid into the toxin-binding site. We demonstrated that resistance to alpha-neurotoxins convergently evolved at least fifteen times across the caecilian tree (three times in Africa, seven times in the Americas, and five times in Asia). Additionally, as several species were shown to possess multiple resistance modifications acting synergistically, caecilians must have undergone at least 20 separate events involving the origin of toxin resistance. On the other hand, resistance in non-caecilian amphibians was found to be limited to five origins. Together, the mutations underlying resistance in caecilians constitute a robust signature of positive selection which strongly correlates with elapid presence through both space (sympatry with caecilian-eating elapids) and time (Cenozoic radiation of elapids). Our study demonstrates the extent of convergent evolution that can be expected when a single widespread predatory adaptation triggers parallel evolutionary arms races at a global scale.

Pediatric Cancer Variant Pathogenicity Information Exchange (PeCanPIE): a cloud-based platform for curating and classifying germline variants.

Variant interpretation in the era of massively parallel sequencing is challenging. Although many resources and guidelines are available to assist with this task, few integrated end-to-end tools exist. Here, we present the Pediatric Cancer Variant Pathogenicity Information Exchange (PeCanPIE), a web- and cloud-based platform for annotation, identification, and classification of variations in known or putative disease genes. Starting from a set of variants in variant call format (VCF), variants are annotated, ranked by putative pathogenicity, and presented for formal classification using a decision-support interface based on published guidelines from the American College of Medical Genetics and Genomics (ACMG). The system can accept files containing millions of variants and handle single-nucleotide variants (SNVs), simple insertions/deletions (indels), multiple-nucleotide variants (MNVs), and complex substitutions. PeCanPIE has been applied to classify variant pathogenicity in cancer predisposition genes in two large-scale investigations involving >4000 pediatric cancer patients and serves as a repository for the expert-reviewed results. PeCanPIE was originally developed for pediatric cancer but can be easily extended for use for nonpediatric cancers and noncancer genetic diseases. Although PeCanPIE's web-based interface was designed to be accessible to non-bioinformaticians, its back-end pipelines may also be run independently on the cloud, facilitating direct integration and broader adoption. PeCanPIE is publicly available and free for research use.

Predicting microbiomes through a deep latent space.

MOTIVATION: Microbial communities influence their environment by modifying the availability of compounds, such as nutrients or chemical elicitors. Knowing the microbial composition of a site is therefore relevant to improve productivity or health. However, sequencing facilities are not always available, or may be prohibitively expensive in some cases. Thus, it would be desirable to computationally predict the microbial composition from more accessible, easily-measured features. RESULTS: Integrating deep learning techniques with microbiome data, we propose an artificial neural network architecture based on heterogeneous autoencoders to condense the long vector of microbial abundance values into a deep latent space representation. Then, we design a model to predict the deep latent space and, consequently, to predict the complete microbial composition using environmental features as input. The performance of our system is examined using the rhizosphere microbiome of Maize. We reconstruct the microbial composition (717 taxa) from the deep latent space (10 values) with high fidelity (>0.9 Pearson correlation). We then successfully predict microbial composition from environmental variables, such as plant age, temperature or precipitation (0.73 Pearson correlation, 0.42 Bray-Curtis). We extend this to predict microbiome composition under hypothetical scenarios, such as future climate change conditions. Finally, via transfer learning, we predict microbial composition in a distinct scenario with only 100 sequences, and distinct environmental features. We propose that our deep latent space may assist microbiome-engineering strategies when technical or financial resources are limited, through predicting current or future microbiome compositions. AVAILABILITY AND IMPLEMENTATION: Software, results and data are available at https://github.com/jorgemf/DeepLatentMicrobiome. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.