Fully automated deep learning-based localization and segmentation of the locus coeruleus in aging and Parkinson's disease using neuromelanin-sensitive MRI.

PURPOSE: Development and performance measurement of a fully automated pipeline that localizes and segments the locus coeruleus in so-called neuromelanin-sensitive magnetic resonance imaging data for the derivation of quantitative biomarkers of neurodegenerative diseases such as Alzheimer's disease and Parkinson's disease. METHODS: We propose a pipeline composed of several 3D-Unet-based convolutional neural networks for iterative multi-scale localization and multi-rater segmentation and non-deep learning-based components for automated biomarker extraction. We trained on the healthy aging cohort and did not carry out any adaption or fine-tuning prior to the application to Parkinson's disease subjects. RESULTS: The localization and segmentation pipeline demonstrated sufficient performance as measured by Euclidean distance (on average around 1.3mm on healthy aging subjects and 2.2mm in Parkinson's disease subjects) and Dice similarity coefficient (overall around [Formula: see text] on healthy aging subjects and [Formula: see text] for subjects with Parkinson's disease) as well as promising agreement with respect to contrast ratios in terms of intraclass correlation coefficient of [Formula: see text] for healthy aging subjects compared to a manual segmentation procedure. Lower values ([Formula: see text]) for Parkinson's disease subjects indicate the need for further investigation and tests before the application to clinical samples. CONCLUSION: These promising results suggest the usability of the proposed algorithm for data of healthy aging subjects and pave the way for further investigations using this approach on different clinical datasets to validate its practical usability more conclusively.

Leveraging biochemical reactions to unravel functional impacts of cancer somatic variants affecting protein interaction interfaces.

Background: Considering protein mutations in their biological context is essential for understanding their functional impact, interpretation of high-dimensional datasets and development of effective targeted therapies in personalized medicine. Methods: We combined the curated knowledge of biochemical reactions from Reactome with the analysis of interaction-mediating 3D interfaces from Mechismo. In addition, we provided a software tool for users to explore and browse the analysis results in a multi-scale perspective starting from pathways and reactions to protein-protein interactions and protein 3D structures. Results: We analyzed somatic mutations from TCGA, revealing several significantly impacted reactions and pathways in specific cancer types. We found examples of genes not yet listed as oncodrivers, whose rare mutations were predicted to affect cancer processes similarly to known oncodrivers. Some identified processes lack any known oncodrivers, which suggests potentially new cancer-related processes (e.g. complement cascade reactions). Furthermore, we found that mutations perturbing certain processes are significantly associated with distinct phenotypes (i.e. survival time) in specific cancer types (e.g. PIK3CA centered pathways in LGG and UCEC cancer types), suggesting the translational potential of our approach for patient stratification. Our analysis also uncovered several druggable processes (e.g. GPCR signalling pathways) containing enriched reactions, providing support for new off-label therapeutic options. Conclusions: In summary, we have established a multi-scale approach to study genetic variants based on protein-protein interaction 3D structures. Our approach is different from previously published studies in its focus on biochemical reactions and can be applied to other data types (e.g. post-translational modifications) collected for many types of disease.

Shared and cell type-specific adaptation strategies of Gag and Env yield high titer bovine foamy virus variants.

During in vitro selection and evolution screens to adapt the tightly cell-associated bovine foamy virus BFV to high titer cell-free transmission, common, cell-type specific and concurrent adaptive changes in Gag and Env, the major players of foamy virus particle assembly and release, were detected. Upon early establishment of cell type-independent pioneering mutations in Env and, subsequently in Gag, a diverse virus pool emerged that was characterized by the occurrence of shared and additional cell type-specific exchanges. At late passages and saturated titers, remarkably homogeneous virus populations characterized by functionally important mutations developed which may be partly due to stochastic evolutionary events that occurred earlier during adaptation. Reverse genetics showed that defined mutations were functionally important for high titer cell-free transmission.

Genetic variants affecting equivalent protein family positions reflect human diversity.

Members of diverse protein families often perform overlapping or redundant functions meaning that different variations within them could reflect differences between individual organisms. We investigated likely functional positions within aligned protein families that contained a significant enrichment of nonsynonymous variants in genomes of healthy individuals. We identified more than a thousand enriched positions across hundreds of family alignments with roles indicative of mammalian individuality, including sensory perception and the immune system. The most significant position is the Arginine from the Olfactory receptor "DRY" motif, which has more variants in healthy individuals than all other positions in the proteome. Odorant binding data suggests that these variants lead to receptor inactivity, and they are mostly mutually exclusive with other loss-of-function (stop/frameshift) variants. Some DRY Arginine variants correlate with smell preferences in sub-populations and all 2,504 humans studied contain a unique spectrum of active and inactive receptors. The many other variant enriched positions, across hundreds of other families might also provide insights into individual differences.

Insights into cancer severity from biomolecular interaction mechanisms.

To attain a deeper understanding of diseases like cancer, it is critical to couple genetics with biomolecular mechanisms. High-throughput sequencing has identified thousands of somatic mutations across dozens of cancers, and there is a pressing need to identify the few that are pathologically relevant. Here we use protein structure and interaction data to interrogate nonsynonymous somatic cancer mutations, identifying a set of 213 molecular interfaces (protein-protein, -small molecule or -nucleic acid) most often perturbed in cancer, highlighting several potentially novel cancer genes. Over half of these interfaces involve protein-small-molecule interactions highlighting their overall importance in cancer. We found distinct differences in the predominance of perturbed interfaces between cancers and histological subtypes and presence or absence of certain interfaces appears to correlate with cancer severity.

Mutagenesis of N-terminal residues of feline foamy virus Gag reveals entirely distinct functions during capsid formation, particle assembly, Gag processing and budding.

BACKGROUND: Foamy viruses (FVs) of the Spumaretrovirinae subfamily are distinct retroviruses, with many features of their molecular biology and replication strategy clearly different from those of the Orthoretroviruses, such as human immunodeficiency, murine leukemia, and human T cell lymphotropic viruses. The FV Gag N-terminal region is responsible for capsid formation and particle budding via interaction with Env. However, the critical residues or motifs in this region and their functional interaction are currently ill-defined, especially in non-primate FVs. RESULTS: Mutagenesis of N-terminal Gag residues of feline FV (FFV) reveals key residues essential for either capsid assembly and/or viral budding via interaction with the FFV Env leader protein (Elp). In an in vitro Gag-Elp interaction screen, Gag mutations abolishing particle assembly also interfered with Elp binding, indicating that Gag assembly is a prerequisite for this highly specific interaction. Gradient sedimentation analyses of cytosolic proteins indicate that wild-type Gag is mostly assembled into virus capsids. Moreover, proteolytic processing of Gag correlates with capsid assembly and is mostly, if not completely, independent from particle budding. In addition, Gag processing correlates with the presence of packaging-competent FFV genomic RNA suggesting that Pol encapsidation via genomic RNA is a prerequisite for Gag processing. Though an appended heterogeneous myristoylation signal rescues Gag particle budding of mutants unable to form capsids or defective in interacting with Elp, it fails to generate infectious particles that co-package Pol, as evidenced by a lack of Gag processing. CONCLUSIONS: Changes in proteolytic Gag processing, intracellular capsid assembly, particle budding and infectivity of defined N-terminal Gag mutants highlight their essential, distinct and only partially overlapping roles during viral assembly and budding. Discussion of these findings will be based on a recent model developed for Gag-Elp interactions in prototype FV.

Genes encoding members of the JAK-STAT pathway or epigenetic regulators are recurrently mutated in T-cell prolymphocytic leukaemia.

T-cell prolymphocytic leukaemia (T-PLL) is an aggressive leukaemia. The primary genetic alteration in T-PLL are the inv(14)(q11q32)/t(14;14)(q11;q32) leading to TRD/TRA-TCL1A fusion, or the t(X;14)(q28;q11) associated with TRD/TRA-MTCP1 fusion. However, additional cooperating abnormalities are necessary for emergence of the full neoplastic phenotype. Though the pattern of secondary chromosomal aberrations is remarkably conserved, targets of the changes are largely unknown. We analysed a cohort of 43 well-characterized T-PLL for hotspot mutations in the genes JAK3, STAT5B and RHOA. Additionally, we selected a subset of 23 T-PLL cases for mutational screening of 54 genes known to be recurrently mutated in T-cell and other haematological neoplasms. Activating mutations in the investigated regions of the JAK3 and STAT5B genes were detected in 30% (13/43) and 21% (8/39) of the cases, respectively, and were mutually exclusive. Further, we identified mutations in the genes encoding the epigenetic regulators EZH2 in 13% (3/23), TET2 in 17% (4/23) and BCOR in 9% (2/23) of the cases. We confirmed that the JAK-STAT pathway is a major mutational target, and identified epigenetic regulators recurrently mutated in T-PLL. These findings complement the mutational spectrum of secondary aberrations in T-PLL and underscore the potential therapeutical relevance of epigenetic regulators in T-PLL.

In Vivo MRI Mapping of Brain Iron Deposition across the Adult Lifespan.

Disruption of iron homeostasis as a consequence of aging is thought to cause iron levels to increase, potentially promoting oxidative cellular damage. Therefore, understanding how this process evolves through the lifespan could offer insights into both the aging process and the development of aging-related neurodegenerative brain diseases. This work aimed to map, in vivo for the first time with an unbiased whole-brain approach, age-related iron changes using quantitative susceptibility mapping (QSM)--a new postprocessed MRI contrast mechanism. To this end, a full QSM standardization routine was devised and a cohort of N = 116 healthy adults (20-79 years of age) was studied. The whole-brain and ROI analyses confirmed that the propensity of brain cells to accumulate excessive iron as a function of aging largely depends on their exact anatomical location. Whereas only patchy signs of iron scavenging were observed in white matter, strong, bilateral, and confluent QSM-age associations were identified in several deep-brain nuclei--chiefly the striatum and midbrain-and across motor, premotor, posterior insular, superior prefrontal, and cerebellar cortices. The validity of QSM as a suitable in vivo imaging technique with which to monitor iron dysregulation in the human brain was demonstrated by confirming age-related increases in several subcortical nuclei that are known to accumulate iron with age. The study indicated that, in addition to these structures, there is a predilection for iron accumulation in the frontal lobes, which when combined with the subcortical findings, suggests that iron accumulation with age predominantly affects brain regions concerned with motor/output functions. SIGNIFICANCE STATEMENT: This study used a whole--brain imaging approach known as quantitative susceptibility mapping (QSM) to provide a novel insight into iron accumulation in the brain across the adult lifespan. Validity of the method was demonstrated by showing concordance with ROI analysis and prior knowledge of iron accumulation in subcortical nuclei. We discovered that, beyond these regions, there is extensive involvement of the frontal lobes that has been missed by past ROI analyses. Broadly speaking, therefore, the motor system selectively accumulates iron with age. The results offer insights into the aging process, but also offer a new approach to studying the role of iron dysregulation in the evolution of age-related neurodegenerative diseases.

DNA methylome analysis in Burkitt and follicular lymphomas identifies differentially methylated regions linked to somatic mutation and transcriptional control.

Although Burkitt lymphomas and follicular lymphomas both have features of germinal center B cells, they are biologically and clinically quite distinct. Here we performed whole-genome bisulfite, genome and transcriptome sequencing in 13 IG-MYC translocation-positive Burkitt lymphoma, nine BCL2 translocation-positive follicular lymphoma and four normal germinal center B cell samples. Comparison of Burkitt and follicular lymphoma samples showed differential methylation of intragenic regions that strongly correlated with expression of associated genes, for example, genes active in germinal center dark-zone and light-zone B cells. Integrative pathway analyses of regions differentially methylated in Burkitt and follicular lymphomas implicated DNA methylation as cooperating with somatic mutation of sphingosine phosphate signaling, as well as the TCF3-ID3 and SWI/SNF complexes, in a large fraction of Burkitt lymphomas. Taken together, our results demonstrate a tight connection between somatic mutation, DNA methylation and transcriptional control in key B cell pathways deregulated differentially in Burkitt lymphoma and other germinal center B cell lymphomas.

The PCBP1 gene encoding poly(rC) binding protein I is recurrently mutated in Burkitt lymphoma.

The genetic hallmark of Burkitt lymphoma is the translocation t(8;14)(q24;q32), or one of its light chain variants, resulting in IG-MYC juxtaposition. However, these translocations alone are insufficient to drive lymphomagenesis, which requires additional genetic changes for malignant transformation. Recent studies of Burkitt lymphoma using next generation sequencing approaches have identified various recurrently mutated genes including ID3, TCF3, CCND3, and TP53. Here, by using similar approaches, we show that PCBP1 is a recurrently mutated gene in Burkitt lymphoma. By whole-genome sequencing, we identified somatic mutations in PCBP1 in 3/17 (18%) Burkitt lymphomas. We confirmed the recurrence of PCBP1 mutations by Sanger sequencing in an independent validation cohort, finding mutations in 3/28 (11%) Burkitt lymphomas and in 6/16 (38%) Burkitt lymphoma cell lines. PCBP1 is an intron-less gene encoding the 356 amino acid poly(rC) binding protein 1, which contains three K-Homology (KH) domains and two nuclear localization signals. The mutations predominantly (10/12, 83%) affect the KH III domain, either by complete domain loss or amino acid changes. Thus, these changes are predicted to alter the various functions of PCBP1, including nuclear trafficking and pre-mRNA splicing. Remarkably, all six primary Burkitt lymphomas with a PCBP1 mutation expressed MUM1/IRF4, which is otherwise detected in around 20-40% of Burkitt lymphomas. We conclude that PCBP1 mutations are recurrent in Burkitt lymphomas and might contribute, in cooperation with other mutations, to its pathogenesis.

Mechismo: predicting the mechanistic impact of mutations and modifications on molecular interactions.

Systematic interrogation of mutation or protein modification data is important to identify sites with functional consequences and to deduce global consequences from large data sets. Mechismo (mechismo.russellab.org) enables simultaneous consideration of thousands of 3D structures and biomolecular interactions to predict rapidly mechanistic consequences for mutations and modifications. As useful functional information often only comes from homologous proteins, we benchmarked the accuracy of predictions as a function of protein/structure sequence similarity, which permits the use of relatively weak sequence similarities with an appropriate confidence measure. For protein-protein, protein-nucleic acid and a subset of protein-chemical interactions, we also developed and benchmarked a measure of whether modifications are likely to enhance or diminish the interactions, which can assist the detection of modifications with specific effects. Analysis of high-throughput sequencing data shows that the approach can identify interesting differences between cancers, and application to proteomics data finds potential mechanistic insights for how post-translational modifications can alter biomolecular interactions.

Recurrent RHOA mutations in pediatric Burkitt lymphoma treated according to the NHL-BFM protocols.

Burkitt lymphoma (BL) is the most frequent B-cell lymphoma in childhood. Genetically, it is characterized by the presence of an IG-MYC translocation which is supposed to be an initiating but not sufficient event in Burkitt lymphomagenesis. In a recent whole-genome sequencing study of four cases, we showed that the gene encoding the ras homolog family member A (RHOA) is recurrently mutated in pediatric BL. Here, we analyzed RHOA by Sanger sequencing in a cohort of 101 pediatric B-cell lymphoma patients treated according to Non-Hodgkin's Lymphoma Berlin-Frankfurt-Münster (NHL-BFM) study protocols. Among the 78 BLs in this series, an additional five had RHOA mutations resulting in a total incidence of 7/82 (8.5%) with c.14G>A (p.R5Q) being present in three cases. Modeling the mutational effect suggests that most of them inactivate the RHOA protein. Thus, deregulation of RHOA by mutation is a recurrent event in Burkitt lymphomagenesis in children.

Recurrent mutation of JAK3 in T-cell prolymphocytic leukemia.

T-cell prolymphocytic leukemia (T-PLL) is an aggressive post-thymic T-cell malignancy characterized by the recurrent inv(14)(q11q32)/t(14;14)(q11;q32) or t(X;14)(q28;q11) leading to activation of either the TCL1 or MTCP1 gene, respectively. However, these primary genetic events are insufficient to drive leukemogenesis. Recently, activating mutations in JAK3 have been identified in other T-cell malignancies. Since JAK3 is essential for T-cell maturation, we analyzed a cohort of 32 T-PLL patients for mutational hot spots in the JAK3 gene using a step-wise screening approach. We identified 14 mutations in 11 of 32 patients (34%). The most frequently detected mutation in our cohort was M511I (seen in 57% of cases) previously described as an activating change in other T-cell malignancies. Three patients carried two mutations in JAK3. In two patients M511I and R657Q were simultaneously detected and in another patient V674F and V678L. In the latter case we could demonstrate that the mutations were on the same allele in cis. Protein modeling and homology analyses of mutations present in other members of the JAK family suggested that these mutations likely activate JAK3, possibly by disrupting the activation loop and the interface between N and C lobes, increasing the accessibility of the catalytic loop. In addition, four of the 21 patients lacking a JAK3 point mutation presented an aberrant karyotype involving the chromosomal band 19p13 harboring the JAK3 locus. The finding of recurrent activating JAK3 mutations in patients with T-PLL could enable the use of JAK3 inhibitors to treat patients with this unfavorable malignancy who otherwise have a very poor prognosis.

Dissecting the genomic complexity underlying medulloblastoma.

Medulloblastoma is an aggressively growing tumour, arising in the cerebellum or medulla/brain stem. It is the most common malignant brain tumour in children, and shows tremendous biological and clinical heterogeneity. Despite recent treatment advances, approximately 40% of children experience tumour recurrence, and 30% will die from their disease. Those who survive often have a significantly reduced quality of life. Four tumour subgroups with distinct clinical, biological and genetic profiles are currently identified. WNT tumours, showing activated wingless pathway signalling, carry a favourable prognosis under current treatment regimens. SHH tumours show hedgehog pathway activation, and have an intermediate prognosis. Group 3 and 4 tumours are molecularly less well characterized, and also present the greatest clinical challenges. The full repertoire of genetic events driving this distinction, however, remains unclear. Here we describe an integrative deep-sequencing analysis of 125 tumour-normal pairs, conducted as part of the International Cancer Genome Consortium (ICGC) PedBrain Tumor Project. Tetraploidy was identified as a frequent early event in Group 3 and 4 tumours, and a positive correlation between patient age and mutation rate was observed. Several recurrent mutations were identified, both in known medulloblastoma-related genes (CTNNB1, PTCH1, MLL2, SMARCA4) and in genes not previously linked to this tumour (DDX3X, CTDNEP1, KDM6A, TBR1), often in subgroup-specific patterns. RNA sequencing confirmed these alterations, and revealed the expression of what are, to our knowledge, the first medulloblastoma fusion genes identified. Chromatin modifiers were frequently altered across all subgroups. These findings enhance our understanding of the genomic complexity and heterogeneity underlying medulloblastoma, and provide several potential targets for new therapeutics, especially for Group 3 and 4 patients.

A novel human NatA Nalpha-terminal acetyltransferase complex: hNaa16p-hNaa10p (hNat2-hArd1).

BACKGROUND: Protein acetylation is among the most common protein modifications. The two major types are post-translational Nepsilon-lysine acetylation catalyzed by KATs (Lysine acetyltransferases, previously named HATs (histone acetyltransferases) and co-translational Nalpha-terminal acetylation catalyzed by NATs (N-terminal acetyltransferases). The major NAT complex in yeast, NatA, is composed of the catalytic subunit Naa10p (N alpha acetyltransferase 10 protein) (Ard1p) and the auxiliary subunit Naa15p (Nat1p). The NatA complex potentially acetylates Ser-, Ala-, Thr-, Gly-, Val- and Cys- N-termini after Met-cleavage. In humans, the homologues hNaa15p (hNat1) and hNaa10p (hArd1) were demonstrated to form a stable ribosome associated NAT complex acetylating NatA type N-termini in vitro and in vivo. RESULTS: We here describe a novel human protein, hNaa16p (hNat2), with 70% sequence identity to hNaa15p (hNat1). The gene encoding hNaa16p originates from an early vertebrate duplication event from the common ancestor of hNAA15 and hNAA16. Immunoprecipitation coupled to mass spectrometry identified both endogenous hNaa15p and hNaa16p as distinct interaction partners of hNaa10p in HEK293 cells, thus demonstrating the presence of both hNaa15p-hNaa10p and hNaa16p-hNaa10p complexes. The hNaa16p-hNaa10p complex acetylates NatA type N-termini in vitro. hNaa16p is ribosome associated, supporting its potential role in cotranslational Nalpha-terminal acetylation. hNAA16 is expressed in a variety of human cell lines, but is generally less abundant as compared to hNAA15. Specific knockdown of hNAA16 induces cell death, suggesting an essential role for hNaa16p in human cells. CONCLUSION: At least two distinct NatA protein Nalpha-terminal acetyltransferases coexist in human cells potentially creating a more complex and flexible system for Nalpha-terminal acetylation as compared to lower eukaryotes.

The molecular architecture of cadherins in native epidermal desmosomes.

Desmosomes are cadherin-based adhesive intercellular junctions, which are present in tissues such as heart and skin. Despite considerable efforts, the molecular interfaces that mediate adhesion remain obscure. Here we apply cryo-electron tomography of vitreous sections from human epidermis to visualize the three-dimensional molecular architecture of desmosomal cadherins at close-to-native conditions. The three-dimensional reconstructions show a regular array of densities at approximately 70 A intervals along the midline, with a curved shape resembling the X-ray structure of C-cadherin, a representative 'classical' cadherin. Model-independent three-dimensional image processing of extracted sub-tomograms reveals the cadherin organization. After fitting the C-cadherin atomic structure into the averaged sub-tomograms, we see a periodic arrangement of a trans W-like and a cis V-like interaction corresponding to molecules from opposing membranes and the same cell membrane, respectively. The resulting model of cadherin organization explains existing two-dimensional data and yields insights into a possible mechanism of cadherin-based cell adhesion.

Characterization of hARD2, a processed hARD1 gene duplicate, encoding a human protein N-alpha-acetyltransferase.

BACKGROUND: Protein acetylation is increasingly recognized as an important mechanism regulating a variety of cellular functions. Several human protein acetyltransferases have been characterized, most of them catalyzing epsilon-acetylation of histones and transcription factors. We recently described the human protein acetyltransferase hARD1 (human Arrest Defective 1). hARD1 interacts with NATH (N-Acetyl Transferase Human) forming a complex expressing protein N-terminal alpha-acetylation activity. RESULTS: We here describe a human protein, hARD2, with 81 % sequence identity to hARD1. The gene encoding hARD2 most likely originates from a eutherian mammal specific retrotransposition event. hARD2 mRNA and protein are expressed in several human cell lines. Immunoprecipitation experiments show that hARD2 protein potentially interacts with NATH, suggesting that hARD2-NATH complexes may be responsible for protein N-alpha-acetylation in human cells. In NB4 cells undergoing retinoic acid mediated differentiation, the level of endogenous hARD1 and NATH protein decreases while the level of hARD2 protein is stable. CONCLUSION: A human protein N-alpha-acetyltransferase is herein described. ARD2 potentially complements the functions of ARD1, adding more flexibility and complexity to protein N-alpha-acetylation in human cells as compared to lower organisms which only have one ARD.