Low Expression of BRCA1 as a Potential Relapse Predictor in B-Cell Acute Lymphoblastic Leukemia.

B-cell acute lymphoblastic leukemia (B-ALL) is the most common childhood hematological malignancy worldwide. Treatment outcomes have improved dramatically in recent years; despite this, relapse is still a problem, and the potential molecular explanation for this remains an important field of study. We performed microarray and single-cell RNA-Seq data mining, and we selected significant data with a P -value<0.05. We validated BRCA1 gene expression by means of quantitative (reverse transcription-polymerase chain reaction.) We performed statistical analysis and considered a P -value<0.05 significant. We identified the overexpression of breast cancer 1, early onset (BRCA1; P -value=2.52 -134 ), by means of microarray analysis. Moreover, the normal distribution of BRCA1 expression in healthy bone marrow. In addition, we confirmed the increases in BRCA1 expression using real-time (reverse transcription-polymerase chain reaction and determined that it was significantly reduced in patients with relapse ( P -values=0.026). Finally, we identified that the expression of the BRCA1 gene could predict early relapse ( P -values=0.01). We determined that low expression of BRCA1 was associated with B-cell acute lymphoblastic leukemia relapse and could be a potential molecular prognostic marker.

The Proliferating Cell Nuclear Antigen (PCNA) Transcript Variants as Potential Relapse Markers in B-Cell Acute Lymphoblastic Leukemia.

Leukemia is the most common childhood malignancy in Mexico, representing more than 50% of all childhood cancers. Although treatment leads to a survival of up to 90% in developing countries, in our country, it is less than 65%. Additionally, ~30% of patients relapse with poor prognosis. Alternative splicing plays an important role in transcriptome diversity and cellular biology. This mechanism promotes an increase in the assortment of proteins with potentially distinct functions from a single gene. The proliferating cell nuclear antigen (PCNA) gene encodes two transcripts for the same protein of 261 amino acids, which is associated with several important cellular processes and with several types of cancer. However, the diversity of the transcript variants expressed in this condition is not clear. Then, we used microarray gene expression to identify changes in the exon expression level of PCNA. The data were validated using RT-PCR and Sanger sequencing, and three additional transcripts (PCNA_V3, PCNA_V4, and PCNA_V5) were identified. Computational analyses were used to determine the potential proteins resulting, their structure, and interactions with PCNA native protein and themselves. Additionally, the PCNA transcript variants were inhibited using specific siRNA, determining that their inhibition contributes to the malignant characteristics in vitro. Finally, we quantified the PCNA transcript variants in acute lymphoblastic leukemia samples and identified their expression in this disease. Based on the clinical characteristics, we determined that PCNA_V2 and PCNA_V4 are expressed at significantly low levels in relapsed B-ALL patients. We conclude that the low expression of PCNA_V2 and PCNA_V4 could be a potential molecular marker of relapse in acute lymphoblastic leukemia patients.

v-myb avian myeloblastosis viral oncogene homolog expression is a potential molecular diagnostic marker for B-cell acute lymphoblastic leukemia.

BACKGROUND: B-cell acute lymphoblastic leukemia (B-ALL) is the most commonly diagnosed childhood malignancy worldwide and is especially common in Mexico. Additionally, the number of cases has increased in recent years. Thus, it is very important to develop molecular strategies to diagnose leukemia. The aim of this study was to investigate MYB expression and to determine its impact on the diagnosis of B-ALL. METHODS: We analyzed the B-ALL gene expression profile by microarray data mining. Bioinformatics analysis was performed to identify the genes that are overexpressed in leukemia. We determined that MYB was highly expressed in leukemia. Then, we validated MYB expression in 70 patients with B-ALL and in 16 healthy controls (HCs) using qRT-PCR. The results were statistically analyzed using the Kolmogorov-Smirnov Z test, Mann-Whitney U test, receiver operating characteristic curves, and the Youden index. RESULTS: The microarrays showed that MYB was overexpressed in B-ALL patients with a fold change of 57.8728 and a P value of 2.56-195 . MYB expression showed great variability among the patients analyzed. However, compared to the HCs, the B-ALL patients had a P value < .0001, an area under the curve of 0.813, and a Youden index of 1.46, indicating the statistical significance. CONCLUSION: MYB expression in B-ALL cells could be a potential molecular marker for childhood leukemia.

Expression of ZNF695 Transcript Variants in Childhood B-Cell Acute Lymphoblastic Leukemia.

B-cell acute lymphoblastic leukemia is the most commonly diagnosed childhood malignancy worldwide; more than 50% of these cases are diagnosed in Mexico. Although the five-year survival rate is >80%, 30% of patients experience relapse with poor prognosis. Cancer-associated gene expression profiles have been identified in several malignancies, and some transcripts have been used to predict disease prognosis. The human transcriptome is incompletely elucidated; moreover, more than 80% of transcripts can be processed via alternative splicing (AS), which increases transcript and protein diversity. The human transcriptome is divided; coding RNA accounts for 2%, and the remaining 98% is noncoding RNA. Noncoding RNA can undergo AS, promoting the diversity of noncoding transcripts. We designed specific primers to amplify previously reported alternative transcript variants of ZNF695 and showed that six ZNF695 transcript variants are co-expressed in cancer cell lines. The amplicons were sequenced and identified. Additionally, we analyzed the expression of these six transcript variants in bone marrow from B-cell acute lymphoblastic leukemia patients and observed that ZNF695 transcript variants one and three were the predominant variants expressed in leukemia. Moreover, our results showed the co-expression of coding and long noncoding RNA. Finally, we observed that long noncoding RNA ZNF695 expression predicted survival rates.

WHOLE SEQUENCING OF THE MITOCHONDRIAL GENOME OF BREAST CANCER TISSUE IN MEXICAN-MESTIZO POSTMENOPAUSAL WOMEN WITH DIFFERENT BODY MASS INDEX.

BACKGROUND: Mitochondrial and oxidative stress has been related to obesity and breast cancer being this cancer more frequent and more aggressive in postmenopausal women with obesity. OBJECTIVE: The objective of this study was to investigate whether Mexican-Mestizo postmenopausal women with breast cancer and obesity present different somatic mutations in the mitochondrial DNA (mtDNA) when compared to women with normal body mass index (BMI). SUBJECTS AND METHODS: We included six Mexican-Mestizo postmenopausal women bearing breast cancer and who underwent mastectomy or breast-conserving surgery. BMI was determined in each case. Patients' genomic DNA was isolated from blood leukocytes and tumor tissue samples. Whole mtDNA sequence was determined by MitoChip v2.0 mitochondrial resequencing array, and data were analyzed using the GeneChip Sequence Analysis Software. Tumor mtDNA sequence was compared with matched leukocyte mtDNA sequence. RESULTS: Three women had a normal BMI and three presented obesity. Overall, we found 64 genetic variants: 53.1% were somatic mutations and 46.9% were polymorphisms; 44.1% were in the non-coding region and 55.9% were in genes that encode for mitochondrial proteins. Among the somatic mutations, 67.7% were in patients with normal BMI and 32.3% in patients with obesity. CONCLUSIONS: We did not find a higher frequency of mitochondrial somatic mutations in postmenopausal women with breast cancer and obesity compared to those with normal BMI. However, results could be due to the small number of women studied.

Whole sequencing of the mitochondrial genome of breast cancer tissue in Mexican-mestizo postmenopausal women with different body mass index

Abstract Background Mitochondrial and oxidative stress has been related to obesity and breast cancer being this cancer more frequent and more aggressive in postmenopausal women with obesity. Objective The objective of this study was to investigate whether Mexican-Mestizo postmenopausal women with breast cancer and obesity present different somatic mutations in the mitochondrial DNA (mtDNA) when compared to women with normal body mass index (BMI). Subjects and Methods We included six Mexican-Mestizo postmenopausal women bearing breast cancer and who underwent mastectomy or breast-conserving surgery. BMI was determined in each case. Patients’ genomic DNA was isolated from blood leukocytes and tumor tissue samples. Whole mtDNA sequence was determined by MitoChip v2.0 mitochondrial resequencing array, and data were analyzed using the GeneChip Sequence Analysis Software. Tumor mtDNA sequence was compared with matched leukocyte mtDNA sequence. Results Three women had a normal BMI and three presented obesity. Overall, we found 64 genetic variants: 53.1% were somatic mutations and 46.9% were polymorphisms; 44.1% were in the non-coding region and 55.9% were in genes that encode for mitochondrial proteins. Among the somatic mutations, 67.7% were in patients with normal BMI and 32.3% in patients with obesity. Conclusions We did not find a higher frequency of mitochondrial somatic mutations in postmenopausal women with breast cancer and obesity compared to those with normal BMI. However, results could be due to the small number of women studied.

Identification and Validation of Novel Reference Genes in Acute Lymphoblastic Leukemia for Droplet Digital PCR.

Droplet digital PCR is the most robust method for absolute nucleic acid quantification. However, RNA is a very versatile molecule and its abundance is tissue-dependent. RNA quantification is dependent on a reference control to estimate the abundance. Additionally, in cancer, many cellular processes are deregulated which consequently affects the gene expression profiles. In this work, we performed microarray data mining of different childhood cancers and healthy controls. We selected four genes that showed no gene expression variations (PSMB6, PGGT1B, UBQLN2 and UQCR2) and four classical reference genes (ACTB, GAPDH, RPL4 and RPS18). Gene expression was validated in 40 acute lymphoblastic leukemia samples by means of droplet digital PCR. We observed that PSMB6, PGGT1B, UBQLN2 and UQCR2 were expressed ~100 times less than ACTB, GAPDH, RPL4 and RPS18. However, we observed excellent correlations among the new reference genes (p < 0.0001). We propose that PSMB6, PGGT1B, UBQLN2 and UQCR2 are housekeeping genes with low expression in childhood cancer.

Centromere-associated protein E expresses a novel mRNA isoform in acute lymphoblastic leukemia.

The alternative splicing plays an important role to generate protein diversity. Recent studies have shown alterations in alternative splicing, resulting in loss, gain or changes of functions in the resulting protein. Specific products of alternative splicing are known to contribute in cancer-related mechanisms, such as angiogenesis, migration, adhesion and cell proliferation, among others. We using high-density microarrays reported a CENP-E as a one of significant transcript expressed and potentially is alternatively spliced in cancer. We focus in validate alternative splicing of CENP-E transcript using RT-PCR and sequencing in different cancer cell lines. We performed RT-PCR using specific primers designed to delimit the non-reported alternative splicing in CENP-E transcript. Our results showed the co-expression of the variant one and two of CENP-E in all cell lines evaluated. We detected more expression of variant one than two. Moreover, we identify an alternative 5'splice site of CENP-E in the exon 38 and was observed in RoVa cell line. Additionally, we characterized alternative skipping from exon 20 (NAT-CENP-E), these alternative splicing was observed in all cell lines evaluated except RoVa. Finally, we corroborate alternative mRNA splicing in leukemia patients using quantitative RT-PCR, in 71.8% of the patients NAT-CENP-E is downregulated and 28.2% is overexpressed.

Duplication of SOX9 associated with 46,XX ovotesticular disorder of sex development.

RESEARCH QUESTION: The purpose of the present study was to investigate whether ten unrelated SRY-negative individuals with this sex differentiation disorder presented a double dose of SOX9 as the cause of their disease. DESIGN: Ten unrelated SRY-negative 46,XX ovotesticular disorder of sexual development (DSD) subjects were molecularly studied. Multiplex-ligation dependent probe amplification (MLPA) and quantitative real-time PCR analysis (qRT-PCR) for SOX9 were performed. RESULTS: The MLPA analysis demonstrated that one patient presented a heterozygous duplication of the entire SOX9 coding region (above 1.3 value of peak ratio), as well as at least a ~ 483 kb upstream duplication. Moreover, no duplication of other SOX9 probes was observed corresponding to the region between -1007 and -1500 kb upstream. A qRT-PCR analysis showed a duplication of at least -581 kb upstream and ~1.63 kb of the coding region that encompasses exon 3. The limits of the duplication were mapped approximately from ~71539762 to 72122741 of Chr17. No molecular abnormalities were found in the remaining nine patients. CONCLUSION: This study is thought to be the first report regarding a duplication of SOX9 that is associated with the presence of 46,XX ovotesticular DSD, encompassing at least -581 kb upstream, and the almost entire coding region of the gene.

Genome-wide mRNA analysis reveals a TUBD1 isoform profile as a potential biomarker for diabetic retinopathy development.

Diabetic retinopathy (DR) affects approximately one third of all diabetic subjects and is the leading cause of blindness in young to middle-aged adults in the developed world. While early diagnosis is crucial for preventing DR-associated visual loss, the identification of accessible biomarkers that could lead to presymptomatic recognition of the disease is of great clinical importance. The aim of this work was to investigate the possible involvement of alternative splicing events in DR development by performing a genome-wide transcriptional profiling comparing blood-derived RNA from DR subjects and from diabetic-non DR controls. A total of 95 RNA samples, 67 from patients with bilateral DR and 28 from diabetic patients without DR after a period of at least 10 years with type 2 DM, were compared in a genome-wide transcriptome analysis using the GeneChip® Human Gene 2.0 ST Array which contains probe sets covering all exons of ∼33,500 coding transcripts of annotated genes. Microarray data analysis followed by RT-PCR and cDNA sequencing identified important differential splicing events in TUBD1 (Tubulin, Delta-1) isoforms between DR and DM samples. Specifically, the co-expression of particular TUBD1 isoforms was significantly associated with NPDR risk (p = 0.039 by Pearson's chi-squared test; OR (CI 95%): 8.1 (1.0-72.7)). Analysis of TUBD1 signal pathways and regulating networks using a MetaCore platform showed that HIF-1, a molecule playing an important role in the pathogenesis of DR, is a direct regulator of TUBD1 expression. In conjunction, our data suggest that TUBD1 mRNA isoform expression profile in peripheral blood could be an accessible biomarker for predicting the risk for diabetic retinopathy development.

Quality Control Usage in High-Density Microarrays Reveals Differential Gene Expression Profiles in Ovarian Cancer.

There are several existing reports of microarray chip use for assessment of altered gene expression in different diseases. In fact, there have been over 1.5 million assays of this kind performed over the last twenty years, which have influenced clinical and translational research studies. The most commonly used DNA microarray platforms are Affymetrix GeneChip and Quality Control Software along with their GeneChip Probe Arrays. These chips are created using several quality controls to confirm the success of each assay, but their actual impact on gene expression profiles had not been previously analyzed until the appearance of several bioinformatics tools for this purpose. We here performed a data mining analysis, in this case specifically focused on ovarian cancer, as well as healthy ovarian tissue and ovarian cell lines, in order to confirm quality control results and associated variation in gene expression profiles. The microarray data used in our research were downloaded from ArrayExpress and Gene Expression Omnibus (GEO) and analyzed with Expression Console Software using RMA, MAS5 and Plier algorithms. The gene expression profiles were obtained using Partek Genomics Suite v6.6 and data were visualized using principal component analysis, heat map, and Venn diagrams. Microarray quality control analysis showed that roughly 40% of the microarray files were false negative, demonstrating over- and under-estimation of expressed genes. Additionally, we confirmed the results performing second analysis using independent samples. About 70% of the significant expressed genes were correlated in both analyses. These results demonstrate the importance of appropriate microarray processing to obtain a reliable gene expression profile.

Data Mining for Identification of Molecular Targets in Ovarian Cancer.

Ovarian cancer is possibly the sixth most common malignancy worldwide, in Mexico representing the fourth leading cause of gynecological cancer death more than 70% being diagnosed at an advanced stage and the survival being very poor. Ovarian tumors are classified according to histological characteristics, epithelial ovarian cancer as the most common (~80%). We here used high-density microarrays and a systems biology approach to identify tissue-associated deregulated genes. Non-malignant ovarian tumors showed a gene expression profile associated with immune mediated inflammatory responses (28 genes), whereas malignant tumors had a gene expression profile related to cell cycle regulation (1,329 genes) and ovarian cell lines to cell cycling and metabolism (1,664 genes).

Exome sequencing identifies a de novo frameshift mutation in the imprinted gene ZDBF2 in a sporadic patient with Nasopalpebral Lipoma-coloboma syndrome.

Nasopalpebral lipoma-coloboma syndrome (NPLCS, OMIM%167730) is an uncommon malformation entity with autosomal dominant inheritance characterized by the combination of nasopalpebral lipoma, colobomas in upper and lower eyelids, telecanthus, and maxillary hypoplasia. To date, no genetic defects have been associated with familial or sporadic NPLCS cases and the etiology of the disease remains unknown. In this work, the results of whole exome sequencing in a sporadic NPLCS patient are presented. Exome sequencing identified a de novo heterozygous frameshift dinucleotide insertion c.6245_6246 insTT (p.His2082fs*67) in ZDBF2 (zinc finger, DBF-type containing 2), a gene located at 2q33.3. This variant was absent in parental DNA, in a set of 300 ethnically matched controls, and in public exome variant databases. This is the first genetic variant identified in a NPLCS patient and evidence supporting the pathogenicity of the identified mutation is discussed. © 2016 Wiley Periodicals, Inc.

Characterization and mRNA expression analysis of a novel ARG1 splicing mutation causing hyperargininemia.

OBJECTIVES: Biallelic mutations in the ARG1 gene result in an uncommon autosomal recessive inborn defect of the urea cycle known as hyperargininemia (OMIM #207800). ARG1 splicing mutations are not reported often, and they are probably related to a more severe phenotype than missense mutations. In this article, we describe the results of molecular studies in a young hyperargininemia patient carrying a novel splicing mutation in ARG1. DESIGN AND METHODS: Molecular analyses included PCR amplification and direct nucleotide sequencing of the ARG1 gene. RT-PCR analysis was performed to investigate the effect of the mutation in mRNA splicing and in the expression of ARG1 isoforms. RESULTS: Mutational analysis identified a novel homozygous ARG1 IVS4-1G>C point mutation in the patient's DNA. Blood leukocyte mRNA was analyzed to demonstrate the splicing defect caused by this mutation. Sequencing of ARG1 RT-PCR products allowed the characterization of a mutated transcript retaining 51-bp from intron 4. In addition, two new, alternatively spliced ARG1 transcripts lacking either exon 4 or exons 4 and 5 were identified in mRNA from the patient and from controls. CONCLUSIONS: Our results expand the mutational spectrum in hyperargininemia patients and indicate that the novel splicing mutation results in an aberrant transcript retaining intronic sequences. Two novel alternatively spliced ARG1 transcripts were also recognized.

Acro-spondylo-pubic dysostosis associated with cataracts, microcephaly, and normal intelligence.

We report on an adult male with normal intelligence who exhibited an unusual combination of microcephaly, dysostoses of limbs, vertebrae, patellae, and pubic bone, camptodactyly of all fingers, and syndactyly of toes, absent nails on thumbs and some fingers, bilateral cataract, cryptorchidism, polythelia, and nipple-like skin pigmentations of shoulders and upper back. We have been unable to find a description of a similar combination of manifestations in literature. The cause of the anomalies remains unknown.

A case of cervical cancer expressed three mRNA variant of Hyaluronan-mediated motility receptor.

Cervical cancer is the second malignancy in Mexico, little is known about the prognostic factors associated with this disease. Several cellular components are important in their transformation and progression. Alternative mRNA splice is an important mechanism for generating protein diversity, nevertheless, in cancer unknown mRNA diversity is expressed. Hyaluronan-mediated motility receptor (HMMR, RHAMM, CD168) is a family member of proteins, hyaluronan acid dependent, and has been associated with different malignant processes such as: angiogenesis, cell invasiveness, proliferation, metastasis and poor outcome in some tumors. In the present study we identified expression of HMMR in cervical cancer by means of RT-PCR and sequencing. Our results indicate co-expression of two HMMR variants in all samples, and one case expressed three alternative HMMR splice transcripts. These results showed the heterogeneity of mRNA transcripts of HMMR that could express in cancer and the expression of HMMR could be marker of malignancy in CC.

Heterogeneity of microRNAs expression in cervical cancer cells: over-expression of miR-196a.

In recent years, the study of microRNAs associated with neoplastic processes has increased. Patterns of microRNA expression in different cell lines and different kinds of tumors have been identified; however, little is known about the alterations in regulatory pathways and genes involved in aberrant set of microRNAs. The identification of these altered microRNAs in several cervical cancer cells and potentially deregulated pathways involved constitute the principal goals of the present study. In the present work, the expression profiles of cellular microRNAs in Cervical Cancer tissues and cell lines were explored using microRNA microarray, Affymetrix. The most over-expressed was miR-196a, which was evaluated by real time PCR, and HOXC8 protein as potential target by immunohistochemistry assay. One hundred and twenty three human microRNAs differentially expressed in the cell tumor, 64 (52%) over-expressed and 59 (48%) under-expressed were observed. Among the microRNAs over-expressed, we focused on miR-196a; at present this microRNA is poorly studied in CC. The expression of this microRNA was evaluated by qRT-PCR, and HOXC8 by immunohistochemistry assay. There is not a specific microRNA expression profile in the CC cells, neither a microRNA related to HPV presence. Furthermore, the miR-196a was over-expressed, while an absence of HOXC8 expression was observed. We suggest that miR-196a could be played as oncomiR in CC.

Identification and expression analysis of a novel intragenic EFNB1 mutation causing craniofrontonasal syndrome.

Craniofrontonasal syndrome (CFNS) is an X-linked disorder caused by mutations in the EFNB1 gene and characterized by distinctive craniofacial and digital malformations. In contrast with most X-linked traits, female patients with CFNS display a more severe phenotype than males. In this report, the clinical, molecular and RNA expression analyses of a female subject with CFNS are described. A novel c.445_449delGAGGG deletion in exon 3 of EFNB1 was demonstrated in this patient. To assess the effect of this novel mutation at the transcript level, the expression of EFNB1 mRNA was studied by quantitative RT-PCR. To our knowledge, this is the first time that an EFNB1 transcript carrying a truncating mutation in exon 3 is demonstrated to undergo degradation by nonsense-mediated mRNA decay. Our results expand the mutational spectrum of CFNS and add to the functional consequences of truncating EFNB1 mutations.

Splice variants of zinc finger protein 695 mRNA associated to ovarian cancer.

BACKGROUND: Studies of alternative mRNA splicing (AS) in health and disease have yet to yield the complete picture of protein diversity and its role in physiology and pathology. Some forms of cancer appear to be associated to certain alternative mRNA splice variants, but their role in the cancer development and outcome is unclear. METHODS: We examined AS profiles by means of whole genome exon expression microarrays (Affymetrix GeneChip 1.0) in ovarian tumors and ovarian cancer-derived cell lines, compared to healthy ovarian tissue. Alternatively spliced genes expressed predominantly in ovarian tumors and cell lines were confirmed by RT-PCR. RESULTS: Among several significantly overexpressed AS genes in malignant ovarian tumors and ovarian cancer cell lines, the most significant one was that of the zinc finger protein ZNF695, with two previously unknown mRNA splice variants identified in ovarian tumors and cell lines. The identity of ZNF695 AS variants was confirmed by cloning and sequencing of the amplicons obtained from ovarian cancer tissue and cell lines. CONCLUSIONS: Alternative ZNF695 mRNA splicing could be a marker of ovarian cancer with possible implications on its pathogenesis.

Prevalence of human papillomavirus in the cervical epithelium of Mexican women: meta-analysis.

BACKGROUND: Human Papillomavirus (HPV) in cervical epithelium has been identified as the main etiological factor in the developing of Cervical Cancer (CC), which has recently become a public health problem in Mexico. This finding has allowed for the development of vaccines that help prevent this infection. In the present study, we aimed to determine the prevalence and HPV type-distribution in Mexican women with CC, high-grade squamous intraepithelial lesion (HSIL), low-grade squamous intraepithelial lesion (LSIL), and Normal cytology (N) to estimate the impact of the HPV vaccines. METHODS: The PubMed database was used to identify and review all articles that reported data on HPV prevalence in CC, precursor lesions, and normal cytology of Mexican women. RESULTS: A total of 8,706 samples of the tissues of Mexican women were stratified according to diagnosis as follows: 499 for CC; 364 for HSIL; 1,425 for LSIL, and 6,418 for N. According to the results, the most prevalent genotypes are the following: HPV16 (63.1%), -18 (8.6%), -58, and -31 (5%) for CC; HPV-16 (28.3%), 58 (12.6%), 18 (7.4%), and 33 (6.5%) for HSIL; HPV-16 (13.1%), 33 (7.4%), 18 (4.2%), and 58 (2.6%) for LSIL, and HPV-16 (3.4%), 33 (2.1%), 18, and 58 (1.2%) for N. CONCLUSIONS: Taken together, genotypes 58 and 31 (10%) are more common than type 18 (8.6%) in CC. Therefore, the inclusion of these two genotypes in a second-generation vaccine would provide optimal prevention of CC in Mexico.