Management of gene promoter mutations in molecular diagnostics

A point mutation in the IAA14 promoter enhances callus formation and regeneration

Callus formation induced by auxin accumulation is considered the first step of in vitro plant regeneration. In Arabidopsis, degradation of the Aux/IAA protein, IAA14, in response to auxin signaling, which activates the AUXIN RESPONSE FACTOR 7 (ARF7) and ARF19 along with a series of downstream transcription factors, also plays a critical role in this process. However, the specific mechanism by which auxin regulates callus formation remains unclear. By screening mutant library in the solitary root 1 (iaa14/slr) Arabidopsis background we obtained the callus formation related 2 (cfr2) mutant. The cfr2 mutant exhibited a stronger capacity for callus formation, as well as lateral root and adventitious root regeneration from leaf explants than wild type (WT) seedlings, but did not recover gravitropism capability. The auxin signal in cfr2 was significantly enhanced, and the expression of some downstream transcription factors was increased. Map-based cloning, whole genome resequencing, and phenotypic complementation experiments showed that the phenotypes observed in the cfr2 mutant were caused by a point mutation in the IAA14 promoter region. This mutation, which is predicted to disrupt the binding of LBD16, LBD19, and LBD30 to the IAA14 promoter, changed the expression pattern of IAA14 in cfr2. Taken together, our results identified a new mutation in the IAA14 promoter region, which affects the expression pattern of IAA14 and in turn its ability to control plant regeneration.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-024-01493-y.

The study of single nucleotide polymorphism in a promoter region of ACE-2 receptor gene in the samples of Iraqi population with COVID-19

BACKGROUND

Angiotensin-converting enzyme 2 (ACE2) is an essential receptor on the host cell's cell membrane. It's interesting to note that the entry point receptor ACE2 protein and the severe acute respiratory syndrome (SARS) coronavirus are correlated. This study aimed to determine the influence of the ACE gene genotype and explore the effects of genetic variation in the promotor region of the ACE-2 gene receptor in SARS COV-2 patients.

METHODS AND RESULTS

The 225 participants were categorized into two groups (75 infected and 150 control) according to the results of Real Time -polymerase chain reaction (RT-PCR), IgM, and IgG, also included two types of samples were collected for diagnosis hematological and molecular study. The hematological and biochemical parameters showed significant differences between the two studied groups according to D. dimer, ferritin, lactate dehydrogenase (LDH), C-reactive protein (CRP), white blood cells (WBC), lymphocyte, packed cell volume (PCV) (P˂0.0001), also red blood cell (RBC) (P = 0.0034). While the results of hemoglobin (HB) and platelet displayed non-significant differences between the two groups (p value 0.6811 and 0.9201 respectively). In addition, the sequencing result in the promotor of the ACE-2 gene detected novel eight polymorphisms and recorded them in NCBI under no. (ON959139).

CONCLUSIONS

The ACE D/D polymorphism associated with increased levels of ACE could represent a genetic risk factor in addition the discovery stems from the prospect that genetic differences could lead to differing responses to COVID-19 therapies.

Unexplored power of CRISPR-Cas9 in neuroscience, a multi-OMICs review

The neuroscience and neurobiology of gene editing to enhance learning and memory is of paramount interest to the scientific community. The advancements of CRISPR system have created avenues to treat neurological disorders by means of versatile modalities varying from expression to suppression of genes and proteins. Neurodegenerative disorders have also been attributed to non-canonical DNA secondary structures by affecting neuron activity through controlling gene expression, nucleosome shape, transcription, translation, replication, and recombination. Changing DNA regulatory elements which could contribute to the fate and function of neurons are thoroughly discussed in this review. This study presents the ability of CRISPR system to boost learning power and memory, treat or cure genetically-based neurological disorders, and alleviate psychiatric diseases by altering the activity and the irritability of the neurons at the synaptic cleft through DNA manipulation, and also, epigenetic modifications using Cas9. We explore and examine how each different OMIC techniques can come useful when altering DNA sequences. Such insight into the underlying relationship between OMICs and cellular behaviors leads us to better neurological and psychiatric therapeutics by intelligently designing and utilizing the CRISPR/Cas9 technology.

Genetic and functional variants of the TBX20 gene promoter in dilated cardiomyopathy

BACKGROUND

Dilated cardiomyopathy (DCM) is a major cause of heart failure and sudden cardiac death. As DCM is a genetically heterogeneous disease, genetic variants of cardiac transcription factor genes may play an important role. Transcription factor TBX20, an indispensable factor in normal heart development, is involved in the regulation of cardiac structure and function. Although the TBX20 gene is associated with the occurrence and development of DCM, the influence of genetic variants of the TBX20 gene promoter region on DCM has not been reported.

METHODS

We conducted a case-control study consisting of 107 DCM patients and 210 healthy controls. Genetic variants within TBX20 gene promoter region were identified using sequencing techniques and were functionally analyzed by dual-luciferase reporting assay. Electrophoretic mobility shift assay (EMSA) was used to investigate DNA-protein interactions.

RESULTS

In this study cohort (n = 317), we identified eight variants within TBX20 gene promoter. One novel DNA sequence variants (DSV) (g.4275G>T) and four single-nucleotide polymorphisms (SNPs) [g.4169G>A (rs1263874255), g.4949C>T (rs1191745927), g.5114G>A (rs112076877), g.5252C>T (rs1356932911)] were identified in DCM patients, but in none of controls. Among them, the DSV (g.4275G>T) and three SNPs [g.4949C>T (rs1191745927), g.5114G>A (rs112076877) and g.5252C>T (rs1356932911)] significantly altered the transcription activity of TBX20 gene promoter by dual-luciferase reporting assay (p < 0.05). Further, EMSA assay indicated that the DSV (g.4275G>T) and three SNPs [g.4949C>T (rs1191745927), g.5114G>A (rs112076877) and g.5252C>T (rs1356932911)] affected the binding of transcription factors.

CONCLUSIONS

These data indicate that the DSV (g.4275G>T) and three SNPs [g.4949C>T (rs1191745927), g.5114G>A (rs112076877) and g.5252C>T (rs1356932911)] increase transcription activity of TBX20 gene promoter in both HEK-293 and neonatal rat cardiomyocytes (NRCMs) cell lines by affecting the binding of transcription factors. But the mechanism remains to be verified in vivo.

Genetic Variants of ISL1 Gene Promoter Identified from Congenital Tetralogy of Fallot Patients Alter Cellular Function Forming Disease Basis

Tetralogy of Fallot (TOF) is the most common cyanotic congenital heart disease in newborns. ISL1 is a master transcription factor in second heart field development, whereas the roles of ISL1 gene promoter variants in TOF patients have not been genetically investigated. Total DNA extraction from 601 human subjects, including 308 TOF patients and 293 healthy controls, and Sanger sequencing were performed. Four variants (including one novel heterozygous variant) within the ISL1 gene promoter were only found in TOF patients. Functional analysis of DNA sequence variants was performed by using the dual-luciferase reporter assay and demonstrated that three of the four variants significantly decreased the transcriptional activity of ISL1 gene promoter in HL-1 cells (p < 0.05). Further, the online JASPAR database and electrophoretic mobility shift assay showed that the three variants affected the binding of transcription factors and altered ISL1 expression levels. In conclusion, the current study for the first time demonstrated that the variants identified from the ISL1 gene promoter region are likely involved in the development of TOF by affecting the transcriptional activity and altering the ISL1 expression level. Therefore, these findings may provide new insights into the molecular etiology and potential therapeutic strategy of TOF.

Promoter–motif extraction from co-regulated genes and their relevance to co-expression using E. coli as a model

Gene expression varies due to the intrinsic stochasticity of transcription or as a reaction to external perturbations that generate cellular mutations. Co-regulation, co-expression and functional similarity of substances have been employed for indoctrinating the process of the transcriptional paradigm. The difficult process of analysing complicated proteomes and biological switches has been made easier by technical improvements, and microarray technology has flourished as a viable platform. Therefore, this research enables Microarray to cluster genes that are co-expressed and co-regulated into specific segments. Copious search algorithms have been employed to ascertain diacritic motifs or a combination of motifs that are performing regular expression, and their relevant information corresponding to the gene patterns is also documented. The associated genes co-expression and relevant cis-elements are further explored by engaging Escherichia coli as a model organism. Various clustering algorithms have also been used to generate classes of genes with similar expression profiles. A promoter database ‘EcoPromDB’ has been developed by referring RegulonDB database; this promoter database is freely available at www.ecopromdb.eminentbio.com and is divided into two sub-groups, depending upon the results of co-expression and co-regulation analyses.

Intronic variants in inborn errors of metabolism: Beyond the exome

Non-coding regions are areas of the genome that do not directly encode protein and were initially thought to be of little biological relevance. However, subsequent identification of pathogenic variants in these regions indicates there are exceptions to this assertion. With the increasing availability of next generation sequencing, variants in non-coding regions are often considered when no causative exonic changes have been identified. There is still a lack of understanding of normal human variation in non-coding areas. As a result, potentially pathogenic non-coding variants are initially classified as variants of uncertain significance or are even overlooked during genomic analysis. In most cases where the phenotype is non-specific, clinical suspicion is not sufficient to warrant further exploration of these changes, partly due to the magnitude of non-coding variants identified. In contrast, inborn errors of metabolism (IEMs) are one group of genetic disorders where there is often high phenotypic specificity. The clinical and biochemical features seen often result in a narrow list of diagnostic possibilities. In this context, there have been numerous cases in which suspicion of a particular IEM led to the discovery of a variant in a non-coding region. We present four patients with IEMs where the molecular aetiology was identified within non-coding regions. Confirmation of the molecular diagnosis is often aided by the clinical and biochemical specificity associated with IEMs. Whilst the clinical severity associated with a non-coding variant can be difficult to predict, obtaining a molecular diagnosis is crucial as it ends diagnostic odysseys and assists in management.

Identification and functional analysis of genetic variants of ISL1 gene promoter in human atrial septal defects

BACKGROUND

Atrial septal defect (ASD) is a common type of congenital heart disease. A gene promoter plays pivotal role in the disease development. This study was designed to investigate the pathological role of variants of the ISL1 gene promoter region in ASD patients.

METHODS

Total DNA extracted from 625 subjects, including 332 ASD patients and 293 healthy controls, was sequenced to identify variants in the promoter region of ISL1 gene. Further functional analyses of the variants were performed with dual luciferase reporter assay and electrophoretic mobility shift assay (EMSA). All possible binding sites of transcription factor affected by the identified variants were predicted using the JASPAR database.

RESULTS

Four variants in the ISL1 gene promoter were found only in patients with ASD by sequencing. Three of the four variants [g.4923 G > C (rs541081886), g.5079 A > G (rs1371835943) and g.5309 G > A (rs116222082)] significantly decreased the transcriptional activities compared with the wild-type ISL1 gene promoter (p < 0.05). The EMSA revealed that these variants [g.4923 G > C (rs541081886), g.5079 A > G (rs1371835943) and g.5309 G > A (rs116222082)] in the ISL1 gene promoter affected the number and affinity of binding sites of transcription factors. Further analysis with the online JASPAR database demonstrated that a cluster of putative binding sites for transcription factors may be altered by these variants.

CONCLUSIONS

These sequence variants identified from the promoter region of ISL1 gene in ASD patients are probably involved in the development of ASD by affecting the transcriptional activity and altering ISL1 levels. Therefore, these findings may provide new insights into the molecular etiology and potential therapeutic strategy of ASD.

The association between PD-1 gene polymorphisms and susceptibility to multiple sclerosis

Programmed cell death 1 (PD-1) is an immune checkpoint and has been reported to be associated with several autoimmune diseases. We aimed to investigate the association between human PD-1 gene (PDCD1) polymorphisms and multiple sclerosis (MS). This case-control study was conducted on 229 MS patients and 246 healthy controls. Genotyping of rs36084323 (PD-1.1 G/A), rs11568821 (PD-1.3 G/A) and rs2227981 (PD-1.5 C/T) polymorphisms was performed by PCR-RFLP technique. The frequency difference of PD-1.1 genotypes and alleles (-536 G/A) between patients and healthy controls was not significant. Regarding PD-1.3, the AA + AG genotype was found to be relatively higher in the control group. Concerning PD-1.5 (+7785 C/T), the frequency of T allele carriers (TT + CT) was relatively higher in MS patients, which was marginally insignificant (p = .07). PD-1 gene polymorphisms may be associated with MS; however, accurate conclusions require further studies with a larger number of samples.

Genetic Variants of CITED2 Gene Promoter in Human Atrial Septal Defects: Case-Control Study and Cellular Functional Verification

Atrial septal defect (ASD) is one of the most common forms of congenital heart disease (CHD). Genetic variants in the coding region of the CITED2 gene are known to be significantly correlated with CHD, but the role of variants in the promoter region of CITED2 is unknown. We investigated variants in the promoter of the CITED2 gene in 625 subjects (332 ASD and 293 healthy controls) through Sanger sequencing. Four variants in the CITED2 gene promoter were found only in eight ASD patients with zero occurrence in the control subjects (one case of g.4078A>C(rs1165649373), one case of g.4240C>A(rs1235857801), four cases of g.4935C>T(rs111470468), two cases of g.5027C>T(rs112831934)). Cellular functional analysis showed that these four variants significantly changed the transcriptional activity of the CITED2 gene promoter in HEK-293 and HL-1 cells. Electrophoretic mobility change assay results and JASPAR database analysis demonstrated that these variants created or destroyed a series of possible transcription factor binding sites, resulting in changes in the expression of CITED2 protein. We conclude that the variants of CITED2 promoter in ASD patients affect the transcriptional activity and are likely involved in the occurrence and development of ASD. These findings provide new perspectives on the pathogenesis and potential therapeutic insights of ASD.

Splicing Analysis of 16 PALB2 ClinVar Variants by Minigene Assays: Identification of Six Likely Pathogenic Variants

PALB2 loss-of-function variants are associated with significant increased risk of breast cancer as well as other types of tumors. Likewise, splicing disruptions are a common mechanism of disease susceptibility. Indeed, we previously showed, by minigene assays, that 35 out of 42 PALB2 variants impaired splicing. Taking advantage of one of these constructs (mgPALB2_ex1-3), we proceeded to analyze other variants at exons 1 to 3 reported at the ClinVar database. Thirty-one variants were bioinformatically analyzed with MaxEntScan and SpliceAI. Then, 16 variants were selected for subsequent RNA assays. We identified a total of 12 spliceogenic variants, 11 of which did not produce any trace of the expected minigene full-length transcript. Interestingly, variant c.49-1G > A mimicked previous outcomes in patient RNA (transcript ∆(E2p6)), supporting the reproducibility of the minigene approach. A total of eight variant-induced transcripts were characterized, three of which (∆(E1q17), ∆(E3p11), and ∆(E3)) were predicted to introduce a premature termination codon and to undergo nonsense-mediated decay, and five (▼(E1q9), ∆(E2p6), ∆(E2), ▼(E3q48)-a, and ▼(E3q48)-b) maintained the reading frame. According to an ACMG/AMP (American College of Medical Genetics and Genomics/Association for Molecular Pathology)-based classification scheme, which integrates mgPALB2 data, six PALB2 variants were classified as pathogenic/likely pathogenic, five as VUS, and five as likely benign. Furthermore, five ±1,2 variants were catalogued as VUS because they produced significant proportions of in-frame transcripts of unknown impact on protein function.

Functional significance of novel variants of the MEF2C gene promoter in congenital ventricular septal defects

Ventricular septal defect (VSD) is the most common congenital heart disease. Although the coding region of MEF2C is highly relevant to cardiac malformations, the role of MEF2C gene promoter variants in VSD patients has not been genetically investigated. We investigated the role of MEF2C gene promoter variants in 400 Han Chinese subjects (200 patients with isolated and sporadic VSD and 200 healthy controls). The promoter region of the MEF2C gene was sequenced that identified 10 variants. Expression vectors encompassing the variants and the firefly luciferase reporter gene plasmid (pGL3-basic) were constructed and subsequently transfected into HEK-293 cells. The luciferase activities were measured by Dual-luciferase reporter assay system. MEF2C gene promoter transcriptional activity was significantly reduced in 4 of the 10 variants in HEK-293 cells (P < 0.05). In addition, the JASPAR database was used to perform bioinformatics analysis, which showed that these variants disrupt the putative binding sites of transcription factors and affected the expression of MEF2C protein. This study for the first time identified the variants in the promoter of the MEF2C gene in Han Chinese population and revealed the role of these variants in the formation of VSD.

Splicing predictions, minigene analyses, and ACMG-AMP clinical classification of 42 germline PALB2 splice-site variants

PALB2 loss-of-function variants confer high risk of developing breast cancer. Here we present a systematic functional analysis of PALB2 splice-site variants detected in approximately 113,000 women in the large-scale sequencing project Breast Cancer After Diagnostic Gene Sequencing (BRIDGES; https://bridges-research.eu/). Eighty-two PALB2 variants at the intron-exon boundaries were analyzed with MaxEntScan. Forty-two variants were selected for the subsequent splicing functional assays. For this purpose, three splicing reporter minigenes comprising exons 1-12 were constructed. The 42 potential spliceogenic variants were introduced into the minigenes by site-directed mutagenesis and assayed in MCF-7/MDA-MB-231 cells. Splicing anomalies were observed in 35 variants, 23 of which showed no traces or minimal amounts of the expected full-length transcripts of each minigene. More than 30 different variant-induced transcripts were characterized, 23 of which were predicted to truncate the PALB2 protein. The pathogenicity of all variants was interpreted according to an in-house adaptation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology (ACMG-AMP) variant classification scheme. Up to 23 variants were classified as pathogenic/likely pathogenic. Remarkably, three ±1,2 variants (c.49-2A>T, c.108+2T>C, and c.211+1G>A) were classified as variants of unknown significance, as they produced significant amounts of either in-frame transcripts of unknown impact on the PALB2 protein function or the minigene full-length transcripts. In conclusion, we have significantly contributed to the ongoing effort of identifying spliceogenic variants in the clinically relevant PALB2 cancer susceptibility gene. Moreover, we suggest some approaches to classify the findings in accordance with the ACMG-AMP rationale. © 2021 The Authors. The Journal of Pathology published by John Wiley & Sons, Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Comparative Analyses of Sperm DNA Methylomes Among Three Commercial Pig Breeds Reveal Vital Hypomethylated Regions Associated With Spermatogenesis and Embryonic Development

Identifying epigenetic changes is essential for an in-depth understanding of phenotypic diversity and pigs as the human medical model for anatomizing complex diseases. Abnormal sperm DNA methylation can lead to male infertility, fetal development failure, and affect the phenotypic traits of offspring. However, the whole genome epigenome map in pig sperm is lacking to date. In this study, we profiled methylation levels of cytosine in three commercial pig breeds, Landrace, Duroc, and Large White using whole-genome bisulfite sequencing (WGBS). The results showed that the correlation of methylation levels between Landrace and Large White pigs was higher. We found that 1,040-1,666 breed-specific hypomethylated regions (HMRs) were associated with embryonic developmental and economically complex traits for each breed. By integrating reduced representation bisulfite sequencing (RRBS) public data of pig testis, 1743 conservated HMRs between sperm and testis were defined, which may play a role in spermatogenesis. In addition, we found that the DNA methylation patterns of human and pig sperm showed high similarity by integrating public data from WGBS and chromatin immunoprecipitation sequencing (ChIP-seq) in other mammals, such as human and mouse. We identified 2,733 conserved HMRs between human and pig involved in organ development and brain-related traits, such as NLGN1 (neuroligin 1) containing a conserved-HMR between human and pig. Our results revealed the similarities and diversity of sperm methylation patterns among three commercial pig breeds and between human and pig. These findings are beneficial for elucidating the mechanism of male fertility, and the changes in commercial traits that undergo strong selection.

Growth, Dietary Intake, and Vitamin D Receptor (VDR) Promoter Genotype in Indonesian School-Age Children

Nutrition has been known as a predominant factor associated with stunting. However, some studies have discovered a genetic contribution in calcium absorption that will affect growth, known as the VDR gene. The aim of this study was to assess the association between VDR gene polymorphism and dietary intake towards height-for-age z-score (HAZ) of elementary school children in Malang District, East Java. This study analyzed the baseline of a randomized trial in East Java, Indonesia. School children aged 8–10 years old (n = 142) were included in this study. Energy, protein, calcium, and vitamin D intakes were obtained using 4-day 24-h dietary recalls. Two SNPs located in the promoter region of VDR gene were selected (rs11568820 and rs4516035) and analyzed using Real-Time PCR. The result showed a significant correlation between energy and protein intake with HAZ of the children (p = 0.030 and p = 0.016, respectively). The association between VDR gene and HAZ was not found (p > 0.05). Adjusted by other factors, protein intake was significantly correlated with HAZ (β = 0.034, 95% CI 0.015–0.052, p < 0.001, adj. R² = 0.089). The children in our study had a favorable VDR gene genotype, however the effect of VDR gene promoter activity might not be revealed due to very low vitamin D and calcium intake to stimulate intestinal calcium absorption which in turn affects HAZ.

Emerging roles of alternative cleavage and polyadenylation (APA) in human disease

In the messenger RNA (mRNA) maturation process, the 3'-end of pre-mRNA is cleaved and a poly(A) sequence is added, this is an important determinant of mRNA stability and its cellular functions. More than 60%-70% of human genes have three or more polyadenylation (APA) sites and can be cleaved at different sites, generating mRNA transcripts of varying lengths. This phenomenon is termed as alternative cleavage and polyadenylation (APA) and it plays role in key biological processes like gene regulation, cell proliferation, senescence, and also in various human diseases. Loss of regulatory microRNA binding sites and interactions with RNA-binding proteins leading to APA are largely investigated in human diseases. However, the functions of the core APA machinery and related factors during disease conditions remain largely unknown. In this review, we discuss the roles of polyadenylation machinery in relation to brain disease, cardiac failure, pulmonary fibrosis, cancer, infectious conditions, and other human diseases. Collectively, we believe this review will be a useful avenue for understanding the emerging role of APA in the pathobiology of various human diseases.

Identification of variants of ISL1 gene promoter and cellular functions in isolated ventricular septal defects

Ventricular septal defects (VSDs) are the most common congenital heart defects (CHDs). Studies have documented that ISL1 has a crucial impact on cardiac growth, but the role of variants in the ISL1 gene promoter in patients with VSD has not been explored. In 400 subjects (200 patients with isolated and sporadic VSDs: 200 healthy controls), we investigated the ISL1 gene promoter variant and performed cellular functional experiments by using the dual-luciferase reporter assay to verify the impact on gene expression. In the ISL1 promoter, five variants were found only in patients with VSD by sequencing. Cellular functional experiments demonstrated that three variants decreased the transcriptional activity of the ISL1 promoter (P < 0.05). Further analysis with the online JASPAR database demonstrated that a cluster of putative binding sites for transcription factors may be altered by these variants, possibly resulting in change of ISL1 protein expression and VSD formation. Our study has, for the first time, identified novel variants in the ISL1 gene promoter region in the Han Chinese patients with isolated and sporadic VSD. In addition, the cellular functional experiments, electrophoretic mobility shift assay, and bioinformatic analysis have demonstrated that these variants significantly alter the expression of the ISL1 gene and affect the binding of transcription factors, likely resulting in VSD. Therefore, this study may provide new insights into the role of the gene promoter region for a better understanding of genetic basis of the formation of CHDs and may promote further investigations on mechanism of the formation of CHDs.

The Role of RNA-Sequencing as a New Genetic Diagnosis Tool

Purpose of Review

Whole exome sequencing (WES) and whole-genome sequencing (WGS) are frontline approaches for the genetic diagnosis of rare diseases. However, WES/WGS fails in up to 75% of cases. Transcriptomics via RNA-sequencing (RNA-Seq) is a novel approach that aims to increase the diagnostic yield in rare diseases.

Recent Findings

Recent publications focus on the success of RNA-Seq for increasing diagnosis rates in WES/WGS-negative patients in up to 36% of cases, across a range of different diseases, sample sizes, and tissue types.

Summary

RNA-Seq is beneficial for aiding prioritisation of causative variants currently not detected or often overlooked by WES/WGS alone. An improvement in diagnostic yields has been demonstrated using multiple source tissues, with muscle and fibroblasts being the most representative, but the more accessible blood still demonstrating diagnostic success, particularly in neuromuscular disorders. The introduction of RNA-Seq to the genetic diagnosis toolbox promises to be a useful complementary tool to WES/WGS for improving genetic diagnosis in patients with rare disease.

Relevance of PD-L1 Non-Coding Polymorphisms on the Prognosis of a Genetically Admixed NSCLC Cohort

Purpose

Although non-small cell lung cancer (NSCLC) remains a deadly disease, new predictive biomarkers have emerged to assist in managing the disease, of which one of the most promising is the programmed death‐ligand 1 (PD-L1). Each, PD-L1 variant seem to modulate the function of immune checkpoints differently and affect response to adjuvant treatment and outcome in NSCLC patients. We thus investigated the influence of these PD-L1 genetic variations in genetically admixed NSCLC tissue samples, and correlated these values with clinicopathological characteristics, including prognosis.

Materials and Methods

We evaluated PD-L1 non-coding genetic variants and protein expression in lung adenocarcinomas (ADC), squamous cell carcinomas (SqCC), and large cell carcinomas (LCC) in silico. Microarray paraffin blocks from 70 samples of ADC (N=33), SqCC (N=24), and LCC (N=13) were used to create PD-L1 multiplex immunofluorescence assays with a Cell Signaling E1L3N clone. Fifteen polymorphisms of the PD-L1 gene were investigated by targeted sequencing and evaluated in silico using dedicated tools.

Results

Although PD-L1 polymorphisms seemed not to interfere with protein expression, PD-L1 expression varied among different histological subtypes, as did clinical outcomes, with the rs4742098A>G, rs4143815G>C, and rs7041009G>A variants being associated with relapse (P=0.01; P=0.05; P=0.02, respectively). The rs7041009 GG genotype showed a significant correlation with younger and alive patients compared to carriers of the A allele (P=0.02 and P<0.01, respectively). The Cox regression model showed that the rs7041009 GG genotype may influence OS (P<0.01) as a co-dependent factor associated with radiotherapy and recurrence in NSCLC patients. Furthermore, the Kaplan–Meier survival curves showed that rs7041009 and rs4742098 might impact PPS in relapsed patients. In silico approaches identified the variants as benign.

Conclusion

PD-L1 non-coding variants play an important role in modulating immune checkpoint function and may be explored as immunotherapy biomarkers. We highlight the rs7041009 variant, which impacts OS and PPS in NSCLC patients.

Genetic analysis of the CITED2 gene promoter in isolated and sporadic congenital ventricular septal defects

Ventricular septal defect (VSD) is the most common congenital heart defect. Previous studies have reported genetic variations in the encoding region of CITED2 highly associated with cardiac malformation but the role of CITED2 gene promoter variations in VSD patients has not yet been explored. We investigated the variation of CITED2 gene promoter and its impacts on gene promoter activity in the DNA of paediatric VSD patients. A total of seven variations were identified by Sanger sequencing in the CITED2 gene promoter region in 400 subjects, including 200 isolated and sporadic VSD patients and 200 healthy controls. Using dual‐luciferase reporter assay, we found four of the 7 variations identified significantly decreased the transcriptional activity of the CITED2 gene promoter in HEK‐293 cells (P < .05). Further, a bioinformatic analysis with the JASPAR databases was performed and a cluster of putative binding sites for transcription factors was created or disrupted by these variations, leading to low expression of CITED2 protein and development of VSD. Our study for the first time demonstrates genetic variations in the CITED2 gene promoter in the Han Chinese population and the role of these variations in the development of VSD, providing new insights into the aetiology of CHD.

Association of XRCC2 rs2040639 with the survival of patients with oral squamous cell carcinoma undergoing concurrent chemoradiotherapy

OBJECTIVE

To investigate the association between the variants of DNA double-strand break repair genes and the clinical outcomes of patients with oral squamous cell carcinoma (OSCC) undergoing concurrent chemoradiotherapy.

METHODS

Five variants of DNA double-strand break repair genes in samples from 319 patients with OSCC were genotyped using the Sequenom iPLEX MassARRAY system. Kaplan-Meier curves and Cox proportional hazards analysis were used to identify the factors associated with patient survival.

RESULTS

The XRCC2 rs2040639 (G3063A) polymorphism in the codominant model was associated with decreased recurrence risk (hazard ratio [HR] = 0.55, 95% confidence interval [CI] = 0.31-0.98; p = 0.042). A marginally significant interaction was observed between XRCC2 rs2040639 and PRKDC rs7003908 in patients carrying the AA and AA genotypes; these patients showed reduced recurrence risk (HR = 0.36, 95% CI = 0.17-0.79; p = 0.010).

CONCLUSION

The A-allele of XRCC2 rs2040639 is a favorable prognostic factor for disease-free survival. Patients with these genotypes may benefit from concurrent chemoradiotherapy. Additional confirmation from studies with larger samples or other ethnic populations is warranted.

Risk Stratification Using a Novel Genetic Classifier Including PLEKHS1 Promoter Mutations for Differentiated Thyroid Cancer with Distant Metastasis

Background: Although most differentiated thyroid carcinomas (DTCs) have indolent behavior, DTCs with distant metastasis have a poor prognosis. However, there are no validated markers that predict the risk of distant metastasis and the prognosis of DTC. We aimed to develop a genetic classifier for predicting the outcomes of DTC patients with distant metastases. Methods: Targeted deep sequencing of 157 cancer-related genes was performed for 61 DTCs with distant metastases. A candidate mutation was validated with independent thyroid cancer samples using digital polymerase chain reaction. Results: The most frequently mutated gene in the 61 DTCs was BRAF (n = 31, 51%), followed by TERT promoter (n = 28, 46%), NRAS (n = 13, 11%), PLEKHS1 promoter (n = 6, 10%), and STK11 (n = 6, 10%) mutations. PLEKHS1 promoter mutations were more common in the radioactive iodine (RAI)-refractory cases (p = 0.003). Losses of 9q and 11q were associated with RAI-refractory disease (p = 0.002) and cancer-specific mortality (p = 0.028), respectively. In multivariate analysis, bone metastasis (adjusted odds ratio [aOR] = 15.17, 95% confidence interval [CI 3.38-68.06], p < 0.001) and at least one mutation in the TERT promoter, the PLEKHS1 promoter, or TP53 (aOR = 7.64 [CI 1.78-32.76], p = 0.006) remained significant factors associated with RAI-refractoriness. In independently collected papillary thyroid carcinomas without initial distant metastasis (n = 75), a PLEKHS1 promoter mutation was only found in one case that developed distant metastasis during the follow-up period. We developed a genetic classifier consisting of BRAF, RAS, the TERT promoter, the PLEKHS1 promoter, and TP53 for categorizing the prognosis of patients with DTC with distant metastasis. In the poor-prognosis group, 61% of the patients were RAI-refractory and death occurred in 21% during the follow-up. In the intermediate-prognosis group, 29% were RAI-refractory, but no death occurred. In the good-prognosis group, all patients were RAI-responsive and no death occurred. Conclusions: Mutations in the PLEKHS1 promoter are a novel genetic marker of aggressive DTC. Our genetic classifier can be useful for predicting RAI-refractory disease and poor prognosis in DTC patients with distant metastases.

A LINE-1 insertion situated in the promoter of IMPG2 is associated with autosomal recessive progressive retinal atrophy in Lhasa Apso dogs

Background

Canine progressive retinal atrophies are a group of hereditary retinal degenerations in dogs characterised by depletion of photoreceptor cells in the retina, which ultimately leads to blindness. PRA in the Lhasa Apso (LA) dog has not previously been clinically characterised or described in the literature, but owners in the UK are advised to have their dog examined through the British Veterinary Association/ Kennel Club/ International Sheep Dog Society (BVA/KC/ISDS) eye scheme annually, and similar schemes that are in operation in other countries. After the exclusion of 25 previously reported canine retinal mutations in LA PRA-affected dogs, we sought to identify the genetic cause of PRA in this breed.

Results

Analysis of whole-exome sequencing data of three PRA-affected LA and three LA without signs of PRA did not identify any exonic or splice site variants, suggesting the causal variant was non-exonic. We subsequently undertook a genome-wide association study (GWAS), which identified a 1.3 Mb disease-associated region on canine chromosome 33, followed by whole-genome sequencing analysis that revealed a long interspersed element-1 (LINE-1) insertion upstream of the IMPG2 gene. IMPG2 has previously been implicated in human retinal disease; however, until now no canine PRAs have been associated with this gene. The identification of this PRA-associated variant has enabled the development of a DNA test for this form of PRA in the breed, here termed PRA4 to distinguish it from other forms of PRA described in other breeds. This test has been used to determine the genotypes of over 900 LA dogs. A large cohort of genotyped dogs was used to estimate the allele frequency as between 0.07–0.1 in the UK LA population.

Conclusions

Through the use of GWAS and subsequent sequencing of a PRA case, we have identified a LINE-1 insertion in the retinal candidate gene IMPG2 that is associated with a form of PRA in the LA dog. Validation of this variant in 447 dogs of 123 breeds determined it was private to LA dogs. We envisage that, over time, the developed DNA test will offer breeders the opportunity to avoid producing dogs affected with this form of PRA.

Whole genome sequencing analysis of high confidence variants of B-cell lymphoma in Canis familiaris

Lymphoma (lymphosarcoma) is the second most frequent cancer in dogs and is clinically comparable to human non-Hodgkin lymphoma. Factors affecting canine lymphoma progression are unknown and complex, but there is evidence that genetic mutations play an important role. We employed Next Gen DNA sequencing of six dogs with multicentric B-cell lymphoma undergoing CHOP chemotherapy to identify genetic variations potentially impacting response. Paired samples from non-neoplastic tissue (blood mononuclear cells) and lymphoma were collected at the time of diagnosis. Cases with progression free survival above the median of 231 days were grouped as 'good' responders and cases below the median were categorized as 'poor' responders. The average number of variants found was 17,138 per case. The variants were filtered to examine those with predicted moderate or high impacts. Many of the genes with variants had human orthologs with links to cancer, but the majority of variants were not previously reported in canine or human lymphoma. Seven genes had variants found in the cancers of at least two 'poor' responders but in no 'good' responders: ATRNL1, BAIAP2L2, ZNF384, ST6GALNAC5, ENSCAFG00000030179 (human ortholog: riboflavin kinase RFK), ENSCAFG00000029320, and ENSCAFG00000007370 (human ortholog: immunoglobin IGKV4-1). Two genes had variants found in the cancers of at least two 'good' responders but in no 'poor' responders: COX18 and ENSCAFG00000030512. ENSCAFG00000030512 has no reported orthologue in any other species. The role of these mutations in the progression of canine lymphoma requires further functional analyses and larger scale study.

PD1 gene polymorphism is associated with a poor prognosis in hepatocellular carcinoma following liver resection, cohort study

BACKGROUND

This study examined whether single nucleotide polymorphism (SNP) in programmed cell death protein (PD)-1 is related to the postoperative prognosis of patients with hepatocellular carcinoma (HCC). The immune checkpoint protein PD-1 is an important inhibitor of T cell responses. SNP in the promoter region of PD-1 -606 G/A has been reported to result in high activation and expression of PD-1 associated with cancer risk.

MATERIALS AND METHODS

We analyzed 321 patients with HCC who underwent hepatectomy between 2010 and 2015. PD-1 SNP was analyzed by polymerase chain reaction, and the prognosis after surgical treatment of patients with HCC was analyzed.

RESULTS

The PD-1 SNP statuses were as follows: 90 AA (28.1%), 163 GA (50.8%), 68 GG (21.2%). The baseline parameters did not statistically differ between the three groups. The overall survival (OS) of patients with the GG genotype was significantly lower than that of those with the other genotypes (P = 0.031). The GG genotype was an independent risk factor for OS (P = 0.009; HR 2.201). There was no significant difference between the GG genotype and other genotypes in recurrent-free survival. The extrahepatic recurrence (EHR) rate of those with the GG genotype was significantly higher than that of those with the other genotypes (P = 0.036). The GG genotype was an independent risk factor for EHR (P = 0.008; HR 2.037).

CONCLUSIONS

The PD-1 SNP GG genotype is associated with poor survival and increased EHR in HCC. Furthermore, the GG genotype is an independent predictive factor for OS and EHR.

Association between PD-1 and PD-L1 Polymorphisms and the Risk of Cancer: A Meta-Analysis of Case-Control Studies

A number of case-control studies regarding the association of the polymorphisms in the programmed cell death 1 (PD-1) and programmed cell death ligand 1 (PD-L1) genes with the risk of cancer have yielded inconsistent findings. Therefore, we have conducted a comprehensive, updated meta-analysis study to identify the impact of PD-1 and PD-L1 polymorphisms on overall cancer susceptibility. The findings revealed that PD-1 rs2227981 and rs11568821 polymorphisms significantly decreased the overall cancer risk (Odds Ratio (OR) = 0.82, 95% CI = 0.68–0.99, p = 0.04, TT vs. CT+CC; OR = 0.79, 95% CI = 0.67–0.94, p = 0.006, AG vs. GG, and OR = 0.82, 95% CI = 0.70–0.96, p = 0.020, AG+AA vs. GG, respectively), while PD-1 rs7421861 polymorphism significantly increased the risk of developing cancer (OR = 1.16, 95% CI = 1.02–1.33, p = 0.03, CT vs. TT). The PD-L1 rs4143815 variant significantly decreased the risk of cancer in homozygous (OR = 0.62, 95% CI = 0.41–0.94, p = 0.02), dominant (OR = 0.70, 95% CI = 0.50–0.97, p = 0.03), recessive (OR = 0.76, 95% CI = 0.60–0.96, p = 0.02), and allele (OR = 0.78, 95% CI = 0.63–0.96, p = 0.02) genetic models. No significant association between rs2227982, rs36084323, rs10204525, and rs2890658 polymorphisms and overall cancer risk has been found. In conclusions, the results of this meta-analysis have revealed an association between PD-1 rs2227981, rs11568821, rs7421861, as well as PD-L1 rs4143815 polymorphisms and overall cancer susceptibility.

Computational Analysis of Protein Structure Changes as a Result of Nondeletion Insertion Mutations in Human β-Globin Gene Suggests Possible Cause of β-Thalassemia

Beta-thalassemia is described as a group of hereditary blood disorders characterized by abnormalities in the synthesis of beta chains of hemoglobin. These anomalies result in different phenotypes ranging from moderate to severe clinical symptoms to no symptoms at all. Most of the defects in hemoglobin arise directly from the mutations in the structural β-globin gene (HBB). Recent advances in computational tools have allowed the study of the relationship between the genotype and phenotype in many diseases including β-thalassemia. Due to high prevalence of β-thalassemia, these analyses have helped to understand the molecular basis of the disease in a better way. In this direction, a relational database, named HbVar, was developed in 2001 by a collective academic effort to provide quality and up-to-date information on the genomic variations leading to hemoglobinopathies and thalassemia. The database recorded details about each variant including the altered sequence, hematological defects, its pathology, and its occurrence along with references. In the present study, an attempt was made to investigate nondeletion mutations in the HBB picked up from HbVar and their effects using the in silico approach. Our study investigated 12 nucleotides insertion mutations in six different altered sequences. These 12 extra nucleotides led to the formation of a loop in the protein structure and did not alter its function. It appears that these mutations act as ‘silent' mutations. However, further in vitro studies are required to reach definitive conclusions.

Systematic analysis of dark and camouflaged genes reveals disease-relevant genes hiding in plain sight

BACKGROUND

The human genome contains "dark" gene regions that cannot be adequately assembled or aligned using standard short-read sequencing technologies, preventing researchers from identifying mutations within these gene regions that may be relevant to human disease. Here, we identify regions with few mappable reads that we call dark by depth, and others that have ambiguous alignment, called camouflaged. We assess how well long-read or linked-read technologies resolve these regions.

RESULTS

Based on standard whole-genome Illumina sequencing data, we identify 36,794 dark regions in 6054 gene bodies from pathways important to human health, development, and reproduction. Of these gene bodies, 8.7% are completely dark and 35.2% are ≥ 5% dark. We identify dark regions that are present in protein-coding exons across 748 genes. Linked-read or long-read sequencing technologies from 10x Genomics, PacBio, and Oxford Nanopore Technologies reduce dark protein-coding regions to approximately 50.5%, 35.6%, and 9.6%, respectively. We present an algorithm to resolve most camouflaged regions and apply it to the Alzheimer's Disease Sequencing Project. We rescue a rare ten-nucleotide frameshift deletion in CR1, a top Alzheimer's disease gene, found in disease cases but not in controls.

CONCLUSIONS

While we could not formally assess the association of the CR1 frameshift mutation with Alzheimer's disease due to insufficient sample-size, we believe it merits investigating in a larger cohort. There remain thousands of potentially important genomic regions overlooked by short-read sequencing that are largely resolved by long-read technologies.

Single nucleotide and structural variants of CYP2D6 gene in Kinh Vietnamese population

CYP2D6 genetic variations could result in alteration of CYP2D6 enzyme activity, leading to dissimilarity among individuals in regard of drug metabolism.This study aims to detect all genetic variants, allele, and genotype frequencies of CYP2D6 gene in 136 unrelated healthy Kinh Vietnamese volunteers. All single nucleotide variants (SNVs) and structural variations (SVs) of CYP2D6 gene were identified by Sanger sequencing and multiplex ligation-dependent probe amplification (MLPA) assay.Totally, 30 SNVs and 9 SVs including a whole gene deletion, 8 hybrid structures, and tandem arrangements were identified. Of the 7 novel SNVs detected, the 3157G>T (R329L) substitution was predicted to be deleterious by PROVEAN; the 3851G>A (W358X) variant resulted in a truncated protein; and the 2988G>A variant located in the intron 6 was predicted to be capable of modifying splicing motif by Human Splicing Finder. We determined 29 different genotypes of CYP2D6 from 136 individuals. The most common alleles were the CYP2D6*10 (43.75%), *1 (18.75%), and tandem arrangement *36-*10 (12.13%).This study provides best information on CYP2D6 polymorphism comprising the newly discovered SNVs, structural variations, and their frequencies in Kinh Vietnamese. These new data would be valuable in view of precise dosing of CYP2D6 metabolized drugs and giving better treatment outcome.

An MTF1 binding site disrupted by a homozygous variant in the promoter of ATP7B likely causes Wilson Disease

Approximately 2% of the human genome accounts for protein-coding genes, yet most known Mendelian disease-causing variants lie in exons or splice sites. Individuals who symptomatically present with monogenic disorders but do not possess function-altering variants in the protein-coding regions of causative genes may harbor variants in the surrounding gene regulatory domains. We present such a case: a male of Afghani descent was clinically diagnosed with Wilson Disease-a disorder of systemic copper buildup-but was found to have no function-altering coding variants in ATP7B (ENST00000242839.4), the typically causative gene. Our analysis revealed the homozygous variant chr13:g.52,586,149T>C (NC_000013.10, hg19) 676 bp into the ATP7B promoter, which disrupts a metal regulatory transcription factor 1 (MTF1) binding site and diminishes expression of ATP7B in response to copper intake, likely resulting in Wilson Disease. Our approach to identify the causative variant can be generalized to systematically discover function-altering non-coding variants underlying disease and motivates evaluation of gene regulatory variants.

The PD-1 rs36084323 A > G polymorphism decrease cancer risk in Asian: A meta-analysis

The rs36084323 A > G polymorphism in programmed cell death-1(PD-1) gene has been reported to be associated with cancer risk. However, the results of previous studies were inconsistent. Therefore, we performed a meta-analysis to identify the potential association, by searching the PubMed, EMBASE, Cochrane Library, and the Chinese CNKI, WANFANG and CBM databases. Data were extracted and odds ratios (ORs) with 95% confidence intervals (CIs) were estimated to assess the strength of the association. A total of 10 relevant studies involving 4445 cancer cases and 5126 controls were recruited. Overall, the results indicated that the PD-1 rs36084323 A > G polymorphism was not statistically associated with cancer risk. However, stratified analysis revealed that there was a statistically reduced cancer risk in Asians(G vs. A, OR = 0.89, 95%CI:0.81-0.97, P = 0.008, I² = 48.8%; GG vs. AA, OR = 0.79, 95% CI:0.66-0.94, P = 0.008, I² = 48.7%; GG/AG vs. AA, OR = 0.87, 95%CI:0.76-0.98, P = 0.017, I² = 34.9%; GG vs. AG/AA, OR = 0.85, 95%CI:0.75-0.97, P = 0.027, I² = 40%) and in the patients with EOC(AG vs. AA, OR = 0.69, 95%CI:0.54-0.90, P = 0.005, I² = 0%; GG/AG vs. AA, OR = 0.67, 95%CI:0.52-0.85, P = 0.001, I² = 0). Meta-regression showed that ethnicity (P = 0.029) but not cancer types (P = 0.792), source of controls (P = 0.207) or ample size (P = 0.585) were the sources of heterogeneity. This meta-analysis demonstrates the PD-1 rs36084323 A > G polymorphism is associated with decreased cancer risk in Asian, and suggests it could potentially serve as a biomarker to screen high-risk individuals. Large-scale and well-designed case-control studies are needed to enrich the evidence of this result.

Promoter Variant Alters Expression of the Autophagic BECN1 Gene: Implications for Clinical Manifestations of Machado-Joseph Disease

Autophagy is especially important in disorders where accumulation of the mutant protein is a hallmark, such as the Machado-Joseph disease/spinocerebellar ataxia type 3 (MJD/SCA3). We analyzed the promoter of the BECN1 gene, whose overexpression has been reported to exert neuroprotective effects in MJD, with the aim of finding variants that could be associated with expression levels of beclin-1 and could be tested as modifiers of onset and disease severity. A fragment encompassing the BECN1 promoter was sequenced in 95 MJD subjects and 120 controls. The impact of the variation detected on transcription factors (TFs) binding affinity was evaluated in silico and inferences concerning levels of expression were confirmed by fluorescence-based quantitative real-time PCR in a subset of 28 MJD subjects and 26 controls. Four previously described (rs60221525, rs116943570, rs34882610, and rs34037822) and one novel (c.-933delG) variants were identified. In silico analysis performed for the most frequent variants-rs60221525 C allele and rs116943570 T allele-predicted an impact of the presence of these alleles on TF binding affinity. BECN1 expression levels were in agreement with the in silico predictions, showing a tendency for decreased levels in samples with the rs60221525 C allele and for increased levels in samples with the rs116943570 T allele. MJD patients carrying the rs60221525 C allele presented a tendency for earlier estimated age at onset. Moreover, patients with the rs60221525 C allele presented a more severe clinical picture, compared to patients without this allele. The analysis of a larger number of patients from different cohorts, currently unavailable, would be required to confirm these results.

Programmed death-1 polymorphisms is associated with risk of esophagogastric junction adenocarcinoma in the Chinese Han population: A case-control study involving 2,740 subjects

Single nucleotide polymorphisms (SNPs) in Programmed cell death 1 (PD-1) gene may contribute to the development of cancer. In this study, we selected PD-1 rs10204525 T>C, rs2227982 A>G, rs36084323 T>C and rs7421861 A>G polymorphisms and designed a hospital-based case-control study to determine the potential relationship between these functional SNPs in PD-1 gene and esophagogastric junction adenocarcinoma (EGJA) risk. A total of 1,063 EGJA patients and 1,677 controls were enrolled from Eastern Chinese Han population. SNPscan^TMgenotyping assay was used to analyze the genotyping of PD-1 polymorphisms. We found that PD-1 rs7421861 A>G polymorphism was associated with the development of EGJA. However, PD-1 rs2227982 A>G polymorphism was a protective factor for EGJA. In addition, PD-1 rs36084323 CC homozygote genotype might be associated with a borderline decreased risk of EGJA. In a subgroup analysis, a decreased risk of EGJA in never drinking and never smoking groups was identified. Haplotype comparison analysis suggested that PD-1 T_rs10204525G_rs2227982C_36084323A_rs7421861 haplotype significantly decreased the risk of EGJA. However, T_rs10204525G_rs2227982C_36084323G_rs7421861 haplotype in PD-1 gene may confer risk to EGJA. In conclusion, our study highlights rs2227982 A>G, rs36084323 T>C and rs7421861 A>G polymorphisms and haplotypes in PD-1 gene, especially within the intron region, are significantly associated with the risk of EGJA. Further case-control studies with larger sample size and detailed gene-environmental data to replicate these findings in different populations are needed to validate our conclusion.

Lysyl oxidase-like-2 mutations and reduced mRNA and protein expression in mid-dermal elastolysis

BACKGROUND

Mid-dermal elastolysis (MDE) is a rare skin condition, characterized by selective loss of elastic fibres in the mid dermis. The pathogenesis of MDE is still unclear.

AIM

To investigate expression of lysyl oxidase-like 2 (LOXL2) in a reasonable sample of patients with MDE and to search for mutations in LOXL2.

METHODS

We investigated archived lesional tissue of 13 patients with MDE and skin tissue samples of 10 sex- and age-matched healthy controls (HCs). Gene and protein expression of LOXL2 was investigated using real-time reverse-transcription PCR and immunohistochemistry. Mutation analysis was performed using the Sanger method.

RESULTS

We observed decreased LOXL2 mRNA expression in lesional skin of patients with MDE (0.48 ± 0.16) compared with healthy skin of the same patients (1.5 ± 0.51) and normal skin of HCs (1.9 ± 0.13). Compared with healthy patient skin (epidermis 2.38 ± 1.6, dermis 1.2 ± 1), LOXL2 protein expression in lesional patient skin (epidermis 1.1 ± 0.7, dermis 0.3 ± 0.45) was significantly decreased (P < 0.04 and P = 0.02, respectively). Mutation analysis of the entire LOXL2 gene could be performed for five patients, all of whom were found to have at least one mutation in the LOXL2 gene. Three of these had a mutation in the promoter region (c.967 G>C, c.1022 C>T, and c.1025 G>A, respectively), and one of them also had a mutation in the splice region of intron 11/exon 12 (IVS11-1 G>A). Of the remaining two patients, one had a mutation in exon 3 (T1391), and the other had a mutation in exon 11 (C663Y).

CONCLUSIONS

Our present data suggest that decreased elastin renewal due to LOXL2 mutations and consecutive reduced LOXL2 expression contribute to the pathogenesis of MDE.

Genetic dissection of the BRCA2 promoter and transcriptional impact of DNA variants

PURPOSE

Promoter mutations may affect transcription and can be associated with human diseases. However, the promoters of the breast cancer (BC) genes are not regularly screened. Our goal was to investigate the BRCA2 promoter in order to study a possible correlation between impaired transcription and disease.

METHODS

The proximal and core promoter of the BRCA2 gene was sequenced in 95 high-risk BC patients. A BRCA2-promoter insert [- 938 to + 312 from the transcription start site (TSS)] was generated and cloned into the firefly luciferase vector pGL4.10. Promoter variants and deletions were introduced by site-directed mutagenesis and quantified by Dual-Luciferase assays and semi-quantitative RT-PCR.

RESULTS

Three different variants were detected in high-risk BC patients: rs3092989, rs206118, and rs563971900. Functional mapping of 13 overlapping deletions revealed four down-regulating segments (TSS positions): -59_-10del/µdel3 (16% of activity of the wild-type construct), -104_-55del/µdel4 (62%), -239_-190del/µdel7 (39%), -464_-415/µdel12 (78%), suggesting the presence therein of putative transcriptional activator motifs. Additionally, six microdeletions rendered luciferase overexpression: +32_+81del/µdel1 (356%), -14_+36del/µdel2 (180%), -194_-145del/µdel6 (154%), -284_-235del/µdel8 (168%), -329_-280del/µdel9 (111%), and -509_-460del/µdel13 (139%), which is indicative of repressor elements. Functional assays of 15 promoter variants (including those detected in patients) showed that ten of them significantly altered expression with seven up-regulating (113-163%) and three down-regulating (rs551887850_G, rs570548398_T, rs55880202_T; 72-83%) SNPs. Eight of them were located in an ENCODE-DNase Hypersensitive Cluster (TSS - 185 to + 105) where most active transcriptional motifs are known to be placed.

CONCLUSIONS

BRCA2 expression is highly sensitive to promoter variations as most of them induced relevant changes. Moreover, we mapped critical regions of the BRCA2 promoter that may constitute potential targets for regulatory variants. Three SNPs moderately decreased luciferase activity, but confirmation of its potential pathogenicity requires further analysis. These data reinforce the need to screen the promoter regions of breast cancer genes with a view to discovering novel deleterious mutations.

Associations between novel genetic variants in the promoter region of MALAT1 and risk of colorectal cancer

The metastasis-associated lung adenocarcinoma transcript 1 (MALAT1), a well-known long non-coding RNA, is involved in pathogenesis and progress of multiple tumors. However, no study has been performed to investigate the relationship between the genetic variants in promoter region of MALAT1 and colorectal cancer risk. In this study, we conducted a two-stage case-control study to evaluate whether MALAT1 genetic variants were associated with colorectal cancer risk. We identified that a single nucleotide polymorphism (SNP) rs1194338 was significantly associated with the decreased colorectal cancer risk with an odds ratio (OR) of 0.70 [95% confidence interval (CI) = 0.49-0.99] in the combined stage. The subsequently stratified analyses showed that the protective effect of rs1194338 was more pronounced in several subgroups. Furthermore, gene expression profiling analysis revealed overexpression of MALAT1 mRNA in colorectal cancer tissue compared with normal controls. Confirmation studies with large sample size and further mechanistic investigations into the function of MALAT1 and its genetic variants are warranted to advance our understanding of their roles in colorectal carcinogenesis, and to aid in the development of novel and targeted therapeutic strategies.

Pathogenic and likely pathogenic genetic alterations and polymorphisms in growth hormone gene (GH1) and growth hormone releasing hormone receptor gene (GHRHR) in a cohort of isolated growth hormone deficient (IGHD) children in Sri Lanka

OBJECTIVE

Genetic alterations in GH1 and GHRHR genes are known to cause isolated growth hormone deficiency (IGHD). Of these, GHRHR codon 72 mutation has been reported to be highly prevalent in the Indian subcontinent, but among Sri Lankans its prevalence was low compared to reports from neighboring countries. The present study was therefore carried out to identify genetic alterations in the GH1 gene and rest of the GHRHR gene in a cohort of Sri Lankan IGHD patients who tested negative for GHRHR codon 72 mutation.

METHODS

Fifty five IGHD children negative for codon 72 (GHRHR) mutation were screened for gross GH1 gene deletion by polymerase chain reaction (PCR) and restriction fragment length polymorphism technique. The coding, intronic and promoter regions of the GH1 gene were sequenced in children who were negative for GH1 deletion (N=53). In a subset (N=40), coding, flanking intronic and promoter regions of the GHRHR gene were screened by single strand conformation polymorphism/sequencing. Identified coding region and intronic variants were subjected to in silico analysis to ascertain pathogenicity. Family members available were screened for the significant variants observed in the index child.

RESULTS

Gross GH1 gene deletions, 6.7kb and 7.0kb were observed in one child each. One novel and 24 reported single nucleotide variants (SNVs) were observed in the GH1 gene and its promoter. These included one reported pathogenic splice site mutation (c.172-2A>T) and one reported likely pathogenic missense mutation (c.406G>T). One large novel deletion of 5875 base pairs that included exon 1, one likely pathogenic novel SNV (c.211G>T) and 18 reported SNVs were observed in the GHRHR gene. Fourteen variants observed were of uncertain significance (8 in GH1 and 6 in GHRHR), twenty three variants were likely benign (11 in GH1 and 12 in GHRHR) and four variants were benign (4 in GH1 and none in GHRHR).

CONCLUSION

In a cohort of IGHD children, six pathogenic or likely pathogenic genetic alterations of either GH1 gene or GHRHR gene were found. These affected a total of six children. Pathogenic status of four of these had been reported in the literature. Novel SNV in the GHRHR gene was predicted to be pathogenic through in silico analysis. The large novel deletion is likely to be pathogenic as it included exon 1 of GHRHR gene. Analysis of other genes will be needed to ascertain the genetic cause of IGHD in the remaining children.

Are programmed cell death 1 gene polymorphisms correlated with susceptibility to rheumatoid arthritis?: A meta-analysis

BACKGROUND

Several studies investigated the relationship between programmed cell death 1 (PDCD1) gene polymorphisms and rheumatoid arthritis (RA) risk, but the results were controversial. To explore whether PDCD1 gene polymorphisms have an effect on RA risk, we conducted this meta-analysis to investigate the relationships between PDCD1 polymorphisms (rs36084323 [PD-1.1 G/A], rs11568821 [PD-1.3 G/A] and rs2227981 [PD-1.5 C/T]) and RA risk under 4 genetic models.

METHODS

PubMed, EMBASE, Web of Science, Cochrane Library China National Knowledge Infrastructure (CNKI), and Chinese Biomedical Literature Database (CBLM) were systematically searched for all eligible case-control studies. The last search was updated on September 10, 2016. Studies were accessed using Newcastle-Ottawa Scale case control study (NOS), and the combined effect size was calculated using STATA software, version 12.0. The pooled odds ratio (OR) with 95% confidence interval (CI) was calculated to assess the association. Heterogeneity analysis and subgroup analysis were also performed. Sensitivity analysis and publication bias were also performed if necessary.

RESULTS

This meta-analysis included 6 studies. The result demonstrated null association between rs36084323 (PD-1.1 G/A) polymorphism and RA susceptibility in all 4 genetic models. With regard to rs11568821 (PD-1.3 G/A), statistically significant association with RA risk was observed under allele model in Caucasians (allele model A vs G, OR = 1.19, 95% CI = 1.03-1.41). There was no significant association between rs2227981 (PD-1.5 C/T) polymorphism and RA risk.

CONCLUSION

The present study suggests that mutant A allele in rs11568821 (PD-1.3 G/A) might increase the susceptibility to RA in Caucasians.

Functional classification of DNA variants by hybrid minigenes: Identification of 30 spliceogenic variants of BRCA2 exons 17 and 18

Mutation screening of the breast cancer genes BRCA1 and BRCA2 identifies a large fraction of variants of uncertain clinical significance (VUS) whose functional and clinical interpretations pose a challenge for genomic medicine. Likewise, an increasing amount of evidence indicates that genetic variants can have deleterious effects on pre-mRNA splicing. Our goal was to investigate the impact on splicing of a set of reported variants of BRCA2 exons 17 and 18 to assess their role in hereditary breast cancer and to identify critical regulatory elements that may constitute hotspots for spliceogenic variants. A splicing reporter minigene with BRCA2 exons 14 to-20 (MGBR2_ex14-20) was constructed in the pSAD vector. Fifty-two candidate variants were selected with splicing prediction programs, introduced in MGBR2_ex14-20 by site-directed mutagenesis and assayed in triplicate in MCF-7 cells. Wild type MGBR2_ex14-20 produced a stable transcript of the expected size (1,806 nucleotides) and structure (V1-[BRCA2_exons_14–20]–V2). Functional mapping by microdeletions revealed essential sequences for exon recognition on the 3’ end of exon 17 (c.7944-7973) and the 5’ end of exon 18 (c.7979-7988, c.7999-8013). Thirty out of the 52 selected variants induced anomalous splicing in minigene assays with >16 different aberrant transcripts, where exon skipping was the most common event. A wide range of splicing motifs were affected including the canonical splice sites (15 variants), novel alternative sites (3 variants), the polypyrimidine tract (3 variants) and enhancers/silencers (9 variants). According to the guidelines of the American College of Medical Genetics and Genomics (ACMG), 20 variants could be classified as pathogenic (c.7806-2A>G, c.7806-1G>A, c.7806-1G>T, c.7806-1_7806-2dup, c.7976+1G>A, c.7977-3_7978del, c.7977-2A>T, c.7977-1G>T, c.7977-1G>C, c.8009C>A, c.8331+1G>T and c.8331+2T>C) or likely pathogenic (c.7806-9T>G, c.7976G>C, c.7976G>A, c.7977-7C>G, c.7985C>G, c.8023A>G, c.8035G>T and c.8331G>A), accounting for 30.8% of all pathogenic/likely pathogenic variants of exons 17–18 at the BRCA Share database. The remaining 8 variants (c.7975A>G, c.7977-6T>G, c.7988A>T, c.7992T>A, c.8007A>G, c.8009C>T, c.8009C>G, and c.8072C>T) induced partial splicing anomalies with important ratios of the full-length transcript (≥70%), so that they remained classified as VUS. Aberrant splicing is therefore especially prevalent in BRCA2 exons 17 and 18 due to the presence of active ESEs involved in exon recognition. Splicing functional assays with minigenes are a valuable strategy for the initial characterization of the splicing outcomes and the subsequent clinical interpretation of variants of any disease-gene, although these results should be checked, whenever possible, against patient RNA.

A significant proportion of disease-causing mutations of inherited disorders impair splicing. Massive sequencing projects of genetic diseases generate thousands of sequence variations that require functional and clinical interpretations. We have shown that splicing reporter minigenes of the breast cancer genes BRCA1 and BRCA2 are useful tools to functionally test DNA variants. In this work, we have constructed a 7-exon BRCA2 minigene (exons 14 to 20) where we mapped critical splicing regulatory sequences and tested 52 selected variants of exons 17 and 18 detected in breast cancer patients. We finely located three DNA segments on both exons that presumably contain splicing enhancer sequences. We observed that a total of 30 variants of any type disrupted the splicing patterns and, given the severity of their outcomes, we classified 20 of them as pathogenic or likely pathogenic. We also showed that a wide range of splicing elements were affected including canonical and novel 5’ and 3’ splice sites, the polypyrimidine tract and enhancer and silencer sequences. We concluded that splicing aberrations are frequent in Hereditary Breast and Ovarian Cancer and that minigenes are valuable tools to functionally classify DNA variants of any human disease gene under the splicing viewpoint.

Single-base substitutions in the CHM promoter as a cause of choroideremia

Although over 150 unique mutations affecting the coding sequence of CHM have been identified in patients with the X-linked chorioretinal disease choroideremia (CHM), no regulatory mutations have been reported, and indeed the promoter has not been defined. Here, we describe two independent families affected by CHM bearing a mutation outside the gene's coding region at position c.-98: C>A and C>T, which segregated with the disease. The male proband of family 1 was found to lack CHM mRNA and its gene product Rab escort protein 1, whereas whole-genome sequencing of an affected male in family 2 excluded the involvement of any other known retinal genes. Both mutations abrogated luciferase activity when inserted into a reporter construct, and by further employing the luciferase reporter system to assay sequences 5' to the gene, we identified the CHM promoter as the region encompassing nucleotides c.-119 to c.-76. These findings suggest that the CHM promoter region should be examined in patients with CHM who lack coding sequence mutations, and reveals, for the first time, features of the gene's regulation.

Autosomal-Dominant Corneal Endothelial Dystrophies CHED1 and PPCD1 Are Allelic Disorders Caused by Non-coding Mutations in the Promoter of OVOL2

Congenital hereditary endothelial dystrophy 1 (CHED1) and posterior polymorphous corneal dystrophy 1 (PPCD1) are autosomal-dominant corneal endothelial dystrophies that have been genetically mapped to overlapping loci on the short arm of chromosome 20. We combined genetic and genomic approaches to identify the cause of disease in extensive pedigrees comprising over 100 affected individuals. After exclusion of pathogenic coding, splice-site, and copy-number variations, a parallel approach using targeted and whole-genome sequencing facilitated the identification of pathogenic variants in a conserved region of the OVOL2 proximal promoter sequence in the index families (c.−339_361dup for CHED1 and c.−370T>C for PPCD1). Direct sequencing of the OVOL2 promoter in other unrelated affected individuals identified two additional mutations within the conserved proximal promoter sequence (c.−274T>G and c.−307T>C). OVOL2 encodes ovo-like zinc finger 2, a C2H2 zinc-finger transcription factor that regulates mesenchymal-to-epithelial transition and acts as a direct transcriptional repressor of the established PPCD-associated gene ZEB1. Interestingly, we did not detect OVOL2 expression in the normal corneal endothelium. Our in vitro data demonstrate that all four mutated OVOL2 promoters exhibited more transcriptional activity than the corresponding wild-type promoter, and we postulate that the mutations identified create cryptic cis-acting regulatory sequence binding sites that drive aberrant OVOL2 expression during endothelial cell development. Our data establish CHED1 and PPCD1 as allelic conditions and show that CHED1 represents the extreme of what can be considered a disease spectrum. They also implicate transcriptional dysregulation of OVOL2 as a common cause of dominantly inherited corneal endothelial dystrophies.

Molecular and cellular analysis of three novel alpha2-globin gene promoter mutations [HBA2: c.-59C>T], [HBA2: c.-81C>A] and [HBA2: c.-91G>A] reveal varying patterns of transcriptional and translational activities

While point mutations affecting the promoter region of β-globin gene are widely described, there are no well characterised reports of any point mutations currently found in the promoter of the α2-globin (HBA2) gene. We present clinical and experimental data for three novel HBA2 gene core and proximal promoter mutations. Using an in vitro system designed to assess the impact of point mutations, the three novel [HBA2:c.-59C>T], [HBA2:c.-81C>A] and [HBA2:c.-91G>A] promoter mutations identified in three unrelated patients were analysed for HBA2 gene transcriptional and translational activities. Following the generation and transfection of expression vectors carrying each mutation, the HBA2 transcription activity of the promoters from each mutant was analysed with quantitative real time-PCR (qReTi-PCR) technique. Immunofluorochemistry (IFC) was used to analyse HBA2 protein synthesis. The analyses showed that [HBA2:c.-59C>T] and [HBA2:c.-91G>A] mutant constructs caused significant reduction in the HBA2 transcription levels by 53.7% (p = 0.0008) and 36.2% (p = 0.004), respectively, resulting in markedly lower HBA2 protein labelling when compared to the wild type as shown with subsequent IFC analysis. Conversely, the [HBA2:c.-81C>A] construct showed no significant changes in either transcription (p = 0.089) or in protein labelling when compared to the wild type. The equal pAmp transcription levels found in each group confirmed that the observed labelling differences were not due to varying transfection efficiencies. This study emphasises the importance of in vitro studies to establish the impact of base substitutions on the level of gene expression, and the value of these studies in clinicopathological correlation so that appropriate advice can be given in genetic counselling.

A novel mutation in the promoter of RARS2 causes pontocerebellar hypoplasia in two siblings

Pontocerebellar hypoplasia (PCH) is characterized by hypoplasia and atrophy of the cerebellum, variable pontine atrophy, microcephaly, severe mental and motor impairments and seizures. Mutations in 11 genes have been reported in 8 out of 10 forms of PCH. Recessive mutations in the mitochondrial arginyl-transfer RNA synthetase gene (RARS2) have been recently associated with PCH type 6, which is characterized by early-onset encephalopathy with signs of oxidative phosphorylation defect. Here we describe the clinical presentation, neuroimaging findings and molecular characterizations of two siblings with a clinical diagnosis of PCH who displayed a novel variant (c.-2A>G) in the 5'-UTR of the RARS2 gene in the homozygous state. This variant was identified through next-generation sequencing testing of a panel of nine genes known to be involved in PCH. Gene expression and functional studies demonstrated that the c.-2A>G sequence change directly leads to a reduced RARS2 messenger RNA expression in the patients by decreasing RARS2 promoter activity, thus providing evidence that mutations in the RARS2 promoter are likely to represent a new causal mechanism of PCH6.

Whole-genome haplotyping approaches and genomic medicine

Genomic information reported as haplotypes rather than genotypes will be increasingly important for personalized medicine. Current technologies generate diploid sequence data that is rarely resolved into its constituent haplotypes. Furthermore, paradigms for thinking about genomic information are based on interpreting genotypes rather than haplotypes. Nevertheless, haplotypes have historically been useful in contexts ranging from population genetics to disease-gene mapping efforts. The main approaches for phasing genomic sequence data are molecular haplotyping, genetic haplotyping, and population-based inference. Long-read sequencing technologies are enabling longer molecular haplotypes, and decreases in the cost of whole-genome sequencing are enabling the sequencing of whole-chromosome genetic haplotypes. Hybrid approaches combining high-throughput short-read assembly with strategic approaches that enable physical or virtual binning of reads into haplotypes are enabling multi-gene haplotypes to be generated from single individuals. These techniques can be further combined with genetic and population approaches. Here, we review advances in whole-genome haplotyping approaches and discuss the importance of haplotypes for genomic medicine. Clinical applications include diagnosis by recognition of compound heterozygosity and by phasing regulatory variation to coding variation. Haplotypes, which are more specific than less complex variants such as single nucleotide variants, also have applications in prognostics and diagnostics, in the analysis of tumors, and in typing tissue for transplantation. Future advances will include technological innovations, the application of standard metrics for evaluating haplotype quality, and the development of databases that link haplotypes to disease.

PMD patient mutations reveal a long-distance intronic interaction that regulates PLP1/DM20 alternative splicing

Alternative splicing of the proteolipid protein 1 gene (PLP1) produces two forms, PLP1 and DM20, due to alternative use of 5' splice sites with the same acceptor site in intron 3. The PLP1 form predominates in central nervous system RNA. Mutations that reduce the ratio of PLP1 to DM20, whether mutant or normal protein is formed, result in the X-linked leukodystrophy Pelizaeus-Merzbacher disease (PMD). We investigated the ability of sequences throughout PLP1 intron 3 to regulate alternative splicing using a splicing minigene construct transfected into the oligodendrocyte cell line, Oli-neu. Our data reveal that the alternative splice of PLP1 is regulated by a long-distance interaction between two highly conserved elements that are separated by 581 bases within the 1071-base intron 3. Further, our data suggest that a base-pairing secondary structure forms between these two elements, and we demonstrate that mutations of either element designed to destabilize the secondary structure decreased the PLP1/DM20 ratio, while swap mutations designed to restore the structure brought the PLP1/DM20 ratio to near normal levels. Sequence analysis of intron 3 in families with clinical symptoms of PMD who did not have coding-region mutations revealed mutations that segregated with disease in three families. We showed that these patient mutations, which potentially destabilize the secondary structure, also reduced the PLP1/DM20 ratio. This is the first report of patient mutations causing disease by disruption of a long-distance intronic interaction controlling alternative splicing. This finding has important implications for molecular diagnostics of PMD.

Expression of PRPF31 and TFPT: regulation in health and retinal disease

PRPF31, a gene located at chromosome 19q13.4, encodes the ubiquitous splicing factor PRPF31. The gene lies in a head-to-head arrangement with TFPT, a poorly characterized gene with a role in cellular apoptosis. Mutations in PRPF31 have been implicated in autosomal dominant retinitis pigmentosa (adRP), a frequent and important cause of blindness worldwide. Disease associated with PRPF31 mutations is unusual, in that there is often non-penetrance of the disease phenotype in affected families, caused by differential expression of PRPF31. This study aimed to characterize the basic promoter elements of PRPF31 and TFPT. Luciferase reporter constructs were made, using genomic DNA from an asymptomatic individual with a heterozygous deletion of the entire putative promoter region. Fragments were tested by the dual-luciferase reporter assay in HeLa and RPE-1 cell lines. A comparison was made between the promoter regions of symptomatic and asymptomatic mutation-carrying individuals. A patient (CAN493) with adRP was identified, harbouring a regulatory region mutation; both alleles were assayed by the dual-luciferase reporter assay. Luciferase assays led to the identification of core promoters for both PRPF31 and TFPT; despite their shared gene architecture, the two genes appear to be controlled by slightly different regulatory regions. One functional polymorphism was identified in the PRPF31 promoter that increased transcriptional activation. The change was not, however, consistent with the observed symptomatic-asymptomatic phenotypes in a family affected by PRPF31-adRP. Analysis of the mutant promoter fragment from CAN493 showed a >50% reduction in promoter activity, suggesting a disease mechanism of functional haploinsufficiency-the first report of this disease mechanism in adRP.