Analysis of genetic variants in myeloproliferative neoplasms using a 22-gene next-generation sequencing panel

Background The Philadelphia (Ph)-negative myeloproliferative neoplasms (MPNs), namely essential thrombocythaemia (ET), polycythaemia vera (PV) and primary myelofibrosis (PMF), are a group of chronic clonal haematopoietic disorders that have the propensity to advance into bone marrow failure or acute myeloid leukaemia; often resulting in fatality. Although driver mutations have been identified in these MPNs, subtype-specific markers of the disease have yet to be discovered. Next-generation sequencing (NGS) technology can potentially improve the clinical management of MPNs by allowing for the simultaneous screening of many disease-associated genes. Methods The performance of a custom, in-house designed 22-gene NGS panel was technically validated using reference standards across two independent replicate runs. The panel was subsequently used to screen a total of 10 clinical MPN samples (ET n = 3, PV n = 3, PMF n = 4). The resulting NGS data was then analysed via a bioinformatics pipeline. Results The custom NGS panel had a detection limit of 1% variant allele frequency (VAF). A total of 20 unique variants with VAFs above 5% (4 of which were putatively novel variants with potential biological significance) and one pathogenic variant with a VAF of between 1 and 5% were identified across all of the clinical MPN samples. All single nucleotide variants with VAFs ≥ 15% were confirmed via Sanger sequencing. Conclusions The high fidelity of the NGS analysis and the identification of known and novel variants in this study cohort support its potential clinical utility in the management of MPNs. However, further optimisation is needed to avoid false negatives in regions with low sequencing coverage, especially for the detection of driver mutations in MPL.

Evaluating the association of TRPA1 gene polymorphisms with pain sensitivity: a protocol for an adaptive recall by genotype study

Background Pain is a complex polygenic trait whose common genetic underpinnings are relatively ill-defined due in part to challenges in measuring pain as a phenotype. Pain sensitivity can be quantified, but this is difficult to perform at the scale required for genome wide association studies (GWAS). Existing GWAS of pain have identified surprisingly few loci involved in nociceptor function which contrasts strongly with rare monogenic pain states. This suggests a lack of resolution with current techniques. We propose an adaptive methodology within a recall-by-genotype (RbG) framework using detailed phenotyping to screen minor alleles in a candidate ‘nociceptor’ gene in an attempt to estimate their genetic contribution to pain. Methods/design Participants of the Avon Longitudinal Study of Parents and Children will be recalled on the basis of genotype at five common non-synonomous SNPs in the ‘nociceptor’ gene transient receptor potential ankylin 1 (TRPA1). Those homozygous for the common alleles at each of the five SNPs will represent a control group. Individuals homozygous for the minor alleles will then be recruited in a series of three sequential test groups. The outcome of a pre-planned early assessment (interim) of the current test group will determine whether to continue recruitment or switch to the next test group. Pain sensitivity will be assessed using quantitative sensory testing (QST) before and after topical application of 10% cinnamaldehyde (a TRPA1 agonist). Discussion The design of this adaptive RbG study offers efficiency in the assessment of associations between genetic variation at TRPA1 and detailed pain phenotypes. The possibility to change the test group in response to preliminary data increases the likelihood to observe smaller effect sizes relative to a conventional multi-armed design, as well as reducing futile testing of participants where an effect is unlikely to be observed. This specific adaptive RbG design aims to uncover the influence of common TRPA1 variants on pain sensation but can be applied to any hypothesis-led genotype study where costly and time intensive investigation is required and / or where there is large uncertainty around the expected effect size. Trial registration: ISRCTN, ISRCTN16294731. Retrospectively registered 25th November 2021.

Identification of the pivotal role of SPP1 in kidney stone disease based on multiple bioinformatics analysis

Background Kidney stone disease (KSD) is a multifactorial disease involving both environmental and genetic factors, whose pathogenesis remains unclear. This study aims to explore the hub genes related to stone formation that could serve as potential therapeutic targets. Methods Based on the GSE73680 dataset with 62 samples, differentially expressed genes (DEGs) between Randall’s plaque (RP) tissues and normal tissues were screened and weighted gene co-expression network analysis (WGCNA) was applied to identify key modules associated with KSD. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis were performed to explore the biological functions. The protein–protein interaction (PPI) network was constructed to identify hub genes. Meanwhile, CIBERSORT and ssGSEA analysis were used to estimate the infiltration level of the immune cells. The correlations between hub genes and immune infiltration levels were also investigated. Finally, the top hub gene was selected for further GSEA analysis. Results A total of 116 DEGs, including 73 up-regulated and 43 down-regulated genes, were screened in the dataset. The red module was identified as the key module correlated with KSD. 53 genes were obtained for functional enrichment analysis by taking the intersection of DEGs and genes in the red module. GO analysis showed that these genes were mainly involved in extracellular matrix organization (ECM) and extracellular structure organization, and others. KEGG analysis revealed that the pathways of aldosterone-regulated sodium reabsorption, cell adhesion molecules, arachidonic acid (AA) metabolism, and ECM-receptor interaction were enriched. Through PPI network construction, 30 hub genes were identified. CIBERSORT analysis revealed a significantly increased proportion of M0 macrophages, while ssGSEA revealed no significant differences. Among these hub genes, SPP1, LCN2, MMP7, MUC1, SCNN1A, CLU, SLP1, LAMC2, and CYSLTR2 were positively correlated with macrophages infiltration. GSEA analysis found that positive regulation of JNK activity was enriched in RP tissues with high SPP1 expression, while negative regulation of IL-1β production was enriched in the low-SPP1 subgroup. Conclusions There are 30 hub genes associated with KSD, among which SPP1 is the top hub gene with the most extensive links with other hub genes. SPP1 might play a pivotal role in the pathogenesis of KSD, which is expected to become a potential therapeutic target, while its interaction with macrophages in KSD needs further investigation.

Association of SMAD7 genetic markers and haplotypes with colorectal cancer risk

Purpose Colorectal cancer (CRC) is one of the common cancers with a high mortality rate worldwide. In Iran, there has been a trend of increased incidence of colorectal cancer in the last three decades that necessitates the early diagnosis. Genetic factors have an influential role in its etiology along with the conventional risk factors such as age, diet, and lifestyle. Results from GWAS have shown significant associations between SMAD7 gene variants and risk of CRC. This study aimed to assess the association of certain polymorphisms as well as haplotypes of this gene and risk of colorectal cancer. Methods and materials This study was designed as a case–control association study. After obtaining ethical approval and informed consent, blood samples from 209 patients with colorectal cancer were collected and DNA was extracted. Four variants: rs4939827, rs34007497, rs8085824 and rs8088297 were genotyped using ARMS-PCR method. Results SMAD7 rs4939827 in the recessive and co-dominant models was associated with colorectal cancer risk [TT/CT + CC: OR = 2.90, 95%CI (1.38–6.09), p = 0.005; CC + TT/CT: OR = 1.66, 95%CI (1.00–2.75), p = 0.01]. Haplotype analysis indicated that some SNP combinations including two for-SNPs haplotypes of T-T-C-C and T-C-C-A were significantly associated with CRC risk. Conclusion Based on the identified association of SMAD7 gene variations and haplotypes with colorectal cancer risk in our population, genetic variations in this gene region may have a role in CRC development. This data may shed light on the genetic predisposition of CRC which involves different pathways including TGF-β.

Clinicopathological investigation of secretory carcinoma cases including a successful treatment outcome using entrectinib for high-grade transformation: a case report

Background Secretory carcinoma (SC) of the salivary gland is a recently described malignant tumor harboring characteristic ETV6-NTRK3 gene fusion. SC generally has a favorable clinical course, and is currently regarded as a low-grade carcinoma. However, a small subset of SCs demonstrates aggressive clinical features with histologically high-grade transformed morphology, the molecular pathogenesis of which has not yet been elucidated. In this study, we performed a clinicopathological and molecular genetic study of patients with SC of the head and neck displaying various clinical characteristics to investigate the differences of pathological and molecular genetics between low-grade and high-grade components of SC. Case presentation Three cases with SC of the head and neck, including a conventional low-grade SC and two high-grade transformed SCs are described. High-grade transformed SCs with histological features such as nuclear polymorphism, distinctive nucleoli and increased mitotic activity developed locoregional recurrence and distant metastasis. Immunohistochemical analysis revealed that low- and high-grade components showed different expression patterns for S-100 protein and mammaglobin, whereas all examined components were positive for p-STAT5. p53-positive cell population was markedly higher in one case with high-grade transformed SC. The proliferative activity of high-grade components was markedly increased, with the Ki-67 labeling index ranging up to 30–32%. A fluorescence in situ hybridization study with an ETV6 (12p13) break apart probe revealed split signals in the nuclei in all 3 cases. A targeted next-generation sequencing-based fusion assay demonstrated that all 6 clinical samples from the 3 patients showed the presence of the ETV6-NTRK3 fusion transcripts. One patient with high-grade transformed SC showed a dramatic clinical response to the pan-TRK inhibitor, entrectinib, for the treatment of locoregional recurrence and pulmonary metastasis. Conclusions High-grade transformed SC showed aggressive clinical and pathological features with increased Ki-67 labeling index. Molecular genetic study of gene rearrangement appears to be beneficial treatment as the presence of ETV6-NTRK3 translocation may represent a therapeutic target in SC, particularly the high-grade transformed type.

Clinical and genetic investigation of ichthyosis in familial and sporadic cases in south of Tunisia: genotype–phenotype correlation

Background Ichthyosis is a heterogeneous group of Mendelian cornification disorders that includes syndromic and non-syndromic forms. Autosomal Recessive Congenital Ichthyosis (ARCI) and Ichthyosis Linearis Circumflexa (ILC) belong to non-syndromic forms. Syndromic ichthyosis is rather a large group of heterogeneous diseases. Overlapping phenotypes and genotypes between these disorders is a major characteristic. Therefore, determining the specific genetic background for each form would be necessary. Methods A total of 11 Tunisian patients with non-syndromic (8 with ARCI and 2 with ILC) and autosomal syndromic ichthyosis (1 patient) were screened by a custom Agilent HaloPlex multi-gene panel and the segregation of causative mutations were analyzed in available family members. Results Clinical and molecular characterization, leading to genotype–phenotype correlation in 11 Tunisian patients was carried out. Overall, we identified 8 mutations in 5 genes. Thus, in patients with ARCI, we identified a novel (c.118T > C in NIPAL4) and 4 already reported mutations (c.534A > C in NIPAL4; c.788G > A and c.1042C > T in TGM1 and c.844C > T in CYP4F22). Yellowish severe keratoderma was found to be associated with NIPAL4 variations and brachydactyly to TGM1 mutations. Two novel variations (c.5898G > C and c.2855A > G in ABCA12) seemed to be features of ILC. Delexon13 in CERS3 was reported in a patient with syndromic ichthyosis. Conclusions Our study further extends the spectrum of mutations involved in ichthyosis as well as clinical features that could help directing genetic investigation.

An intron SNP rs2069837 in IL-6 is associated with osteonecrosis of the femoral head development

Background Genetic polymorphisms play a crucial role in the development of osteonecrosis of the femoral head (ONFH). This study mainly explored the association of IL-6 variants and ONFH susceptibility among the Chinese Han population. Methods Two variants (rs2069837, and rs13306435) in the IL-6 gene were identified and genotyped from 566 patients with ONFH and 566 healthy controls. The associations between IL-6 polymorphisms and ONFH susceptibility were assessed using odds ratio (OR) and 95% confidence interval (95% CI) via logistic regression. The potential function of these two variants was predicted by the HaploReg online database. Results The results of the overall analysis revealed that IL-6 rs2069837 was correlated with decreased risk of ONFH among the Chinese Han population (p < 0.05). In stratified analysis, rs2069837 also reduced the susceptibility to ONFH in older people (> 51 years), males, nonsmokers, and nondrinkers (p < 0.05). However, no associations between rs13306435 and ONFH susceptibility were observed (p > 0.05). Conclusions To sum up, we suggested that rs2069837 G>A polymorphism in the IL-6 gene was significantly associated with a decreased risk of ONFH among the Chinese Hans. These findings underscored the crucial role of IL-6 rs2069837 in the occurrence of ONFH.

The frequency of the known mitochondrial variants associated with drug-induced toxicity in a Korean population

Background Few studies have annotated the whole mitochondrial DNA (mtDNA) genome associated with drug responses in Asian populations. This study aimed to characterize mtDNA genetic profiles, especially the distribution and frequency of well-known genetic biomarkers associated with diseases and drug-induced toxicity in a Korean population. Method Whole mitochondrial genome was sequenced for 118 Korean subjects by using a next-generation sequencing approach. The bioinformatic pipeline was constructed for variant calling, haplogroup classification and annotation of mitochondrial mutation. Results A total of 681 variants was identified among all subjects. The MT-TRNP gene and displacement loop showed the highest numbers of variants (113 and 74 variants, respectively). The m.16189T > C allele, which is known to reduce the mtDNA copy number in human cells was detected in 25.4% of subjects. The variants (m.2706A > G, m.3010A > G, and m.1095T > C), which are associated with drug-induced toxicity, were observed with the frequency of 99.15%, 30.51%, and 0.08%, respectively. The m.2150T > A, a genotype associated with highly disruptive effects on mitochondrial ribosomes, was identified in five subjects. The D and M groups were the most dominant groups with the frequency of 34.74% and 16.1%, respectively. Conclusions Our finding was consistent with Korean Genome Project and well reflected the unique profile of mitochondrial haplogroup distribution. It was the first study to annotate the whole mitochondrial genome with drug-induced toxicity to predict the ADRs event in clinical implementation for Korean subjects. This approach could be extended for further study for validation of the potential ethnic-specific mitochondrial genetic biomarkers in the Korean population.

Novel deletion of exon 3 in TYR gene causing Oculocutaneous albinism 1B in an Indian family along with intellectual disability associated with chromosomal copy number variations

Background Oculocutaneous albinism (OCA) is an autosomal recessive disorder characterized by hypo-pigmentation of skin, hair, and eyes. The OCA clinical presentation is due to a deficiency of melanin biosynthesis. Intellectual disability (ID) in OCA cases is a rare clinical presentation and appropriate diagnosis of ID is challenging through clinical examination. We report an Indian family with a rare co-inheritance of OCA1B and ID due to a novel TYR gene variant and chromosomal copy number variations. Methods We have done a study on three siblings (2 males and 1 female) of a family where all of them presented with hypopigmented skin, hair and eyes. The male children and their father was affected with ID. Targeted exome sequencing and multiplex ligation-dependent probe amplification analysis were carried out to identify the OCA1B and ID associated genomic changes. Further Array-CGH was performed using SurePrint G3 Human CGH + SNP, 8*60 K array. Results A rare homozygous deletion of exon 3 in TYR gene causing OCA1B was identified in all three children. The parents were found to be heterozygous carriers. The Array-CGH analysis revealed paternally inherited heterozygous deletion(1.9 MB) of 15q11.1-> 15q11.2 region in all three children. Additionally, paternally inherited heterozygous deletion(2.6 MB)of 10q23.2-> 10q23.31 region was identified in the first male child; this may be associated with ID as the father and the child both presented with ID. While the 2nd male child had a denovo duplication of 13q31.1-> 13q31.3 chromosomal region. Conclusion A rare homozygous TYR gene exon 3 deletion in the present study is the cause of OCA1B in all three children, and the additional copy number variations are associated with the ID. The study highlights the importance of combinational genetic approaches for diagnosing two different co-inherited disorders (OCA and ID). Hence, OCA cases with additional clinical presentation need to be studied in-depth forthe appropriate management of the disease.

A genomic mutation spectrum of collecting duct carcinoma in the Chinese population

Background Renal collecting duct carcinoma (CDC) is a rare and lethal subtype of renal cell carcinoma (RCC). The genomic profile of the Chinese population with CDC remains unclear. In addition, clinical treatments are contradictory. In this study, we aimed to identify the genomic mutation spectrum of CDC in the Chinese population. Methods Whole-exome sequencing was performed using the Illumina Novaseq™ 6000 platform. MuTect2 detects single-nucleotide variants (SNVs) and small scale insertions/deletions (INDELs). The identified mutations were annotated with ANNOVAR and validated by Sanger sequencing. Control-FREEC was used to detect copy number variation (CNV), and GISTIC was applied to detect frequently mutated altered regions. These data were compared with associated The Cancer Genome Atlas cohorts. Results Ten normal-matched CDC patients were included. The mean tumour mutation burden was 1.37 Mut/Mb. Six new recurrent somatic mutated genes were identified, including RBM14, MTUS1, GAK, DST, RNF213 and XIRP2 (20% and 2 of 10, respectively), and validated by Sanger sequencing. In terms of common mutated genes, SETD2 was altered in both CDC and other RCC subtypes but not in bladder urothelial carcinoma (BLCA); CDKN2A was a driver gene in both CDC (SNV: 10%, 1 of 10) and BLCA but not in other RCC subtypes. Next, 29 amplifications and 6 deletions of recurrent focal somatic CNVs were identified by GISTIC2.0, which displayed differences from kidney renal clear cell carcinoma (KIRC), kidney renal papillary cell carcinoma (KIRP) and BLCA cohorts. Of note, CDKN2A (CNV alteration: 30%, 3 of 10) and CDKN2A-AS1 were the only overlapping genes of these four cohorts. Importantly, the CDKN2A mutation in our cohort differed from previous studies in urinary carcinomas. Moreover, CDKN2A-altered cases had significantly worse overall survival than wild-type cases in both KIRC and KIRP cohorts. In addition, the most frequently altered genomic pathway of our CDC cohort was the CDKN2A-mediated p53/RB1 pathway. Conclusions Our study offers the first genomic spectrum of the Chinese population with CDC, which differs from that of the Western population. The altered CDKN2A-mediated p53/RB1 pathway might provide new insight into potential therapeutic targets for CDC patients.