The BCAR1 locus in carotid intima-media thickness and atherosclerosis

Carotid intima media thickness (IMT) is a marker of subclinical atherosclerosis that can predict cardiovascular events over traditional risk factors. This thesis focused on a chromosome 16 locus associated with IMT and coronary artery disease, in order to identify functional variation affecting protein structure or gene expression, and investigate the role of novel genes in atherosclerosis. Bioinformatics tools were used to identify variants in strong LD with the lead SNP and to filter these to a shortlist of potential regulatory variants. Electrophoretic mobility shift assays on these 6 SNPs detected allele-specific protein binding to the lead SNP rs4888378 in CFDP1, and implicated the protein FOXA. Luciferase reporter assays showed a 35-92% decrease in gene expression with the A allele. Expression-QTL analysis confirmed associations of the protective allele of rs4888378 with higher expression of BCAR1 in vascular tissues. Genotyping and analysis of the lead SNP in the PLIC cohort of 2144 healthy men and women suggested a sex-specific effect of the SNP on IMT progression. Meta-analysis of five cohort studies supported a protective effect of the A allele on common-carotid IMT in women only. As illustrated by this locus, functional variation often lies in enhancers far from its target promoter. Circular chromosome conformation capture (4C) was explored to investigate interactions between enhancers and promoters at this locus, with the aim of capturing regions interacting with rs4888378 and the BCAR1 promoter. Analysis of exome sequencing data identified a SNP in BCAR1 (coding amino acid change P76S), which was genotyped in two cohorts (IMPROVE and PLIC). Associations with lower plaque and IMT were found (wild-type proline form), but differed between the cohorts and require further validation. Wild-type and variant vectors were expressed in cells to assess the variant’s effect on protein function and pathways involved in plaque, showing cells with wild-type protein to migrate more quickly. These analyses on regulatory and protein-coding mechanisms implicate BCAR1 in atherosclerosis and provide new pathways for analysis in the understanding of cardiovascular disease

High-temporal resolution fluvial sediment source fingerprinting with uncertainty: a Bayesian approach

This contribution addresses two developing areas of sediment fingerprinting research. Specifically, how to improve the temporal resolution of source apportionment estimates whilst minimizing analytical costs and, secondly, how to consistently quantify all perceived uncertainties associated with the sediment mixing model procedure. This first matter is tackled by using direct X-ray fluorescence spectroscopy (XRFS) and diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS) analyses of suspended particulate matter (SPM) covered filter papers in conjunction with automatic water samplers. This method enables SPM geochemistry to be quickly, accurately, inexpensively and non-destructively monitored at high-temporal resolution throughout the progression of numerous precipitation events. We then employed a Bayesian mixing model procedure to provide full characterization of spatial geochemical variability, instrument precision and residual error to yield a realistic and coherent assessment of the uncertainties associated with source apportionment estimates. Applying these methods to SPM data from the River Wensum catchment, UK, we have been able to apportion, with uncertainty, sediment contributions from eroding arable topsoils, damaged road verges and combined subsurface channel bank and agricultural field drain sources at 60- and 120-minute resolution for the duration of five precipitation events. The results presented here demonstrate how combining Bayesian mixing models with the direct spectroscopic analysis of SPM-covered filter papers can produce high-temporal resolution source apportionment estimates that can assist with the appropriate targeting of sediment pollution mitigation measures at a catchment level

SMAD6 variants in craniosynostosis : genotype and phenotype evaluation

PURPOSE: Enrichment of heterozygous missense and truncating SMAD6 variants was previously reported in nonsyndromic sagittal and metopic synostosis, and interaction of SMAD6 variants with a common polymorphism near BMP2 (rs1884302) was proposed to contribute to inconsistent penetrance. We determined the occurrence of SMAD6 variants in all types of craniosynostosis, evaluated the impact of different missense variants on SMAD6 function, and tested independently whether rs1884302 genotype significantly modifies the phenotype. METHODS: We performed resequencing of SMAD6 in 795 unsolved patients with any type of craniosynostosis and genotyped rs1884302 in SMAD6-positive individuals and relatives. We examined the inhibitory activity and stability of SMAD6 missense variants. RESULTS: We found 18 (2.3%) different rare damaging SMAD6 variants, with the highest prevalence in metopic synostosis (5.8%) and an 18.3-fold enrichment of loss-of-function variants comparedwith gnomAD data (P < 10-7). Combined with eight additional variants, ≥20/26 were transmitted from an unaffected parent but rs1884302 genotype did not predict phenotype. CONCLUSION: Pathogenic SMAD6 variants substantially increase the risk of both nonsyndromic and syndromic presentations of craniosynostosis, especially metopic synostosis. Functional analysis is important to evaluate missense variants. Genotyping of rs1884302 is not clinically useful. Mechanisms to explain the remarkable diversity of phenotypes associated with SMAD6 variants remain obscure

Identifying low density lipoprotein cholesterol associated variants in the Annexin A2 (ANXA2) gene.

BACKGROUND AND AIMS: Annexin-A2 (AnxA2) is an endogenous inhibitor of proprotein convertase subtilisin/kexin type-9 (PCSK9). The repeat-one (R1) domain of AnxA2 binds to PCSK9, blocking its ability to promote degradation of low-density lipoprotein cholesterol-receptors (LDL-R) and thereby regulate low-density lipoprotein cholesterol (LDL-C) levels. Here we identify variants in ANXA2 influencing LDL-C levels and we determine the molecular mechanisms of their effects. RESULTS: The ANXA2 single nucleotide polymorphism (SNP) genotype-phenotype association was examined using the Second-Northwick-Park Heart Study (NPHSII) (n∼2700) and the UCL-LSHTM-Edinburgh-Bristol (UCLEB) consortium (n∼14,600). The ANXA2-R1 domain coding-SNP rs17845226 (V98L) associated with LDL-C, homozygotes for the minor allele having ≈18.8% higher levels of LDL-C (p = 0.004), and higher risk of coronary heart disease (CHD) (p = 0.04). The SNP is in modest linkage disequilibrium (r(2) > 0.5) with two intergenic SNPs, rs17191344 and rs11633032. Both SNPs showed allele-specific protein binding, and the minor alleles caused significant reduction in reporter gene expression (≈18%, p < 0.001). In the expression quantitative trait loci (eQTL) study, minor allele homozygotes have significantly lower levels of ANXA2-mRNA expression (p = 1.36 × 10(-05)). CONCLUSIONS: Both rs11633032 and rs17191344 SNPs are functional variants, where the minor alleles create repressor-binding protein sites for transcription factors that contribute to reduced ANXA2 gene expression. Lower AnxA2 levels could increase plasma levels of PCSK9 and thus increase LDL-C levels and risk of CHD. This supports, for the first time in humans, previous observations in mouse models that changes in the levels of AnxA2 directly influence plasma LDL-C levels, and thus implicate this protein as a potential therapeutic target for LDL-C lowering

The inner junction protein CFAP20 functions in motile and non-motile cilia and is critical for vision

Motile and non-motile cilia are associated with mutually-exclusive genetic disorders. Motile cilia propel sperm or extracellular fluids, and their dysfunction causes primary ciliary dyskinesia. Non-motile cilia serve as sensory/signalling antennae on most cell types, and their disruption causes single-organ ciliopathies such as retinopathies or multi-system syndromes. CFAP20 is a ciliopathy candidate known to modulate motile cilia in unicellular eukaryotes. We demonstrate that in zebrafish, cfap20 is required for motile cilia function, and in C. elegans, CFAP-20 maintains the structural integrity of non-motile cilia inner junctions, influencing sensory-dependent signalling and development. Human patients and zebrafish with CFAP20 mutations both exhibit retinal dystrophy. Hence, CFAP20 functions within a structural/functional hub centered on the inner junction that is shared between motile and non-motile cilia, and is distinct from other ciliopathy-associated domains or macromolecular complexes. Our findings suggest an uncharacterised pathomechanism for retinal dystrophy, and potentially for motile and non-motile ciliopathies in general.</p

Identification of 4 novel human ocular coloboma genes ANK3, BMPR1B, PDGFRA, and CDH4 through evolutionary conserved vertebrate gene analysis

Purpose: Ocular coloboma arises from genetic or environmental perturbations that inhibit optic fissure (OF) fusion during early eye development. Despite high genetic heterogeneity, 70% to 85% of patients remain molecularly undiagnosed. In this study, we have identified new potential causative genes using cross-species comparative meta-analysis. Methods: Evolutionarily conserved differentially expressed genes were identified through in silico analysis, with in situ hybridization, gene knockdown, and rescue performed to confirm spatiotemporal gene expression and phenotype. Interrogation of the 100,000 Genomes Project for putative pathogenic variants was performed. Results: Nine conserved differentially expressed genes between zebrafish and mouse were identified. Expression of zebrafish ank3a, bmpr1ba/b, cdh4, and pdgfaa was localized to the OF, periocular mesenchyme cells, or ciliary marginal zone, regions traversed by the OF. Knockdown of ank3, bmpr1b, and pdgfaa revealed a coloboma and/or microphthalmia phenotype. Novel pathogenic variants in ANK3, BMPR1B, PDGFRA, and CDH4 were identified in 8 unrelated coloboma families. We showed BMPR1B rescued the knockdown phenotype but variant messenger RNAs failed, providing evidence of pathogenicity. Conclusion: We show the utility of cross-species meta-analysis to identify several novel coloboma disease-causing genes. There is a potential to increase the diagnostic yield for new and unsolved patients while adding to our understanding of the genetic basis of OF morphogenesis

Assessing the digenic model in rare disorders using population sequencing data

An important fraction of patients with rare disorders remains with no clear genetic diagnostic, even after whole-exome or whole-genome sequencing, posing a difficulty in giving adequate treatment and genetic counseling. The analysis of genomic data in rare disorders mostly considers the presence of single gene variants in coding regions that follow a concrete monogenic mode of inheritance. A digenic inheritance, with variants in two functionally-related genes in the same individual, is a plausible alternative that might explain the genetic basis of the disease in some cases. In this case, digenic disease combinations should be absent or underrepresented in healthy individuals. We develop a framework to evaluate the significance of digenic combinations and test its statistical power in different scenarios. We suggest that this approach will be relevant with the advent of new sequencing efforts including hundreds of thousands of samples

Author Correction: Nuclear-mitochondrial DNA segments resemble paternally inherited mitochondrial DNA in humans

An amendment to this paper has been published and can be accessed via a link at the top of the paper

Clinical, genetic, epidemiologic, evolutionary, and functional delineation of TSPEAR-related autosomal recessive ectodermal dysplasia 14

TSPEAR variants cause autosomal recessive ectodermal dysplasia (ARED) 14. The function of TSPEAR is unknown. The clinical features, the mutation spectrum, and the underlying mechanisms of ARED14 are poorly understood. Combining data from new and previously published individuals established that ARED14 is primarily characterized by dental anomalies such as conical tooth cusps and hypodontia, like those seen in individuals with WNT10A-related odontoonychodermal dysplasia. AlphaFold-predicted structure-based analysis showed that most of the pathogenic TSPEAR missense variants likely destabilize the β-propeller of the protein. Analysis of 100000 Genomes Project (100KGP) data revealed multiple founder TSPEAR variants across different populations. Mutational and recombination clock analyses demonstrated that non-Finnish European founder variants likely originated around the end of the last ice age, a period of major climatic transition. Analysis of gnomAD data showed that the non-Finnish European population TSPEAR gene-carrier rate is ∼1/140, making it one of the commonest AREDs. Phylogenetic and AlphaFold structural analyses showed that TSPEAR is an ortholog of drosophila Closca, an extracellular matrix-dependent signaling regulator. We, therefore, hypothesized that TSPEAR could have a role in enamel knot, a structure that coordinates patterning of developing tooth cusps. Analysis of mouse single-cell RNA sequencing (scRNA-seq) data revealed highly restricted expression of Tspear in clusters representing enamel knots. A tspeara−/−;tspearb−/− double-knockout zebrafish model recapitulated the clinical features of ARED14 and fin regeneration abnormalities of wnt10a knockout fish, thus suggesting interaction between tspear and wnt10a. In summary, we provide insights into the role of TSPEAR in ectodermal development and the evolutionary history, epidemiology, mechanisms, and consequences of its loss of function variants

Biallelic variants in KARS1 are associated with neurodevelopmental disorders and hearing loss recapitulated by the knockout zebrafish

Purpose: Pathogenic variants in Lysyl-tRNA synthetase 1 (KARS1) have increasingly been recognized as a cause of early-onset complex neurological phenotypes. To advance the timely diagnosis of KARS1-related disorders, we sought to delineate its phenotype and generate a disease model to understand its function in vivo. Methods: Through international collaboration, we identified 22 affected individuals from 16 unrelated families harboring biallelic likely pathogenic or pathogenic in KARS1 variants. Sequencing approaches ranged from disease-specific panels to genome sequencing. We generated loss-of-function alleles in zebrafish. Results: We identify ten new and four known biallelic missense variants in KARS1 presenting with a moderate-to-severe developmental delay, progressive neurological and neurosensory abnormalities, and variable white matter involvement. We describe novel KARS1-associated signs such as autism, hyperactive behavior, pontine hypoplasia, and cerebellar atrophy with prevalent vermian involvement. Loss of kars1 leads to upregulation of p53, tissue-specific apoptosis, and downregulation of neurodevelopmental related genes, recapitulating key tissue-specific disease phenotypes of patients. Inhibition of p53 rescued several defects of kars1−/− knockouts. Conclusion: Our work delineates the clinical spectrum associated with KARS1 defects and provides a novel animal model for KARS1-related human diseases revealing p53 signaling components as potential therapeutic targets