Genome-wide association studies and polygenic risk score phenome-wide association studies across complex phenotypes in the human phenotype project.

BACKGROUND: Genome-wide association studies (GWASs) associate phenotypes and genetic variants across a study cohort. GWASs require large-scale cohorts with both phenotype and genetic sequencing data, limiting studied phenotypes. The Human Phenotype Project is a longitudinal study that has measured a wide range of clinical and biomolecular features from a self-assignment cohort over 5 years. The phenotypes collected are quantitative traits, providing higher-resolution insights into the genetics of complex phenotypes. METHODS: We present the results of GWASs and polygenic risk score phenome-wide association studies with 729 clinical phenotypes and 4,043 molecular features from the Human Phenotype Project. This includes clinical traits that have not been previously associated with genetics, including measures from continuous sleep monitoring, continuous glucose monitoring, liver ultrasound, hormonal status, and fundus imaging. FINDINGS: In GWAS of 8,706 individuals, we found significant associations between 169 clinical traits and 1,184 single-nucleotide polymorphisms. We found genes associated with both glycemic control and mental disorders, and we quantify the strength of genetic signals in serum metabolites. In polygenic risk score phenome-wide association studies for clinical traits, we found 16,047 significant associations. CONCLUSIONS: The entire set of findings, which we disseminate publicly, provides newfound resolution into the genetic architecture of complex human phenotypes. FUNDING: E.S. is supported by the Minerva foundation with funding from the Federal German Ministry for Education and Research and by the European Research Council and the Israel Science Foundation.

mRNA splicing is modulated by intronic microRNAs.

Splicing of transcripts is catalyzed by the spliceosome, a mega-complex consisting of hundreds of proteins and five snRNAs, which employs direct interactions. When U1 snRNA forms high-affinity binding, namely more than eight base pairs, with the 5'SS, the result is usually a suppressing effect on the splicing activity. This likely occurs due to the inefficient unwinding of U1/5'SS base-pairing or other regulatory obstructions. Here, we show in vitro and in patient-derived cell lines that pre-microRNAs can modulate the splicing reaction by interacting with U1 snRNA. This leads to reduced binding affinity to the 5'SS, and hence promotes the inclusion of exons containing 5'SS, despite sequence-based high affinity to U1. Application of the mechanism resulted in correction of the splicing defect in the disease-causing VCAN gene from an individual with Wagner syndrome. This pre-miRNA/U1 interaction can regulate the expression of alternatively spliced exons, thus extending the scope of mechanisms regulating splicing.

Phenotypic suppression of acral peeling skin syndrome in a patient with autosomal recessive congenital ichthyosis.

Autosomal recessive congenital ichthyosis (ARCI) manifests with generalized scaling often associated with generalized erythema. Mutations in at least 13 different genes have been reported to cause ARCI. Acral peeling skin syndrome (APSS) is a rare autosomal recessive disorder manifesting with peeling over the distal limbs and dorsal surfaces of hands and feet. APSS is mostly due to mutations in TGM5, encoding transglutaminase 5. Both ARCI and APSS are fully penetrant genetic traits. Here, we describe a consanguineous family in which one patient with mild ARCI was found to carry a homozygous mutation in ALOXE3 (c.1238G > A; p.Gly413Asp). The patient was also found to carry a known pathogenic homozygous mutation in TGM5 (c.1335G > C; p.Lys445Asn) but did not display acral peeling skin. Her uncle carried the same homozygous mutation in TGM5 but carried the ALOXE3 mutation in a heterozygous state and showed clinical features typical of APSS. Taken collectively, these observations suggested that the ALOXE3 mutation suppresses the clinical expression of the TGM5 variant. We hypothesized that ALOXE3 deficiency may affect the expression of a protein capable of compensating for the lack of TGM5 expression. Downregulation of ALOXE3 in primary human keratinocytes resulted in increased levels of corneodesmosin, which plays a critical role in the maintenance of cell-cell adhesion in the upper epidermal layers. Accordingly, ectopic corneodesmosin expression rescued the cell-cell adhesion defect caused by TGM5 deficiency in keratinocytes as ascertained by the dispase dissociation assay. The present data thus provide evidence for phenotypic suppression in a human hereditary skin disorder.

Comparison of breast cancer metastasis models reveals a possible mechanism of tumor aggressiveness.

In breast cancer patients, the lungs are among the first sites of cancer metastasis, and in nearly one quarter of metastatic patients, the exclusive first event. Two common mouse models mimic breast cancer lung colonization and distal metastasis: an orthotopic model and intravenous (IV) cell injections. Gene expression analysis of pulmonary lesions from these two methods demonstrated high inter-model resemblance. However, microRNA (miRNA) expression profiles were not compared. In this study, we compared the overall miRNA expression profiles (miRNome) of the orthotopic and IV breast cancer metastasis models and identified significant miRNome changes between the two models. Overexpression of the most significant candidate, miR-96 or downregulation of its validated gene-target, ABCE1 reduced cancer cells 2D/3D cell movement and proliferation in vitro, and abated tumor growth and metastasis formation in vivo. Human data analysis further strengthened miR-96/ABCE1 role in breast cancer tumor aggression. Taken together, our results indicate that IV- and orthotopic models differ by their miRNome. Specifically in our study, breast cancer aggressiveness was dictated by miR-96 regulating ABCE1. Overall, miRNome analysis of various metastatic cancer models may lead to the identification of candidate genes critical to metastasis development.

Loss of Protocadherin-12 Leads to Diencephalic-Mesencephalic Junction Dysplasia Syndrome.

OBJECTIVE: To identify causes of the autosomal-recessive malformation, diencephalic-mesencephalic junction dysplasia (DMJD) syndrome. METHODS: Eight families with DMJD were studied by whole-exome or targeted sequencing, with detailed clinical and radiological characterization. Patient-derived induced pluripotent stem cells were derived into neural precursor and endothelial cells to study gene expression. RESULTS: All patients showed biallelic mutations in the nonclustered protocadherin-12 (PCDH12) gene. The characteristic clinical presentation included progressive microcephaly, craniofacial dysmorphism, psychomotor disability, epilepsy, and axial hypotonia with variable appendicular spasticity. Brain imaging showed brainstem malformations and with frequent thinned corpus callosum with punctate brain calcifications, reflecting expression of PCDH12 in neural and endothelial cells. These cells showed lack of PCDH12 expression and impaired neurite outgrowth. INTERPRETATION: DMJD patients have biallelic mutations in PCDH12 and lack of protein expression. These patients present with characteristic microcephaly and abnormalities of white matter tracts. Such pathogenic variants predict a poor outcome as a result of brainstem malformation and evidence of white matter tract defects, and should be added to the phenotypic spectrum associated with PCDH12-related conditions. Ann Neurol 2018;84:646-655.

Analysis of microRNAs in familial Mediterranean fever.

OBJECTIVES: Although Familial Mediterranean fever (FMF) is categorized as autosomal recessive, frequent exceptions to this model exist and therefore we aimed to search epigenetic modifications in this disease. METHODS: Ten M694V homozygous FMF patients (the most severe phenotype) were recruited for this study. Patients with inflammatory flare were excluded. Total RNA was extracted from peripheral blood, and microRNA expression profiled using NanoString nCounter technology. These patients were compared to 10 healthy age- and sex-matched controls. RESULTS: Seven hundred nighty-eight mature human miRNAs were probed, 103 of which had expression levels above the negative control probes. Seven miRNAs showed significant differences in expression in samples from FMF patients compared to healthy controls: four miRNAs were upregulated (miR-144-3p, miR-21-5p, miR-4454, and miR-451a), and three were downregulated (miR-107, let-7d-5p, and miR-148b-3p). CONCLUSION: In this pilot study, we identified epigenetic modifications in clinically quiescent FMF patients. More studies are required for exploration of their contribution to FMF pathogenesis and their potential role as clinical biomarkers.

Differential analysis of mutations in the Jewish population and their implications for diseases.

Sequencing large cohorts of ethnically homogeneous individuals yields genetic insights with implications for the entire population rather than a single individual. In order to evaluate the genetic basis of certain diseases encountered at high frequency in the Ashkenazi Jewish population (AJP), as well as to improve variant annotation among the AJP, we examined the entire exome, focusing on specific genes with known clinical implications in 128 Ashkenazi Jews and compared these data to other non-Jewish populations (European, African, South Asian and East Asian). We targeted American College of Medical Genetics incidental finding recommended genes and the Catalogue of Somatic Mutations in Cancer (COSMIC) germline cancer-related genes. We identified previously known disease-causing variants and discovered potentially deleterious variants in known disease-causing genes that are population specific or substantially more prevalent in the AJP, such as in the ATP and HGFAC genes associated with colorectal cancer and pancreatic cancer, respectively. Additionally, we tested the advantage of utilizing the database of the AJP when assigning pathogenicity to rare variants of independent whole-exome sequencing data of 49 Ashkenazi Jew early-onset breast cancer (BC) patients. Importantly, population-based filtering using our AJP database enabled a reduction in the number of potential causal variants in the BC cohort by 36%. Taken together, population-specific sequencing of the AJP offers valuable, clinically applicable information and improves AJP filter annotation.

X-linked elliptocytosis with impaired growth is related to mutated AMMECR1.

In this study, we report a family with X-linked recessive syndrome caused by mutated AMMECR1 and characterized by elliptocytosis with or without anemia, midface hypoplasia, proportionate short stature and hearing loss. Recently, mutations in AMMECR1 were reported in two maternal half-brothers, presenting with nephrocalcinosis, midface hypoplasia and, in one of the siblings, deafness and elliptocytosis. AMMECR1 gene is localized in the critical region of contiguous deletion syndrome on Xq22.3 implicated in Alport syndrome, mental retardation, midface hypoplasia, and elliptocytosis (AMME complex). Interestingly, alternative splicing of exon 2, the same exon harboring the truncating mutation, was observed in the proband and in his unaffected mother. Alternative splicing of this exon is predicted to lead to an in-frame deletion. We provide further evidence that mutated AMMECR1 gene is responsible for this clinically recognizable X-linked condition with variable expressivity.

Mutations in TSPEAR, Encoding a Regulator of Notch Signaling, Affect Tooth and Hair Follicle Morphogenesis.

Despite recent advances in our understanding of the pathogenesis of ectodermal dysplasias (EDs), the molecular basis of many of these disorders remains unknown. In the present study, we aimed at elucidating the genetic basis of a new form of ED featuring facial dysmorphism, scalp hypotrichosis and hypodontia. Using whole exome sequencing, we identified 2 frameshift and 2 missense mutations in TSPEAR segregating with the disease phenotype in 3 families. TSPEAR encodes the thrombospondin-type laminin G domain and EAR repeats (TSPEAR) protein, whose function is poorly understood. TSPEAR knock-down resulted in altered expression of genes known to be regulated by NOTCH and to be involved in murine hair and tooth development. Pathway analysis confirmed that down-regulation of TSPEAR in keratinocytes is likely to affect Notch signaling. Accordingly, using a luciferase-based reporter assay, we showed that TSPEAR knock-down is associated with decreased Notch signaling. In addition, NOTCH1 protein expression was reduced in patient scalp skin. Moreover, TSPEAR silencing in mouse hair follicle organ cultures was found to induce apoptosis in follicular epithelial cells, resulting in decreased hair bulb diameter. Collectively, these observations indicate that TSPEAR plays a critical, previously unrecognized role in human tooth and hair follicle morphogenesis through regulation of the Notch signaling pathway.

Local microRNA delivery targets Palladin and prevents metastatic breast cancer.

Metastasis is the primary cause for mortality in breast cancer. MicroRNAs, gene expression master regulators, constitute an attractive candidate to control metastasis. Here we show that breast cancer metastasis can be prevented by miR-96 or miR-182 treatment, and decipher the mechanism of action. We found that miR-96/miR-182 downregulate Palladin protein levels, thereby reducing breast cancer cell migration and invasion. A common SNP, rs1071738, at the miR-96/miR-182-binding site within the Palladin 3'-UTR abolishes miRNA:mRNA binding, thus diminishing Palladin regulation by these miRNAs. Regulation is successfully restored by applying complimentary miRNAs. A hydrogel-embedded, gold-nanoparticle-based delivery vehicle provides efficient local, selective, and sustained release of miR-96/miR-182, markedly suppressing metastasis in a breast cancer mouse model. Combined delivery of the miRNAs with a chemotherapy drug, cisplatin, enables significant primary tumour shrinkage and metastasis prevention. Our data corroborate the role of miRNAs in metastasis, and suggest miR-96/miR-182 delivery as a potential anti-metastatic drug.

Somatic Mosaicism for a "Lethal" GJB2 Mutation Results in a Patterned Form of Spiny Hyperkeratosis without Eccrine Involvement.

BACKGROUND: Spiny hyperkeratosis refers to a rare clinical phenotype characterized by nonfollicular keratotic projections and sometimes associated with other acquired and inherited conditions. We describe a case of congenital patterned spiny hyperkeratosis. METHODS: To identify the cause of this disorder, we used a combination of whole exome sequencing, direct sequencing and TaqMan assay. RESULTS: We found that the peculiar clinical features displayed by the patient are due to somatic mosaicism for a heterozygous mutation in the GJB2 gene. CONCLUSION: Because histopathologic examination of two independent biopsies did not reveal porokeratotic eccrine ostial and dermal duct nevus (PEODDN), previously reported to result from somatic mutations in GJB2, it appears that mutations in this gene can cause nevoid spiny hyperkeratosis in the context of PEODDN or as an isolated finding.

Whole-exome sequencing in individuals with multiple cardiovascular risk factors and normal coronary arteries.

OBJECTIVES: Most studies on the genes involved in coronary artery disease (CAD) targeted individuals with angiographically or clinically proven CAD. Focusing on high-risk individuals with normal coronary arteries (NCA) may offer novel insights into the pathogenesis of CAD. We aimed to identify genes putatively protective for development of CAD. METHODS: Pooled whole-exome sequencing (WES) was performed on 17 patients with multiple cardiovascular risk factors and NCA and on 17 controls with multivessel CAD. Rare NCA-unique sequence variants were subsequently individually validated using the Fluidigm platform in 100 additional CAD controls and 100 general population controls. RESULTS: In total, 555 100 variants were detected in at least one WES pool in the study group and in none of the control WES pools. For second phase validation, we focused on rare, nonsynonymous variants, resulting in a total of 144 variants in 40 genes, of which 96 were selected for subsequent genotyping. Validation phase genotyping resulted in 19 variants in 16 genes that were found in the NCA group and in none of the CAD controls. The SPTBN5, NID2, and ADAMTSL4 genes harbored sequence variants in more than one CAD-protected patient and none of the 117 CAD controls. CONCLUSION: Applying WES technology and focusing on individuals seemingly protected from developing CAD successfully identified 19 variants that may offer protection from CAD by undetermined mechanisms. Studying the genetics of high-risk individuals apparently protected from CAD may provide novel insights into the pathogenesis of CAD.

The risk for developing cancer in Israeli ATM, BLM, and FANCC heterozygous mutation carriers.

Cancer risks in heterozygous mutation carriers of the ATM, BLM, and FANCC genes are controversial. To shed light on this issue, cancer rates were evaluated by cross referencing asymptomatic Israeli heterozygous mutation carriers in the ATM, BLM, and FANCC genes with cancer diagnoses registered at the Israeli National Cancer Registry (INCR). Comparison of observed to expected Standardized Incidence Rates (SIR) was performed. Overall, 474 individuals participated in the study: 378 females; 25 Arab and 31 Jewish ATM carriers, 152 BLM carriers, and 170 FANCC carriers (all Ashkenazim). Age range at genotyping was 19-53 years (mean + SD 30.6 + 5 years). In addition, 96 males were included; 5, 34, and 57 ATM, BLM, and FANCC mutation carriers, respectively. Over 5-16 years from genotyping (4721 person/years), 15 new cancers were diagnosed in mutation carriers: 5 breast, 4 cervical, 3 melanomas, and one each bone sarcoma, pancreatic, and colorectal cancer. No single cancer diagnosis was more prevalent then expected in all groups combined or per gene analyzed. Specifically breast cancer SIR was 0.02-0.77. We conclude that Israeli ATM, BLM, and FANCC heterozygous mutation carriers are not at an increased risk for developing cancer.

Whole-Genome Sequencing Analysis from the Chikungunya Virus Caribbean Outbreak Reveals Novel Evolutionary Genomic Elements.

BACKGROUND: Chikungunya virus (CHIKV), an alphavirus and member of the Togaviridae family, is capable of causing severe febrile disease in humans. In December of 2013 the Asian Lineage of CHIKV spread from the Old World to the Americas, spreading rapidly throughout the New World. Given this new emergence in naïve populations we studied the viral genetic diversity present in infected individuals to understand how CHIKV may have evolved during this continuing outbreak. METHODOLOGY/PRINCIPLE FINDINGS: We used deep-sequencing technologies coupled with well-established bioinformatics pipelines to characterize the minority variants and diversity present in CHIKV infected individuals from Guadeloupe and Martinique, two islands in the center of the epidemic. We observed changes in the consensus sequence as well as a diverse range of minority variants present at various levels in the population. Furthermore, we found that overall diversity was dramatically reduced after single passages in cell lines. Finally, we constructed an infectious clone from this outbreak and identified a novel 3' untranslated region (UTR) structure, not previously found in nature, that led to increased replication in insect cells. CONCLUSIONS/SIGNIFICANCE: Here we preformed an intrahost quasispecies analysis of the new CHIKV outbreak in the Caribbean. We identified novel variants present in infected individuals, as well as a new 3'UTR structure, suggesting that CHIKV has rapidly evolved in a short period of time once it entered this naïve population. These studies highlight the need to continue viral diversity surveillance over time as this epidemic evolves in order to understand the evolutionary potential of CHIKV.

Actionable clinical decisions based on comprehensive genomic evaluation in asymptomatic adults.

Whole-exome sequencing (WES) arises as a new approach in diagnosing individuals affected by multigenic and complex phenotypes. Herein, we aim to examine whether WES is useful in screening asymptomatic individuals for actionable interventions, which has not yet been established. Twenty-five healthy adults underwent WES, bioinformatics, and manual curation of their exomes. Six participants (24%) harbored significant, management-changing variants in cancer predisposition genes, American College of Medical Genetics, and genomics reportable cardiac diseases and pharmacogenomic biomarkers that have led to clinical recommendations and interventions. Furthermore, more than 80% of the participants (21) carried 1-3 genetic variants with an associated clinical guideline for an altered drug dosing or administration based on the FDA's table of pharmacogenomics. These results support WES potential not only to answer specific diagnostic questions presented by the relevant personal and/or family history but also to uncover clinically important genetic findings unrelated to the primary indication for sequencing.

Cancer risks in Jewish male BRCA1 and BRCA2 mutation carriers.

Cancer risks and tumor types in male BRCA1 and BRCA2 mutation carriers are still unsettled. Cancer risks in men who were found to harbor a BRCA1 (n = 150) or a BRCA2 (n = 88) mutation or both (n = 2) were assessed by cross referencing with data on cancer occurrence in the Israeli National Cancer Registry. Incidence rates in mutation carriers were compared with men who were counseled, genotyped, and found not to harbor the familial mutation (true negative n = 122), and with standardized incidence rates (SIRs). Of 210 cancer-free individuals at initial counseling, 11 cancers were diagnosed after a mean follow-up of 5.06 ± 4.1 years (1064 person/years) compared with 1/122 in a BRCA true-negative man. The SIR for all BRCA1/2 mutation carriers compared with the rates in the general population were elevated for pancreatic cancer [2.97 (95 % CI 1.83-4.29)] and breast cancer [16.44 (95 % CI 9.65-26.24)]. For prostate cancer these rates were 0.59 (95 % CI 0.4-0.84). Jewish BRCA1/2 mutation carriers are at an increased risk for breast and pancreatic, but not prostate cancer. These cancer risks and the consequent recommendations, if validated, should be transmitted to carriers at test result disclosure.

Amerindian-specific regions under positive selection harbour new lipid variants in Latinos.

Dyslipidemia and obesity are especially prevalent in populations with Amerindian backgrounds, such as Mexican-Americans, which predispose these populations to cardiovascular disease. Here we design an approach, known as the cross-population allele screen (CPAS), which we conduct prior to a genome-wide association study (GWAS) in 19,273 Europeans and Mexicans, in order to identify Amerindian risk genes in Mexicans. Utilizing CPAS to restrict the GWAS input variants to only those differing in frequency between the two populations, we identify novel Amerindian lipid genes, receptor-related orphan receptor alpha (RORA) and salt-inducible kinase 3 (SIK3), and three loci previously unassociated with dyslipidemia or obesity. We also detect lipoprotein lipase (LPL) and apolipoprotein A5 (APOA5) harbouring specific Amerindian signatures of risk variants and haplotypes. Notably, we observe that SIK3 and one novel lipid locus underwent positive selection in Mexicans. Furthermore, after a high-fat meal, the SIK3 risk variant carriers display high triglyceride levels. These findings suggest that Amerindian-specific genetic architecture leads to a higher incidence of dyslipidemia and obesity in modern Mexicans.

The WWOX gene modulates high-density lipoprotein and lipid metabolism.

BACKGROUND: Low levels of high-density lipoprotein (HDL) cholesterol constitutes a major risk factor for atherosclerosis. Recent studies from our group reported a genetic association between the WW domain-containing oxidoreductase (WWOX) gene and HDL cholesterol levels. Here, through next-generation resequencing, in vivo functional studies and gene microarray analyses, we investigated the role of WWOX in HDL and lipid metabolism. METHODS AND RESULTS: Using next-generation resequencing of the WWOX region, we first identified 8 variants significantly associated and perfectly segregating with the low-HDL trait in 2 multigenerational French Canadian dyslipidemic families. To understand in vivo functions of WWOX, we used liver-specific Wwox(hep-/-) and total Wwox(-/-) mice models, where we found decreased ApoA-I and Abca1 levels in hepatic tissues. Analyses of lipoprotein profiles in Wwox(-/-), but not Wwox(hep-/-) littermates, also showed marked reductions in serum HDL cholesterol concentrations, concordant with the low-HDL findings observed in families. We next obtained evidence of a sex-specific effect in female Wwox(hep-/-) mice, where microarray analyses revealed an increase in plasma triglycerides and altered lipid metabolic pathways. We further identified a significant reduction in ApoA-I and Lpl and an upregulation in Fas, Angptl4, and Lipg, suggesting that the effects of Wwox involve multiple pathways, including cholesterol homeostasis, ApoA-I/ABCA1 pathway, and fatty acid biosynthesis/triglyceride metabolism. CONCLUSIONS: Our data indicate that WWOX disruption alters HDL and lipoprotein metabolism through several mechanisms and may account for the low-HDL phenotype observed in families expressing the WWOX variants. These findings thus describe a novel gene involved in cellular lipid homeostasis, which effects may impact atherosclerotic disease development.

Genomic study in Mexicans identifies a new locus for triglycerides and refines European lipid loci.

BACKGROUND: The Mexican population and others with Amerindian heritage exhibit a substantial predisposition to dyslipidemias and coronary heart disease. Yet, these populations remain underinvestigated by genomic studies, and to date, no genome-wide association (GWA) studies have been reported for lipids in these rapidly expanding populations. METHODS AND FINDINGS: We performed a two-stage GWA study for hypertriglyceridemia and low high-density lipoprotein cholesterol (HDL-C) in Mexicans (n=4361), and identified a novel Mexican-specific genome-wide significant locus for serum triglycerides (TGs) near the Niemann-Pick type C1 protein gene (p=2.43×10(-08)). Furthermore, three European loci for TGs (APOA5, GCKR and LPL), and four loci for HDL-C (ABCA1, CETP, LIPC and LOC55908) reached genome-wide significance in Mexicans. We used cross-ethnic mapping to narrow three European TG GWA loci, APOA5, MLXIPL, and CILP2 that were wide and contained multiple candidate variants in the European scan. At the APOA5 locus, this reduced the most likely susceptibility variants to one, rs964184. Importantly, our functional analysis demonstrated a direct link between rs964184 and postprandial serum apoAV protein levels, supporting rs964184 as the causative variant underlying the European and Mexican GWA signal. Overall, 52 of the 100 reported associations from European lipid GWA meta-analysis generalised to Mexicans. However, in 82 of the 100 European GWA loci, a different variant other than the European lead/best-proxy variant had the strongest regional evidence of association in Mexicans. CONCLUSIONS: This first Mexican GWA study of lipids identified a novel GWA locus for high TG levels; used the interpopulation heterogeneity to significantly restrict three previously known European GWA signals, and surveyed whether the European lipid GWA SNPs extend to the Mexican population.