Search

Background Sequence-based identification is one of the most effective methods for species-level identification of nontuberculous mycobacteria (NTM). However, it is time-consuming because of the bioinformatics processes involved, including sequence trimming, consensus sequence generation, and public database searches. We developed a simple and fully automated software that enabled species-level identification of NTM from trace files, SnackNTM (https://github.com/Young-gonKim/SnackNTM). Methods JAVA programing language was used for software development. The SnackNTM diagnostic algorithm utilized 16S rRNA gene sequences, according to the Clinical & Laboratory Standards Institute guidelines, and an rpoB gene region was adjunctively utilized to narrow down the species. The software performance was validated using trace files of 234 clinical cases, comprising 217 consecutive cases and 17 additionally selected cases of unique species. Results SnackNTM could analyze multiple cases at once, and all the bioinformatics processes required for sequence-based NTM identification were automatically performed with a single mouse click. SnackNTM successfully identified 95.9% (208/217) of consecutive clinical cases, and the results showed 99.0% (206/208) agreement with manual classification results. SnackNTM successfully identified all 17 cases of unique species. In a processing time comparison test, the analysis and reporting of 30 cases, which took 150 minutes manually, took only 40 minutes with SnackNTM. Conclusions SnackNTM is expected to reduce the workload for NTM identification, especially in clinical laboratories that process large numbers of cases.

Background Identifying the causal pathogen of encephalitis remains a clinical challenge. A 50-year-old man without a history of neurological disease was referred to our department for the evaluation of an intracranial lesion observed on brain magnetic resonance imaging (MRI) scans, and the pathology results suggested protozoal infection. We identified the species responsible for encephalitis using thymine-adenine (TA) cloning, suitable for routine clinical practice. Methods We extracted DNA from a paraffin-embedded brain biopsy sample and performed TA cloning using two universal eukaryotic primers targeting the V4-5 and V9 regions of the 18S rRNA gene. The recombinant plasmids were extracted, and the inserted amplicons were identified by Sanger sequencing and a homology search of sequences in the National Center for Biotechnology Information Basic Local Alignment Search Tool. Results The infection was confirmed to be caused by the free-living amoeba Balamuthia mandrillaris. Two of 41 colonies recombinant with 18S V4-5 primers and 35 of 63 colonies recombinant with the 18S V9 primer contained B. mandrillaris genes; all other colonies contained human genes. Pathogen-specific PCR ruled out Entamoeba histolytica, Naegleria fowleri, Acanthamoeba spp., and Toxoplasma gondii infections. Conclusions This is the first report of B. mandrillaris-induced encephalitis in Korea based on molecular identification. TA cloning with the 18S rRNA gene is a feasible and affordable diagnostic tool for the detection of infectious agents of unknown etiology.

Emerging roles of circular RNAs (circRNAs) in various biological processes have advanced our knowledge of transcriptional and posttranscriptional gene regulation. To date, no research has been conducted to explore their roles in the rice–Xanthomonas oryzae pv. oryzae interaction. Therefore, we identified 3,517 circRNAs from rice leaves infected with the highly virulent X. oryzae pv. oryzae strain PXO99A by using rRNA depleted RNA sequencing technique coupled with the CIRI2 and CIRCexplorer2 pipeline. Characterization analyses showed that these circRNAs were distributed across the whole genome of rice, and most circRNAs arose from exons (85.13%), ranged from 200 to 1,000 bp, and were with a noncanonical GT/AG (including CT/AC equivalent) splicing signal. Functional annotation and enrichment analysis of the host genes that produced the differentially expressed circRNAs (DEcircRNAs) suggested that these identified circRNAs might play an important role in reprogramming rice responses to PXO99A invasion, mainly by mediating photorespiration and chloroplast, peroxisome, and diterpenoid biosynthesis. Moreover, 31 DEcircRNAs were predicted to act as microRNA decoys in rice. The expression profile of four DEcircRNAs were validated by quantitative real-time PCR with divergent primers, and the back-splicing sites of seven DEcircRNAs were verified by PCR analysis and Sanger sequencing. Collectively, these results inferred a potential functional role of circRNAs in the regulation of rice immunity and provide novel clues about the molecular mechanisms of rice–PXO99A interaction.

BACKGROUND: Somatic variations in rearranged during transfection (RET) proto-oncogene acts to influence Thyroid cancer (TC) in a low penetrance manner, but their effects tend to vary between different populations. OBJECTIVE: This case-control study was aimed to evaluate effect of RET G691S, S904S and L769L single nucleotide polymorphisms (SNPs) on the risk for differentiated thyroid carcinoma (DTC). METHODS: A total of 180 patients and 220 controls were genotyped by Polymerase chain reaction – restriction fragment length polymorphism (PCR-RFLP). Di-Deoxy Sanger sequencing was performed on 100 samples with variations and 20 wild samples for each amplified exon. In addition, In Silico tools were used to evaluate structural and functional impact of individual SNPs in disease progression. RESULTS: In RET G691S/L769L/S904S SNPs, frequency of variant genotypes in DTC cases was 61.1%, 54.4% and 76.6% as compared to 45.9%, 43.6% and 89.09% in controls respectively (P⩽ 0.05). In Silico analysis revealed that different protein formed due to G691S substitution decreases the stability of 3D structure of protein. The RET G691S and L769L SNP followed “Dominant” but RET S904S SNP confirmed an “Additive” mode of inheritance. CONCLUSION: RET G691S/L769L/S904S SNPs are significantly associated with DTC with G691S SNP declining the stability of final protein product.

Synovial sarcoma is characterized by variable epithelial differentiation and specific SS18-SSX gene fusions. The diagnosis is primarily based on phenotype, but fusion gene detection is increasingly being considered indispensable, with SS18 break-apart fluorescence in situ hybridization (FISH) being favored in many laboratories. However, SS18 FISH assay produces negative or atypical results in a minority of cases, leaving uncertainties in diagnosis and management. Here, we analyzed this challenging subset of SS18 FISH-negative/atypical synovial sarcoma using RNA sequencing and monoclonal antibodies that recognize SS18-SSX and the SSX C-terminus. Among 99 synovial sarcoma cases that were previously subjected to SS18 break-apart FISH, eight cases were reported as negative and three cases were indeterminate, owing to atypical signal patterns. Three of these 11 tumors (two monophasic and one biphasic) harbored novel EWSR1-SSX1 fusions, were negative for SS18-SSX staining, and were positive for SSX C-terminus staining. One monophasic tumor harbored a novel MN1-SSX1 fusion, and showed negative SS18-SSX expression and positive SSX C-terminus staining. Another monophasic tumor carried an SS18L1-SSX1 fusion, and was weakly positive for SS18-SSX, while SMARCB1 expression was reduced. The presence of these novel and/or rare fusions was confirmed using RT-PCR and Sanger sequencing. EWSR1-SSX1 was further validated by EWSR1 FISH assay. The remaining six tumors (five monophasic and one biphasic) showed strong SS18-SSX expression, and RNA sequencing successfully performed in three cases identified canonical SS18-SSX2 fusions. Based on a DNA methylation-based unsupervised clustering, the tumors with EWSR1-SSX1 and SS18L1-SSX1 clustered with synovial sarcoma, while the MN1-SSX1-positive tumor was not co-clustered despite classic histology and immunoprofile. In summary, we discovered novel and rare SSX1 fusions to non-SS18 genes in synovial sarcoma. The expanded genetic landscape carries significant diagnostic implications and advances our understanding of the oncogenic mechanism.

Recurrent respiratory papillomatosis (RRP) is a rare and chronic disease affecting the upper airway with papillomatous lesions caused by the human papillomavirus (HPV) infection, especially HPV-6 and/or HPV-11 types. Little is known about the genetic and epigenetic-drivers in RRP pathophysiology. For this purpose, we analyzed 27 papillomatous lesions from RRP patients to evaluate somatic mutations and methylation status in CDKN2A (p14ARF/p16INK4A) and TP53, which are key tumor suppressor genes for the cell cycle control. Sanger sequencing analysis revealed one somatic mutation in TP53 (c.733_734insA) and four mutations in CDKN2A (c.-30G>T, c.29_30insA, c.69delT, and c.300C>A). These mutations were observed in 10 patients, 6 of which carried double mutation. Furthermore, 50% (5/10) of these patients carrying somatic mutations had RRP severity, representing 62.5% (5/8) of the severity cases in this study, albeit no significant association was found between somatic mutations and disease severity. Methylation-specific PCR (MSP) assays revealed p14ARF promoter hypermethylation in 100% of cases, followed by TP53 (96.3%), and p16INK4A (55.6%), suggesting the influence of HPV in the DNA methylation machinery. In conclusion, somatic mutations were not common events identified in RRP patients. However, epigenetic modulation by high methylation rates, particularly for the p14ARF/TP53 pathway, seems to be in the course of RRP development.

Autosomal primary microcephaly (MCPH) is a heterogenetic disorder that affects brain's cerebral cortex size and leads to a reduction in the cranial vault. Along with the hallmark feature of reduced head circumference, microcephalic patients also exhibit a variable degree of intellectual disability as well. Genetic studies have reported 28 MCPH genes, most of which produce microtubule-associated proteins and are involved in cell division. Herein this study, 14 patients from seven Pashtun origin Pakistani families of primary microcephaly were analyzed. Mutation analysis was performed through targeted Sanger DNA sequencing on the basis of phenotype-linked genetic makeup. Genetic analysis in one family found a novel pathogenic DNA change in the abnormal spindle microtubule assembly (ASPM) gene (NM_018136.4:c.3871dupGA), while the rest of the families revealed recurrent nonsense mutation c.3978GA (p.Trp1326*) in the same gene. The novel reported frameshift insertion presumably truncates the protein p.(Lys1291Glyfs*14) and deletes the N-terminus domains. Identification of novel ASPM-truncating mutation expands the mutational spectrum of the ASPM gene, while mapping of recurrent mutation c.3978GA (p.Trp1326*) will aid in establishing its founder effect in the Khyber Pakhtunkhwa (KPK) inhabitant population of Pakistan and should be suggestively screened for premarital counseling of MCPH susceptible families. Most of the recruited families are related to first-degree consanguinity. Hence, all the family elders were counseled to avoid intrafamilial marriages.

Cenani-Lenz syndrome (CLS) is a rare autosomal-recessive congenital disorder affecting development of distal limbs. It is characterized mainly by syndactyly and/or oligodactyly, renal anomalies, and characteristic facial features. Mutations in the LRP4 gene, located on human chromosome 11p11.2-q13.1, causes the CLS. The gene LRP4 encodes a low-density lipoprotein receptor-related protein-4, which mediates SOST-dependent inhibition of bone formation and Wnt signaling. In the study, presented here, three families of Pakistani origin, segregating CLS in the autosomal recessive manner were clinically and genetically characterized. In two families (A and B), microsatellite-based homozygosity mapping followed by Sanger sequencing identified a novel homozygous missense variant [NM_002334.3: c.295G>C; p.(Asp99His)] in the LRP4 gene. In the third family C, exome sequencing revealed a second novel homozygous missense variant [NM_002334.3: c.1633C>T; p.(Arg545Trp)] in the same gene. To determine the functional relevance of these variants, we tested their ability to inhibit canonical WNT signaling in a luciferase assay. Wild type LRP4 was able to inhibit LRP6-dependent WNT signaling robustly. The two mutants p.(Asp99His) and p.(Arg545Trp) inhibited WNT signaling less effectively, suggesting they reduced LRP4 function.

Background The mitochondrial tRNA translation optimization 1 (MTO1) gene, which is closely related to defective mitochondrial oxidative phosphorylation, is an evolutionarily conserved protein expressed in high energy-demanding tissues and is associated with complex oxidative phosphorylation deficiency type 10 (COXPD10) in humans. Related cases and studies are still scarce and have not been reported in the Chinese region. Materials and methods Detailed clinical assessment was applied to the patient. Based on next-generation sequencing technology, we performed whole-exome sequencing of the patient and the parents. Sanger sequencing was used for validation. Bioinformatics software and protein simulations were used to predict the pathogenicity of the variants. Results The patient was diagnosed with a possible association with mitochondrial disease according to the clinical manifestations and physical examination. A novel frameshift mutation c.344delA (p. Asn115Thrfs*11) and a novel point mutation c.1055C > T (p. Thr352Met) in the MTO1 gene were identified. They were found to cause abnormal changes in amino acids and the protein by biochemical tools, indicating it may be pathogenic. Conclusion We present two novel and possibly pathogenic variants in the MTO1 gene in a Chinese Han family.

European bison are susceptible to a range of pathogens which may influence their health, and hence, to ensure their protection, it is essential to provide effective monitoring of potential exposure. This study presents the first molecular confirmation of Sarcocystis cruzi infection in European bison based on PCR amplification of the cytochrome c oxidase subunit I (cox1) gene. A sample of heart tissue taken from one fifteen-year-old European bison cow was examined by light microscopy for the presence of heart sarcocysts. The genomic DNA isolated from any identified sarcocysts was subjected to PCR to amplify cox1 gene sequences, and the obtained amplicons were sequenced by Sanger dideoxy sequencing. Two partial cox1 sequences were obtained; they were identified as S. cruzi and deposited in the GenBank™ database under the accession numbers MW490605 and MW490606. BLAST analysis found them to demonstrate the closest similarity to S. levinei (MH255771-MH255779 and KU247874-KU247884), sharing an identity of 93.14–93.8 %. This is the first report to identify sarcocysts isolated from heart tissue of infected European bison living in the Białowieża forest to species level using cox1 analysis. Our findings confirm that the European bison is a natural intermediate host for S. cruzi. As such, coordinators of future conservation programmes should consider the impact of these diseases on reintroduced animals., Graphical abstract Image 1, Highlights • Molecular identification targeting cox1 gene confirmed S. cruzi infection in European bison. • The pathogenic species of Sarcocystis in endangered bison from Białowieża forest. • Bison bonasus as a natural intermediate host for Sarcocystis cruzi.

Citrin deficiency is one of the most common inborn errors of metabolism in East Asians, which may manifest as neonatal cholestasis, failure to thrive and dyslipidaemia, or recurrent hyperammonaemic encephalopathy. Its molecular diagnosis requires confirmation of the presence of biallelic pathogenic variants in SLC25A13 gene by sequencing, and analysis for a common insertion IVS16ins3kb. However, patients with compatible biochemical features but only one monoallelic pathogenic variant have remained a diagnostic challenge. Here we report the development, validation and application of a multiplex ligation-dependent probe amplification (MLPA) assay using an in-house oligonucleotide probemix and a customised Coffalyer.NET worksheet for detection of exonic copy number variations in SLC25A13. With this MLPA assay, we successfully identified the presence of a heterozygous exonic deletion in SLC25A13 in three of 15 (20%) unrelated individuals with only one monoallelic pathogenic variant detected using conventional methods. Three exonic deletions, two novel involving exon 14 and one reported involving exon 5, were subsequently confirmed with Sanger sequencing. In summary, we developed, evaluated, and demonstrated the clinical utility of an in-house MLPA assay to look for exonic deletions in SLC25A13 in patients with citrin deficiency. With the discovery of novel deletions, MLPA should be considered a test of choice for molecular diagnosis of citrin deficiency when the sequencing result is inconclusive.

Background Three genes clustered together on chromosome 12 comprise a group of hydroxycarboxylic acid receptors (HCARs): HCAR1, HCAR2, and HCAR3. These paralogous genes encode different G-protein coupled receptors responsible for detecting glycolytic metabolites and controlling fatty acid oxidation. Though better known for regulating lipid metabolism in adipocytes, more recently, HCARs have been functionally associated with breast cancer proliferation/survival; HCAR2 has been described as a tumor suppressor and HCAR1 and HCAR3 as oncogenes. Thus, we sought to identify germline variants in HCAR1, HCAR2, and HCAR3 that could potentially be associated with breast cancer risk. Methods Two different cohorts of breast cancer cases were investigated, the Alabama Hereditary Cancer Cohort and The Cancer Genome Atlas, which were analyzed through nested PCRs/Sanger sequencing and whole-exome sequencing, respectively. All datasets were screened for rare, non-synonymous coding variants. Results Variants were identified in both breast cancer cohorts, some of which appeared to be associated with breast cancer BC risk, including HCAR1 c.58C > G (p.P20A), HCAR2 c.424C > T (p.R142W), HCAR2 c.517_518delinsAC (p.G173T), HCAR2 c.1036A > G (p.M346V), HCAR2 c.1086_1090del (p.P363Nfs*26), HCAR3 c.560G > A (p.R187Q), and HCAR3 c.1117delC (p.Q373Kfs*82). Additionally, HCAR2 c.515C > T (p.S172L), a previously identified loss-of-function variant, was identified. Conclusions Due to the important role of HCARs in breast cancer, it is vital to understand how these genetic variants play a role in breast cancer risk and proliferation and their consequences on treatment strategies. Additional studies will be needed to validate these findings. Nevertheless, the identification of these potentially pathogenic variants supports the need to investigate their functional consequences.

Rotaviruses by virtue of its segmented genome generate numerous genotypes. G1P[8] is the most common genotype reported globally. We intend to identify the evolutionary differences among G1P[8] strains from the study with vaccine strains. Stool samples collected from children

Labeo rohita, one of the Indian major carps, is the most popular culture species in Indian subcontinent due to its consumer preference and delicacy. A selective breeding program for harvest body weight has resulted in an average genetic gain of 17% per generation. Transcriptome resource for this species is scanty. Here, we have characterized the liver and muscle transcriptomes of rohu using Roche 454 GS-FLX next generation sequencing platform. In total, 1.2 million reads were generated, de novo assembly and clustering resulted in 4171 transcripts. Out of these, 4171 had significant blast hit against NCBI nr database, and 2130 transcripts were successfully annotated. In total, 289 SSRs were identified with an identification rate of 5.8%, and dinucleotide repeat motifs were observed to be the most abundant SSRs. Further, 2231 putative SNPs were identified with high confidence. Validation of eight putative SNPs using Sanger sequencing resulted in 100% true SNPs. Significant allelic imbalance of M1, M4 and M5 loci between growth selected and control individual were observed. Furthermore, 13 transcription factors were identified in the present study belonging to six different transcription factor families. The present study demonstrated the utility of RNAseq to develop genomics resources in non-model fish species, and the marker resources developed would support the genetic improvement program of this species.

Primary Sjogren’s syndrome (pSS) is a complex systemic autoimmune disease, which is difficult to accurately diagnose due to symptom diversity in patients, especially at earlier stages. We tried to find potential single nucleotide polymorphisms (SNPs), copy number variations (CNVs) and related signaling pathways. Genomic DNA was extracted from peripheral blood of 12 individuals (7 individuals from 3 pSS pedigrees and 5 sporadic cases) for whole-exome sequencing (WES) analysis. SNPs and CNVs were identified, followed by functional annotation of genes with SNPs and CNVs. Gene expression profile (involving 64 normal controls and 166 cases) was downloaded from the Gene Expression Omnibus database (GEO) dataset for differentially expression analysis. Sanger sequencing and in vitro validation was used to validate the identified SNPs and differentially expressed genes, respectively. A total of 5 SNPs were identified in both pedigrees and sporadic cases, such as FES, PPM1J, and TRAPPC9. A total of 3402 and 19 CNVs were identified in pedigrees and sporadic cases, respectively. Fifty-one differentially expressed genes were associated with immunity, such as BATF3, LAP3, BATF2, PARP9, and IL15RA. AMPK signaling pathway and cell adhesion molecules (CAMs) were the most significantly enriched signaling pathways of identified SNPs. Identified CNVs were associated with systemic lupus erythematosus, mineral absorption, and HTLV-I infection. IL2-STAT5 signaling, interferon-gamma response, and interferon-alpha response were significantly enriched immune related signaling pathways of identified differentially expressed genes. In conclusion, our study found some potential SNPs, CNVs, and related signaling pathways, which could be useful in understanding the pathological mechanism of pSS.

This study aims to identify the underlying genetic cause of a Chinese patient with Leber congenital amaurosis (LCA). Detailed clinical data and family history were collected. A medical diagnostic panel sequencing covering 4450 genes was conducted. Two candidate disease-causing mutations detected in CEP290 were then validated with Sanger sequencing and bioinformatic analysis. Reverse transcription polymerase chain reaction (RT-PCR) and cDNA sequencing were performed to understand the effect of the novel CEP290 mutation on CEP290 mRNA splicing. A five-month-old LCA patient with both parents was enrolled. Medical diagnostic panel sequencing revealed that the patient is a compound heterozygote for two potentially pathogenic CEP290 mutations. Among them, c.1666dupA (p.I556NfsX20) was previously reported and has a significant association with LCA phenotype. A novel CEP290 mutation (c.3310–1_3313delGCTTA) was confirmed in both the patient and her father. RT-PCR and cDNA sequencing confirmed that the novel mutation affects the CEP290 mRNA splicing and results in a complete skipping of exon 29. The frameshift resulted in an early stop codon and truncation of the mutant protein by 1371 amino acid residues (p.L1104fsX6). Our report provided a new mutation to the spectrum of CEP290-related diseases. The data suggested that the c.3310–1_3313delGCTTA mutation affects the CEP290 mRNA splicing and the CEP290 protein function. This valuable information is important for future studies on the mRNA splicing of CEP290 and the pathogenesis of CEP290-related diseases.

Purpose To date, many genes have been associated with congenital hypothyroidism (CH). Our aim was to identify the mutational spectrum of 23 causative genes in Turkish patients with permanent CH, including thyroid dysgenesis (TD) and dyshormonogenesis (TDH) cases. Methods A total of 134 patients with permanent CH (130 primary, 4 central) were included. To identify the genetic etiology, we screened 23 candidate genes associated with CH by next-generation sequencing. For confirmation and to detect the status of the specific familial variant in relatives, Sanger sequencing was also performed. Results Possible pathogenic variants were found in 5.2% of patients with TD and in 64.0% of the patients with normal-sized thyroid or goiter. In all patients, variants were most frequently found in TSHR, followed by TPO and TG. The same homozygous TSHB variant (c.162 + 5G > A) was identified in four patients with central CH. In addition, we detected novel variants in the TSHR, TG, SLC26A7, FOXE1, and DUOX2. Conclusion Genetic causes were determined in the majority of CH patients with TDH, however, despite advances in genetics, we were unable to identify the genetic etiology of most CH patients with TD, suggesting the effect of unknown genes or environmental factors. The previous studies and our findings suggest that TSHR and TPO mutations is the main genetic defect of CH in the Turkish population.

Introduction: Congenital stationary night blindness (CSNB) is a rare, largely nonprogressive, inherited retinal disorder that can be clinically classified on the basis of fundus and electroretinogram abnormalities. Methods: We analyzed four large consanguineous families from the Southern Punjab region of Pakistan including multiple individuals affected with CSNB. Exome sequencing was performed in probands of all four families; Sanger sequencing was performed in additional members to test co-segregation of the variants identified. Results: We identified two novel and likely pathogenic variants in two pedigrees, namely, NM_002905.4:c.668A>C (p.Gln223Pro) in RDH5 and NM_022567.2:c.908del (p.Gly303ValfsTer45) in NYX. In the two other families, the variants NM_002905.4:c.319G>C (p.Gly107Arg) in RDH5 and NM_000541.5:c.874C>T (p.Arg292Ter) in SAG were identified. These latter mutations have been reported previously, but not in the Pakistani population. Conclusions: Our findings expand the mutational spectrum of CSNB, in particular within the population of Southern Punjab.

Cleidocranial dysplasia (CCD) is an autosomal dominant hereditary disease associated with the gene RUNX2. Disease-specific induced pluripotent stem cells (iPSCs) have emerged as a useful resource to further study human hereditary diseases such as CCD. In this study, we identified a novel CCD-specific RUNX2 mutation and established iPSCs with this mutation. Biopsies were obtained from familial CCD patients and mutation analyses were performed through Sanger sequencing and next generation sequencing. CCD-specific human iPSCs (CCD-hiPSCs) were established and maintained under completely defined serum, feeder, and integration-free condition using a non-integrating replication-defective Sendai virus vector. We identified the novel mutation RUNX2_c.371C>G and successfully established CCD-hiPSCs. The CCD-hiPSCs inherited the same mutation, possessed pluripotency, and showed the ability to differentiate the three germ layers. We concluded that RUNX2_c.371C>G was likely pathogenic because our results, derived from next generation sequencing, are supported by actual clinical evidence, familial tracing, and genetic data. Thus, we concluded that hiPSCs with a novel CCD-specific RUNX2 mutation are viable as a resource for future studies on CCD.

PURPOSE: Genetic etiology of idiopathic male infertility is enigmatic owing to involvement of multiple gene regulatory networks in spermatogenesis process. Any change in optimal function of the transcription factors involved in this process owing to polymorphisms/mutations may increase the risk of infertility. We investigated polymorphisms/mutations of spermatogenic transcription regulators TAF7 and RFX2 and analysed their association with incidence of azoospermia among the men from West Bengal, India. METHODS: Genotyping was carried by Sanger’s dideoxy sequencing of 130 azoospermic men who were detected negative in Y chromosome microdeletion screening and 140 healthy controls. Association study was done by suitable statistical methods. In silico analysis was performed to infer the intuitive damaging effects of detected variants at transcripts and protein level. RESULTS: We found significant association of TAF7 C16T (MW827584 G > A), RFX2 562delT (MZ560629delA), rs11547633 A > C, rs17606721 A > G, MW827583 C > T, and MZ379836 C > T variants with the incidence of azoospermia. In silico analysis predicted that the variants either alter the natural splice junctions of the transcript or cause probable damage in the structure of proteins of respective genes. CONCLUSION: Polymorphisms/mutations of TAF7 and RFX2 genes increase risk of male infertility in Bengali population. The novel variants may be used as markers for male infertility screening in ART practise. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s10815-021-02352-5.

Introduction GATA1 mutation plays an important role in initiating transient abnormal myelopoiesis (TAM) and in the clonal evolution towards acute megakaryoblastic leukaemia (AMKL) associated with Down syndrome (DS). This study aimed to develop and validate the clinical utility of a complementary DNA (cDNA) analysis in parallel with the conventional genomic DNA (gDNA) Sanger sequencing (Ss), as an initial screening test for GATA1 mutations. Methods GATA1 mutations were evaluated using both gDNA and cDNA in 14 DS patients using Ss and fragment analysis (FA), respectively. Results The detection sensitivity of conventional gDNA sequencing was limited in low blast percentage TAM (LBP-TAM); however, cDNA-based Ss readily detected all the pathognomonic GATA1 mutations. The cDNA-based FA readily detected GATA1 frameshift mutation with a reliable sensitivity ranging from 0.005% to 0.01% of clonal cells. Conclusions GATA1 mutations are heterogeneous; therefore, we would like to propose a dual cDNA and gDNA analysis as a standard diagnostic approach, especially for LBP-TAM. cDNA-based FA promises an excellent sensitivity for detecting frameshift GATA1 mutations in the longitudinal clonal evolution towards AMKL without using a patient specific primer.

Pseudoachondroplasia (PSACH) is an autosomal dominant disorder characterized by rhizomelic short-limbed skeletal dysplasia. The primary clinical and radiographic features include disproportionate dwarfism, joint laxity and hyperextensibility, exaggerated lumbar lordosis, and late ossification of the epiphyses. Identification of disease-causing variants in heterozygous state in COMP establishes the molecular diagnosis of PSACH. We examined 11 families with clinical features suggestive of PSACH. In nine families, we used Sanger sequencing of exons 8-19 of COMP (NM_000095.2) and in two families exome sequencing was used for confirming the diagnosis. We identified 10 de novo variants, including five known variants (c.925G>A, c.976G>A, c.1201G>T, c.1417_1419del, and c.1511G>A) and five variants (c.874T>C, c.1201G>C, c.1309G>A, c.1416_1421delCGACAA, and c.1445A>T) which are not reported outside Indian ethnicity. We hereby report the largest series of individuals with molecular diagnosis of PSACH from India and reiterate the well-known genotype-phenotype corelation in PSACH.

Background: Congenital cataracts are the most common cause of visual impairment worldwide. Inherited cataract is a clinically and genetically heterogeneous disease. Here we report disease-causing variants in a novel gene, CYP21A2, causing autosomal dominant posterior polar cataract. Variants in this gene are known to cause autosomal recessive congenital adrenal hyperplasia (CAH). Methods: Using whole-exome sequencing (WES), we have identified disease-causing sequence variants in two families of British and Irish origin, and in two isolated cases of Asian-Indian and British origin. Bioinformatics analysis confirmed these variants as rare with damaging pathogenicity scores. Segregation was tested within the families using direct Sanger sequencing. Results: A nonsense variant NM_000500.9 c.955 C > T; p.Q319* was identified in CYP21A2 in two families with posterior polar cataract and in an isolated case with unspecified congenital cataract phenotype. This is the same variant previously linked to CAH and identified as Q318* in the literature. We have also identified a rare missense variant NM_000500.9 c.770 T > C; p.M257T in an isolated case with unspecified congenital cataract phenotype. Conclusion: This is the first report of separate sequence variants in CYP21A2 associated with congenital cataract. Our findings extend the genetic basis for congenital cataract and add to the phenotypic spectrum of CYP21A2 variants and particularly the CAH associated Q318* variant. CYP21A2 has a significant role in mineralo- and gluco-corticoid biosynthesis. These findings suggest that CYP21A2 may be important for extra-adrenal biosynthesis of aldosterone and cortisol in the eye lens.

Emerging evidence has demonstrated that circular RNAs (circRNAs) take part in the initiation and development of pancreatic ductal adenocarcinoma (PDA), a deadly neoplasm with an extremely low 5-year survival rate. Reprogrammed glucose metabolism is a key feature of tumour development, including PDA. In this research, we evaluated the role of circRNAs in reprogrammed glucose metabolism in PDA. RNA sequencing under various glucose incubation circumstances was performed. A new circMYOF was identified. Sanger sequencing and RNase R treatment confirmed its circular RNA characteristics. Real-time PCR indicated that it was highly expressed in PDA clinical specimens and cell lines. Gain-of- and loss-of-function assays showed that circMYOF induced progression in PDA. Mechanistically, RNA pull-down and luciferase reporter experiments elucidated that circMYOF, as a competing endogenous RNA for miR-4739, facilitated glycolysis via the VEGFA/PI3K/AKT pathway. Taken together, our findings indicate that circMYOF may work as a desirable biomarker and therapeutic target for PDA patients.

Knowledge of the genetic landscape of a specific population group is vital for population-specific diagnosis and treatment of familial breast cancer. Although BRCA-related diagnostic testing has long been implemented in South Africa, the genotyping approach previously failed for the SA Indian population as it was based on other SA population groups. Because this population is uniquely admixed, the lack of population-specific data resulted in the implementation of comprehensive mutation screens for BRCA1/2. A total of 223 female patients were screened for clinically actionable variants. High-resolution melting analysis (HRMA) was used to screen 88 patients for DNA alterations in the coding and splice site boundaries of BRCA1 exons 2-9, BRCA1 exons 11-23, BRCA2 exons 2-9 and BRCA2 exons 12-27. The protein truncation test (PTT) was used to screen the three larger exons (BRCA1 exon 10 and BRCA2 exons 10 and 11) for protein termination changes. Multiplex ligation-dependent probe amplification (MLPA) was used to determine the presence of larger indels and possible copy number differences. Next Generation Sequencing (NGS) was performed on the remaining 135 samples. All potential variants were confirmed by performing Sanger DNA sequencing. The search revealed 28 different pathogenic heterozygotic variants, together with nine variants of unknown significance (VUS). The results suggested that the SA Indian population represents a different genetic admixture compared to that of mainland India, as only five pathogenic variants corresponded to those reported for mainland India. Familial breast cancer testing for SA Indian patients should therefore be performed as comprehensively as possible as the pathogenic variants seem to be family- rather than population-specific. Furthermore, predictive testing of family members will contribute to relieve the financial burden on the country's healthcare system, as increased surveillance and appropriate management could prevent disease.

In more than one-third of primary immunodeficiency (PID) patients, extensive genetic analysis including whole-exome sequencing (WES) fails to identify the genetic defect. Whole-genome sequencing (WGS) is able to detect variants missed by other genomics platforms, enabling the molecular diagnosis of otherwise unresolved cases. Here, we report two siblings, offspring of consanguineous parents, who experienced similar severe events encompassing early onset of colitis, lymphoproliferation, and hypogammaglobulinemia, typical of lipopolysaccharide-responsive and beige-like anchor (LRBA) or cytotoxic T lymphocyte antigen 4 (CTLA4) deficiencies. Gene-panel sequencing, comparative genomic hybridization (CGH) array, and WES failed to reveal a genetic aberration in relevant genes. WGS of these patients detected a 12.3 kb homozygous tandem duplication that was absent in control cohorts and is predicted to disrupt the reading frame of the LRBA gene. The variant was validated by PCR and Sanger sequencing, demonstrating the presence of the junction between the reference and the tandem-duplicated sequence. Droplet digital PCR (ddPCR) further confirmed the copy number in the unaffected parents (CN = 3, heterozygous) and affected siblings (CN = 4, homozygous), confirming the expected segregation pattern. In cases of suspected inherited immunodeficiency, WGS may reveal a mutation when other methods such as microarray and WES analysis failed to detect an aberration.

Progressive familial intrahepatic cholestasis type 1 (PFIC1) results from biallelic pathogenic variants in ATP8B1. This study sought second pathogenic variants in ATP8B1 by whole-genome sequencing (WGS) in four unrelated low γ-glutamyl transpeptidase cholestasis patients in whom clinical suspicion of PFIC1 was high and gene-panel or Sanger sequencing had identified only one pathogenic variant in ATP8B1. Sanger sequencing confirmed WGS findings and determined the origin of each variant. Novel nonrecurrent structural variants in three patients (patient 1 to patient 3) were identified in trans: g.55396652_55403080del (6427-bp deletion), g.55335906_55346620dup (10,715-bp duplication), and g.55362063_55364293dup (2231-bp duplication). One synonymous variant in patient 4 was recognized in trans (c.1029GA, p. Thr343Thr) and demonstrated as deleterious. In conclusion, WGS improves genetic diagnostic yield in PFIC1. These findings expand the gene-variant spectrum associated with familiar intrahepatic cholestasis 1 (FIC1) disease and for the first time report tandem duplication in ATP8B1 associated with cholestasis.

The potato/tomato psyllid Bactericera cockerelli (Hemiptera: Triozidae) is a pest of Solanaceae plants and a vector of the pathogenic bacterium ‘Candidatus Liberibacter solanacearum’, which is associated with zebra chip disease in potato. This disease is controlled through insecticide treatments, and more environmentally friendly management options are desirable. The objective of this study was to identify viruses present in potato psyllid populations that might be used as biocontrol agents for this insect pest. A new picorna-like virus, tentatively named “Bactericera cockerelli picorna-like virus” (BcPLV), was discovered in B. cockerelli populations maintained in greenhouses, through the use of high-throughput sequencing data and subsequent confirmation by RT-PCR and Sanger sequencing. BcPLV has a positive-sense 9,939-nt RNA genome encoding a single 2,947-aa polyprotein and is related to the Diaphorina citri picorna-like virus (DcPLV) found in Asian citrus psyllid Diaphorina citri populations. Based on their genome organization and the phylogeny of their RNA-dependent RNA polymerase domains, BcPLV and DcPLV together are proposed to comprise a new genus, provisionally named "Psylloidivirus", within the family Iflaviridae.

To compare single-molecule real-time technology (SMRT) and conventional genetic diagnostic technology of rare types of thalassemia mutations, and to analyze the molecular characteristics and phenotypes of rare thalassemia gene variants, we used 434 cases with positive hematology screening as the cohort, then used SMRT technology and conventional gene diagnosis technology [(Gap-PCR, multiple ligation probe amplification technology (MLPA), PCR-reverse dot blot (RDB)] for thalassemia gene screening. Among the 434 enrolled cases, conventional technology identified 318 patients with variants (73.27%) and 116 patients without variants (26.73%), SMRT identified 361 patients with variants (83.18%), and 73 patients without variants (16.82%). The positive detection rate of SMRT was 9.91% higher than conventional technology. Combination of the two methods identified 485 positive alleles among 49 types of variant. The genotypes of 354 cases were concordant between the two methods, while 80 cases were discordant. Among the 80 cases, 76 cases had variants only identified in SMRT method, 3 cases had variants only identified in conventional method, and 1 false positive result by the traditional PCR detection technology. Except the three variants in HS40 and HBG1-HBG2 loci, which was beyond the design of SMRT method in this study, all the other discordant variants identified by SMRT were validated by further Sanger sequencing or MLPA. The hematological phenotypic parameters of 80 discordant cases were also analyzed. SMRT technology increased the positive detection rate of thalassemia genes, and detected rare thalassemia cases with variable phenotypes, which had great significance for clinical thalassemia gene screening.

Congenital absence of the vas deferens (CAVD), a congenital malformation of the male reproductive system, causes obstructive azoospermia and male infertility. Currently, the cystic fibrosis transmembrane conductance regulator (CFTR) has been recognized as the main pathogenic gene in CAVD, with some other genes, such as adhesion G-protein-coupled receptor G2 (ADGRG2), solute carrier family 9 isoform 3 (SLC9A3), sodium channel epithelial 1 subunit beta (SCNN1B), and carbonic anhydrase 12 (CA12), being candidate genes in the pathogenesis of CAVD. However, the frequency and spectrum of these mutations, as well as the pathogenic mechanisms of CAVD, have not been fully investigated. Here, we sequenced all genes with potentially pathogenic mutations using next-generation sequencing and verified all identified variants by Sanger sequencing. Further bioinformatic analysis was performed to predict the pathogenicity of mutations. We described the distribution of the p.V470M, poly-T, and TG-repeat CFTR polymorphisms and identified novel missense mutations in the CFTR and SLC9A3 genes, respectively. Taken together, we identified mutations in the CFTR, ADGRG2, SLC9A3, SCNN1B, and CA12 genes in 22 patients with CAVD, thus broadening the genetic spectrum of Chinese patients with CAVD.

To develop an effective genome editing tool for blueberry breeding, CRISPR-Cas9 and CRISPR-Cas12a were evaluated for their editing efficiencies of a marker gene, beta-glucuronidase (gusA), which was previously introduced into two blueberry cultivars each a single-copy transgene. Four expression vectors were built, with CRISPR-Cas9 and CRISPR-Cas12a each driven by a 35S promoter or AtUbi promoter. Each vector contained two editing sites in the gusA. These four vectors were respectively transformed into the leaf explants of transgenic gusA blueberry and the resulting transgenic calli were induced under hygromycin selection. GUS staining showed that some small proportions of the hygromycin-resistant calli had non-GUS stained sectors, suggesting some possible occurrences of gusA editing. We sequenced GUS amplicons spanning the two editing sites in three blueberry tissues and found about 5.5% amplicons having editing features from the calli transformed with the 35S-Cas9 vector. Further, we conducted a second round of shoot regeneration from leaf explants derived from the initial Cas9- and Cas12a-containing calli (T0) and analyzed amplicons of the target editing region. Of the newly induced shoots, 15.5% for the 35S-Cas9 and 5.3% for the AtUbi-Cas9 showed non-GUS staining, whereas all of the shoots containing the Cas12a vectors showed blue staining. Sanger sequencing confirmed the editing-induced mutations in two representative non-GUS staining lines. Clearly, the second round of regeneration had enriched editing events and enhanced the production of edited shoots. The results and protocol described will be helpful to facilitating high-precision breeding of blueberries using CRISPR Cas technologies. A second round of regeneration enriched editing events and enhanced the production of edited blueberry shoots. The new protocol described facilitates high-precision breeding of blueberries using CRISPR Cas technologies.

Mutations in the SLC12A3 gene have been reported to cause Gitelman syndrome (GS). This study aimed to investigate the genetic mutations and clinical features of patients with GS. Four pedigrees (4 GS patients and 14 family members) were enrolled. The symptoms, laboratory results, management, and genotypes were analyzed. Genomic DNA was screened for gene variations using Sanger sequencing. DNA sequences were compared with reference sequences. The effects of the mutations were predicted using prediction tools (Mutation Taster, PolyPhen-2, SIFT, and PROVEAN). Genetic analysis revealed six genetic variants of SLC12A3, including three novel heterozygous mutations (c.2T > C, c.1609C > T, c.3055G > A) and three previously characterized mutations (c.1456G > A, c.2542G > A, c.1077C > G). These mutations were predicted to exert a damaging effect based on predictive in silico tools. GS patients had low blood pressure and low levels of serum K+, serum Mg2+, and 24-h urinary Ca2+ but high levels of 24-h urinary K+. These clinical manifestations and genotypes were consistent with the diagnostic criteria of GS. The study described the phenotypes and genotypes of 4 pedigrees involving GS patients, demonstrating the importance of SLC12A3 gene screening for GS.

The gold standard for malaria diagnosis is evidence of parasitological confirmation but the traditional method by light microscopy and the routinely used rapid diagnostic tests (RDTs) have limitations. Molecular assays are known to have higher sensitivity and specificity but there are indications that they may also be compromised by genetic variability of the target sequence. The aim of this study therefore, was to evaluate the DNA sequence profile of the diagnostic target of the P. falciparum 18S rRNA PCR assay in field isolates from North-Central Nigeria. Blood samples were collected from 324 children presenting with acute febrile illness suspected to be uncomplicated malaria. Light microscopy and 18S rRNA PCR assay were employed to determine the presence of P. falciparum parasites. Sequence profile of the diagnostic target was evaluated by Sanger sequencing of the PCR products on ABI PRISM® 3100 DNA sequencer (PE Applied Biosystems). Of the 324 children enrolled into this study, 134 (41.4%) were positive for P. falciparum by microscopy while 218 (67.3%) were positive by PCR. The sensitivity of microscopy was 61.47%(95% CI= 57.88% - 69.64%) using the PCR assay as reference standard. The degree of agreement between microscopy and PCR as measured by Cohen's kappa was moderate (κ = 0.502, 95% CI = 0.463 - 0.715).Sequence analysis showed that the DNA target of the P. falciparum 18S rRNA PCR from the field isolates were highly conserved. Only one A>T single nucleotide polymorphism was found within the target sequence among the isolates in this study. This study showed that the DNA target sequence of the18S rRNA PCR assay is highly conserved in field isolates in the study region suggesting little or no impact of selective pressure acting on the locus and has implications for the enhanced sensitivity of the molecular assay.

The pX genetic region of bovine leukemia virus (BLV) includes four genes with overlapping reading frames that code for the Tax, Rex, R3, and G4 proteins. These proteins are involved in the regulation of transcriptional and post-transcriptional viral expression, as well as having oncogenic potential. Our goal was to investigate the pathogenicity of the pX region of BLV genotype 1 in terms of lymphocytosis, lymphomas, and proviral DNA load. We screened 724 serological samples from mixed-age Holstein Friesian cattle from six states in Mexico. Peripheral blood leukocytes (PBLs) were isolated from whole blood with anticoagulant, and genomic DNA was extracted from the PBLs using a commercial kit. Then, a set of primers that hybridize in conserved regions of the BLV pX region were used, which allowed for PCR standardization to detect proviral DNA in infected cells. Positive amplicons were sequenced using the Sanger method, resulting in 1156-nucleotide-long final sequences that included the four pX region genes. The experimental group consisted of 30 animals. Twelve of these had lymphocytosis, six had lymphoma, and 12 were apparently healthy cattle without any signs of lymphocytosis or lymphoma. The presence of lymphoma was detected in six bovine tumor tissues using histopathology, and the presence of BLV was detected by in situ hybridization. Phylogenetic analysis demonstrated that the 30 sequences were associated with genotype 1, and the genetic distance between the sequences ranged from 0.2% to 2.09%. We identified two sequences in the G4 gene: one with a three-nucleotide deletion resulting in the loss of a leucine (AGU_7488L, in a cow with lymphocytosis), and one with a nine-nucleotide deletion resulting in the loss of leucine, proline, and leucine (AGU_18A, in a cow without lymphocytosis). Analysis of the PX region indicated that positive selection had occurred in the G4, rex, and R3 genes, and we found no difference in proviral DNA load between the studied groups. We were unable to establish an association between variations in the pX region and the development of lymphocytosis, lymphoma, asymptomatic status, or proviral DNA load in BLV-infected cattle.

Japanese star anise (Illicium anisatum L., JSA) is seriously damaged by a ringspot disease in Japan. Herein, to determine the causal agent using high-throughput sequencing, we discovered viral RNAs associated with JSA ringspot disease. We then determined the complete or near-complete nucleotide sequences of these RNAs using Sanger sequencing and RACE. The complementary strand of viral RNAs 1, 2, 3, 4, and 5 encoded a single protein, which shared sequence identity with P1 (RNA-dependent RNA polymerase), P2 (glycoprotein precursor), P3 (nucleocapsid protein), P4 (movement protein), and a protein with unknown function of emaraviruses (genus Emaravirus), respectively; however, the highest amino acid sequence identity for the P1–P5 proteins between JSARaV and known emaraviruses was 41.9%, 30.0%, 30.1%, 52.2%, and 38.0%, respectively, all of which were lower than the species demarcation criterion. Furthermore, RNA segments harbored conserved 12-nt terminal sequences at the 5′- and 3′-termini, and a high complementarity of approximately 20 nt in 5′- and 3′-terminal sequences. Transmission electron microscopy confirmed the presence of virus-like particles. JSA ringspot disease was found to be transmitted by an eriophyid mite (subclass Acari, superfamily Eriophyoidea) that belongs to the family Diptilomiopidae. Taken together, these results identified the virus responsible for the ringspot disease of JSA as a new member of the genus, Emaravirus, which we named as the Japanese star anise ringspot-associated virus (JSARaV). Moreover, this is the first report noting that eriophyid mites of the family Diptilomiopidae are capable of transmitting emaravirus.

Here, we report the complete genome sequence of a novel polerovirus, "Plantago asiatica virus A" (PlaVA), detected in Plantago asiatica using high-throughput RNA sequencing and validated by Sanger sequencing. The complete PlaVA genome contains 5,881 nucleotides and has seven open reading frames (ORF0-5 and ORF3a) encoding putative proteins (P0-5 and P3a, respectively) in an arrangement that is similar to that of typical Polerovirus members. Pairwise sequence comparisons revealed that P0 to P5 encoded by PlaVA had the highest sequence identity (25.48%-79.21%) to the corresponding proteins of previously reported poleroviruses. A phylogenetic analysis using the PlaVA P1-2 and P3 amino acid sequences and those of members of the family Solemoviridae (formerly Luteoviridae) indicated that although PlaVA belongs to the genus Polerovirus, it does not represent a known species. Consequently, PlaVA should be considered a member of a new species within the genus Polerovirus.

ETO2 is a nuclear co-repressor, which plays a critical role in the regulation of the cell cycle, self-renewal capacity, and differentiation of hematopoietic progenitor cells. We identified novel fusion transcripts involving ETO2 and CTCF by RNA-seq in a multiple relapsed AML case. The CTCF-ETO2 and ETO2-CTCF chimeric genes were validated by RT-PCR and Sanger sequencing. In addition, both transcripts apparently promoted cell proliferation via JAK/STAT3 pathway that is sensitive to STAT3 inhibitors. The novel fusions may have prognostic value and pathogenic mechanisms in acute myeloid leukemia.

Acinetobacter baumannii (A. baumannii) has developed several resistance mechanisms. The bacteria have been reported as origin of multiple outbreaks. This study aims to investigate the use of efflux pumps and quinolone resistance-associated genotypic mutations as mechanisms of resistance in A. baumannii isolates at a tertiary hospital. A total number of 103 A. baumannii isolates were investigated after identification and antimicrobial susceptibility testing by VITEK2 followed by PCR amplification of blaOXA-51. Conventional PCR amplification of the AdeABC efflux pump (adeB, adeS, and adeR) and quinolone (parC and gyrA) resistance genes were performed, followed by quantitative real-time PCR of AdeABC efflux pump genes. Phenotypic evaluation of efflux pump expression was performed by determining the difference between the MIC of tigecycline before and after exposure to an efflux pump inhibitor. The Sanger sequencing method was used to sequence the parC and gyrA amplicons. A phylogenetic tree was drawn using MEGA 4.0 to evaluate evolutionary relatedness of the strains. All the collected isolates were blaOXA-51-positive. High resistance to almost all the tested antibiotics was observed. Efflux pump was found in 75% of isolates as a mechanism of resistance. The study detected parC gene mutation in 60% and gyrA gene mutation in 85%, while 37% of isolates had mutations on both genes. A minimal evolutionary distance between the isolates was reported. The use of the AdeABC efflux pump system as an active mechanism of resistance combined with point mutation mainly in gyrA was shown to contribute to broaden the resistance spectrum of A. baumannii isolates.

BackgroundRecurrent preimplantation embryo developmental arrest (RPEA) is the most common cause of assisted reproductive technology treatment failure associated with identified genetic abnormalities. Variants in known maternal genes can only account for 20%–30% of these cases. The underlying genetic causes for the other affected individuals remain unknown.MethodsWhole exome sequencing was performed for 100 independent infertile females that experienced RPEA. Functional characterisations of the identified candidate disease-causative variants were validated by Sanger sequencing, bioinformatics and in vitro functional analyses, and single-cell RNA sequencing of zygotes.ResultsBiallelic variants in ZFP36L2 were associated with RPEA and the recurrent variant (p.Ser308_Ser310del) prevented maternal mRNA decay in zygotes and HeLa cells.ConclusionThese findings emphasise the relevance of the relationship between maternal mRNA decay and human preimplantation embryo development and highlight a novel gene potentially responsible for RPEA, which may facilitate genetic diagnoses.

Familial periodic paralyses (PPs) are inherited disorders of skeletal muscle characterized by recurrent episodes of flaccid muscle weakness. PPs are classified as hypokalemic (HypoPP), normokalemic (NormoPP), or hyperkalemic (HyperPP) according to the potassium level during the paralytic attacks. HypoPP is an autosomal dominant disease caused by mutations in the CACNA1S gene, encoding for Cav1.1 channel (HypoPP-1), or SCN4A gene, encoding for Nav1.4 channel (HypoPP-2). In the present study, we included 60 patients with a clinical diagnosis of HypoPP. Fifty-one (85%) patients were tested using the direct sequencing (Sanger method) of all reported HypoPP mutations in CACNA1S and SCN4A genes; the remaining 9 (15%) patients were analyzed through a next-generation sequencing (NGS) panel, including the whole CACNA1S and SCN4A genes, plus other genes rarely associated to PPs. Fifty patients resulted mutated: 38 (76%) cases showed p.R528H and p.R1239G/H CACNA1S mutations and 12 (24%) displayed p.R669H, p.R672C/H, p.R1132G/Q, and p.R1135H SCN4A mutations. Forty-one mutated cases were identified among the 51 patients managed with Sanger sequencing, while all the 9 cases directly analyzed with the NGS panel showed mutations in the hotspot regions of SCN4A and CACNA1S. Ten out of the 51 patients unresolved through the Sanger sequencing were further analyzed with the NGS panel, without the detection of any mutation. Hence, our data suggest that in HypoPP patients, the extension of genetic analysis from the hotspot regions using the Sanger method to the NGS sequencing of the entire CACNA1S and SCN4A genes does not lead to the identification of new pathological mutations.

The province of Misiones is considered a region with a high mortality rate due to cervical cancer (CC). To gain insight into this problem, we explored the association between genetic variation in the E6 and E7 oncogenes of HPV16 and the risk of CC. We studied 160 women with cytological diagnoses of negative for intraepithelial lesion or malignity, low-grade squamous intraepithelial lesion, and high-grade squamous intraepithelial lesion/CC and a positive test for HPV16 infection. The genetic characterization of E6 and E7 genes was undertaken through PCR amplification and direct Sanger sequencing. Phylogenetic classification was conducted using Bayesian methods. To estimate the odds ratio (OR) for an association between genetic variants in the E6 and E7 genes and the risk of CC, we used ordinal logistic regression adjusted by age. The final data set comprised 112 samples. Diagnostic single-nucleotide polymorphisms (SNPs) and phylogenetic trees confirmed the presence of Lineage A (95.5%) and D (4.5%) in the samples. For the E6 gene, we identified eleven different sequences, with the most common ones being Lineage A E6 350G (58.9%) and E6 350T (37.5%). The E6 350G was associated with progression to HSIL/CC, with an OR of 19.41 (4.95-76.10). The E7 gene was more conserved than E6, probably due to the functional constraints of this small protein. Our results confirmed the association of the E6 350G SNP with a higher risk of developing CC. These data will contribute to understanding the biological bases of CC incidence in this region.

Paula Triana-Fonseca,1 Juan Fernando Parada-Márquez,1 Claudia T Silva-Aldana,2 Daniela Zambrano-Arenas,1 Laura Lucia Arias-Gomez,2 Natalia Morales-Fonseca,3 Esteban Medina-Méndez,2 Carlos M Restrepo,4 Daniel Felipe Silgado-Guzmán,2 Dora Janeth Fonseca-Mendoza4 1Facultad de Medicina, Universidad del Bosque, Bogotá, DC, Colombia; 2Department of Molecular Diagnosis, Genética Molecular de Colombia SAS, Bogotá, DC, Colombia; 3Department of Medicine, Instituto Colombiano de Neurociencias, Bogotá, DC, Colombia; 4Center for Research in Genetics and Genomics – CIGGUR, GENIUROS Research Group, School of Medicine and Health Sciences, Universidad Del Rosario, Bogotá, DC, ColombiaCorrespondence: Dora Janeth Fonseca-Mendoza; Daniel Felipe Silgado Guzmán Email dora.fonseca@urosario.edu.co; dansilgado@hotmail.comBackground: Duchenne and Becker muscular dystrophies (DMD/BMD) are the most common human dystrophinopathies with recessive X-linked inheritance. Dystrophin gene deletions and duplications are the most common mutations, followed by point mutations. The aim of this study is to characterize the mutational profile of the dystrophin gene in Colombian patients with DMD/BMD.Material and Methods: Mutational profiling was determined in 69 affected patients using Sanger sequencing, next-generation sequencing (NGS) and/or multiplex ligation dependent-probes amplification (MLPA). Genetic variants were classified according to molecular consequence and new variants were determined through database and literature analysis.Results: Mutational profile in affected patients revealed that large deletions/duplications analyzed by MLPA accounted for 72.5% of all genetic variations. By using Sanger sequencing or NGS, we identified point mutations in 15.9% and small deletions in 11.6% of the patients. New mutations were found, most of them were point mutations or small deletions (10.1%).Conclusion: Our results described the genetic profile of the dystrophin gene in Colombian patients with DMD and contribute to efforts to identify molecular variants in Latin American populations. For our population, 18.8% of cases could be treated with FDA or MDA approved molecular therapies based on specific mutations. These data contribute to the establishment of appropriate genetic counseling and potential treatment.Keywords: Duchenne–Becker muscular dystrophy, DMD, MLPA, next-generation sequencing, target molecular therapy, exon skipping, mutation

Recessive variants in GJB2 are the most common genetic cause of sensorineural hearing impairment. However, in many patients, only one variant in the GJB2 coding region is identified using conventional sequencing strategy (eg, Sanger sequencing), resulting in nonconfirmative diagnosis. Conceivably, there might be other unidentified pathogenic variants in the noncoding region of GJB2 or other deafness-causing genes in these patients. To address this, a next-generation sequencing-based diagnostic panel targeting the entire GJB2 gene and the coding regions of 158 other known deafness-causing genes was designed and applied to 95 patients with nonsyndromic sensorineural hearing impairment (including 81 Han Taiwanese and 14 Mongolian patients) in whom only a single GJB2 variant had been detected using conventional Sanger sequencing. The panel confirmed the genetic diagnosis in 24 patients (25.3%). Twenty-two of them had causative variants in several deafness-causing genes other than GJB2, including MYO15A, MYO7A, TECTA, POU4F3, KCNQ4, SLC26A4, OTOF, MT-RNR1, MITF, WFS1, and USH2A. The other two patients had causative variants in GJB2, including a Taiwanese patient with a mosaic maternal uniparental disomy c.235delC variant (approximately 69% mosaicism) and a Mongolian patient with compound heterozygous c.35dupG and c.35delG variants, which occurred at the same site. This study demonstrates the utility of next-generation sequencing in clarifying the genetic diagnosis of hearing-impaired patients with nonconfirmative GJB2 genotypes on conventional genetic examinations.

Despite the immense genetic heterogeneity of B-lymphoblastic leukemia [or precursor B-cell acute lymphoblastic leukemia (B-ALL)], RNA sequencing (RNA-Seq) could comprehensively interrogate its genetic drivers, assigning a specific molecular subtype in >90% of patients. However, study groups have only started to use RNA-Seq. For broader clinical use, technical, quality control, and appropriate performance validation are needed. We describe the development and validation of an RNA-Seq workflow for subtype classification, TPMT/NUDT15/TP53 variant discovery, and immunoglobulin heavy chain (IGH) disease clone identification for Malaysia-Singapore acute lymphoblastic leukemia (ALL) 2020. We validated this workflow in 377 patients in our preceding Malaysia-Singapore ALL 2003/Malaysia-Singapore ALL 2010 studies and proposed the quality control measures for RNA quality, library size, sequencing, and data analysis using the International Organization for Standardization 15189 quality and competence standard for medical laboratories. Compared with conventional methods, we achieved >95% accuracy in oncogene fusion identification, digital karyotyping, and TPMT and NUDT15 variant discovery. We found seven pathogenic TP53 mutations, confirmed with Sanger sequencing, which conferred a poorer outcome. Applying this workflow prospectively to the first 21 patients in Malaysia-Singapore ALL 2020, we identified the genetic drivers and IGH disease clones in >90% of patients with concordant TPMT, NUDT15, and TP53 variants using PCR-based methods. The median turnaround time was 12 days, which was clinically actionable. In conclusion, RNA-Seq workflow could be used clinically in management of B-cell ALL patients.

Genetic variants of SARS-CoV-2 have been emerging and circulating in many places across the world. Rapid detection of these variants is essential since their dissemination can impact transmission rates, diagnostic procedures, disease severity, response to vaccines or patient management. Sanger sequencing has been used as the preferred approach for variant detection among circulating human immunodeficiency and measles virus genotypes. Using primers to amplify a fragment of the SARS-CoV-2 genome encoding part of the Spike protein, we showed that Sanger sequencing allowed us to rapidly detect the introduction and spread of three distinct SARS-CoV-2 variants in two major Brazilian cities. In both cities, after the predominance of variants closely related to the virus first identified in China, the emergence of the P.2 variant was quickly followed by the identification of the P1 variant, which became dominant in less than one month after it was first detected.

Background Hereditary angioedema (HAE) is a potentially fatal disorder resulting in recurrent attacks of severe swelling. It may be associated with a genetic deficiency of functional C1 inhibitor or with normal C1 inhibitor (HAEnCI). In families with HAEnCI, HAE-linked mutations in the F12, PLG, KNG1, ANGPT1, or MYOF genes have been identified. In many families with HAEnCI the genetic cause of the disease is currently unknown. Objective The aim of this study was to identify a novel disease-linked mutation for HAEnCI. Methods The study methods comprised whole exome sequencing, Sanger sequencing analysis, pedigree analysis, bioinformatic analysis of the mutation, and biochemical analysis of parameters of the kallikrein-kinin (contact) system. Results By performing whole exome sequencing on a multigenerational family with HAEnCI we were able to identify the heparan sulfate (HS)-glucosamine 3-O-sulfotransferase 6 (HS3ST6) mutation c.430A>T (p.Thr144Ser) in all 3 affected family members who were sequenced. This gene encodes HS-glucosamine 3-O-sulfotransferase 6 (3-OST-6), which is involved in the last step of HS biosynthesis. The p.Thr144Ser mutation is likely to affect the interaction between 2 β-sheets stabilizing the active center of the 3-OST-6 protein. Conclusions We conclude that mutant 3-OST-6 fails to transfer sulfo groups to the 3-OH position of HS, resulting in incomplete HS biosynthesis. This likely affects cell surface interactions of key players in angioedema formation and is a novel mechanism for disease development.

Protein S (PS) deficiency is a risk factor for venous thromboembolism (VTE) and can be caused by variants of the gene encoding PS (PROS1). This study aimed to evaluate the clinical value of molecular analysis of the PROS1 gene in PS-deficient participants. We performed Sanger sequencing of the coding region of the PROS1 gene and multiplex ligation-dependent probe amplification to exclude large structural rearrangements. Free PS was measured by a particle-enhanced immunoassay, while PS activity was assessed by a clotting method.A total of 87 PS-deficient participants and family members were included. In 22 index participants, we identified 13 PROS1 coding variants. Five variants were novel. In 21 index participants, no coding sequence variants or structural rearrangements were identified. The free PS level was lower in index participants carrying a PROS1 variant compared with index participants with no variant (0.51 [0.32–0.61] vs. 0.62 [0.57–0.73] × 103 IU/L; p In conclusion, we report 13 PROS1 coding variants including five novel variants. PS levels differ by PROS1 variant and the frequency of VTE was higher when a coding PROS1 variant was present. Hence, molecular analysis of the PROS1 gene may add clinical value in the diagnostic work-up of PS deficiency.

The association of taste genetics and the oral microbiome in autoimmune diseases such as rheumatoid arthritis (RA) has not been reported. We explored a novel oral mucosal innate immune pathway involving the bitter taste G protein-coupled receptor T2R38. This case–control study aimed to evaluate whether T2R38 polymorphisms associate with the buccal microbial composition in RA. Genomic DNA was obtained from buccal swabs of 35 RA patients and 64 non-RA controls. TAS2R38 genotypes were determined by Sanger sequencing. The buccal microbiome was assessed by Illumina MiSeq sequencing of the V4-16S rRNA gene. Bacterial community differences were analyzed with alpha and beta diversity measures. Linear discriminant analysis effect size identified taxa discriminating between RA versus non-RA and across TAS2R38 genotypes. TAS2R38 genotype frequency was similar between RA and non-RA controls (PAV/PAV, PAV/AVI, AVI/AVI: RA 42.9%, 45.7%, 11.4% versus controls 32.8%, 48.4%, 18.8%, chi-square (2, N = 99) = 2.1, p = 0.35). The relative abundance of Porphyromonas, among others, differed between RA and non-RA controls. The relative abundance of several bacterial species also differed across TAS2R38 genotypes. These findings suggest an association between T2R38 polymorphisms and RA buccal microbial composition. However, further research is needed to understand the impact of T2R38 in oral health and RA development.

INTRODUCTION Hemophilia B is associated with molecular heterogeneity, with more than 1200 unique variants in the F9 gene. We hereby describe the mutational spectrum of severe hemophilia B patients presenting in a tertiary-care center in India. METHOD DNA was extracted from peripheral blood samples of 35 diagnosed severe hemophilia B patients belonging to 32 families, and were subjected to Sanger sequencing. Determination of the effect of novel variants on the protein structure and correlation between genotype and phenotype was attempted using in-silico tools. RESULTS Twenty-seven different mutations were detected in 30 probands, including 20 known and 7 novel variants. Also, we found one suspected case of whole gene deletion. The serine peptidase domain harbored most of the variants (48.1%). Inhibitory antibodies were found in two patients. CONCLUSIONS This study provides a comprehensive mutational spectrum and mutation screening strategy by Sanger sequencing of F9 gene in severe hemophilia B patients, in a resource-constraint setting.