Childhood neuroblastoma has a remarkable variability in outcome. Age at diagnosis is one of the most important prognostic factors, with children less than 1 year old having favorable outcomes. Here we study single-cell and single-nuclei transcriptomes of neuroblastoma with different clinical risk groups and stages, including healthy adrenal gland. We compare tumor cell populations with embryonic mouse sympatho-adrenal derivatives, and post-natal human adrenal gland. We provide evidence that low and high-risk neuroblastoma have different cell identities, representing two disease entities. Low-risk neuroblastoma presents a transcriptome that resembles sympatho- and chromaffin cells, whereas malignant cells enriched in high-risk neuroblastoma resembles a subtype of TRKB+ cholinergic progenitor population identified in human post-natal gland. Analyses of these populations reveal different gene expression programs for worst and better survival in correlation with age at diagnosis. Our findings reveal two cellular identities and a composition of human neuroblastoma tumors reflecting clinical heterogeneity and outcome.

The petrous bone is considered the most efficient source of endogenous DNA across skeletal tissues in ancient DNA (aDNA) research as well as in forensic work. Recently, ancient DNA in auditory ossicle bones was shown to be comparably well-preserved as in the petrous, although no attempt was made to distinguish among the three ossicle bones. In this study, we compare aDNA profiles across matched ossicle- and petrous-derived sequencing libraries prepared from 29 human skeletons from Neolithic Anatolia and Medieval Iberia. We find that the stapes and incus provide higher human endogenous aDNA than the petrous bone, with >2×; higher median rates of endogenous aDNA recovery, while the malleus performs similarly to the petrous. Human aDNA fragments retrieved from the stapes were 8% longer than those from the petrous, whereas postmortem damage, clonality and contamination rates were comparable among the studied bone types. These observations are corroborated by data from non-matched ossicle or petrous libraries from 81 individuals from the same contexts, with the highest endogenous aDNA content observed in the stapes. Despite being the smallest bone in the human skeleton, the stapes, along with the incus, may be among the most optimal aDNA sources yet identified.

Reduced fitness through genetic drift and inbreeding is a major threat to small and isolated populations. Although previous studies have generally used genetically verified pedigrees to document effects of inbreeding and gene flow, these often fail to capture the whole inbreeding history of the species. By assembling a draft arctic fox (Vulpes lagopus) genome and resequencing complete genomes of 23 additional foxes born before and after a well-documented immigration event in Scandinavia, we here look into the genomic consequences of inbreeding and genetic rescue. We found a difference in genome-wide diversity, with 18% higher heterozygosity and 81% lower FROH in immigrant F1 compared to native individuals. However, more distant descendants of immigrants (F2, F3) did not show the same pattern. We also found that foxes with lower inbreeding had higher probability to survive their first year of life. Our results demonstrate the important link between genetic variation and fitness as well as the transient nature of genetic rescue. Moreover, our results have implications in conservation biology as they demonstrate that inbreeding depression can effectively be detected in the wild by a genomic approach.

The structure of 5' untranslated regions (5' UTRs) of bacterial mRNAs often determines the fate of the transcripts. Using a dimethyl sulfate mutational profiling with sequencing (DMS-MaPseq) approach, we developed a protocol to generate sequence libraries to determine the base-pairing status of adenines and cytosines in the 5' UTRs of bacterial mRNAs. Our method increases the sequencing depth of the 5' UTRs and allows detection of changes in their structures by sequencing libraries of moderate sizes. For complete details on the use and execution of this protocol, please refer to Ignatov et al. (2020).

Craniosynostosis (CS), the premature closure of one or more cranial sutures, occurs both as part of a syndrome or in isolation (nonsyndromic form). Here, we have studied the prevalence and spectrum of genetic alterations associated with coronal suture closure in 100 Scandinavian patients treated at a single craniofacial unit. All patients were phenotypically assessed and analyzed with a custom-designed 63 gene NGS-panel. Most cases (78%) were syndromic forms of CS. Pathogenic and likely pathogenic variants explaining the phenotype were found in 80% of the families with syndromic CS and in 14% of those with nonsyndromic CS. Sixty-five percent of the families had mutations in the CS core genes FGFR2, TWIST1, FGFR3, TCF12, EFNB1, FGFR1, and POR. Five novel pathogenic/likely pathogenic variants in TWIST1, TCF12, and EFNB1 were identified. We also found novel variants in SPECC1L, IGF1R, and CYP26B1 with a possible modulator phenotypic effect. Our findings demonstrate that NGS targeted sequencing is a powerful tool to detect pathogenic mutations in patients with coronal CS and further emphasize the importance of thorough assessment of the patient's phenotype for reliable interpretation of the molecular findings. This is particularly important in patients with complex phenotypes and rare forms of CS.

An important part of the Human Proteome Project is to characterize the protein complement of the genome with antibody-based profiling. Within the framework of this effort, a new version 12 of the Human Protein Atlas ( www.proteinatlas.org ) has been launched, including transcriptomics data for 27 tissues and 44 cell lines to complement the protein expression data from antibody-based profiling. Besides the extensive addition of transcriptomics data, the Human Protein Atlas now contains antibody-based protein profiles for 82% of the 20 329 putative protein-coding genes. The comprehensive data resulting from RNA-seq analysis and antibody-based profiling performed within the Human Protein Atlas as well as information from UniProt were used to generate evidence summary scores for each of the 20 329 genes, of which 94% now have experimental evidence at least at transcript level. The evidence scores for all individual genes are displayed with regards to both RNA- and antibody-based protein profiles, including chromosome-centric visualizations. An analysis of the human chromosome 19 shows that ∼43% of the genes are expressed at the transcript level in all 27 tissues analyzed, suggesting a "house-keeping" function, while 12% of the genes show a more tissue-specific pattern with enriched expression in one of the analyzed tissues only.

Morabine grasshoppers in the Vandiemenella viatica species group, which show karyotype diversity, have been studied for their ecological distribution and speciation in relation to their genetic and chromosomal diversity. They are good models for studying sex chromosome evolution as "old" and newly emerged sex chromosomes co-exist within the group. Here we present a reference genome for the viatica19 chromosomal race, that possesses the ancestral karyotype within the group. Using PacBio HiFi and Hi-C sequencing, we generated a chromosome-level assembly of 4.09 Gb in span, scaffold N50 of 429 Mb, and complete BUSCO score of 98.1%, containing 10 pseudo-chromosomes. We provide Illumina datasets of males and females, used to identify the X chromosome. The assembly contains 19,034 predicted protein-coding genes, and a total of 75.21% of repetitive DNA sequences. By leveraging HiFi reads, we mapped the genome-wide distribution of methylated bases (5mC and 6 mA). This comprehensive assembly offers a robust reference for morabine grasshoppers and supports further research into speciation and sex chromosome diversification within the group and its related species.