Research Article - (2020) Volume 9, Issue 2

Researching Novel Variants in Endometriosis using Next Generation Sequencing Variant Analysis
Shambo Datta1 and Manisha Rana2*
 
1Department of Chemistry, University of Lyon, Lyon, France
2Department of Forensic Science, Amity University, Uttar Pradesh, India
 
*Correspondence: Manisha Rana, Department of Forensic Science, Amity University, Uttar Pradesh, India, Tel: 011-24363012, Email:

Received: 17-Aug-2020 Published: 07-Aug-2020, DOI: 10.35248/2329-6682.20.9.153

Abstract

Endometriosis is characterized as the presence of ectopic endometrial tissue outside of the uterine, most generally in the ovaries and peritoneum. It is an illness that is impacted by various elements. It is additionally a typical gynaecological confusion and influences roughly 10-15% of all women of regenerative age. Later molecular and pathological examinations demonstrate that endometriosis may fill in as an antecedent of ovarian malignant growth (endometriosis associated ovarian disease, EAOC), especially endometrioid furthermore, clear cell ovarian malignant growths. Albeit histological and epidemiological investigations have shown that endometriosis has a malignant potential, the molecular component that underlies the harmful change of endometriosis is as yet questionable, and the exact component of carcinogenesis must be completely illustrated. At present, the advancement and improvement of another sequencing innovation, next-generation sequencing (NGS), has been progressively significant in malignant growth genomics examine. Lately, NGS has likewise been used in clinical oncology to propel the customized treatment of malignancy. Also, the affectability, speed, and cost make NGS a profoundly alluring stage contrasted with other sequencing modalities. Thus, NGS may lead to the recognizable proof of driver mutations and fundamental pathways related with EAOC. Our sole motivation behind the study was to decipher new variants if any and report any unreported variants identified with genes. We have performed a variant analysis investigation with the assistance of Next Generation Sequencing GALAXY device accessible on the web.

Keywords

Endometriosis; Next Generation Sequencing; Variant analysis; Novel variants

Introduction

Endometriosis (E) is benign gynaecological condition, debilitating, estrogens-subordinate, progesterone-safe, inflammatory issue related with pelvic pain and infertility, with endometrial (uterine covering)-like tissue present outside the uterus (Giudice). By retrograde menstruation, endometrial tissue cells are transplanted to the pelvis (Sampson) where they set up a blood supply, react to cyclic hormones, develop, attack encompassing structures, progress toward becoming innervated (Berkley, et al.; Tokushige, et al.), and inspire a nearby inflammatory reaction and scarring (Giudice LC) [1,2].

Endometriosis influences 5%–10% of regenerative age women (Eskenazi and Warner) [3] and half of women with pelvic pain as well as infertility (>100 million women around the world) (Meuleman et al.) [4] and is a noteworthy reason for inability and bargained personal satisfaction (Sasson and Taylor; Anglesio, et al.) [5-7]. Pelvic, lower stomach and back pain, and urinary and gastrointestinal indications make diagnosis challenging, on the grounds that numerous indications are nonspecific or are related with different disorders (Giudice) [8]. Pelvic inflammation and nerve invasion result in pain (Berkley, et al.; Tokushige, et al.) [9,10], and infertility is expected to ovulatory dysfunction, poor egg quality, unusual (progesterone-safe) uterine endometrium, and bargained embryo implantation (Giudice; Bulun SE) [1,-11]. The meaning of endometriosis is histological and requires the distinguishing proof of the presence of endometrial organ and stroma-like tissue outside the uterus (Sourial) [12]. A few hypotheses have shown that the histogenesis of endometriosis is that emanating streams retrograde through the lumen of the fallopian tubes into the pelvic-peritoneal depressions at feminine cycle (Robboy and Bean; Robboy, et al.) [13,14].

Moreover, it can create distant foci through expansion, connection, and intrusion of endometrial glandular epithelial tissue to distant organs (Somigliana, et al.) [15]. The most normally influenced parts of the body incorporate the ovaries, fallopian tubes, bladder, rectosigmoid colon, and myometrium (Giudice; Pavone and Lyttle) [1,16]. Another hypothesis, the coelomic metaplasia hypothesis, recommends that endometriosis emerges from the metaplasia of cells that line the instinctive and stomach peritoneum following hormonal, ecological, or irresistible incitement (Overton, et al.) [17]. A later hypothesis underpins stem/ progenitor cells and bone marrow-determined immature microorganisms in the pathogenesis of endometriosis (Sasson and Taylor) [5]. Notwithstanding, Anglesio. recognized substantial malignant growth driver mutations in the glandular epithelium of deep infiltrating endometriosis sores, and the authors recommended that the undifferentiated cell related hypothesis requires extra investigations to affirm the rational of a speculation (Anglesio, et al.) [7]. Additionally, (Noe, et al.) recognized 19 mutations enhanced in epithelial however not in stromal sores utilizing bead advanced PCR innovation [18]. The authors proposed another theory that epithelial and stromal segments in creating endometriotic sores co-create from independent ancestors. As right on time as 1925, Sampson proposed a potential relationship among's endometriosis and malignant change (Sampson) [19]. Czernobilsky and Morris portrayed a "middle stage" in the harmful change alluded to as "atypical endometriosis"; it is as of now characterized by the level of dysplastic histologic atypia (Czernobilsky and Morris) [20]. Quite, endometriosis is viewed as a potential pre-intrusive sore and is as of now named a tumor-like sore under the World Health Organization (WHO) histologic arrangement of ovarian tumors. Lately, Tsai et al. (Tai et al.) [21] demonstrated that patients with pelvic incendiary ailment had a three-fold increment in the danger of creating endometriosis dependent on the National Health Insurance Research Database (NHIRD) of Taiwan. The hidden mechanism of endometrisis might be related with three unique procedures:

Endometriosis pieces move from the uterus through the fallopian tubes amid retrograde feminine cycle, spreading these endometriosis sections to the peritoneal depression and embedding on the serosal surface. Metaplasia of the coelom and Vascular and lymphatic metastatic spread (Sasson and Taylor; Anglesio; Bulun; Sampson; Sampson; Figueira) [5,6,11,19,22].

Risk factors and Etiology of endometriosis

Huge hazard factors for the development of endometriosis incorporate conditions that increase the odds of retrograde menstruation and hereditary/genetic factors. Hazard factors for endometriosis incorporate early menarche, nulliparity, broken uterine bleeding, variant estrogen levels (Darrow; Signorello, et al.; Cramer, et al.; Candiani, et al.), and low weight record (Signorello, et al.) [23-26]. Factors, for example, sufficient exercise might be precaution against development of endometriosis (Kvaskoff) [27]. It is realized that the occurrence of endometriosis in women with first-degree relatives who likewise have the ailment might be up to multiple times higher than that of the all inclusive population (Matalliotakis; Treloar) [28,29]. There is probably going to be a multifactorial hereditary inclination for endometriosis, and genome-wide association studies (GWAS) have shown singlenucleotide polymorphism (SNP) profiles which may expand the danger of endometriosis in people (Rahmioglu) [30]. In 2012, Nyholt et al (Nyholt, et al. 2012] [31] distinguished 18 genomic areas harboring 38 putative endometriosis-related SNPs in a GWAS including 4,604 instances of endometriosis.

Among the huge aberrations distinguished were SNPs related with the WNT4 gene, known to be critical in reproductive tract differentiation and advancement in mammalian females (Jaaskelainen; Vainio et al.) [32,33] just as steroidigenesis (Boyer A et al.) [34], VEZT, appeared to be down regulated in gastric diseases (Guo X et al.) [35], and GREB1, an estrogen-managed gene not yet been elucidated. The various elements detailed in the pathogenesis of endometriosis-related ovarian malignancy are outlined in (Figure 1).

gene-technology-pathogenesis

Figure 1: Pathogenesis of EAOC.

Development of EAOC from endometriosis

The idea that endometriosis is the forerunner lesion of some ovarian malignant growth subtypes has been upheld by various lines of examination. The affiliation was noted by pathological techniques, however epidemiological, and hereditary examinations have been important (Sampson, et al.; Sampson, et al.; Vercellini.; LaGrenade and Silverberg, et al.; Fukunaga, et al.; Pearce, et al.; Jiang, et al.; Scott; Lu, et al.; Prowse, et al.; McMeekin, et al.; Sainz de la Cuesta, et al.) [73,79,80-88]. Jiang et al depicted a portion of the principal contemplates recommending a molecular basis connecting endometriosis with cancer development in 1998. They exhibited a similar loss of heterozygosity (LOH) occasions in endometriosis lesion and contiguous endometrioid ovarian malignant growths in 82% of cases inspected (n=11) (Jiang, et al.) [83]. Comparable proof was accounted for by Prowse et al in 2006, who exhibited normal LOH occasions in both endometrioid and clear cell OCs and their related endometriosis lesion, including both nearby and contralateral endometriosis (Prowse, et al.) [86]. Moreover, LOH bringing about PTEN loss might be an early driver occasion in the beginning of in EAOC from endometriosis (Worley, et al.; Sato, et al.) [89,90]. Throughout the most recent 7 years, sequencing and immunohistochemical research have given corroborative proof that changes found in endometriosis-related malignant growths are found in adjoining endometriosis. These sequencing examines unmistakably exhibit a clonal connection among benign and malignant partners affirming that the malignant growths have actuality emerged from the endometriotic lesions (Stamp, et al.; Anglesio, et al.; Wiegand, et al.; Chene, et al.) [9194]. Somatic mutations and other genomic deviations are found in endometriosis that have been embroiled in the advancement of cancer. Mutations in TP53 (Bischoff, et al.; Sainz de la Cuesta, et al.) [95,96] KRAS (Anglesio, et al.; Vestergaard, et al.) [7,97], PTEN (Sato, et al.), PIK3CA (Laudanski, et al.; Yamamoto, et al.) [98,99], and ARID1A gene locales (Anglesio, et al.) have been portrayed. Loss of expression of mismatch repair proteins (Grassi, et al.) [100], microsatellite precariousness (Fuseya, et al.) [101], and tissue-explicit gene copy number changes (Yang, et al. 2013; Mafra, et al.) [102,103], may likewise be found in endometriosis sores. LOH in endometriosis at known oncogenic loci is additionally habitually observed (Sato, et al.; Ali-Fehmi, et al.; Xu, et al.; Obata and Hoshiai, et al.; Thomas and Campbell, et al.; Jiang, et al.; Silveira, et al.) [83,90,104-108]. SNPs that are related with oncogenic change (seen in GWAS datasets) have been recognized in instances of endometriosis (Nyholt, et al.; Albertsen, et al.; Uno, et al.; Painter, et al.) [31,38-40]. In 2015, a meta-investigation detailed by Lee et al including more than 15,000 ovarian disease patients, assessed the 38 putative endometriosis- related SNPs distinguished by Nyholt in 2012 (Nyholt, et al.). Eight of these were related with critical hazard for ovarian malignancy (rs7515106, rs7521902, rs742356, rs4858692, rs1603995, rs4241991, rs6907340, and rs10777670) (Lee, et al.) [109]. Likewise in 2015, Lu et al exhibited shared hereditary hazard among endometriosis and epithelial ovarian malignancy, especially clear- cell and endometrioid histotypes utilizing genome wide affiliation (GWAS) datasets (Lu, et al.) [85].

ARID1A is a tumor silencer gene that was observed to be transformed in an extensive number of EAOC (Wiegand, et al.) [93]. Examiners were initially eager to find that up to 42–61% of CCC and 21–33% EnOC show loss of the comparing ARID1A gene protein articulation (BAF250a) on IHC (Stamp, et al.; Wiegand, et al.; Yamamoto, et al.) [91,93,110]. ARID1A manages essential celularl capacities (expansion and genomic stability) as a tumor silencer gene; along these lines, it was believed that it may play a role in the change of endometriosis to malignancy (Wu, et al.) [111]. In 2015, Anglesio et al showed that clear cell ovarian carcinomas imparted numerous transformations to related simultaneous endometriosis sores, incorporating mutations in ARID1A. Shared transformations in PIK3CA were additionally distinguished among endometriosis and clear-cell sores, an occasion happening in early movement components in other malignant growth types (Anglesio, et al.). This investigation unmistakably exhibited depicted transformations in coterminous endometriosis shared by EAOC, and even some distant sores contained the equivalent (PIK3CA and ARID1A) transformations. Studies looking at BAF250a expression by IHC demonstrate that in simply over half of the announced instances of EAOC, loss of BAF250a expression is seen most of the time (67–80%) in regions of coterminous endometriosis or atypical endometriosis, and that lost Baf250a protein expression appeared to be an early molecular occasion in the advancement of Baf250a-negative EAOC (Stamp, et al.; Chene, et al.; Nishikimi, et al.) [91,94,112]. Strangely, ARID1A transformations are not adequate all alone to cause malignancy (Guan, et al.) [113]. In help of this perception, Borrelli et al portrayed halfway loss of BAF250a in ordinary endometrium without disease (Borrelli, et al.) [114]. An imperative examination lately announced that that 65% of malignancy causing genomic variations are arbitrary DNA repair anomalies (Tomasetti and Vogelstein, et al.) [115]. Bringing this data into context, one can infer that BAF250a loss in endometriosis could speak to an EAOC antecedent sore; nonetheless, ARID1A transformations are neither a fundamental driver transformation nor a critical determinant of the malignant phenotype. The presence of transformations in endometriosis is an indication of more extensive genomic interruption prompting the advancement of EAOC. Figure 2 demonstrates a schematic of the foundation and development of endometriosis sores to EAOC. Investigations have been finished looking at patient results in EAOC dependent on the presence or absence of BAF250a expression. In view of the accessible proof, it still can't seem to be resolved concerning whether there are contrasts in visualization or treatment results identified with BAF250a loss in EAOC (Katagiri, et al.; Lowery, et al.) [116,117]. There are couple of recognizable proteomic changes in a board of proteins assessed by reverse phase protein array (RPPA) recommending that BAF250a loss does not characterize a particular proteomic signature (Wiegand, et al.) [118]. Moreover, the presence or absence of an endometriosis antecedent sore in EAOC has not been related with change in overall disease result (Minlikeeva, et al.) [119].

gene-technology-heliotrope

Figure 2: Erythema in heliotrope.

Endometriosis as neoplasm

Deep infiltrating endometriosis is an intriguing uncommon subtype of endometriosis which was lately exposed to genomic assessment. Deep endometriosis has a penchant to locally attack encompassing structures (entrail, bladder, ureter) yet seldom metastasises. Anglesio et al showed the presence of somatic mutation occasions in 79% of 24 cases, with 26% of all cases screened harbouring measurably somatic mutation in known malignant growth driver genes, for example, KRAS, PIK3CA, ARID1A, and PPP2R1A. In the investigation of a littler subset of tests, mutations in KRAS observed to be available in the epithelial part of endometriosis sores were missing in the stroma. Moreover, one patient was found to have the equivalent KRAS transformation in three spatially unmistakable endometriosis sores. While these molecular occasions are usually found in EAOCs, this investigation showed their essence in deep infiltrating endometriosis. While customarily oncogenic driver transformations (like KRAS) were available in a quarter of tests, they didn't seem to demonstrate the probability of the sore to advance into a gynaecologic malignant growth nor have all the earmarks of being required for the improvement of the deep-infiltrating sores.

This recommends extra or distinctive molecular components might be having an effect on everything in the improvement of endometriosis, and future research utilizing an expansive cluster of molecular advances (epigenetic, grafting deviations, complex chromosomal adjustments, transcriptome, proteome and post- translational changes) to examine the functional science of endometriosis is justified. Novel molecular innovations may likewise help clarify the biology of clonally indistinguishable sores in a similar patient. At long last, the bizarre presence of endometriosis in lymph hubs has been portrayed, with a few cases indicating BAF250a loss (Borrelli, et al.) [114]. Consequently, one may expect that these extremely irregular cases are molecularly particular as they copy locally metastatic malignant growths. Maybe even the deep-infiltrating subtype of endometriosis, which shows unequivocal intrusion of encompassing tissues, might be more fittingly considered a neoplasm than a benign condition. Better comprehension of the molecular pathology of this disease may give helpful procedures to analyze and treat complex cases, with the objective of decreasing morbidity and ailment inconveniences like infertility.

Advanced technologies are revolutionizing the aspects of the athogenesis of endometriosis

The next-generation sequencing (NGS) stage will significantly affect disease diagnosis, management and treatment and anticipating result and reaction (Meldrum, et al.) [120]. NGS innovation is a plausible and solid strategy with that might be utilized to identify novel and uncommon somatic mutations. Also, NGS has been effectively utilized to distinguish germline and somatic mutations in a different of malignancies, including gynecological cancer (Evans, and Matuloni) [121], and it can go about as a diagnostic technique and aiding the customized treatment of malignant growth (Valtcheva, et al.) [122]. What's more, NGS innovation substantially affects precision medication and hazard assessment, including early diagnosis, prognosis, and optimization of treatment choice (Morash, et al.; Fountzilas and Tsimberidou) [123,124]. By performing genomic screening by means of NGS innovation, it is conceivable to distinguish whether a patient has previous hereditary conditions that would make them progressively susceptible to creating malignancy in their lifetime (Meldrum, et al.) [120]. In the ongoing years, NGS has been used to describe genomic alterations in EAOC. A few investigations had shown the utility of NGS in recognizing driver mutations in EAOC patients utilizing whole genome sequencing and target sequencing (Wiegand, et al.; Er, et al.) [93,125]. In our past investigation, ultra-deep (>1000×) target sequencing was performed on 409 cancer related genes to distinguish pathogenic changes related with EAOC, and hopeful genes prescient of threatening change were recognized (Zondervan, et al.) [78]. In light of these discoveries, the recognized driver mutations for benign to premalignant sores could be focuses to control the early diagnosis and avoidance of EAOC. As recently examined, endometriosis is a confusion in which the endometriotic tissue is outside the uterus, andit is commonly thought to be a benign sickness. Also, we realized that NGS or ultra-deep sequencing empowers the revelation of novel sequence variants. (Li, et al.) [126] recently demonstrated that hereditary changes in cyto-skeletal and chromatin re-modelling proteins assume a critical job in the pathogenesis of endometriosis utilizing whole-exome sequencing. Lately, exome sequencing likewise yielded promising discoveries that sores in deep infiltrating endometriosis, which are related with for all intents and purposes no danger of malignant transformation, harbor substantial malignant growth driver mutations (Anglesio, et al.). In spite of the fact that endometriosis is viewed as a benign issue, the consequences of NGS innovation recommend another point of view, that the glandular epithelium of deep infiltrating endometriosis injuries harbor understood malignant related somatic transformations. Suda K et al. (Suda, et al.) [127] distinguished numerous malignant related somatic transformations in epithelial cells from ovarian endometriosis and ordinary endometrium utilizing whole exome sequencing. They affirmed that KRAS and PIK3CA were the most oftentimes transformed genes in endometriotic and ordinary uterine endometrial epithelium tests utilizing target-gene sequencing. They additionally showed that clonal extension of epithelial cells with malignant related somatic transformations prompts the advancement of endometriosis. These discoveries reinforce the past hypothesis that the root of endometriosis happens at the genomic level. Lately, Lac, et al.) [128] distinguished physical somatic driver transformations in incisional endometriosis and profound invading endometriosis utilizing an overly sensitive malignant growth hotspot sequencing board, incorporating hotspot changes in KRAS, ERBB2, PIK3CA and CTNNB1. Taken together, NGS innovation may enable us to grow our insight into the pathogenesis of endometriosis and subvert the traditional hypothesis. These examinations have involved endometriosis as a potential premalignant issue and have demonstrated it might give chances to diagnostics and treatments sooner rather than later. In any case, the impact and job of malignant related transformations in the pathogenesis of endometriosis must be completely clarified.

Materials

The variant analysis was performed on study accession PRJNA326570 where sample SRR3711510 and SRR3711512 were considered as a control sample for rest all 8 samples (SRR3711641, SRR3711642, SRR3711644, SRR3711645, SRR3711646, SRR3711647, SRR3711648 and SRR3711649). For NGS data analysis the library layout was Illumina sequenced. In the Illumina platform, the raw reads produced by the sequencing machine are shown in FASTQ, viewed as the standard design configuration of sequencing reads. The prepared library for the sample is a singleend library.

Methods

Next-generation sequencing is an incredible asset for recognizing uncommon and de novo variations, disease mapping, and evaluating expression levels. For the investigation, NGS reads are first adjusted to a reference genome, and afterward exposed to variant calling after fundamental quality control strategies. The alignment is pivotal for variant calling precision, and BWA is a broadly utilized aligner with great execution. Galaxy system is a web open application for high-throughput genomics, uncovering well known third-party data sources and standard bioinformatics investigation bundles in an incorporated and steady structure, intended to help scholar clients performing reproducible examinations. There is a free open site (http://usegalaxy.org). To import information, we utilized the ENA (European nucleotide document) governs by EMBL (website https://www.ebi.ac.uk/ena). When the file is uploaded from the ENA FASTQ Groomer (Galaxy Tool Version 1.1.1) is performed. It changes over between different FASTQ quality organizations. FASTQ Groomer is open-source toolset was executed in Python and has been coordinated into the online data examination platform Galaxy (Goss, et al.; Nichols et al.) [129,130]. After grooming of the data quality check is done using FASTQC tool. FastQC Read Quality reports (Galaxy Tool Version 0.72) gives quality control keeps an eye on raw sequence data originating from high throughput sequencing pipelines. The report incorporates synopsis charts and tables in an H. T. M. L based configuration. These outcomes got from QC investigations give us adequate data concerning whether the data has any issues or not before continuing forward. For above samples the quality was not good enough to perform mapping, therefore before mapping trimming is performed using TRIMMOMATIC (Galaxy Version 0.36.5) methods. BOWTIE2 (Galaxy Tool Version 2.3.4.2) is utilized to list reference genome which works at rapid and memory proficient way. Bowtie2 is utilized for short read alignment. What makes bowtie2 fascinating is the utilization of almost no RAM with precision and unobtrusive execution in ordering the alignment (Langmead and Salzberg) [131]. The alignment results yield in SAM format (Li, et al.) [132] after mapping to remove PCR duplicates RmDup tool (Galaxy Tool Version 2.0.1) is used and after removing duplicates quality is checked before proceeding. Mpileup (Galaxy Tool Version 2.1.4) reports variants for one or various B.A.M documents. Alignments records are gathered giving one log document (content organization) and other Variant Calling record (V.C.F format) which will give data like probability genotype, position on reads, mapping quality (Blankenber, et al.; Blankenberg, Daniel, et al.; Giardine, Belinda, et al.; Goecks, et al.; Sherry, et al.; Team The Galaxy] [133-136]. Varscan for variants (Galaxy Tool Version 2.4.2) performs variant location for enormously parallel sequencing data, for example, exome, W.G.S, and transcriptome information. It calls variants from M Pileup dataset and produces a Variant Calling File (V.C.F) (Andrew) [137]. Finally, we used wANNOVAR to perform regional and functional annotations. Variant calls are then clarified utilizing Annovar (Wang, et al.) [138]. The comment incorporates the utilization of databases, for example, ClinVar, Exac, dbSNP, and dbNSFP.

Results and Discussion

The novel variants acquired from the outcomes as appeared in the Table 1 which were not recently observed associated with Endometriosis. We found an aggregate of 24 new variations from Shenzhen Second Hospital (Shenzhen, Guangdong, China) endometrium sample. Among which, the majority of the variants got were non-synonymous SNVs, aside from them just a single of the variant (ADRA1B) indicated stop-gain SNP. This could be then additionally considered upon for their jobs in different disease or can be contrasted with different samples for same disease (Table 2).

Sample No. Novel Variants Obtained Type Mutations in Exonic Functions Chromosome Location
SRR3711641 ATP6V0D2 Nonsynonymous SNV Chr8: 86150277
VPS13B Nonsynonymous SNV Chr8: 99859386
GLG1 Nonsynonymous SNV Chr16: 74493027
SRR3711642 LTBP3 Nonsynonymous SNV Chr11: 65540877
CLCN7 Nonsynonymous SNV Chr16: 1474949
MNX1 Nonsynonymous SNV Chr7: 1.57E+08
ADRA1B Stopgain Mutation Chr5: 1.6E+08
SRR3711644 UROD Nonsynonymous SNV Chr1:  45015374
CLCN7 Nonsynonymous SNV Chr16: 1474949
RABGEF1 Nonsynonymous SNV Chr7: 66805250
SRR3711645 ASCL2 Nonsynonymous SNV Chr11: 2269837
DSCAML1 Nonsynonymous SNV Chr11: 1.17E+08
ZNF274 Nonsynonymous SNV Chr19: 58211630
PLXNA1 Nonsynonymous SNV Chr3: 1.27E+08
SEMA6A Nonsynonymous SNV Chr5: 1.16E+08
RABGEF1 Nonsynonymous SNV Chr7: 66805250
SRR3711646 LPR5 Nonsynonymous SNV Chr11: 68410020
MFAP3L Nonsynonymous SNV Chr4: 1.7E+08
GPR22 Nonsynonymous SNV Chr7: 1.07E+08
ATP6V0D2 Nonsynonymous SNV Chr8: 86150277
SRR3711647 DENND3 Nonsynonymous SNV Chr8: 1.41E+08
SRR3711647 USP7 Nonsynonymous SNV Chr16: 8904506
L3MBTL1 Nonsynonymous SNV Chr20: 43534909
MFAP3L Nonsynonymous SNV Chr4: 1.7E+08
ADRA1B Stopgain Mutation Chr5: 1.6E+08
PRICKLE4 Nonsynonymous SNV Chr6: 41786956
SRR3711649 DSCAML1 Nonsynonymous SNV Chr11: 1.17E+08
ARHGAP40 Nonsynonymous SNV Chr20: 38637795
CEBPB Nonsynonymous SNV Chr20: 50191512
SEMA3A Nonsynonymous SNV Chr7: 83961468
VPS13B Nonsynonymous SNV Chr8: 99859386

Table 1: Novel Variants with chromosome location and SNP.

Genes Mutation Frequencies (From intogen)
VPS13B 5.65%
PLXNA1 3.91%
DSCAML1 2.61%
GLG1, LRP5 1.74%
RABGEF1, DENND3, PRICKLE4, SEMA3A, SEMA6A, ADRA1B 1.30%
CLCN7, USP7, MFAP3L 0.87%
L3MBTL1, UROD, GPR22, ATP6V0D2, ZNF274, ARHGAP40 0.43%
LTBP3, MNX1, CEBPB 0%
ASCL2 No Data

Table 2: The frequency for mutation of the novel variants for entometriosis was discovered utilizing Intogen Database (https://www.intogen.org/seek).

The variants involvement in endometriosis was affirmed utilizing Driver: A database for malignancy driver gene (driverdb.tms. cmu.edu.tw/ddbv2/index.php). Four of the variants acquired ARHGAP40, UROD, MNX1 and MFAP3L were appeared to have some association in endometriosis as saw on Driver. The remaining genes are novel and once in a while connected with endometriosis. Majority of the variants obtained showed relativeness in other diseases, apart from endometriosis. The molecular genetics of some of the novel variants is discussed:

LRP5: Gong et al. (2001) demonstrated that LRP5 influences bone mass gathering during development and recognized changes in the LRP5 gene (e.g., 603506.0001) that develop autosomal recessive osteoporosis-pseudoglioma disorder (OPPG; 259770) [139]. They found that obligate bearers of mutant LRP5 gene had decreased bone mass when contrasted with age and sexual orientation coordinated controls. Little et al. (2002) recognized a gly171-to-val transformation in the LRP5 gene (G171V; 603506.0013) that outcomes in an autosomal prevailing high bone mass attribute (see 601884) [140]. Boyden, et al. (2002) found the equivalent LRP5 transformation in a family with autosomal dominant [141], high bone density related with square jaw and torus palatinus. Guo, et al. (2006) genotyped 1,873 Caucasian people from 405 family units for SNPs and haplotypes of the LRP5 gene and found that the regular allele A for SNP4 (rs4988300) and the minor allele G for SNP6 (rs634008) were essentially connected with obesity and body mass index (BMI) [142]. Critical affiliations were additionally seen between the regular haplotype A-G-G-G in block 2 (intron 1) with obesity, BMI, and fat mass (p under 0.001, p under 0.001, and p=0.003, individually). Guo et al. (2006) inferred that intronic variations of the LRP5 gene are particularly connected with weight. In affected people from 4 irrelevant families with polycystic liver disease-4 with or without kidney cysts (PCLD4; 617875), Cnossen et al. (2014) recognized 4 diverse heterozygous missense transformations in the in the LRP5 gene (603506.0035- 603506.0038) [143]. Two transformations influenced the intracellular domain, and 2 influenced the extracellular domian. The transformation in the main family was found by whole exome sequencing and affirmed by Sanger sequencing; the 3 different transformations were found by direct sequencing of the LRP5 gene in a cohort of 150 probands with cystic liver disease. The transformations isolated with the turmoil in the families, with some proof for age-subordinate deficient penetrance. None of the patients conveying transformations had proof of clinical highlights of other LRP5-related disease, including bone density or ocular abnormalities.

CLCN7: In light of the closeness between the phenotype of patients with childish harmful osteopetrosis (see OPTB4; 611490) which create serious osteopetrosis and retinal degeneration, Kornak et al. (2001) hunt down transformations in the human CLCN7 gene in 12 patients with juvenile osteopetrosis [144]. They recognized compound heterozygosity for a nonsense (Q555X; 602727.0001) and a missense (R762Q; 602727.0002) change in the CLCN7 quality in 1 persistent with the illness who had early visual hindrance. No retinal histology was accessible. Blair et al. (2004) developed CD14 cells from control and 4 osteopetrotic human subjects within the sight of bone and analysed their osteoclastic separation in vitro [145]. The osteopetrotic cells indicated absconds in acid transport, natural framework evacuation, and cell fusion with inadequate connection compared with the ordinary cells. Genotype investigation demonstrated that cells from 2 patients compound heterozygous for TCIRG1 (604592) transformations had acid transport defects, though cells from 1 patient compound heterozygous for CLCN7 transformation had natural framework evacuation defects. The cells with a connection defect were from a patient who needed TCIRG1 and CLCN7 transformations. In affected people from 12 disconnected families with autosomal prevailing osteopetrosis-2 (OPTA2; 166600), Cleiren et al. (2001) distinguished heterozygosity for 7 unique transformations in the CLCN7 gene (see, e.g., 602727.0004 and 602727.0005) [146]. Examination of microsatellite markers showed that the changes emerged autonomously in every family. Among these families was the Danish family that Van Hul et al. (1997) at first connected to chromosome 1p21. Also, Cleiren et al. (2001) distinguished 1 patient with the extreme autosomal recessive puerile type of osteopetrosis (OPTB4) who was homozygous for a CLCN7 missense transformation (L766P; 602727.0003), for which her asymptomatic guardians were heterozygous [146].

UROD: In the UROD cDNA from a patient with familial porphyria cutanea tarda (PCT; 176100), Garey et al. (1989) showed a heterozygous gly281-to-val substitution (G281V; 613521.0001). The change was not distinguished in affected people from 7 other PCT families with an autosomal dominant pattern of legacy. In a Tunisian family with hepatoerythropoietic porphyria (HEP; see 176100), de Verneuil et al. identified homozygosity for a G281E change (613521.0002) in the UROD gene product [147].

SEMA3A: In 2 sibs and their dad with Kallmann disorder (HH16; 614897), Young et al. distinguished heterozygosity for a 213-kb cancellation in the SEMA3A gene (603961.0001). Sequencing of the nondeleted SEMA3A allele and of 12 known HH-related gene in affected individuals from the family did not reveal some other transformations. Youthful et al. reasoned that SEMA3A play a job in anosmic hypogonadotropic hypogonadism.

USP7: In a 13-year-old young lady with formative deferral, hypotonia, and seizures, Hao et al. distinguished a once more heterozygous c.429C-G transversion in the USP7 gene, bringing about a tyr143-to-ter (Y143X) substitution and anticipated to result in haplo insufficiency. Direct utilitarian investigations of the variation and investigations of patient cells were not performed [148]. Be that as it may, in vitro knockdown of USP7 in cells brought about a diminishing in TRIM27 (602165) protein levels and impeded endosomal protein reusing with diminished F-actin collection. Hao et al. announced 6 random kids with variable neuro developmental issue related with de novo heterozygous micro deletions of chromosome 16p13.2 and 1 patient with a new heterozygous truncating variation in the USP7 gene (602519) on chromosome 16p13.2. All had formative postponement and scholarly incapacity, and 5 were determined to have chemical imbalance range issue. Extra regular highlights included seizures (5 patients), cryptorchidism or micro penis (in 4 of 5 guys), hypotonia (4 patients), and aggressive conduct (4 patients). Different highlights included gentle nonspecific dysmorphic highlights and poor or missing speech with speech apraxia. Practically all patients were in a specialized curriculum.

LTBP3: In affected individuals from a consanguineous Pakistani family with specific tooth agenesis and short stature (DASS; 601216), Noor et al. distinguished a homozygous nonsense transformation in the LTBP3 gene (Y744X; 602090.0001). Two affected guys were analyzed in detail [149]. The phenotype was described by absence of a large number of the perpetual teeth, just as obvious expanded bone density in the spine and skull base. The discoveries proposed an essential job for LTBP3-intervened transcription being developed of the axial skeleton. In a mother and her 2 children who indicated highlights reliable with mellow geleophysic dysplasia (GPHYSD3; 617809), McInerney-Leo et al. recognized heterozygosity for a missense transformation in the LTBP3 gene (S696C; 602090.0008) [150]. In 2 inconsequential young men determined to have geleophysic dysplasia, who kicked the bucket in early youth from respiratory failure, McInerney- Leo et al. recognized heterozygosity for a stop-loss transformation (602090.0009) and a splice site transformation (602090.0010) in LTBP3, respectively.

VPS13B: In a 33-year-elderly person who showed the typical facial gestalt of Cohen disorder and had neutropenia and retinopathy, yet who did not show truncal stoutness or mental impediment, Gueneau et al. distinguished compound heterozygosity for 2 splice site transformations in the VPS13Bgene (607817.0014; 607817.0015) [151]. The authors proposed that a dose impact of remaining typical VPS13B protein may clarify the deficient phenotype in this patient. In 2 Lebanese siblings with Cohen disorder and the extra highlights of cutis verticis gyrata and sensorineural deafness, initially announced by Megarbane et al. as an unmistakable disorder, Megarbane et al. recognized a homozygous grafting transformation in the VPS13B gene (607817.0016) [152].

MNX1: In 2 predominantly acquired sacral agenesis families, Lynch et al. discovered linkage to 7q36 markers. Ross et al. refined the sub chromosomal confinement in a few extra inherited sacral agenesis families and recognized causative transformations in the MNX1 gene (142994.0001-142994.0006) [153]. In affected individuals from a 3-age family isolating Currarino disorder, Urioste et al. identified a frameshift transformation in the MNX1 gene (142994.0009). Malignant mutation of a presacral teratoma was seen in the 22-year-old proband, and presacral teratomas were found in 6 other relatives, including the 3 asymptomatic people. Of 9 influenced individuals, just 2 showed the total set of three. In affected individuals from a 4-age family with Currarino disorder, Wang et al. (2006) recognized heterozygosity for a nonsense transformation in the MNX1 gene (142994.0010) [154].

ADRA1B: The distal end of 5q, 5q31.1-qter, contains the genes for 2 adrenergic receptors, ADRB2 (109690) and ADRA1B, and the dopamine receptor type 1A gene (DRD1A; 126449). Krushkal et al. utilized an effective conflicting sib-pair ascertainment plan to examine the effect of this area of the genome on variety in systolic blood pressure in youthful Caucasians [155]. They quantified 8 exceedingly polymorphic markers crossing this positional applicant gene rich district in 427 people from 55 3-age families containing 69 conflicting sib-pair, and determined multipoint character by plunge probabilities. The after effects of hereditary linkage and affiliation tests showed that the district between markers D5S2093 and D5S462 was altogether connected to at least 1 polymorphic genes influencing inter individual variety in systolic blood pressure. Since the ADRA1B and DRD1A genes are found near these markers, the information recommended that hereditary variety in 1 or both of these G protein-coupled receptors, which partake in the control of vascular tone, assumes an essential job in affecting inter individual variety in systolic blood pressure levels (Table 3).

Gene Name Expression Associated Cancer (Mutation Frequency)
VPS13B Ubiquitous expression in endometrium (RPKM 3.0) Cutaneous melanoma (12.47%)
PLXNA1 Ubiquitous expression in lung (RPKM 7.9) Lung squamous cell carcinoma (5.71%)
DSCAML1 Biased expression in brain (RPKM 3.3) Cutaneous melanoma (9.76%)
CDK11A Ubiquitous expression in bone marrow (RPKM 22.4) Cutaneous Melanoma (1.08%)
GLG1 Ubiquitous expression in ovary (RPKM 26.0) Small cell lung carcinoma (8.07%)
LRP5 Ubiquitous expression in fat (RPKM 20.9) Cutaneous melanoma (6.78%)
RABGEF1 Ubiquitous expression in bone marrow (RPKM 17.2) Bladder carcinoma(2.04%)
DENND3 Broad expression in bone marrow (RPKM 17.0) Cutaneous melanoma (8.40%)
PRICKLE4 Ubiquitous expression in spleen (RPKM 16.5) Stomach adenocarcinoma (1.86%)
SEMA3A Broad expression in placenta (RPKM 2.5) Bladder carcinoma (4.08%)
SEMA6A Broad expression in adrenal (RPKM 13.5) Lung adenocarcinoma (2.56%)
ADRA1B Biased expression in spleen (RPKM 2.0) Stomach adenocarcinoma (1.86%)
CLCN7 Ubiquitous expression in spleen (RPKM 16.9) Non-small cell lung carcinoma (3.23%)
USP7 Ubiquitous expression in testis (RPKM 31.0) Stomach adenocarcinoma (3.73%)
MFAP3L Broad expression in kidney (RPKM 7.1) Stomach adenocarcinoma (1.86%)
L3MBTL1 Broad expression in testis (RPKM 4.3) Cutaneous melanoma (2.17%)
UROD Ubiquitous expression in bone marrow (RPKM 73.9) Bladder carcinoma (1.08%)
GPR22 Biased expression in heart (RPKM 6.9) Lung adenocarcinoma (1.08%)
ATP6V0D2 Biased expression in kidney (RPKM 22.7) Small cell lung carcinoma (2.90%)
ZNF274 Ubiquitous expression in thyroid (RPKM 9.7) Cutaneous melanoma (2.44%)
ARHGAP40 Biased expression in skin (RPKM 12.2) Acute myeloid leukemia (0.51%)
LTBP3 Ubiquitous expression in ovary (RPKM 27.2) Cutaneous melanoma(3.25%)
MNX1 Biased expression in colon (RPKM 3.3) Lung squamous cell carcinoma(1.15%)
CEBPB No Data Bladder carcinoma(1.02%)
ASCL2 Broad expression in colon (RPKM 4.3) No Data

Table 3: Novel variants with their expression (taken from gene database of NCBI- https://www.ncbi.nlm.nih.gov/gene) and their highest cancer associated mutation frequency (taken from IntOgen- https://www.intogen.org/seek).

Conclusion

The investigation can additionally expand in discovering the job of such novel genes in interaction and metabolic pathways and can additionally be contemplated for DNA-protein interaction investigation to help novel research particularly towards its molecular relationship or interaction of the powerful gene products.

REFERENCES

Citation: Datta S, Rana M (2020) Researching Novel Variants in Endometriosis Using Next Generation Sequencing Variant Analysis. 9:153. DOI: 10.35248/2329-6682.20.9.153.

Copyright: © 2020 Datta S, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.