Exome Sequencing to Identify Novel Variants Associated with Secondary Amenorrhea and Premature Ovarian Insufficiency (POI) in Saudi Women.

BACKGROUND AND OBJECTIVES: Post-pubertal disappearance of menstrual cycles (secondary amenorrhea) associated with premature follicular depletion is a heterogeneous condition. Patients with this disease have low levels of gonadal hormones and high levels of gonadotropins. It is one of the causes of female infertility and a strong genetic component is attributed as an underlying cause of this condition. Although variants in several genes have been associated with the condition, the cause of the disease remains undetermined in the vast majority of cases. Methodology and Materials: Ten Saudi married women experiencing secondary amenorrhea were referred to a center for genetics and inherited diseases for molecular investigation. A family-based study design was used. Intensive clinical examinations, including pelvic ultra-sonography (U/S) and biochemical evaluations, were carried out. Karyotypes were normal in all cases and polycystic ovarian syndrome (PCOS) was excluded by using Rotterdam consensus criteria. Patients' DNA samples were whole-exome sequenced (WES). Bidirectional Sanger sequencing was then utilized to validate the identified candidate variants. The pathogenicity of detected variants was predicted using several types of bioinformatics software. RESULTS: Most of the patients have a normal uterus with poor ovarian reserves. Exome sequence data analysis identified candidate variants in genes associated with POI in 60% of cases. Novel variants were identified in HS6ST1, MEIOB, GDF9, and BNC1 in POI-associated genes. Moreover, a homozygous variant was also identified in the MMRN1 gene. Interestingly, mutations in MMRN1 have never been associated with any human disease. The variants identified in this study were not present in 125 healthy Saudi individuals. CONCLUSIONS: WES is a powerful tool to identify the underlying variants in genetically heterogeneous diseases like secondary amenorrhea and POI. In this study, we identified six novel variants and expanded the genotype continuum of POI. Unravelling the genetic landscape of POI will help in genetic counselling, management, and early intervention.

Association of Homozygous PROP1 Mutation in a Saudi Family with Combined Pituitary Hormone Deficiency.

Background and Objectives: Combined pituitary hormone deficiency (CPHD) is a rare heterogeneous disease. It is characterized by the deficiency of growth hormone (GH) and shortage of at least one or more other hormones of the pituitary gland including thyroid-stimulating hormone (TSH), luteinizing hormone (LH), follicle-stimulating hormone (FSH), and prolactin. Rare pathogenic variants in nearly 30 genes have been identified as an underlying cause of CPHD pathogenicity. Among these genes, paired-like homeobox 1 (PROP1) has been reported to be the most common cause of CPHD. Materials and Methods: In the present study, we investigated a large family of Saudi origin with three adult sisters suffering from short stature in combination of secondary amenorrhea. Results: Whole-exome sequencing followed by Sanger sequencing shows a homozygous missense variant (NM_006261.5; c.211C > T; p.R71C) in the PROP1 gene segregating with the disease phenotype within the family. In silico analysis studies show that this variant is highly conserved among several orthologues and is predicted as likely pathogenic using various bioinformatics tools. Conclusions: Our finding presents the first Saudi familial case of autosomal recessive form of CPHD caused by the PROP1 variant.

Investigating the pathogenic SNPs in BLM helicase and their biological consequences by computational approach.

The BLM helicase protein plays a vital role in DNA replication and the maintenance of genomic integrity. Variation in the BLM helicase gene resulted in defects in the DNA repair mechanism and was reported to be associated with Bloom syndrome (BS) and cancer. Despite extensive investigation of helicase proteins in humans, no attempt has previously been made to comprehensively analyse the single nucleotide polymorphism (SNPs) of the BLM gene. In this study, a comprehensive analysis of SNPs on the BLM gene was performed to identify, characterize and validate the pathogenic SNPs using computational approaches. We obtained SNP data from the dbSNP database version 150 and mapped these data to the genomic coordinates of the "NM_000057.3" transcript expressing BLM helicase (P54132). There were 607 SNPs mapped to missense, 29 SNPs mapped to nonsense, and 19 SNPs mapped to 3'-UTR regions. Initially, we used many consensus tools of SIFT, PROVEAN, Condel, and PolyPhen-2, which together increased the accuracy of prediction and identified 18 highly pathogenic non-synonymous SNPs (nsSNPs) out of 607 SNPs. Subsequently, these 18 high-confidence pathogenic nsSNPs were analysed for BLM protein stability, structure-function relationships and disease associations using various bioinformatics tools. These 18 mutants of the BLM protein along with the native protein were further investigated using molecular dynamics simulations to examine the structural consequences of the mutations, which might reveal their malfunction and contribution to disease. In addition, 28 SNPs were predicted as "stop gained" nonsense SNPs and one SNP was predicted as "start lost". Two SNPs in the 3'UTR were found to abolish miRNA binding and thus may enhance the expression of BLM. Interestingly, we found that BLM mRNA overexpression is associated with different types of cancers. Further investigation showed that the dysregulation of BLM is associated with poor overall survival (OS) for lung and gastric cancer patients and hence led to the conclusion that BLM has the potential to be used as an important prognostic marker for the detection of lung and gastric cancer.

Homozygosity mapping and whole exome sequencing provide exact diagnosis of Cohen syndrome in a Saudi family.

BACKGROUND: Cohen syndrome (CS) is a rare multi-system autosomal recessive disorder with a high prevalence in the Finnish population. Clinical features of Finnish-type CS are homogeneous, however, in non-Finnish populations, CS diagnosis is challenging due to broad phenotypic variability. METHODS: We studied a consanguineous family having three affected individuals with clinical features of severe intellectual disability and global developmental delay. Clinical diagnosis of the phenotype could not be established based on the features. Therefore, whole genome SNP genotyping and whole exome sequencing (WES) were performed on DNA samples from affected and unaffected family members. RESULTS: Homozygosity mapping identified a shared loss of heterozygosity region on chromosome 8q22.1-q22.3 and WES data analysis revealed an insertion-deletion (indel) mutation (c.11519_11521delCAAinsT) in the VPS13B gene. The indel is predicted to cause a frameshift resulting in a premature termination of the VPS13B protein (NP_060360.3:p.Pro3840Leufs*2). CONCLUSION: VPS13B encodes a giant transmembrane protein called vacuolar protein sorting 13 homolog B. VPS13B is known to play a role in the glycosylation of Golgi proteins and in endosomal-lysosomal trafficking. Moreover, it is thought to function in vesicle mediated transport and sorting of proteins within the cell. The mechanism by which abnormalities of the VPS13B protein lead to the phenotype of CS is currently unknown. Here, in this study, we successfully established a clinical diagnosis of CS cases from a family using genomic analyses. Clinical re-examination of the patients revealed characteristic ocular abnormalities.

KMT2C, a histone methyltransferase, is mutated in a family segregating non-syndromic primary failure of tooth eruption.

Primary failure of tooth eruption (PFE) is a rare odontogenic defect and is characterized by failure of eruption of one or more permanent teeth. The aim of the study is to identify the genetic defect in a family with seven affected individuals segregating autosomal dominant non-syndromic PFE. Whole genome single-nucleotide polymorphism (SNP) genotyping was performed. SNP genotypes were analysed by DominantMapper and multiple shared haplotypes were detected on different chromosomes. Four individuals, including three affected, were exome sequenced. Variants were annotated and data were analysed while considering candidate chromosomal regions. Initial analysis of variants obtained by whole exome sequencing identified damaging variants in C15orf40, EPB41L4A, TMEM232, KMT2C, and FBXW10 genes. Sanger sequencing of all family members confirmed segregation of splice acceptor site variant (c.1013-2 A > G) in the KMT2C gene with the phenotype. KMT2C is considered as a potential candidate gene based on segregation analysis, the absence of variant in the variation databases, the presence of variant in the shared identical by descent (IBD) region and in silico pathogenicity prediction. KMT2C is a histone methyltransferase and recently the role of another member of this family (KMT2D) has been implicated in tooth development. Moreover, protein structures of KMT2C and KMT2D are highly similar. In conclusion, we have identified that the KMT2C gene mutation causes familial non-syndromic PFE. These findings suggest the involvement of KMT2C in the physiological eruption of permanent teeth.

Whole exome sequencing identified a novel single base pair insertion mutation in the EYS gene in a six generation family with retinitis pigmentosa.

Retinitis pigmentosa (RP) is a group of inherited progressive retinal dystrophies (RD) and is characterized by photoreceptor degeneration. RP is clinically and genetically heterogeneous disorder. More than 70 genes are known and, thus, identification of causative genes and mutations in known genes is challenging. This study was designed to identify the underlying genetic defect in a large extended Saudi family with multiple RP affected members. Fundus photography, Optical Coherence Tomography (OCT) and visual field perimetry were performed for affected individuals. Whole exome sequencing was used to detect the underlying genetic defect in a large family with 12 affected individuals showing autosomal recessive isolated RP. WES data analysis identified a novel insertion mutation in the EYS (eyes shut homolog) gene (c.910_911insT; p.Trp304LeufsTer8). Sanger sequencing validates the variant discovered through exome in all 12 affected individuals and showed that this mutation is segregating with RP phenotype in an autosomal recessive manner in 51 individuals of the family tested here. Our study expands the mutation spectrum of EYS gene in RP patients and extends the body of evidence that supports the importance of EYS gene in eye development.

Identification of a class of human cancer germline genes with transcriptional silencing refractory to the hypomethylating drug 5-aza-2'-deoxycytidine.

Bona fide germline genes have expression restricted to the germ cells of the gonads. Testis-specific germline development-associated genes can become activated in cancer cells and can potentially drive the oncogenic process and serve as therapeutic/biomarker targets; such germline genes are referred to as cancer/testis genes. Many cancer/testis genes are silenced via hypermethylation of CpG islands in their associated transcriptional control regions and become activated upon treatment with DNA hypomethylating agents; such hypomethylation-induced activation of cancer/testis genes provides a potential combination approach to augment immunotherapeutics. Thus, understanding cancer/testis gene regulation is of increasing clinical importance. Previously studied cancer/testis gene activation has focused on X chromosome encoded cancer/testis genes. Here we find that a sub-set of non-X encoded cancer/testis genes are silenced in non-germline cells via a mechanism that is refractory to epigenetic dysregulation, including treatment with the hypomethylating agent 5-aza-2'-deoxycytidine and the histone deacetylase inhibitor tricostatin A. These findings formally indicate that there is a sub-group of the clinically important cancer/testis genes that are unlikely to be activated in clinical therapeutic approaches using hypomethylating agents and it indicates a unique transcriptional silencing mechanism for germline genes in non-germline cells that might provide a target mechanism for new clinical therapies.

Meta-analysis of clinical data using human meiotic genes identifies a novel cohort of highly restricted cancer-specific marker genes.

Identifying cancer-specific biomarkers represents an ongoing challenge to the development of novel cancer diagnostic, prognostic and therapeutic strategies. Cancer/testis (CT) genes are an important gene family with expression tightly restricted to the testis in normal individuals but which can also be activated in cancers. Here we develop a pipeline to identify new CT genes. We analysed and validated expression profiles of human meiotic genes in normal and cancerous tissue followed by meta-analyses of clinical data sets from a range of tumour types resulting in the identification of a large cohort of highly specific cancer biomarker genes, including the recombination hot spot activator PRDM9 and the meiotic cohesin genes SMC1beta and RAD21L. These genes not only provide excellent cancer biomarkers for diagnostics and prognostics, but may serve as oncogenes and have excellent drug targeting potential.