TUBB1 mutations cause thyroid dysgenesis associated with abnormal platelet physiology

The genetic causes of congenital hypothyroidism due to thyroid dysgenesis (TD) remain largely unknown. We identified three novel TUBB1 gene mutations that co‐segregated with TD in three distinct families leading to 1.1% of TUBB1 mutations in TD study cohort. TUBB1 (Tubulin, Beta 1 Class VI) encodes for a member of the β‐tubulin protein family. TUBB1 gene is expressed in the developing and adult thyroid in humans and mice. All three TUBB1 mutations lead to non‐functional α/β‐tubulin dimers that cannot be incorporated into microtubules. In mice, Tubb1 knock‐out disrupted microtubule integrity by preventing β1‐tubulin incorporation and impaired thyroid migration and thyroid hormone secretion. In addition, TUBB1 mutations caused the formation of macroplatelets and hyperaggregation of human platelets after stimulation by low doses of agonists. Our data highlight unexpected roles for β1‐tubulin in thyroid development and in platelet physiology. Finally, these findings expand the spectrum of the rare paediatric diseases related to mutations in tubulin‐coding genes and provide new insights into the genetic background and mechanisms involved in congenital hypothyroidism and thyroid dysgenesis.

Mutations in GREB1L Cause Bilateral Kidney Agenesis in Humans and Mice

Congenital anomalies of the kidney and urinary tract (CAKUT) constitute a major cause of chronic kidney disease in children and 20% of prenatally detected anomalies. CAKUT encompass a spectrum of developmental kidney defects, including renal agenesis, hypoplasia, and cystic and non-cystic dysplasia. More than 50 genes have been reported as mutated in CAKUT-affected case subjects. However, the pathophysiological mechanisms leading to bilateral kidney agenesis (BKA) remain largely elusive. Whole-exome or targeted exome sequencing of 183 unrelated familial and/or severe CAKUT-affected case subjects, including 54 fetuses with BKA, led to the identification of 16 heterozygous variants in GREB1L (growth regulation by estrogen in breast cancer 1-like), a gene reported as a target of retinoic acid signaling. Four loss-of-function and 12 damaging missense variants, 14 being absent from GnomAD, were identified. Twelve of them were present in familial or simplex BKA-affected case subjects. Female BKA-affected fetuses also displayed uterus agenesis. We demonstrated a significant association between GREB1L variants and BKA. By in situ hybridization, we showed expression of Greb1l in the nephrogenic zone in developing mouse kidney. We generated a Greb1l knock-out mouse model by CRISPR-Cas9. Analysis at E13.5 revealed lack of kidneys and genital tract anomalies in male and female Greb1l^-/- embryos and a slight decrease in ureteric bud branching in Greb1l^+/- embryos. We showed that Greb1l invalidation in mIMCD3 cells affected tubulomorphogenesis in 3D-collagen culture, a phenotype rescued by expression of the wild-type human protein. This demonstrates that GREB1L plays a major role in early metanephros and genital development in mice and humans.

Targeted Exome Sequencing Identifies PBX1 as Involved in Monogenic Congenital Anomalies of the Kidney and Urinary Tract

Congenital anomalies of the kidney and urinary tract (CAKUT) occur in three to six of 1000 live births, represent about 20% of the prenatally detected anomalies, and constitute the main cause of CKD in children. These disorders are phenotypically and genetically heterogeneous. Monogenic causes of CAKUT in humans and mice have been identified. However, despite high-throughput sequencing studies, the cause of the disease remains unknown in most patients, and several studies support more complex inheritance and the role of environmental factors and/or epigenetics in the pathophysiology of CAKUT. Here, we report the targeted exome sequencing of 330 genes, including genes known to be involved in CAKUT and candidate genes, in a cohort of 204 unrelated patients with CAKUT; 45% of the patients were severe fetal cases. We identified pathogenic mutations in 36 of 204 (17.6%) patients. These mutations included five de novo heterozygous loss of function mutations/deletions in the PBX homeobox 1 gene (PBX1), a gene known to have a crucial role in kidney development. In contrast, the frequency of SOX17 and DSTYK variants recently reported as pathogenic in CAKUT did not indicate causality. These findings suggest that PBX1 is involved in monogenic CAKUT in humans and call into question the role of some gene variants recently reported as pathogenic in CAKUT. Targeted exome sequencing also proved to be an efficient and cost-effective strategy to identify pathogenic mutations and deletions in known CAKUT genes.

Antitumour activity of an inhibitor of miR-34a in liver cancer with β-catenin-mutations

OBJECTIVE

Hepatocellular carcinoma (HCC) is the most prevalent primary tumour of the liver. About a third of these tumours presents activating mutations of the β-catenin gene. The molecular pathogenesis of HCC has been elucidated, but mortality remains high, and new therapeutic approaches, including treatments based on microRNAs, are required. We aimed to identify candidate microRNAs, regulated by β-catenin, potentially involved in liver tumorigenesis.

DESIGN

We used a mouse model, in which β-catenin signalling was overactivated exclusively in the liver by the tamoxifen-inducible and Cre-Lox-mediated inactivation of the Apc gene. This model develops tumours with properties similar to human HCC.

RESULTS

We found that miR-34a was regulated by β-catenin, and significantly induced by the overactivation of β-catenin signalling in mouse tumours and in patients with HCC. An inhibitor of miR-34a (locked nucleic acid, LNA-34a) exerted antiproliferative activity in primary cultures of hepatocyte. This inhibition of proliferation was associated with a decrease in cyclin D1 levels, orchestrated principally by HNF-4α, a target of miR-34a considered to act as a tumour suppressor in the liver. In vivo, LNA-34a approximately halved progression rates for tumours displaying β-catenin activation together with an activation of caspases 2 and 3.

CONCLUSIONS

This work demonstrates the key oncogenic role of miR-34a in liver tumours with β-catenin gene mutations. We suggest that patients diagnosed with HCC with β-catenin mutations could be treated with an inhibitor of miR-34a. The potential value of this strategy lies in the modulation of the tumour suppressor HNF-4α, which targets cyclin D1, and the induction of a proapoptotic programme.

Integrin alpha 8 recessive mutations are responsible for bilateral renal agenesis in humans

Renal hypodysplasia (RHD) is a heterogeneous condition encompassing a spectrum of kidney development defects including renal agenesis, hypoplasia, and (cystic) dysplasia. Heterozygous mutations of several genes have been identified as genetic causes of RHD with various severity. However, these genes and mutations are not associated with bilateral renal agenesis, except for RET mutations, which could be involved in a few cases. The pathophysiological mechanisms leading to total absence of kidney development thus remain largely elusive. By using a whole-exome sequencing approach in families with several fetuses with bilateral renal agenesis, we identified recessive mutations in the integrin α8-encoding gene ITGA8 in two families. Itga8 homozygous knockout in mice is known to result in absence of kidney development. We provide evidence of a damaging effect of the human ITGA8 mutations. These results demonstrate that mutations of ITGA8 are a genetic cause of bilateral renal agenesis and that, at least in some cases, bilateral renal agenesis is an autosomal-recessive disease.

Mutations in endothelin 1 cause recessive auriculocondylar syndrome and dominant isolated question-mark ears

Auriculocondylar syndrome (ACS) is a rare craniofacial disorder with mandibular hypoplasia and question-mark ears (QMEs) as major features. QMEs, consisting of a specific defect at the lobe-helix junction, can also occur as an isolated anomaly. Studies in animal models have indicated the essential role of endothelin 1 (EDN1) signaling through the endothelin receptor type A (EDNRA) in patterning the mandibular portion of the first pharyngeal arch. Mutations in the genes coding for phospholipase C, beta 4 (PLCB4) and guanine nucleotide binding protein (G protein), alpha inhibiting activity polypeptide 3 (GNAI3), predicted to function as signal transducers downstream of EDNRA, have recently been reported in ACS. By whole-exome sequencing (WES), we identified a homozygous substitution in a furin cleavage site of the EDN1 proprotein in ACS-affected siblings born to consanguineous parents. WES of two cases with vertical transmission of isolated QMEs revealed a stop mutation in EDN1 in one family and a missense substitution of a highly conserved residue in the mature EDN1 peptide in the other. Targeted sequencing of EDN1 in an ACS individual with related parents identified a fourth, homozygous mutation falling close to the site of cleavage by endothelin-converting enzyme. The different modes of inheritance suggest that the degree of residual EDN1 activity differs depending on the mutation. These findings provide further support for the hypothesis that ACS and QMEs are uniquely caused by disruption of the EDN1-EDNRA signaling pathway.

Converging evidence for an association of ATP2B2 allelic variants with autism in male subjects

BACKGROUND

Autism is a severe developmental disorder, with strong genetic underpinnings. Previous genome-wide scans unveiled a linkage region spanning 3.5 Mb, located on human chromosome 3p25. This region encompasses the ATP2B2 gene, encoding the plasma membrane calcium-transporting ATPase 2 (PMCA2), which extrudes calcium (Ca2+) from the cytosol into the extracellular space. Multiple lines of evidence support excessive intracellular Ca2+ signaling in autism spectrum disorder (ASD), making ATP2B2 an attractive candidate gene.

METHODS

We performed a family-based association study in an exploratory sample of 277 autism genetic resource exchange families and in a replication sample including 406 families primarily recruited in Italy.

RESULTS

Several markers were significantly associated with ASD in the exploratory sample, and the same risk alleles at single nucleotide polymorphisms rs3774180, rs2278556, and rs241509 were found associated with ASD in the replication sample after correction for multiple testing. In both samples, the association was present in male subjects only. Markers associated with autism are all comprised within a single block of strong linkage disequilibrium spanning several exons, and the "risk" allele seems to follow a recessive mode of transmission.

CONCLUSIONS

These results provide converging evidence for an association between ATP2B2 gene variants and autism in male subjects, spurring interest into the identification of functional variants, most likely involved in the homeostasis of Ca2+ signaling. Additional support comes from a recent genome-wide association study by the Autism Genome Project, which highlights the same linkage disequilibrium region of the gene.

RET and GDNF mutations are rare in fetuses with renal agenesis or other severe kidney development defects

BACKGROUND

The RET/GDNF signalling pathway plays a crucial role during development of the kidneys and the enteric nervous system. In humans, RET activating mutations cause multiple endocrine neoplasia, whereas inactivating mutations are responsible for Hirschsprung disease. RET mutations have also been reported in fetuses with renal agenesis, based on analysis of a small series of samples.

OBJECTIVE AND METHODS

To characterise better the involvement of RET and GDNF in kidney development defects, a series of 105 fetuses with bilateral defects, including renal agenesis, severe hypodysplasia or multicystic dysplastic kidney, was studied. RET and GDNF coding sequences, evolutionary conserved non-coding regions (ECRs) in promoters, 3'UTRs, and RET intron 1 were analysed. Copy number variations at these loci were also investigated.

RESULTS

The study identified: (1) a low frequency (<7%) of potential mutations in the RET coding sequence, with inheritance from the healthy father for four of them; (2) no GDNF mutation; (3) similar allele frequencies in patients and controls for most single nucleotide polymorphism variants, except for RET intron 1 variant rs2506012 that was significantly more frequent in affected fetuses than in controls (6% vs 2%, p=0.01); (4) distribution of the few rare RET variants unidentified in controls into the various 5'-ECRs; (5) absence of copy number variations.

CONCLUSION

These results suggest that genomic alteration of RET or GDNF is not a major mechanism leading to renal agenesis and other severe kidney development defects. Analysis of a larger series of patients will be necessary to validate the association of the RET intron 1 variant rs2506012 with renal development defects.

Evaluating diagnostic accuracy of genetic profiles in affected offspring families

Diagnostic accuracy of a genetic test involving multiple disease genes is evaluated using sensitivity and specificity. For estimation, data from both affected and unaffected subjects are required. For early onset diseases, such as autism spectrum disorder (ASD), only data from families with affected offspring are available. To enable estimation of specificity when no data for unaffected offspring are available (single affected offspring, SAO, data), we combine the pseudocontrol method of Cordell and Clayton (Am. J. Hum. Genet. 2002; 70:124-141) with the approach of DeLong et al. (Biometrics 1985; 41:947-958) in a logistic regression model for disease outcome with a risk score (RS) constructed from genotype information as prognostic variable. The area under the receiver operating characteristic curve (AUC) is then computed using the non-parametric Mann-Whitney method. Extensive simulation studies show that, analogous to other approaches utilizing pseudocontrols, the resulting estimates of AUC using SAO data are slightly conservative when compared with the estimates computed using the full population-based data. The method is illustrated using data from a study of ASD.

Robust physical methods that enrich genomic regions identical by descent for linkage studies: confirmation of a locus for osteogenesis imperfecta

BACKGROUND

The monogenic disease osteogenesis imperfecta (OI) is due to single mutations in either of the collagen genes ColA1 or ColA2, but within the same family a given mutation is accompanied by a wide range of disease severity. Although this phenotypic variability implies the existence of modifier gene variants, genome wide scanning of DNA from OI patients has not been reported. Promising genome wide marker-independent physical methods for identifying disease-related loci have lacked robustness for widespread applicability. Therefore we sought to improve these methods and demonstrate their performance to identify known and novel loci relevant to OI.

RESULTS

We have improved methods for enriching regions of identity-by-descent (IBD) shared between related, afflicted individuals. The extent of enrichment exceeds 10- to 50-fold for some loci. The efficiency of the new process is shown by confirmation of the identification of the Col1A2 locus in osteogenesis imperfecta patients from Amish families. Moreover the analysis revealed additional candidate linkage loci that may harbour modifier genes for OI; a locus on chromosome 1q includes COX-2, a gene implicated in osteogenesis.

CONCLUSION

Technology for physical enrichment of IBD loci is now robust and applicable for finding genes for monogenic diseases and genes for complex diseases. The data support the further investigation of genetic loci other than collagen gene loci to identify genes affecting the clinical expression of osteogenesis imperfecta. The discrimination of IBD mapping will be enhanced when the IBD enrichment procedure is coupled with deep resequencing.

Convergent evidence identifying MAP/microtubule affinity-regulating kinase 1 (MARK1) as a susceptibility gene for autism

Autism spectrum disorders (ASDs) are common, heritable, but genetically heterogeneous neurodevelopmental conditions. We recently defined a susceptibility locus for ASDs on chromosome 1q41-q42. High-resolution single-nucleotide polymorphisms (126 SNPs) genotyping across the chromosome 1q41-q42 region, followed by a MARK1 (microtubule affinity-regulating kinase 1)-tagged-SNP association study in 276 families with autism from the Autism Genetic Research Exchange, showed that several SNPs within the MARK1 gene were significantly associated with ASDs by transmission disequilibrium tests. Haplotype rs12740310*C-rs3737296*G-rs12410279*A was overtransmitted (P(corrected)= 0.0016), with a relative risk for autism of 1.8 in homozygous carriers. Furthermore, ASD-associated SNP rs12410279 modulates the level of transcription of MARK1. We found that MARK1 was overexpressed in the prefrontal cortex (BA46) but not in cerebellar granule cells, on postmortem brain tissues from patients. MARK1 displayed an accelerated evolution along the lineage leading to humans, suggesting possible involvement of this gene in cognition. MARK1 encodes a kinase-regulating microtubule-dependent transport in axons and dendrites. Both overexpression and silencing of MARK1 resulted in significantly shorter dendrite length in mouse neocortical neurons and modified dendritic transport speed. As expected for a gene encoding a key polarity determinant Par-1 protein kinase, MARK1 is involved in axon-dendrite specification. Thus, MARK1 overexpression in humans may be responsible for subtle changes in dendritic functioning.

Association of autism with polymorphisms in the paired-like homeodomain transcription factor 1 (PITX1) on chromosome 5q31: a candidate gene analysis

Background

Autism is a complex, heterogeneous, behaviorally-defined disorder characterized by disruptions of the nervous system and of other systems such as the pituitary-hypothalamic axis. In a previous genome wide screen, we reported linkage of autism with a 1.2 Megabase interval on chromosome 5q31. For the current study, we hypothesized that 3 of the genes in this region could be involved in the development of autism: 1) paired-like homeodomain transcription factor 1 (PITX1), which is a key regulator of hormones within the pituitary-hypothalamic axis, 2) neurogenin 1, a transcription factor involved in neurogenesis, and 3) histone family member Y (H2AFY), which is involved in X-chromosome inactivation in females and could explain the 4:1 male:female gender distortion present in autism.

Methods

A total of 276 families from the Autism Genetic Resource Exchange (AGRE) repository composed of 1086 individuals including 530 affected children were included in the study. Single nucleotide polymorphisms tagging the three candidate genes were genotyped on the initial linkage sample of 116 families. A second step of analysis was performed using tightly linked SNPs covering the PITX1 gene. Association was evaluated using the FBAT software version 1.7.3 for single SNP analysis and the HBAT command from the same package for haplotype analysis respectively.

Results

Association between SNPs and autism was only detected for PITX1. Haplotype analysis within PITX1 showed evidence for overtransmission of the A-C haplotype of markers rs11959298 – rs6596189 (p = 0.0004). Individuals homozygous or heterozygous for the A-C haplotype risk allele were 2.54 and 1.59 fold more likely to be autistic than individuals who were not carrying the allele, respectively.

Conclusion

Strong and consistent association was observed between a 2 SNPs within PITX1 and autism. Our data suggest that PITX1, a key regulator of hormones within the pituitary-hypothalamic axis, may be implicated in the etiology of autism.

Estimating haplotype relative risks in complex disease from unphased SNPs data in families using a likelihood adjusted for ascertainment

The understanding of complex diseases and insights to improve their medical management may be achieved through the deduction of how specific haplotypes may play a joint effect to change relative risk information. In this paper we describe an ascertainment adjusted likelihood-based method to estimate haplotype relative risks using pooled family data coming from association and/or linkage studies that were used to identify specific haplotypes. Haplotype-based analysis tends to require a large amount of parameters to capture all the information that leads to efficiency problems. An adaptation of the Stochastic Expectation Maximization algorithm is used for haplotypes inference from genotypic data and to reduce the number of nuisance parameters for risk estimation. Using different simulations, we show that this method provides unbiased relative risk estimates even in case of departure from Hardy-Weinberg equilibrium.

Haplotypes in the gene encoding protein kinase c-beta (PRKCB1) on chromosome 16 are associated with autism

Autism is a developmental disorder characterized by impairments in social interaction and communication associated with repetitive patterns of interest or behavior. Autism is highly influenced by genetic factors. Genome-wide linkage and candidate gene association approaches have been used to try and identify autism genes. A few loci have repeatedly been reported linked to autism. Several groups reported evidence for linkage to a region on chromosome 16p. We have applied a direct physical identity-by-descent (IBD) mapping approach to perform a high-density (0.85 megabases) genome-wide linkage scan in 116 families from the AGRE collection. Our results confirm linkage to a region on chromosome 16p with autism. High-resolution single-nucleotide polymorphism (SNP) genotyping and analysis of this region show that haplotypes in the protein kinase c-beta gene are strongly associated with autism. An independent replication of the association in a second set of 167 trio families with autism confirmed our initial findings. Overall, our data provide evidence that the PRKCB1 gene on chromosome 16p may be involved in the etiology of autism.

On confidence intervals for genotype relative risks and attributable risks from case parent trio designs for candidate-gene studies

Scherag et al. [Hum Hered 2002;54:210-217] recently proposed point estimates and asymptotic as well as exact confidence intervals for genotype relative risks (GRRs) and the attributable risk (AR) in case parent trio designs using single nucleotide polymorphism (SNP) data. The aim of this study was the investigation of coverage probabilities and bias in estimates if the marker locus is not identical to the disease locus. Using a variety of parameter constellations, including marker allele frequencies identical to and different from the SNP at the disease locus, we performed an analytical study to quantify the bias and a Monte-Carlo simulation study for quantifying both bias and coverage probabilities. No bias was observed if marker and trait locus coincided. Two parameters had a strong impact on coverage probabilities of confidence intervals and bias in point estimates if they did not coincide: the linkage disequilibrium (LD) parameter delta and the allele frequency at the marker SNP. If marker allele frequencies were different from the allele frequencies at the functional SNP, substantial biases occurred. Further, if delta between the marker and the disease locus was lower than the maximum possible delta, estimates were also biased. In general, biases were towards the null hypothesis for both GRRs and AR. If one GRR was not increased, as e.g. in a recessive genetic model, biases away from the null could be observed. If both GRRs were in identical directions and if both were substantially larger than 1, the bias always was towards the null. When applying point estimates and confidence intervals for GRRs and AR in candidate gene studies, great care is needed. Effect estimates are substantially biased towards the null if either the allele frequencies at the marker SNP and the true disease locus are different or if the LD between the marker SNP and the disease locus is not at its maximum. A bias away from the null occurs only in uncommon study situations; it is small and can therefore be ignored for applications.

The art of pedigree drawing: algorithmic aspects

MOTIVATION

Giving a meaningful representation of a pedigree is not obvious when it includes consanguinity loops, individuals with multiple mates or several related families.

RESULTS

We show that finding a perfectly meaningful representation of a pedigree is equivalent to the interval graph sandwich problem and we propose an algorithm for drawing pedigrees.

Sib-pair linkage analysis of alcohol dependence taking into account covariates and age-of-onset variability: evaluation of the residual approach

The construction of residuals allowed us to consider a phenotype related to alcohol dependence (ALDX1) adjusted for age, gender and current smoking status. Using the Haseman and Elston regression test, we compared genome scan results for detection of loci for the unadjusted and adjusted ALDX1 phenotypes in the Collaborative Study on the Genetics of Alcoholism (COGA) data. More markers were significant at the 1% level for the binary than for the residual analyses. Overall, our residual analyses do not show substantial improvement in their ability to detect linkage for alcohol dependence in the COGA data. This may be due to the relatively small sample size of discordant sib pairs.