Comparative genomics and genome-wide SNPs of endangered Eld's deer provide breeder selection for inbreeding avoidance.

Eld's deer, a conserved wildlife species of Thailand, is facing inbreeding depression, particularly in the captive Siamese Eld's deer (SED) subspecies. In this study, we constructed genomes of a male SED and a male Burmese Eld's deer (BED), and used genome-wide single nucleotide polymorphisms to evaluate the genetic purity and the inbreeding status of 35 SED and 49 BED with limited pedigree information. The results show that these subspecies diverged approximately 1.26 million years ago. All SED were found to be purebred. A low proportion of admixed SED genetic material was observed in some BED individuals. Six potential breeders from male SED with no genetic relation to any female SED and three purebred male BED with no relation to more than 10 purebred female BED were identified. This study provides valuable insights about Eld's deer populations and appropriate breeder selection in efforts to repopulate this endangered species while avoiding inbreeding.

Exome sequencing as first-tier genetic testing in infantile-onset pharmacoresistant epilepsy: diagnostic yield and treatment impact.

Pharmacoresistant epilepsy presenting during infancy poses both diagnostic and therapeutic challenges. We aim to identify diagnostic yield and treatment implications of exome sequencing (ES) as first-tier genetic testing for infantile-onset pharmacoresistant epilepsy. From June 2016 to December 2020, we enrolled patients with infantile-onset (age ≤ 12 months) pharmacoresistant epilepsy. 103 unrelated patients underwent ES. Clinical characteristics and changes in management due to the molecular diagnosis were studied. 42% (43/103) had epilepsy onset within the first month of life. After ES as first-tier genetic testing, 62% (64/103) of the cases were solved. Two partially solved cases (2%; 2/103) with heterozygous variants identified in ALDH7A1 known to cause autosomal recessive pyridoxine dependent epilepsy underwent genome sequencing (GS). Two novel large deletions in ALDH7A1 were detected in both cases. ES identified 66 pathogenic and likely pathogenic single nucleotide variants (SNVs) in 27 genes. 19 variants have not been previously reported. GS identified two additional copy number variations (CNVs). The most common disease-causing genes are SCN1A (13%; 13/103) and KCNQ2 (8%; 8/103). Eight percent (8/103) of the patients had treatable disorders and specific treatments were provided resulting in seizure freedom. Pyridoxine dependent epilepsy was the most common treatable epilepsy (6%; 6/103). Furthermore, 35% (36/103) had genetic defects which guided gene-specific treatments. Altogether, the diagnostic yield is 64%. Molecular diagnoses change management in 43% of the cases. This study substantiates the use of next generation sequencing (NGS) as the first-tier genetic investigation in infantile-onset pharmacoresistant epilepsy.

A LILRB1 variant with a decreased ability to phosphorylate SHP-1 leads to autoimmune diseases.

Inborn errors of immunity are known to cause not only immunodeficiencies and allergies but also autoimmunity. Leukocyte immunoglobulin-like receptor B1 (LILRB1) is a receptor on leukocytes playing a role in regulating immune responses. No phenotypes have been reported to be caused by germline mutations in LILRB1. We aimed to identify the causative variant in a three-generation family with nine members suffering from one of the three autoimmune diseases-Graves' disease, Hashimoto's thyroiditis, or systemic lupus erythematosus. Whole-genome linkage study revealed a locus on chromosome 19q13.4 with the maximum LOD score of 2.71. Whole-exome sequencing identified a heterozygous missense variant, c.479G > A (p. G160E) in LILRB1, located within the chromosomal-linked region, in all nine affected members. The variant has never been previously reported. Jurkat cells transfected with the mutant LILRB1, compared with those with the wild-type LILRB1, showed decreased phosphorylation of both LILRB1 and its downstream protein, SHP-1. Flow cytometry was used to study immunophenotype and revealed that LILRB1 was significantly lower on the surface of activated regulatory T lymphocytes (Treg) cells of patients. Single-cell RNA sequencing showed substantially increased M1-like monocytes in peripheral blood mononuclear cells of affected individuals. This study, for the first time, implicates LILRB1 as a new disease gene for autoimmunity.

Host genetic factors of COVID-19 susceptibility and disease severity in a Thai population.

Host genetic factors have been shown to play a role in SARs-CoV-2 infection in diverse populations. However, the genetic landscape differs among various ethnicities; therefore, we explored the host genetic factors associated with COVID-19 disease susceptivity and disease severity in a Thai population. We recruited and genotyped 212 unrelated COVID-19 Thai patients and 36 controls using AxiomTM Human Genotyping SARs-COV-2 array, including 847,384 single nucleotide polymorphisms related to SARs-COV-2 pathogenesis, immune response, and related comorbidity No SNPs passed the genome-wide significance threshold of p value <1 × 10-8. However, with a threshold of p value <1 × 10-5, a locus on chromosome 5q32 was found to have a suggestive association with COVID-19 disease susceptibility (p value 6.9 × 10-6; Q-Q plot λ = 0.805, odds ratio 0.02). Notably, IL17B is a gene located in this linkage disequilibrium block and is previously shown to play a part in inflammation and pneumonia. Additionally, a suggestive locus on chromosome 12q22, harboring EEA1 and LOC643339, was associated with COVID-19 disease severity (p value 1.3 × 10-6 - 4.4 × 10-6, Q-Q plot λ = 0.997, odds ratio 0.28-0.31). EEA1 is involved in viral entry into cells, while LOC643339 is a long non-coding RNA. In summary, our study suggested loci on chromosomes 5q32 and 12q22 to be linked to COVID-19 disease susceptibility and disease severity, respectively. The small sample size of this study may lessen the likelihood that the association found is real, but it could still be true. Further study with a larger cohort is required to confirm these findings.

Age-related dental phenotypes and tooth characteristics of FAM83H-associated hypocalcified amelogenesis imperfecta.

OBJECTIVES: Autosomal-dominant hypocalcified amelogenesis imperfecta (ADHCAI) shows phenotypic heterogeneity. Our aim was to characterise the ADHCAI phenotypes, tooth properties and genotypes. METHODS: Three unrelated ADHCAI probands and seven additional affected members of the three families were recruited. Mutations were identified by exome and Sanger sequencing, and haplotypes by SNP array. Tooth colour, roughness, density, nanohardness, minerals and ultrastructure were investigated. RESULTS: Ten participants were heterozygous for the FAM83H mutation c.1387C>T (p.Gln463*). All shared a 3.43 Mbp region on chromosome 8q24.3 encompassing the FAM83H variant, indicating a common ancestry. The c.1387C>T was estimated to be 23.8 generations or 600 years. The FAM83H enamel had higher roughness and lower lightness, density, nanohardness, and calcium and phosphorus levels than controls. Blunted enamel rods, wide interrod spaces and disorganised dentinoenamel junctions were observed. Evaluating the patients with the same mutation and reviewing others with different mutations in FAM83H revealed that the FAM83H heterogeneous phenotypes are age-influenced. Tooth colour and surface texture change with ageing. CONCLUSIONS: FAM83H enamel demonstrated decreased lightness, density, hardness, calcium, phosphorus and defective ultrastructure. We have identified that the phenotypic variation in FAM83H-associated ADHCAI is age-related. Awareness of the correlation between age and clinical features of FAM83H-ADHCAI can help dentists make an accurate diagnosis.

Metagenomic analysis of viral genes integrated in whole genome sequencing data of Thai patients with Brugada syndrome.

Brugada syndrome (BS) is an autosomal dominant inheritance cardiac arrhythmia disorder associated with sudden death in young adults. Thailand has the highest prevalence of BS worldwide, and over 60% of patients with BS still have unclear disease etiology. Here, we performeda new viral metagenome analysis pipeline called VIRIN and validated it with whole genome sequencing (WGS) data of HeLa cell lines and hepatocellular carcinoma. Then the VIRIN pipelinewas applied to identify viral integration positions from unmapped WGS data of Thai males, including 100 BS patients (case) and 100 controls. Even though the sample preparation had noviral enrichment step, we can identify several virus genes from our analysis pipeline. The predominance of human endogenous retrovirus K (HERV-K) viruses was found in both cases andcontrols by blastn and blastx analysis. This study is the first report on the full-length HERV-Kassembled genomes in the Thai population. Furthermore, the HERV-K integration breakpointpositions were validated and compared between the case and control datasets. Interestingly,Brugada cases contained HERV-K integration breakpoints at promoters five times more oftenthan controls. Overall, the highlight of this study is the BS-specific HERV-K breakpoint positionsthat were found at the gene coding region "NBPF11" (n = 9), "NBPF12" (n = 8) and longnon-coding RNA (lncRNA) "PCAT14" (n = 4) region. The genes and the lncRNA have been reported to be associated with congenital heart and arterial diseases. These findings provide another aspect of the BS etiology associated with viral genome integrations within the humangenome.

The Thai reference exome (T-REx) variant database.

To maximize the potential of genomics in medicine, it is essential to establish databases of genomic variants for ethno-geographic groups that can be used for filtering and prioritizing candidate pathogenic variants. Populations with non-European ancestry are poorly represented among current genomic variant databases. Here, we report the first high-density survey of genomic variants for the Thai population, the Thai Reference Exome (T-REx) variant database. T-REx comprises exome sequencing data of 1092 unrelated Thai individuals. The targeted exome regions common among four capture platforms cover 30.04 Mbp on autosomes and chromosome X. 345 681 short variants (18.27% of which are novel) and 34 907 copy number variations were found. Principal component analysis on 38 469 single nucleotide variants present worldwide showed that the Thai population is most genetically similar to East and Southeast Asian populations. Moreover, unsupervised clustering revealed six Thai subpopulations consistent with the evidence of gene flow from neighboring populations. The prevalence of common pathogenic variants in T-REx was investigated in detail, which revealed subpopulation-specific patterns, in particular variants associated with erythrocyte disorders such as the HbE variant in HBB and the Viangchan variant in G6PD. T-REx serves as a pivotal addition to the current databases for genomic medicine.

Severe neonatal haemolytic anaemia caused by compound heterozygous KLF1 mutations: report of four families and literature review.

Mutations in the KLF1 gene, which encodes a transcription factor playing a role in erythropoiesis, have recently been demonstrated to be a rare cause of hereditary haemolytic anaemia. We described the genotypic and phenotypic spectra of four unrelated families with compound heterozygous class 2/class 3 KLF1 mutations. All patients had p.G176RfsX179 on one allele and either p.A298P, p.R301H or p.G335R on the other allele. All presented on the first day of life with severe haemolytic anaemia with abnormal red blood cell morphology, markedly increased nucleated red blood cells and hyperbilirubinaemia. Three patients later became transfusion-dependent. All parents with heterozygous KLF1 mutation without co-inherited thalassaemia had normal to borderline mean corpuscular volume (MCV) and normal to slightly elevated Hb F. Fifteen previously reported cases of biallelic KLF1 mutations were identified from a literature review. All except one presented with severe haemolytic anaemia in the neonatal period. Our finding substantiates that compound heterozygous KLF1 mutations are associated with severe neonatal haemolytic anaemia and expands the haematologic phenotypic spectrum. In carriers, the previously suggested findings of low MCV, high Hb A2 and high Hb F are inconsistent; thus this necessitates molecular studies for the identification of carriers.

Expanding the genotypic spectrum of PYCR2 and a common ancestry in Thai patients with hypomyelinating leukodystrophy 10.

PYCR2 pathogenic variants lead to an autosomal recessive hypomyelinating leukodystrophy 10 (HLD10), characterized by global developmental delay, microcephaly, facial dysmorphism, movement disorder, and hypomyelination. This study identified the first two unrelated Thai patients with HLD10. Patient 1 harbored the novel compound heterozygous variants, c.257T>G (p.Val86Gly) and c.400G>A (p.Val134Met), whereas patient 2 possessed the homozygous variant, c.400G>A (p.Val134Met), in PYCR2. Haplotype analysis revealed that the two families' members shared a 2.3 Mb region covering the c.400G>A variant, indicating a common ancestry. The variant was estimated to age 1450 years ago. Since the c.400G>A was detected in three out of four mutant alleles and with a common ancestry, this variant might be common in Thai patients. We also reviewed the phenotype and genotype of all 35 previously reported PYCR2 patients and found that majorities of cases were homozygous with a consanguineous family history, except patient 1 and another reported case who were compound heterozygous. All patients had microcephaly and developmental delay. Hypotonia and peripheral spasticity were common. Hypomyelination or delayed myelination was a typical radiographic feature. Here, we report the first two Thai patients with HLD10 with the novel PYCR2 variants expanding the genotypic spectrum and suggest that the c.400G>A might be a common mutation in Thai patients.

Rapid exome sequencing as the first-tier investigation for diagnosis of acutely and severely ill children and adults in Thailand.

The use of rapid DNA sequencing technology in severely ill children in developed countries can accurately identify diagnoses and positively impact patient outcomes. This study sought to evaluate the outcome of Thai children and adults with unknown etiologies of critical illnesses with the deployment of rapid whole exome sequencing (rWES) in Thailand. We recruited 54 unrelated patients from 11 hospitals throughout Thailand. The median age was 3 months (range, 2 days-55 years) including 47 children and 7 adults with 52% males. The median time from obtaining blood samples to issuing the rWES report was 12 days (range, 5-27 days). A molecular diagnosis was established in 25 patients (46%), resulting in a change in clinical management for 24 patients (44%) resulting in improved clinical outcomes in 16 patients (30%). Four out of seven adult patients (57%) received the molecular diagnosis which led to a change in management. The 25 diagnoses comprised 23 different diseases. Of the 34 identified variants, 15 had never been previously reported. This study suggests that use of rWES as a first-tier investigation tool can provide tremendous benefits in critically ill patients with unknown etiology across age groups in Thailand.

Prenatal Sonographic Features of CHARGE Syndrome.

CHARGE syndrome is a rare autosomal dominant disorder, associated with coloboma (C), heart defects (H), choanal atresia (A), retardation of growth and/or central nervous system (R), genitourinary anomalies (G) and ear abnormalities (E). Prenatal diagnosis of the syndrome is very rare but may be suspected when a combination of such abnormalities is identified. We describe a prenatally suspected case of CHARGE syndrome due to unique findings of cardiac defects (DORV) in combination with minor clues, including a structurally malformed ear with persistent non-response to an acoustic stimulation (which has never been prenatally described elsewhere), renal malrotation and growth restriction. Postnatal diagnosis was made based on confirmation of the prenatal findings and additional specific findings of bilateral coloboma, choanal atresia and ear canal stenosis. Finally, molecular genetic testing by whole exome sequencing of the neonate and her parents revealed a novel de novo heterozygous frameshift c.3506_3509dup variant in the CHD7 gene, confirming the clinical diagnosis of CHARGE syndrome. In conclusion, we describe unique prenatal features of CHARGE syndrome. Educationally, this is one of the rare examples of CHARGE syndrome, comprising all of the six specific anomalies as originally described; it is also supported by the identification of a specific genetic mutation. The identified genetic variant has never been previously reported, thereby expanding the mutational spectrum of CHD7. Finally, this case can inspire prenatal sonographers to increase awareness of subtle or minor abnormalities as genetic sonomarkers.

Double heterozygous variants in FBN1 and FBN2 in a Thai woman with Marfan and Beals syndromes.

A phenotype of an individual is resulted from an interaction among variants in several genes. Advanced molecular technologies allow us to identify more patients with mutations in more than one genes. Here, we studied a Thai woman with combined clinical features of Marfan (MFS) and Beals (BS) syndromes including frontal bossing, enophthalmos, myopia, the crumpled appearance to the top of the pinnae, midface hypoplasia, high arched palate, dermal stretch marks, aortic enlargement, mitral valve prolapse and regurgitation, aortic root dilatation, and progressive scoliosis. The aortic root enlargement was progressive to a diameter of 7.2 cm requiring an aortic root replacement at the age of 8 years. At her last visit when she was 19 years old, she had moderate aortic regurgitation. Exome sequencing revealed that she carried the c.3159C > G (p.Cys1053Trp) in exon 26 of FBN1 and c.2638G > A (p. Gly880Ser) in exon 20 of FBN2. The variant in FBN1 was de novo, while that in FBN2 was inherited from her unaffected mother. Both genes encode for fibrillins, which are essential for elastic fibers and can form the heterotypic microfibrils. Two defective fibrillins may synergistically worsen cardiovascular manifestations seen in our patient. In this study, we identified the fourth patient with both MFS and BS, carrying mutations in both FBN1 and FBN2.

A case of GABRA5-related developmental and epileptic encephalopathy with response to a combination of antiepileptic drugs and a GABAering agent.

BACKGROUND: GABAA receptors are ligand-gated chloride channels that regulate inhibitory neurotransmission in the central nervous system. Recently, monoallelic de novo mutations in GABRA5 resulting in altered inhibitory synapses were found in three patients with developmental and epileptic encephalopathy. Patient description: We report on a four-year and six-month-old girl with epilepsy and global developmental delay. Serial head growth measurement revealed postnatal onset microcephaly. RESULTS: Magnetic resonance imaging (MRI) of the brain was normal at the age of eight months and subsequently showed a decrease in white matter volume and thin corpus callosum at the age of 3 years. Using whole-genome sequencing, we identified the fourth patient harboring a de novo missense mutation in GABRA5. Interestingly, the c.880G > C (p.V294F) affects the same position found in two of the three previously reported patients. CONCLUSION: This study suggests that the nucleotide c.880G is a mutation hotspot. Our patient also demonstrated significant seizure reduction after benzodiazepine. To our knowledge, this is the first case describing the favorable outcome of a GABAergic agent in seizure control for GABRA5-related developmental and epileptic encephalopathy.

A patient with combined pituitary hormone deficiency and osteogenesis imperfecta associated with mutations in LHX4 and COL1A2.

Genetic disorders have been shown to co-occur in individual patient. A Thai boy with features of osteogenesis imperfecta (OI) and combined pituitary hormone deficiency (CPHD) was identified. The causative mutations were investigated by whole exome and Sanger sequencing. Pathogenicity and pathomechanism of the variants were studied by luciferase assay. The proband was found to harbor a novel de novo heterozygous missense mutation, c.1531G > T (p.G511C), in COL1A2 leading to OI and a heterozygous missense variant, c.364C > T (p.R122W), in LHX4. The LHX4 p.R122W has never been reported to cause CPHD. The variant was predicted to be deleterious and found in the highly conserved LIM2 domain of LHX4. The luciferase assays revealed that the p.R122W was unable to activate POU1F1, GH1, and TSHB promoters, validating its pathogenic effect in CPHD. Moreover, the variant did not alter the function of wild-type LHX4, indicating its hypomorphic pathomechanism. In conclusion, the novel de novo heterozygous p.G511C mutation in COL1A2 and the heterozygous pathogenic p.R122W mutation in LHX4 were demonstrated in a patient with OI and CPHD. This study proposes that the mutations in two different genes should be sought in the patients with clinical features unable to be explained by a mutation in one gene.

A regulatory variant of CHRM3 is associated with cannabis-induced hallucinations in European Americans.

Cannabis, the most widely used illicit drug, can induce hallucinations. Our understanding of the biology of cannabis-induced hallucinations (Ca-HL) is limited. We used the Semi-Structured Assessment for Drug Dependence and Alcoholism (SSADDA) to identify cannabis-induced hallucinations (Ca-HL) among long-term cannabis users (used cannabis ≥1 year and ≥100 times). A genome-wide association study (GWAS) was conducted by analyzing European Americans (EAs) and African Americans (AAs) in Yale-Penn 1 and 2 cohorts individually, then meta-analyzing the two cohorts within population. In the meta-analysis of Yale-Penn EAs (n = 1917), one genome-wide significant (GWS) signal emerged at the CHRM3 locus, represented by rs115455482 (P = 1.66 × 10-10), rs74722579 (P = 2.81 × 10-9), and rs1938228 (P = 1.57 × 10-8); signals were GWS in Yale-Penn 1 EAs (n = 1092) and nominally significant in Yale-Penn 2 EAs (n = 825). Two SNPs, rs115455482 and rs74722579, were available from the Collaborative Study on the Genetics of Alcoholism data (COGA; 3630 long-term cannabis users). The signals did not replicate, but when meta-analyzing Yale-Penn and COGA EAs, the two SNPs' association signals were increased (meta-P-values 1.32 × 10-10 and 2.60 × 10-9, respectively; n = 4291). There were no significant findings in AAs, but in the AA meta-analysis (n = 3624), nominal significance was seen for rs74722579. The rs115455482*T risk allele was associated with lower CHRM3 expression in the thalamus. CHRM3 was co-expressed with three psychosis risk genes (GABAG2, CHRNA4, and HRH3) in the thalamus and other human brain tissues and mouse GABAergic neurons. This work provides strong evidence for the association of CHRM3 with Ca-HL and provides insight into the potential involvement of thalamus for this trait.

Whole exome sequencing revealed mutations in FBXL4, UNC80, and ADK in Thai patients with severe intellectual disabilities.

Intellectual disabilities (ID) are etiologically heterogeneous. Advanced molecular techniques could be helpful in identification of the underlying genetic defects. We aimed to characterize clinical and molecular features of three Thai patients with ID. Patient 1 had ID, hypotonia and lactic acidosis. Patient 2 had ID and growth failure. Patient 3 had ID, seizure, diarrhea and hypoglycemia. Whole exome sequencing found that Patient 1 was homozygous for a nonsense, c.1303C>T (p.Arg435Ter), mutation in FBXL4, a gene responsible for encephalomyopathic mitochondrial DNA depletion syndrome-13 (MTDPS13). Patient 2 was compound heterozygous for two novel mutations, c.3226C>T (p.Arg1076Ter) and c.3205C>T (p.Arg1069Ter), in UNC80, a known gene of infantile hypotonia with psychomotor retardation and characteristic facies-2 (IHPRF2). Patient 3 was homozygous for a novel missense, c.427T>C (p.Cys143Arg), mutation in ADK, a known gene of adenosine kinase deficiency leading to hypermethioninemia. This study expands the mutational spectra of ID genes.

Hypomethylation of intragenic LINE-1 represses transcription in cancer cells through AGO2.

In human cancers, the methylation of long interspersed nuclear element -1 (LINE-1 or L1) retrotransposons is reduced. This occurs within the context of genome wide hypomethylation, and although it is common, its role is poorly understood. L1s are widely distributed both inside and outside of genes, intragenic and intergenic, respectively. Interestingly, the insertion of active full-length L1 sequences into host gene introns disrupts gene expression. Here, we evaluated if intragenic L1 hypomethylation influences their host gene expression in cancer. First, we extracted data from L1base (http://l1base.molgen.mpg.de), a database containing putatively active L1 insertions, and compared intragenic and intergenic L1 characters. We found that intragenic L1 sequences have been conserved across evolutionary time with respect to transcriptional activity and CpG dinucleotide sites for mammalian DNA methylation. Then, we compared regulated mRNA levels of cells from two different experiments available from Gene Expression Omnibus (GEO), a database repository of high throughput gene expression data, (http://www.ncbi.nlm.nih.gov/geo) by chi-square. The odds ratio of down-regulated genes between demethylated normal bronchial epithelium and lung cancer was high (p<1E(-27); OR = 3.14; 95% CI = 2.54-3.88), suggesting cancer genome wide hypomethylation down-regulating gene expression. Comprehensive analysis between L1 locations and gene expression showed that expression of genes containing L1s had a significantly higher likelihood to be repressed in cancer and hypomethylated normal cells. In contrast, many mRNAs derived from genes containing L1s are elevated in Argonaute 2 (AGO2 or EIF2C2)-depleted cells. Hypomethylated L1s increase L1 mRNA levels. Finally, we found that AGO2 targets intronic L1 pre-mRNA complexes and represses cancer genes. These findings represent one of the mechanisms of cancer genome wide hypomethylation altering gene expression. Hypomethylated intragenic L1s are a nuclear siRNA mediated cis-regulatory element that can repress genes. This epigenetic regulation of retrotransposons likely influences many aspects of genomic biology.

Detection of LINE-1s hypomethylation in oral rinses of oral squamous cell carcinoma patients.

This study aimed to (i) investigate long interspersed nuclear element-1 (LINE-1) methylation levels of oral squamous cell carcinomas (OSCCs), the major type of oral malignancies; and (ii) investigate whether the hypomethylation of LINE-1s can be detected in oral rinses of OSCC patients. The combined bisulfite restriction analysis polymerase chain reaction (PCR) of LINE-1s (COBRALINE-1) was used. We found that tissues from OSCC specimens had lower methylation levels of LINE-1s than cells collected from the oral rinses of normal volunteers. Interestingly, cells collected from oral rinses of OSCC patients also revealed hypomethylated LINE-1s at the same level as OSCC tissues. There was no difference in the level of hypomethylation due to stages, locations, histological grades, and history of betel chewing, smoking and/or alcohol consumption. In conclusion, OSCCs possessed global hypomethylation and this alteration could be detected from oral rinses of OSCC patients by a simple PCR technique, COBRALINE-1. Therefore, COBRALINE-1 of oral rinses may be applied for non-invasive detection of oral malignancies.

LINE-1 methylation patterns of different loci in normal and cancerous cells.

This study evaluated methylation patterns of long interspersed nuclear element-1 (LINE-1) sequences from 17 loci in several cell types, including squamous cell cancer cell lines, normal oral epithelium (NOE), white blood cells and head and neck squamous cell cancers (HNSCC). Although sequences of each LINE-1 are homologous, LINE-1 methylation levels at each locus are different. Moreover, some loci demonstrate the different methylation levels between normal tissue types. Interestingly, in some chromosomal regions, wider ranges of LINE-1 methylation levels were observed. In cancerous cells, the methylation levels of most LINE-1 loci demonstrated a positive correlation with each other and with the genome-wide levels. Therefore, the loss of genome-wide methylation in cancerous cells occurs as a generalized process. However, different LINE-1 loci showed different incidences of HNSCC hypomethylation, which is a lower methylation level than NOE. Additionally, we report a closer direct association between two LINE-1s in different EPHA3 introns. Finally, hypermethylation of some LINE-1s can be found sporadically in cancer. In conclusion, even though the global hypomethylation process that occurs in cancerous cells can generally deplete LINE-1 methylation levels, LINE-1 methylation can be influenced differentially depending on where the particular sequences are located in the genome.