Bi-allelic variants in HMGCR cause an autosomal-recessive progressive limb-girdle muscular dystrophy.

Statins are a mainstay intervention for cardiovascular disease prevention, yet their use can cause rare severe myopathy. HMG-CoA reductase, an essential enzyme in the mevalonate pathway, is the target of statins. We identified nine individuals from five unrelated families with unexplained limb-girdle like muscular dystrophy and bi-allelic variants in HMGCR via clinical and research exome sequencing. The clinical features resembled other genetic causes of muscular dystrophy with incidental high CPK levels (>1,000 U/L), proximal muscle weakness, variable age of onset, and progression leading to impaired ambulation. Muscle biopsies in most affected individuals showed non-specific dystrophic changes with non-diagnostic immunohistochemistry. Molecular modeling analyses revealed variants to be destabilizing and affecting protein oligomerization. Protein activity studies using three variants (p.Asp623Asn, p.Tyr792Cys, and p.Arg443Gln) identified in affected individuals confirmed decreased enzymatic activity and reduced protein stability. In summary, we showed that individuals with bi-allelic amorphic (i.e., null and/or hypomorphic) variants in HMGCR display phenotypes that resemble non-genetic causes of myopathy involving this reductase. This study expands our knowledge regarding the mechanisms leading to muscular dystrophy through dysregulation of the mevalonate pathway, autoimmune myopathy, and statin-induced myopathy.

Cancer-associated polybromo-1 bromodomain 4 missense variants variably impact bromodomain ligand binding and cell growth suppression.

The polybromo, brahma-related gene 1-associated factors (PBAF) chromatin remodeling complex subunit polybromo-1 (PBRM1) contains six bromodomains that recognize and bind acetylated lysine residues on histone tails and other nuclear proteins. PBRM1 bromodomains thus provide a link between epigenetic posttranslational modifications and PBAF modulation of chromatin accessibility and transcription. As a putative tumor suppressor in several cancers, PBRM1 protein expression is often abrogated by truncations and deletions. However, ∼33% of PBRM1 mutations in cancer are missense and cluster within its bromodomains. Such mutations may generate full-length PBRM1 variant proteins with undetermined structural and functional characteristics. Here, we employed computational, biophysical, and cellular assays to interrogate the effects of PBRM1 bromodomain missense variants on bromodomain stability and function. Since mutations in the fourth bromodomain of PBRM1 (PBRM1-BD4) comprise nearly 20% of all cancer-associated PBRM1 missense mutations, we focused our analysis on PBRM1-BD4 missense protein variants. Selecting 16 potentially deleterious PBRM1-BD4 missense protein variants for further study based on high residue mutational frequency and/or conservation, we show that cancer-associated PBRM1-BD4 missense variants exhibit varied bromodomain stability and ability to bind acetylated histones. Our results demonstrate the effectiveness of identifying the unique impacts of individual PBRM1-BD4 missense variants on protein structure and function, based on affected residue location within the bromodomain. This knowledge provides a foundation for drawing correlations between specific cancer-associated PBRM1 missense variants and distinct alterations in PBRM1 function, informing future cancer personalized medicine approaches.

Advanced computational analysis of CD40LG variants in atypical X-linked hyper-IgM syndrome.

X-Linked Hyper-IgM Syndrome is caused by pathogenic variants in CD40LG. Three patients with atypical clinical and immunological features were identified with variants in CD40LG requiring further characterization. Flow cytometry was used to evaluate CD40L protein expression and binding capacity to a surrogate receptor, CD40-muIg. Though functional anomalies were observed, there was still a lack of clarity regarding the underlying mechanism. We developed structural models for wild-type and the three variants of CD40L protein observed in these patients (p. Lys143Asn, Leu225Ser and Met36Arg) to evaluate structural alterations by molecular mechanic calculations, and assess protein movement by molecular dynamic simulations. These studies demonstrate that functional analysis of variants of unknown significance in CD40LG can be supplemented by advanced computational analysis in atypical clinical contexts. These studies in combination identify the deleterious effects of these variants and potential mechanisms for protein dysfunction.

SASH3 variants cause a novel form of X-linked combined immunodeficiency with immune dysregulation.

Sterile alpha motif (SAM) and Src homology-3 (SH3) domain-containing 3 (SASH3), also called SH3-containing lymphocyte protein (SLY1), is a putative adaptor protein that is postulated to play an important role in the organization of signaling complexes and propagation of signal transduction cascades in lymphocytes. The SASH3 gene is located on the X-chromosome. Here, we identified 3 novel SASH3 deleterious variants in 4 unrelated male patients with a history of combined immunodeficiency and immune dysregulation that manifested as recurrent sinopulmonary, cutaneous, and mucosal infections and refractory autoimmune cytopenias. Patients exhibited CD4+ T-cell lymphopenia, decreased T-cell proliferation, cell cycle progression, and increased T-cell apoptosis in response to mitogens. In vitro T-cell differentiation of CD34+ cells and molecular signatures of rearrangements at the T-cell receptor α (TRA) locus were indicative of impaired thymocyte survival. These patients also manifested neutropenia and B-cell and natural killer (NK)-cell lymphopenia. Lentivirus-mediated transfer of the SASH3 complementary DNA-corrected protein expression, in vitro proliferation, and signaling in SASH3-deficient Jurkat and patient-derived T cells. These findings define a new type of X-linked combined immunodeficiency in humans that recapitulates many of the abnormalities reported in mice with Sly1-/- and Sly1Δ/Δ mutations, highlighting an important role of SASH3 in human lymphocyte function and survival.

IL6Myc mouse is an immunocompetent model for the development of aggressive multiple myeloma.

Multiple Myeloma (MM) is a plasma cell neoplasm originating in the bone marrow and is the second most common blood cancer in the United States. One challenge in understanding the pathogenesis of MM and improving treatment is a lack of immunocompetent mouse models. We previously developed the IL6Myc mouse that generates plasmacytomas at 100% penetrance that phenotypically resemble aggressive MM. Using comprehensive genomic analysis, we found that the IL6Myc tumors resemble aggressive MM by RNA and protein expression. We also found that IL6Myc tumors accumulated fusions and missense mutations in genes that overlap significantly with human myeloma, indicating that the mouse is good model for studying disease etiology. Lastly, we derived cell lines from IL6Myc tumors that express cell surface markers typical of MM and readily engraft into mice, home to the bone marrow, and induce osteolytic disease. The cell lines may be useful in developing immunotherapies directed against BAFF-R and TACI, though not BCMA, and may also be a good model for studying dexamethasone resistance. These data indicate that the IL6Myc model is useful for studying development of aggressive MM and for developing new treatments against such forms of the disease.

Computational modeling reveals key molecular properties and dynamic behavior of disruptor of telomeric silencing 1-like (DOT1L) and partnering complexes involved in leukemogenesis.

Disruptor of telomeric silencing 1-like (DOT1L) is the only non-SET domain histone lysine methyltransferase (KMT) and writer of H3K79 methylation on nucleosomes marked by H2B ubiquitination. DOT1L has elicited significant attention because of its interaction or fusion with members of the AF protein family in blood cell biology and leukemogenic transformation. Here, our goal was to extend previous structural information by performing a robust molecular dynamic study of DOT1L and its leukemogenic partners combined with mutational analysis. We show that statically and dynamically, D161, G163, E186, and F223 make frequent time-dependent interactions with SAM, while additional residues T139, K187, and N241 interact with SAM only under dynamics. Dynamics models reveal DOT1L, SAM, and H4 moving as one and show that more than twice the number of DOT1L residues interacts with these partners, relative to the static structure. Mutational analyses indicate that six of these residues are intolerant to substitution. We describe the dynamic behavior of DOT1L interacting with AF10 and AF9. Studies on the dynamics of a heterotrimeric complex of DOT1L1-AF10 illuminated describe coordinated motions that impact the relative position of the DOT1L HMT domain to the nucleosome. The molecular motions of the DOT1L-AF9 complex are less extensive and highly dynamic, resembling a swivel-like mechanics. Through molecular dynamics and mutational analysis, we extend the knowledge previous provided by static measurements. These results are important to consider when describing the biochemical properties of DOT1L, under normal and in disease conditions, as well as for the development of novel therapeutic agents.

Structural bioinformatics enhances mechanistic interpretation of genomic variation, demonstrated through the analyses of 935 distinct RAS family mutations.

MOTIVATION: Protein-coding genetic alterations are frequently observed in Clinical Genetics, but the high yield of variants of uncertain significance remains a limitation in decision making. RAS-family GTPases are cancer drivers, but only 54 variants, across all family members, fall within well-known hotspots. However, extensive sequencing has identified 881 non-hotspot variants for which significance remains to be investigated. RESULTS: Here, we evaluate 935 missense variants from seven RAS genes, observed in cancer, RASopathies and the healthy adult population. We characterized hotspot variants, previously studied experimentally, using 63 sequence- and 3D structure-based scores, chosen by their breadth of biophysical properties. Applying scores that display best correlation with experimental measures, we report new valuable mechanistic inferences for both hot-spot and non-hotspot variants. Moreover, we demonstrate that 3D scores have little-to-no correlation with those based on DNA sequence, which are commonly used in Clinical Genetics. Thus, combined, these new knowledge bear significant relevance. AVAILABILITY AND IMPLEMENTATION: All genomic and 3D scores, and markdown for generating figures, are provided in our supplemental data. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

NADPH oxidase 4 contributes to TRPV4-mediated endothelium-dependent vasodilation in human arterioles by regulating protein phosphorylation of TRPV4 channels.

Impaired endothelium-dependent vasodilation has been suggested to be a key component of coronary microvascular dysfunction (CMD). A better understanding of endothelial pathways involved in vasodilation in human arterioles may provide new insight into the mechanisms of CMD. The goal of this study is to investigate the role of TRPV4, NOX4, and their interaction in human arterioles and examine the underlying mechanisms. Arterioles were freshly isolated from adipose and heart tissues obtained from 71 patients without coronary artery disease, and vascular reactivity was studied by videomicroscopy. In human adipose arterioles (HAA), ACh-induced dilation was significantly reduced by TRPV4 inhibitor HC067047 and by NOX 1/4 inhibitor GKT137831, but GKT137831 did not further affect the dilation in the presence of TRPV4 inhibitors. GKT137831 also inhibited TRPV4 agonist GSK1016790A-induced dilation in HAA and human coronary arterioles (HCA). NOX4 transcripts and proteins were detected in endothelial cells of HAA and HCA. Using fura-2 imaging, GKT137831 significantly reduced GSK1016790A-induced Ca2+ influx in the primary culture of endothelial cells and TRPV4-WT-overexpressing human coronary artery endothelial cells (HCAEC). However, GKT137831 did not affect TRPV4-mediated Ca2+ influx in non-phosphorylatable TRPV4-S823A/S824A-overexpressing HCAEC. In addition, treatment of HCAEC with GKT137831 decreased the phosphorylation level of Ser824 in TRPV4. Finally, proximity ligation assay (PLA) revealed co-localization of NOX4 and TRPV4 proteins. In conclusion, both TRPV4 and NOX4 contribute to ACh-induced dilation in human arterioles from patients without coronary artery disease. NOX4 increases TRPV4 phosphorylation in endothelial cells, which in turn enhances TRPV4-mediated Ca2+ entry and subsequent endothelium-dependent dilation in human arterioles.

Somatic PIK3R1 variation as a cause of vascular malformations and overgrowth.

PURPOSE: Somatic activating variants in the PI3K-AKT pathway cause vascular malformations with and without overgrowth. We previously reported an individual with capillary and lymphatic malformation harboring a pathogenic somatic variant in PIK3R1, which encodes three PI3K complex regulatory subunits. Here, we investigate PIK3R1 in a large cohort with vascular anomalies and identify an additional 16 individuals with somatic mosaic variants in PIK3R1. METHODS: Affected tissue from individuals with vascular lesions and overgrowth recruited from a multisite collaborative network was studied. Next-generation sequencing targeting coding regions of cell-signaling and cancer-associated genes was performed followed by assessment of variant pathogenicity. RESULTS: The phenotypic and variant spectrum associated with somatic variation in PIK3R1 is reported herein. Variants occurred in the inter-SH2 or N-terminal SH2 domains of all three PIK3R1 protein products. Phenotypic features overlapped those of the PIK3CA-related overgrowth spectrum (PROS). These overlapping features included mixed vascular malformations, sandal toe gap deformity with macrodactyly, lymphatic malformations, venous ectasias, and overgrowth of soft tissue or bone. CONCLUSION: Somatic PIK3R1 variants sharing attributes with cancer-associated variants cause complex vascular malformations and overgrowth. The PIK3R1-associated phenotypic spectrum overlaps with PROS. These data extend understanding of the diverse phenotypic spectrum attributable to genetic variation in the PI3K-AKT pathway.

Impact of integrated translational research on clinical exome sequencing.

PURPOSE: Exome sequencing often identifies pathogenic genetic variants in patients with undiagnosed diseases. Nevertheless, frequent findings of variants of uncertain significance necessitate additional efforts to establish causality before reaching a conclusive diagnosis. To provide comprehensive genomic testing to patients with undiagnosed disease, we established an Individualized Medicine Clinic, which offered clinical exome testing and included a Translational Omics Program (TOP) that provided variant curation, research activities, or research exome sequencing. METHODS: From 2012 to 2018, 1101 unselected patients with undiagnosed diseases received exome testing. Outcomes were reviewed to assess impact of the TOP and patient characteristics on diagnostic rates through descriptive and multivariate analyses. RESULTS: The overall diagnostic yield was 24.9% (274 of 1101 patients), with 174 (15.8% of 1101) diagnosed on the basis of clinical exome sequencing alone. Four hundred twenty-three patients with nondiagnostic or without access to clinical exome sequencing were evaluated by the TOP, with 100 (9% of 1101) patients receiving a diagnosis, accounting for 36.5% of the diagnostic yield. The identification of a genetic diagnosis was influenced by the age at time of testing and the disease phenotype of the patient. CONCLUSION: Integration of translational research activities into clinical practice of a tertiary medical center can significantly increase the diagnostic yield of patients with undiagnosed disease.

Recurrent MSC E116K mutations in ALK-negative anaplastic large cell lymphoma.

Anaplastic large cell lymphomas (ALCLs) represent a relatively common group of T-cell non-Hodgkin lymphomas (T-NHLs) that are unified by similar pathologic features but demonstrate marked genetic heterogeneity. ALCLs are broadly classified as being anaplastic lymphoma kinase (ALK)+ or ALK-, based on the presence or absence of ALK rearrangements. Exome sequencing of 62 T-NHLs identified a previously unreported recurrent mutation in the musculin gene, MSC E116K, exclusively in ALK- ALCLs. Additional sequencing for a total of 238 T-NHLs confirmed the specificity of MSC E116K for ALK- ALCL and further demonstrated that 14 of 15 mutated cases (93%) had coexisting DUSP22 rearrangements. Musculin is a basic helix-loop-helix (bHLH) transcription factor that heterodimerizes with other bHLH proteins to regulate lymphocyte development. The E116K mutation localized to the DNA binding domain of musculin and permitted formation of musculin-bHLH heterodimers but prevented their binding to authentic target sequence. Functional analysis showed MSCE116K acted in a dominant-negative fashion, reversing wild-type musculin-induced repression of MYC and cell cycle inhibition. Chromatin immunoprecipitation-sequencing and transcriptome analysis identified the cell cycle regulatory gene E2F2 as a direct transcriptional target of musculin. MSCE116K reversed E2F2-induced cell cycle arrest and promoted expression of the CD30-IRF4-MYC axis, whereas its expression was reciprocally induced by binding of IRF4 to the MSC promoter. Finally, ALCL cells expressing MSC E116K were preferentially targeted by the BET inhibitor JQ1. These findings identify a novel recurrent MSC mutation as a key driver of the CD30-IRF4-MYC axis and cell cycle progression in a unique subset of ALCLs.

Mitochondrial Metabolic Reprogramming by CD36 Signaling Drives Macrophage Inflammatory Responses.

RATIONALE: A hallmark of chronic inflammatory disorders is persistence of proinflammatory macrophages in diseased tissues. In atherosclerosis, this is associated with dyslipidemia and oxidative stress, but mechanisms linking these phenomena to macrophage activation remain incompletely understood. OBJECTIVE: To investigate mechanisms linking dyslipidemia, oxidative stress, and macrophage activation through modulation of immunometabolism and to explore therapeutic potential targeting specific metabolic pathways. METHODS AND RESULTS: Using a combination of biochemical, immunologic, and ex vivo cell metabolic studies, we report that CD36 mediates a mitochondrial metabolic switch from oxidative phosphorylation to superoxide production in response to its ligand, oxidized LDL (low-density lipoprotein). Mitochondrial-specific inhibition of superoxide inhibited oxidized LDL-induced NF-κB (nuclear factor-κB) activation and inflammatory cytokine generation. RNA sequencing, flow cytometry, 3H-labeled palmitic acid uptake, lipidomic analysis, confocal and electron microscopy imaging, and functional energetics revealed that oxidized LDL upregulated effectors of long-chain fatty acid uptake and mitochondrial import, while downregulating fatty acid oxidation and inhibiting ATP5A (ATP synthase F1 subunit alpha)-an electron transport chain component. The combined effect is long-chain fatty acid accumulation, alteration of mitochondrial structure and function, repurposing of the electron transport chain to superoxide production, and NF-κB activation. Apoe null mice challenged with high-fat diet showed similar metabolic changes in circulating Ly6C+ monocytes and peritoneal macrophages, along with increased CD36 expression. Moreover, mitochondrial reactive oxygen species were positively correlated with CD36 expression in aortic lesional macrophages. CONCLUSIONS: These findings reveal that oxidized LDL/CD36 signaling in macrophages links dysregulated fatty acid metabolism to oxidative stress from the mitochondria, which drives chronic inflammation. Thus, targeting to CD36 and its downstream effectors may serve as potential new strategies against chronic inflammatory diseases such as atherosclerosis.

Germline evaluation of patients undergoing tumor genomic profiling: An academic cancer center's experience with implementing a germline review protocol.

Tumor genomic profiling (TGP) has the potential to identify germline variants in addition to its primary use of informing cancer treatment based on genetic alterations within the tumor. However, there are no formal consensus guidelines to identify patients who would be eligible for genetic counseling (GC) and germline testing (GT) testing in patients undergoing TGP. The purpose of this study is to describe an institutionally developed Germline Review Protocol (GRP) to evaluate adult cancer patient cases already undergoing TGP to determine GC referral eligibility. We report on our retrospective experience implementing this protocol into practice wherein 172 patients out of 638 patients reviewed (27%) were recommended for a GC referral over a 17-month time period. Of those 172 patients recommended for a GC referral, only 34 patients (20%) completed GC and GT. Among patients who received GT, 15 (44%) were positive for at least one pathogenic or likely pathogenic (P/LP) variant, seven patients (21%) were negative and 12 patients (35%) had at least 1 variant of uncertain significance (VUS). The primary reason GC and GT was not completed was because the patient moved to hospice care or was deceased. This is one of the first studies outlining the process and results of a formalized institutional protocol to facilitate patient referrals for GC and GT based on TGP results.

Modeling post-translational modifications and cancer-associated mutations that impact the heterochromatin protein 1α-importin α heterodimers.

Heterochromatin protein 1α (HP1α) is a protein that mediates cancer-associated processes in the cell nucleus. Proteomic experiments, reported here, demonstrate that HP1α complexes with importin α (IMPα), a protein necessary for its nuclear transport. This data is congruent with Simple Linear Motif (SLiM) analyses that identify an IMPα-binding motif within the linker that joins the two globular domains of this protein. Using molecular modeling and dynamics simulations, we develop a model of the IMPα-HP1α complex and investigate the impact of phosphorylation and genomic variants on their interaction. We demonstrate that phosphorylation of the HP1α linker likely regulates its association with IMPα, which has implications for HP1α access to the nucleus, where it functions. Cancer-associated genomic variants do not abolish the interaction of HP1α but instead lead to rearrangements where the variant proteins maintain interaction with IMPα, but with less specificity. Combined, this new mechanistic insight bears biochemical, cell biological, and biomedical relevance.

Genetic variations in human papillomavirus and cervical cancer outcomes.

Cervical cancer is driven by persistent infection of human papillomavirus (HPV), which is influenced by HPV type and intratypic variants, yet the impact of HPV type and intratypic variants on patient outcomes is far less understood. Here, we examined the association of cervical cancer stage and survival with HPV type, clade, lineage, and intratypic variants within the HPV E6 locus. Of 1,028 HPV-positive cases recruited through the CerGE study, 301 were in-situ and 727 were invasive cervical cancer (ICC), with an average post-diagnosis follow-up of 4.8 years. HPV sequencing was performed using tumor-isolated DNA to assign HPV type, HPV 16 lineage, clade, and intratypic variants within the HPV 16 E6 locus, of which nonsynonomous variants were functionally annotated by molecular modeling. HPV 18-related types were more prevalent in ICC compared to in-situ disease and associated with significantly worse recurrence-free survival (RFS) compared to HPV 16-related types. The HPV 16 Asian American lineage D3 and Asian lineage A4 associated more frequently with ICC than with in situ disease and women with an intratypic HPV 16 lineage B exhibited a trend toward worse RFS than those with A, C, or D lineages. Participants with intratypic E6 variants predicted to stabilize the E6-E6AP-p53 complex had worse RFS. Variants within the highly immunogenic HPV 16 E6 region (E14-I34) were enriched in ICC compared to in-situ lesions but were not associated with survival. Collectively, our results suggest that cervical cancer outcome is associated with HPV variants that affect virus-host interactions.

Molecular characterization of known and novel ACVR1 variants in phenotypes of aberrant ossification.

Diffuse idiopathic skeletal hyperostosis (DISH) is a disorder principally characterized by calcification and ossification of spinal ligaments and entheses. Fibrodysplasia ossificans progressiva (FOP) is a rare autosomal dominant disabling disorder characterized by progressive ossification of skeletal muscle, fascia, tendons, and ligaments. These conditions manifest phenotypic overlap in the ossification of tendons and ligaments. We describe herein a patient with DISH, exhibiting heterotopic ossification of the posterior longitudinal ligament where clinical whole exome sequencing identified a variant within ACVR1, a gene implicated in FOP. This variant, p.K400E, is a novel variant, not identified previously, and occurs in a highly conserved region across orthologs. We used sequence-based predicative algorithms, molecular modeling, and molecular dynamics simulations, to test the potential for p.K400E to alter the structure and dynamics of ACVR1. We applied the same modeling and simulation methods to established FOP variants, to identify the detailed effects that they have on the ACVR1 protein, as well as to act as positive controls against which the effects of p.K400E could be evaluated. Our in silico molecular analyses support p.K400E as altering the behavior of ACVR1. In addition, functional testing to measure the effect of this variant on BMP-pSMAD 1/5/8 target genes was carried out which revealed this variant to cause increased ID1 and Msx2 expression compared with the wild-type receptor. This analysis supports the potential for the variant of uncertain significance to contribute to the patient's phenotype.

A Novel Kleefstra Syndrome-associated Variant That Affects the Conserved TPLX Motif within the Ankyrin Repeat of EHMT1 Leads to Abnormal Protein Folding.

Kleefstra syndrome (KS) (Mendelian Inheritance in Man (MIM) no. 610253), also known as 9q34 deletion syndrome, is an autosomal dominant disorder caused by haploinsufficiency of euchromatic histone methyltransferase-1 (EHMT1). The clinical phenotype of KS includes moderate to severe intellectual disability with absent speech, hypotonia, brachycephaly, congenital heart defects, and dysmorphic facial features with hypertelorism, synophrys, macroglossia, protruding tongue, and prognathism. Only a few cases of de novo missense mutations in EHMT1 giving rise to KS have been described. However, some EHMT1 variants have been described in individuals presenting with autism spectrum disorder or mild intellectual disability, suggesting that the phenotypic spectrum resulting from EHMT1 alterations may be quite broad. In this report, we describe two unrelated patients with complex medical histories consistent with KS in whom next generation sequencing identified the same novel c.2426C>T (p.P809L) missense variant in EHMT1 To examine the functional significance of this novel variant, we performed molecular dynamics simulations of the wild type and p.P809L variant, which predicted that the latter would have a propensity to misfold, leading to abnormal histone mark binding. Recombinant EHMT1 p.P809L was also studied using far UV circular dichroism spectroscopy and intrinsic protein fluorescence. These functional studies confirmed the model-based hypotheses and provided evidence for protein misfolding and aberrant target recognition as the underlying pathogenic mechanism for this novel KS-associated variant. This is the first report to suggest that missense variants in EHMT1 that lead to protein misfolding and disrupted histone mark binding can lead to KS.

Correction to: the impact of pharmacokinetic gene profiles across human cancers.

Correction to: BMC Cancer (2018) 18:577 DOI https://doi.org/10.1186/s12885-018-4345-2.