Genetic predisposition to porto-sinusoidal vascular disorder: A functional genomic-based, multigenerational family study.

BACKGROUND AND AIMS: Porto-sinusoidal vascular disorder (PSVD) is a group of liver vascular diseases featuring lesions encompassing the portal venules and sinusoids unaccompanied by cirrhosis, irrespective of the presence/absence of portal hypertension. It can occur secondary to coagulation disorders or insult by toxic agents. However, the cause of PSVD remains unknown in most cases. Hereditary cases of PSVD are exceptionally rare, but they are of particular interest and may unveil genetic alterations and molecular mechanisms associated with the disease. APPROACH AND RESULTS: We performed genome sequencing of four patients and two healthy individuals of a large multigenerational Lebanese family with PSVD and identified a heterozygous deleterious variant (c.547C>T, p.R183W) of FCH and double SH3 domains 1 ( FCHSD1 ), an uncharacterized gene, in patients. This variant segregated with the disease, and its pattern of inheritance was suggestive of autosomal dominant with variable expressivity. RNA structural modelling of human FCHSD1 suggests that the C-to-T substitution at position 547, corresponding to FCHSD1R183W , may increase both messenger RNA (mRNA) and protein stability and its interaction with MTOR-associated protein, LST8 homolog, a key protein of the mechanistic target of rapamycin (mTOR pathway). These predictions were substantiated by biochemical analyses, which showed that FCHSD1R183W induced high FCHSD1 mRNA stability, overexpression of FCHSD1 protein, and an increase in mTORC1 activation. This human FCHSD1 variant was introduced into mice through CRISPR/Cas9 genome editing. Nine out of the 15 mice carrying the human FCHSD1R183W variant mimicked the phenotype of human PSVD, including splenomegaly and enlarged portal vein. CONCLUSIONS: Aberrant FCHSD1 structure and function leads to mTOR pathway overactivation and may cause PSVD.

Transcriptome Profile Identifies Actin as an Essential Regulator of Cardiac Myosin Binding Protein C3 Hypertrophic Cardiomyopathy in a Zebrafish Model.

Variants in cardiac myosin-binding protein C (cMyBP-C) are the leading cause of inherited hypertrophic cardiomyopathy (HCM), demonstrating the key role that cMyBP-C plays in the heart's contractile machinery. To investigate the c-MYBPC3 HCM-related cardiac impairment, we generated a zebrafish mypbc3-knockout model. These knockout zebrafish displayed significant morphological heart alterations related to a significant decrease in ventricular and atrial diameters at systolic and diastolic states at the larval stages. Immunofluorescence staining revealed significant hyperplasia in the mutant's total cardiac and ventricular cardiomyocytes. Although cardiac contractility was similar to the wild-type control, the ejection fraction was significantly increased in the mypbc3 mutants. At later stages of larval development, the mutants demonstrated an early cardiac phenotype of myocardium remodeling, concurrent cardiomyocyte hyperplasia, and increased ejection fraction as critical processes in HCM initiation to counteract the increased ventricular myocardial wall stress. The examination of zebrafish adults showed a thickened ventricular cardiac wall with reduced heart rate, swimming speed, and endurance ability in both the mypbc3 heterozygous and homozygous groups. Furthermore, heart transcriptome profiling showed a significant downregulation of the actin-filament-based process, indicating an impaired actin cytoskeleton organization as the main dysregulating factor associated with the early ventricular cardiac hypertrophy in the zebrafish mypbc3 HCM model.

Genetic predisposition to cancer across people of different ancestries in Qatar: a population-based, cohort study.

BACKGROUND: Disparities in the genetic risk of cancer among various ancestry groups and populations remain poorly defined. This challenge is even more acute for Middle Eastern populations, where the paucity of genomic data could affect the clinical potential of cancer genetic risk profiling. We used data from the phase 1 cohort of the Qatar Genome Programme to investigate genetic variation in cancer-susceptibility genes in the Qatari population. METHODS: The Qatar Genome Programme generated high-coverage genome sequencing on DNA samples collected from 6142 native Qataris, stratified into six distinct ancestry groups: general Arab, Persian, Arabian Peninsula, Admixture Arab, African, and South Asian. In this population-based, cohort study, we evaluated the performance of polygenic risk scores for the most common cancers in Qatar (breast, prostate, and colorectal cancers). Polygenic risk scores were trained in The Cancer Genome Atlas (TCGA) dataset, and their distributions were subsequently applied to the six different genetic ancestry groups of the Qatari population. Rare deleterious variants within 1218 cancer susceptibility genes were analysed, and their clinical pathogenicity was assessed by ClinVar and the CharGer computational tools. FINDINGS: The cohort included in this study was recruited by the Qatar Biobank between Dec 11, 2012, and June 9, 2016. The initial dataset comprised 6218 cohort participants, and whole genome sequencing quality control filtering led to a final dataset of 6142 samples. Polygenic risk score analyses of the most common cancers in Qatar showed significant differences between the six ancestry groups (p<0·0001). Qataris with Arabian Peninsula ancestry showed the lowest polygenic risk score mean for colorectal cancer (-0·41), and those of African ancestry showed the highest average for prostate cancer (0·85). Cancer-gene rare variant analysis identified 76 Qataris (1·2% of 6142 individuals in the Qatar Genome Programme cohort) carrying ClinVar pathogenic or likely pathogenic variants in clinically actionable cancer genes. Variant analysis using CharGer identified 195 individuals carriers (3·17% of the cohort). Breast cancer pathogenic variants were over-represented in Qataris of Persian origin (22 [56·4%] of 39 BRCA1/BRCA2 variant carriers) and completely absent in those of Arabian Peninsula origin. INTERPRETATION: We observed a high degree of heterogeneity for cancer predisposition genes and polygenic risk scores across ancestries in this population from Qatar. Stratification systems could be considered for the implementation of national cancer preventive medicine programmes. FUNDING: Qatar Foundation.

The QChip1 knowledgebase and microarray for precision medicine in Qatar.

Risk genes for Mendelian (single-gene) disorders (SGDs) are consistent across populations, but pathogenic risk variants that cause SGDs are typically population-private. The goal was to develop "QChip1," an inexpensive genotyping microarray to comprehensively screen newborns, couples, and patients for SGD risk variants in Qatar, a small nation on the Arabian Peninsula with a high degree of consanguinity. Over 108 variants in 8445 Qatari were identified for inclusion in a genotyping array containing 165,695 probes for 83,542 known and potentially pathogenic variants in 3438 SGDs. QChip1 had a concordance with whole-genome sequencing of 99.1%. Testing of QChip1 with 2707 Qatari genomes identified 32,674 risk variants, an average of 134 pathogenic alleles per Qatari genome. The most common pathogenic variants were those causing homocystinuria (1.12% risk allele frequency), and Stargardt disease (2.07%). The majority (85%) of Qatari SGD pathogenic variants were not present in Western populations such as European American, South Asian American, and African American in New York City and European and Afro-Caribbean in Puerto Rico; and only 50% were observed in a broad collection of data across the Greater Middle East including Kuwait, Iran, and United Arab Emirates. This study demonstrates the feasibility of developing accurate screening tools to identify SGD risk variants in understudied populations, and the need for ancestry-specific SGD screening tools.

Genome sequencing unveils mutational landscape of the familial Mediterranean fever: Potential implications of IL33/ST2 signalling.

Familial Mediterranean fever (FMF) is the most common auto-inflammatory disease. It is transmitted as autosomal recessive trait with mutations in MEditerranean FeVer (MEFV) gene. Despite a typical clinical expression, many patients have either a single or no mutation in MEFV. The current work is aimed to revisit the genetic landscape of FMF disease using high-coverage whole genome sequencing. In atypical patients (carrying a single or no mutation in MEFV), we revealed many rare variants in genes associated with auto-inflammatory disorders, and more interestingly, we discovered a novel variant ( a 2.1-Kb deletion) in exon 11 of IL1RL1 gene, present only in patients. To validate and screen this patient-specific variant, a tandem of allele-specific PCR and quantitative real-time PCR was performed in 184 FMF patients and 218 healthy controls and we demonstrated that the novel deletion was absent in controls and was present in more than 19% of patients. This study sheds more light on the mutational landscape of FMF. Our discovery of a disease-specific variant in IL1RL1 gene may constitute a novel genetic marker for FMF. This finding suggesting a potential role of the IL33/ST2 signalling in the disease pathogenicity highlights a new paradigm in FMF pathophysiology.

Differential responsiveness to BRAF inhibitors of melanoma cell lines BRAF V600E-mutated.

BACKGROUND: Most mutations in melanoma affect one critical amino acid on BRAF gene, resulting in the V600E substitution. Patient management is often based on the use of specific inhibitors targeting this mutation. METHODS: DNA and RNA mutation status was assessed in 15 melanoma cell lines by Sanger sequencing and RNA-seq. We tested the cell lines responsiveness to BRAF inhibitors (vemurafenib and PLX4720, BRAF-specific and sorafenib, BRAF non-specific). Cell proliferation was assessed by MTT colorimetric assay. BRAF V600E RNA expression was assessed by qPCR. Expression level of phosphorylated-ERK protein was assessed by Western Blotting as marker of BRAF activation. RESULTS: Three cell lines were discordant in the mutation detection (BRAF V600E at DNA level/Sanger sequencing and BRAF WT on RNA-seq). We initially postulated that those cell lines may express only the WT allele at the RNA level although mutated at the DNA level. A more careful analysis showed that they express low level of BRAF RNA and the expression may be in favor of the WT allele. We tested whether the discordant cell lines responded differently to BRAF-specific inhibitors. Their proliferation rate decreased after treatment with vemurafenib and PLX4720 but was not affected by sorafenib, suggesting a BRAF V600E biological behavior. Yet, responsiveness to the BRAF specific inhibitors was lower as compared to the control. Western Blot analysis revealed a decreased expression of p-ERK protein in the BRAF V600E control cell line and in the discordant cell lines upon treatment with BRAF-specific inhibitors. The discordant cell lines showed a lower responsiveness to BRAF inhibitors when compared to the BRAF V600E control cell line. The results obtained from the inhibition experiment and molecular analyses were also confirmed in three additional cell lines. CONCLUSION: Cell lines carrying V600E mutation at the DNA level may respond differently to BRAF targeted treatment potentially due to a lower V600E RNA expression.

Keratinocytes Derived from Patient-Specific Induced Pluripotent Stem Cells Recapitulate the Genetic Signature of Psoriasis Disease.

Psoriasis is characterized by hyperproliferation and defective differentiation of keratinocytes (KCs). Patients with psoriasis are at a high risk of developing diabetes and cardiovascular diseases. The debate on the genetic origin of psoriasis pathogenesis remains unresolved due to lack of suitable in vitro human models mimicking the disease phenotypes. In this study, we provide the first human induced pluripotent stem cell (iPSC) model for psoriasis carrying the genetic signature of the patients. iPSCs were generated from patients with psoriasis (PsO-iPSCs) and healthy donors (Ctr-iPSCs) and were efficiently differentiated into mature KCs. RNA sequencing of KCs derived from Ctr-iPSCs and PsO-iPSCs identified 361 commonly upregulated and 412 commonly downregulated genes. KCs derived from PsO-iPSCs showed dysregulated transcripts associated with psoriasis and KC differentiation, such as HLA-C, KLF4, chemokines, type I interferon-inducible genes, solute carrier family, IVL, DSG1, and HLA-DQA1, as well as transcripts associated with insulin resistance, such as IRS2, GDF15, GLUT10, and GLUT14. Our data suggest that the KC abnormalities are the main driver triggering psoriasis pathology and highlights the substantial contribution of genetic predisposition in the development of psoriasis and insulin resistance.

Monoallelic expression in melanoma.

BACKGROUND: Monoallelic expression (MAE) is a frequent genomic phenomenon in normal tissues, however its role in cancer is yet to be fully understood. MAE is defined as the expression of a gene that is restricted to one allele in the presence of a diploid heterozygous genome. Constitutive MAE occurs for imprinted genes, odorant receptors and random X inactivation. Several studies in normal tissues have showed MAE in approximately 5-20% of the cases. However, little information exists on the MAE rate in cancer. In this study we assessed the presence and rate of MAE in melanoma. The genetic basis of melanoma has been studied in depth over the past decades, leading to the identification of mutations/genetic alterations responsible for melanoma development. METHODS: To examine the role of MAE in melanoma we used 15 melanoma cell lines and compared their RNA-seq data with genotyping data obtained by the parental TIL (tumor infiltrating lymphocytes). Genotyping was performed using the Illumina HumanOmni1 beadchip. The RNA-seq library preparation and sequencing was performed using the Illumina TruSeq Stranded Total RNA Human Kit and subsequently sequenced using a HiSeq 2500 according to manufacturer's guidelines. By comparing genotyping data with RNA-seq data, we identified SNPs in which DNA genotypes were heterozygous and corresponding RNA genotypes were homozygous. All homozygous DNA genotypes were removed prior to the analysis. To confirm the validity to detect MAE, we examined heterozygous DNA genotypes from X chromosome of female samples as well as for imprinted and olfactory receptor genes and confirmed MAE. RESULTS: MAE was detected in all 15 cell lines although to a different rate. When looking at the B-allele frequencies we found a preferential pattern of complete monoallelic expression rather then differential monoallelic expression across the 15 melanoma cell lines. As some samples showed high differences in the homozygous and heterozygous call rate, we looked at the single chromosomes and showed that MAE may be explained by underlying large copy number imbalances in some instances. Interestingly these regions included genes known to play a role in melanoma initiation and progression. Nevertheless, some chromosome regions showed MAE without CN imbalances suggesting that additional mechanisms (including epigenetic silencing) may explain MAE in melanoma. CONCLUSION: The biological implications of MAE are yet to be realized. Nevertheless, our findings suggest that MAE is a common phenomenon in melanoma cell lines. Further analyses are currently being undertaken to evaluate whether MAE is gene/pathway specific and to understand whether MAE can be employed by cancers to achieve a more aggressive phenotype.

Osteolytic cancer cells induce vascular/axon guidance processes in the bone/bone marrow stroma.

Prostate and breast cancers frequently metastasize to bone. The physiological bone homeostasis is perturbed once cancer cells proliferate at the bone metastatic site. Tumors are complex structures consisting of cancer cells and numerous stroma cells. In this study, we show that osteolytic cancer cells (PC-3 and MDA-MB231) induce transcriptome changes in the bone/bone marrow microenvironment (stroma). This stroma transcriptome differs from the previously reported stroma transcriptome of osteoinductive cancer cells (VCaP). While the biological process "angiogenesis/vasculogenesis" is enriched in both transcriptomes, the "vascular/axon guidance" process is a unique process that characterizes the osteolytic stroma. In osteolytic bone metastasis, angiogenesis is denoted by vessel morphology and marker expression specific for arteries/arterioles. Interestingly, intra-tumoral neurite-like structures were in proximity to arteries. Additionally, we found that increased numbers of mesenchymal stem cells and vascular smooth muscle cells, expressing osteolytic cytokines and inhibitors of bone formation, contribute to the osteolytic bone phenotype. Osteoinductive and osteolytic cancer cells induce different types of vessels, representing functionally different hematopoietic stem cell niches. This finding suggests different growth requirements of osteolytic and osteoinductive cancer cells and the need for a differential anti-angiogenic strategy to inhibit tumor growth in osteolytic and osteoblastic bone metastasis.

Accelerating next generation sequencing data analysis with system level optimizations.

Next generation sequencing (NGS) data analysis is highly compute intensive. In-memory computing, vectorization, bulk data transfer, CPU frequency scaling are some of the hardware features in the modern computing architectures. To get the best execution time and utilize these hardware features, it is necessary to tune the system level parameters before running the application. We studied the GATK-HaplotypeCaller which is part of common NGS workflows, that consume more than 43% of the total execution time. Multiple GATK 3.x versions were benchmarked and the execution time of HaplotypeCaller was optimized by various system level parameters which included: (i) tuning the parallel garbage collection and kernel shared memory to simulate in-memory computing, (ii) architecture-specific tuning in the PairHMM library for vectorization, (iii) including Java 1.8 features through GATK source code compilation and building a runtime environment for parallel sorting and bulk data transfer (iv) the default 'on-demand' mode of CPU frequency is over-clocked by using 'performance-mode' to accelerate the Java multi-threads. As a result, the HaplotypeCaller execution time was reduced by 82.66% in GATK 3.3 and 42.61% in GATK 3.7. Overall, the execution time of NGS pipeline was reduced to 70.60% and 34.14% for GATK 3.3 and GATK 3.7 respectively.

Diagnosis implications of the whole genome sequencing in a large Lebanese family with hyaline fibromatosis syndrome.

BACKGROUND: Hyaline fibromatosis syndrome (HFS) is a recently introduced alternative term for two disorders that were previously known as juvenile hyaline fibromatosis (JHF) and infantile systemic hyalinosis (ISH). These two variants are secondary to mutations in the anthrax toxin receptor 2 gene (ANTXR2) located on chromosome 4q21. The main clinical features of both entities include papular and/or nodular skin lesions, gingival hyperplasia, joint contractures and osteolytic bone lesions that appear in the first few years of life, and the syndrome typically progresses with the appearance of new lesions. METHODS: We describe five Lebanese patients from one family, aged between 28 and 58 years, and presenting with nodular and papular skin lesions, gingival hyperplasia, joint contractures and bone lesions. Because of the particular clinical features and the absence of a clinical diagnosis, Whole Genome Sequencing (WGS) was carried out on DNA samples from the proband and his parents. RESULTS: A mutation in ANTXR2 (p. Gly116Val) that yielded a diagnosis of HFS was noted. CONCLUSIONS: The main goal of this paper is to add to the knowledge related to the clinical and radiographic aspects of HFS in adulthood and to show the importance of Next-Generation Sequencing (NGS) techniques in resolving such puzzling cases.

Temple-Baraitser Syndrome and Zimmermann-Laband Syndrome: one clinical entity?

BACKGROUND: KCNH1 encodes a voltage-gated potassium channel that is predominantly expressed in the central nervous system. Mutations in this gene were recently found to be responsible for Temple-Baraitser Syndrome (TMBTS) and Zimmermann-Laband syndrome (ZLS). METHODS: Here, we report a new case of TMBTS diagnosed in a Lebanese child. Whole genome sequencing was carried out on DNA samples of the proband and his parents to identify mutations associated with this disease. Sanger sequencing was performed to confirm the presence of detected variants. RESULTS: Whole genome sequencing revealed three missense mutations in TMBTS patient: c.1042G > A in KCNH1, c.2131 T > C in STK36, and c.726C > A in ZNF517. According to all predictors, mutation in KCNH1 is damaging de novo mutation that results in substitution of Glycine by Arginine, i.e., p.(Gly348Arg). This mutation was already reported in a patient with ZLS that could affect the connecting loop between helices S4-S5 of KCNH1 with a gain of function effect. CONCLUSIONS: Our findings demonstrate that KCNH1 mutations cause TMBTS and expand the mutational spectrum of KCNH1 in TMBTS. In addition, all cases of TMBTS were reviewed and compared to ZLS. We suggest that the two syndromes are a continuum and that the variability in the phenotypes is the result of the involvement of genetic modifiers.

Semantic Assembly and Annotation of Draft RNAseq Transcripts without a Reference Genome.

Transcriptomes are one of the first sources of high-throughput genomic data that have benefitted from the introduction of Next-Gen Sequencing. As sequencing technology becomes more accessible, transcriptome sequencing is applicable to multiple organisms for which genome sequences are unavailable. Currently all methods for de novo assembly are based on the concept of matching the nucleotide context overlapping between short fragments-reads. However, even short reads may still contain biologically relevant information which can be used as hints in guiding the assembly process. We propose a computational workflow for the reconstruction and functional annotation of expressed gene transcripts that does not require a reference genome sequence and can be tolerant to low coverage, high error rates and other issues that often lead to poor results of de novo assembly in studies of non-model organisms. We start with either raw sequences or the output of a context-based de novo transcriptome assembly. Instead of mapping reads to a reference genome or creating a completely unsupervised clustering of reads, we assemble the unknown transcriptome using nearest homologs from a public database as seeds. We consider even distant relations, indirectly linking protein-coding fragments to entire gene families in multiple distantly related genomes. The intended application of the proposed method is an additional step of semantic (based on relations between protein-coding fragments) scaffolding following traditional (i.e. based on sequence overlap) de novo assembly. The method we developed was effective in analysis of the jellyfish Cyanea capillata transcriptome and may be applicable in other studies of gene expression in species lacking a high quality reference genome sequence. Our algorithms are implemented in C and designed for parallel computation using a high-performance computer. The software is available free of charge via an open source license.

Obesity susceptibility loci in Qataris, a highly consanguineous Arabian population.

OBJECTIVES: In Qataris, a population characterized by a small size and a high rate of consanguinity, between two-thirds to three-quarters of adults are overweight or obese. We investigated the relevance of 23 obesity-related loci in the Qatari population. METHODS: Eight-hundred-four individuals assessed to be third generation Qataris were included in the study and assigned to 3 groups according to their body mass index (BMI): 190 lean (BMI < 25 kg/m(2)); 131 overweight (25 kg/m(2) ≤ BMI < 30 kg/m(2)) and 483 obese (BMI ≥ 30 kg/m(2)). Genomic DNA was isolated from peripheral blood and genotyped by TaqMan. RESULTS: Two loci significantly associated with obesity in Qataris: the TFAP2B variation (rs987237) (A allele versus G allele: chi-square = 10.3; P = 0.0013) and GNPDA2 variation (rs10938397) (A allele versus G allele: chi-square = 6.15; P = 0.013). The TFAP2B GG genotype negatively associated with obesity (OR = 0.21; P = 0.0031). Conversely, the GNDPA2 GG homozygous genotype associated with higher risk of obesity in subjects of age < 32 years (P = 0.0358). CONCLUSION: We showed a different genetic profile associated with obesity in the Qatari population compared to Western populations. Studying the genetic background of Qataris is of primary importance as the etiology of a given disease might be population-specific.

The molecular signature of the stroma response in prostate cancer-induced osteoblastic bone metastasis highlights expansion of hematopoietic and prostate epithelial stem cell niches.

The reciprocal interaction between cancer cells and the tissue-specific stroma is critical for primary and metastatic tumor growth progression. Prostate cancer cells colonize preferentially bone (osteotropism), where they alter the physiological balance between osteoblast-mediated bone formation and osteoclast-mediated bone resorption, and elicit prevalently an osteoblastic response (osteoinduction). The molecular cues provided by osteoblasts for the survival and growth of bone metastatic prostate cancer cells are largely unknown. We exploited the sufficient divergence between human and mouse RNA sequences together with redefinition of highly species-specific gene arrays by computer-aided and experimental exclusion of cross-hybridizing oligonucleotide probes. This strategy allowed the dissection of the stroma (mouse) from the cancer cell (human) transcriptome in bone metastasis xenograft models of human osteoinductive prostate cancer cells (VCaP and C4-2B). As a result, we generated the osteoblastic bone metastasis-associated stroma transcriptome (OB-BMST). Subtraction of genes shared by inflammation, wound healing and desmoplastic responses, and by the tissue type-independent stroma responses to a variety of non-osteotropic and osteotropic primary cancers generated a curated gene signature ("Core" OB-BMST) putatively representing the bone marrow/bone-specific stroma response to prostate cancer-induced, osteoblastic bone metastasis. The expression pattern of three representative Core OB-BMST genes (PTN, EPHA3 and FSCN1) seems to confirm the bone specificity of this response. A robust induction of genes involved in osteogenesis and angiogenesis dominates both the OB-BMST and Core OB-BMST. This translates in an amplification of hematopoietic and, remarkably, prostate epithelial stem cell niche components that may function as a self-reinforcing bone metastatic niche providing a growth support specific for osteoinductive prostate cancer cells. The induction of this combinatorial stem cell niche is a novel mechanism that may also explain cancer cell osteotropism and local interference with hematopoiesis (myelophthisis). Accordingly, these stem cell niche components may represent innovative therapeutic targets and/or serum biomarkers in osteoblastic bone metastasis.