Search

Gephyrin (GPHN) regulates the clustering of postsynaptic components at inhibitory synapses and is involved in pathophysiology of neuropsychiatric disorders. Here, we uncover an extensive diversity of GPHN transcripts that are tightly controlled by splicing during mouse and human brain development. Proteomic analysis reveals at least a hundred isoforms of GPHN incorporated at inhibitory Glycine and GABA-A receptors containing synapses. They exhibit different localization and postsynaptic clustering properties, and altering the expression level of one isoform is sufficient to affect the number, size, and density of inhibitory synapses in cerebellar Purkinje cells. Furthermore, we discovered that splicing defects reported in neuropsychiatric disorders are carried by multiple alternative GPHN transcripts, demonstrating the need for a thorough analysis of the GPHN transcriptome in patients. Overall, we show that alternative splicing of GPHN is an important genetic variation to consider in neurological diseases and a determinant of the diversity of postsynaptic inhibitory synapses.

Indirect somatic genetic rescue (SGR) of a germline mutation is thought to be rare in inherited Mendelian disorders. Here, we establish that acquired mutations in the EIF6 gene are a frequent mechanism of SGR in Shwachman-Diamond syndrome (SDS), a leukemia predisposition disorder caused by a germline defect in ribosome assembly. Biallelic mutations in the SBDS or EFL1 genes in SDS impair release of the anti-association factor eIF6 from the 60S ribosomal subunit, a key step in the translational activation of ribosomes. Here, we identify diverse mosaic somatic genetic events (point mutations, interstitial deletion, reciprocal chromosomal translocation) in SDS hematopoietic cells that reduce eIF6 expression or disrupt its interaction with the 60S subunit, thereby conferring a selective advantage over non-modified cells. SDS-related somatic EIF6 missense mutations that reduce eIF6 dosage or eIF6 binding to the 60S subunit suppress the defects in ribosome assembly and protein synthesis across multiple SBDS-deficient species including yeast, Dictyostelium and Drosophila. Our data suggest that SGR is a universal phenomenon that may influence the clinical evolution of diverse Mendelian disorders and support eIF6 suppressor mimics as a therapeutic strategy in SDS., Shwachman-Diamond syndrome (SDS) is a leukemia predisposition disorder that is caused by defective release of eIF6 during ribosome assembly. Here the authors show that acquired somatic EIF6 mutations are frequent in the hematopoietic cells from individuals with SDS and provide a selective advantage over non-modified cells.

PURPOSE: Inborn errors of IFN-γ immunity underlie Mendelian susceptibility to mycobacterial disease (MSMD). Autosomal recessive complete IL-12Rβ1 deficiency is the most frequent genetic etiology of MSMD. Only two of the 84 known mutations are copy number variations (CNVs), identified in two of the 213 IL-12Rβ1-deficient patients and two of the 164 kindreds reported. These two CNVs are large deletions found in the heterozygous or homozygous state. We searched for novel families with IL-12Rβ1 deficiency due to CNVs. METHODS: We studied six MSMD patients from five unrelated kindreds displaying adverse reactions to BCG vaccination. Three of the patients also presented systemic salmonellosis, two had mucocutaneous candidiasis, and one had disseminated histoplasmosis. We searched for CNVs and other variation by IL12RB1-targeted next-generation sequencing (NGS). RESULTS: We identified six new IL-12Rβ1-deficient patients with a complete loss of IL-12Rβ1 expression on PHA-activated T cells and/or EBV-transformed B cells. The cells of these patients did not respond to IL-12 and IL-23. Five different CNVs encompassing IL12RB1 (four deletions and one duplication) were identified in these patients by NGS coverage analysis, in either the homozygous state (n=1) or in trans (n=4) with a single nucleotide variation (n=3) or a small indel (n=1). Seven of the nine mutations are novel. Interestingly, four of the five CNVs were predicted to be driven by nearby Alu elements, as well as the two previously reported large deletions. The IL12RB1 locus is actually enriched in Alu elements (44.7%), when compared with the rest of the genome (10.5%). CONCLUSION: The IL12RB1 locus is Alu-enriched and therefore prone to rearrangements at various positions. CNVs should be considered in the genetic diagnosis of IL-12Rβ1 deficiency.

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.

Background While mitochondrial DNA (mtDNA) copy number is strictly regulated during differentiation and according to cell type, very little is known regarding the mechanism which accurately controls mtDNA copy number in human. Exon 2 of the human POLG gene, encoding the catalytic subunit of the mitochondrial-specific DNA polymerase gamma, contains a CpG island, highly conserved in mice and human. Changes of DNA methylation at the POLG locus have been shown to modulate mtDNA copy number during cell differentiation in both mouse and human. Methods We have investigated the epigenetic modification of the POLG gene, by assessing the methylation level of its exon 2 using deep-Next Generation Sequencing analysis of bisulfite-treated DNA. Analysis were performed on various tissues at either postnatal or prenatal stages, on samples from carriers of mtDNA mutations, patients carrying two loss-of-function POLG mutations and controls. Results Very high methylation levels at POLG exon 2 were found (94±3%) and no variation was observed according to either developmental stage or tissue of origin, except for sperm samples for which lower methylation levels were found (80%). This high level of methylation was neither correlated with the presence of mtDNA mutations (94±1% of methylated alleles), nor with biallelic POLG mutations (93%±2%), even in tissues where a mtDNA depletion had been observed. Conclusions This study suggests that, at variance with mouse and un/de-differentiated human cells, differentiated human cells control mtDNA levels irrespective of POLG methylation. The factors which actually control the mtDNA levels in such cell types remain to be identified.

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.

A loss-of-function mutation in tetratricopeptide repeat domain 7A (TTC7A) is a recently identified cause of human intestinal and immune disorders. However, clues to related underlying molecular dysfunctions remain elusive. It is now shown based on the study of TTC7A-deficient and wild-type cells that TTC7A is an essential nuclear protein. It binds to chromatin, preferentially at actively transcribed regions. Its depletion results in broad range of epigenomic changes at proximal and distal transcriptional regulatory elements and in altered control of the transcriptional program. Loss of WT_TTC7A induces general decrease in chromatin compaction, unbalanced cellular distribution of histones, higher nucleosome accessibility to nuclease digestion along with genome instability, and reduced cell viability. Our observations characterize for the first time unreported functions for TTC7A in the nucleus that exert a critical role in chromatin organization and gene regulation to safeguard healthy immune and intestinal status.

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.

International audience; 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.

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.

International audience; 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.

Background During meiosis, the X and Y chromosomes are transcriptionally silenced. The persistence of repressive chromatin marks on the sex chromatin after meiosis initially led to the assumption that XY gene silencing persists to some extent in spermatids. Considering the many reports of XY-linked genes expressed and needed in the post-meiotic phase of mouse spermatogenesis, it is still unclear whether or not the mouse sex chromatin is a repressive or permissive environment, after meiosis. Results To determine the transcriptional and chromatin state of the sex chromosomes after meiosis, we re-analyzed ten ChIP-Seq datasets performed on mouse round spermatids and four RNA-seq datasets from male germ cells purified at different stages of spermatogenesis. For this, we used the last version of the genome (mm10/GRCm38) and included reads that map to several genomic locations in order to properly interpret the high proportion of sex chromosome-encoded multicopy genes. Our study shows that coverage of active epigenetic marks H3K4me3 and Kcr is similar on the sex chromosomes and on autosomes. The post-meiotic sex chromatin nevertheless differs from autosomal chromatin in its enrichment in H3K9me3 and its depletion in H3K27me3 and H4 acetylation. We also identified a posttranslational modification, H3K27ac, which specifically accumulates on the Y chromosome. In parallel, we found that the X and Y chromosomes are enriched in genes expressed post-meiotically and display a higher proportion of spermatid-specific genes compared to autosomes. Finally, we observed that portions of chromosome 14 and of the sex chromosomes share specific features, such as enrichment in H3K9me3 and the presence of multicopy genes that are specifically expressed in round spermatids, suggesting that parts of chromosome 14 are under the same evolutionary constraints than the sex chromosomes. Conclusions Based on our expression and epigenomic studies, we conclude that, after meiosis, the mouse sex chromosomes are no longer silenced but are nevertheless regulated differently than autosomes and accumulate different chromatin marks. We propose that post-meiotic selective constraints are at the basis of the enrichment of spermatid-specific genes and of the peculiar chromatin composition of the sex chromosomes and of parts of chromosome 14. Electronic supplementary material The online version of this article (doi:10.1186/s13072-016-0099-8) contains supplementary material, which is available to authorized users.

While mitochondrial DNA (mtDNA) copy number is strictly regulated during differentiation and according to cell type, very little is known regarding the mechanism which accurately controls mtDNA copy number in human. Exon 2 of the human POLG gene, encoding the catalytic subunit of the mitochondrial-specific DNA polymerase gamma, contains a CpG island, highly conserved in mice and human. Changes of DNA methylation at the POLG locus have been shown to modulate mtDNA copy number during cell differentiation in both mouse and human.We have investigated the epigenetic modification of the POLG gene, by assessing the methylation level of its exon 2 using deep-Next Generation Sequencing analysis of bisulfite-treated DNA. Analysis were performed on various tissues at either postnatal or prenatal stages, on samples from carriers of mtDNA mutations, patients carrying two loss-of-function POLG mutations and controls.Very high methylation levels at POLG exon 2 were found (94±3%) and no variation was observed according to either developmental stage or tissue of origin, except for sperm samples for which lower methylation levels were found (80%). This high level of methylation was neither correlated with the presence of mtDNA mutations (94±1% of methylated alleles), nor with biallelic POLG mutations (93%±2%), even in tissues where a mtDNA depletion had been observed.This study suggests that, at variance with mouse and un/de-differentiated human cells, differentiated human cells control mtDNA levels irrespective of POLG methylation. The factors which actually control the mtDNA levels in such cell types remain to be identified.

International audience; Objective: Hepatocellular carcinoma (HCC) is the most prevalent primary tumor of the liver. About a third of these tumors 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 signaling was overactivated exclusively in the liver, by the tamoxifen-inducible and Cre-Lox-mediated inactivation of the Apc gene. This model develops tumors with properties similar to human HCC. Results: We found that miR-34a was regulated by β-catenin, and significantly induced by the overactivation of β-catenin signaling in mouse tumors and in HCC patients. An inhibitor of miR-34a (LNA-34a) exerted anti-proliferative 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 tumor suppressor in the liver. In vivo, LNA-34a approximately halved progression rates for tumors displaying β-catenin activation together with an activation of caspases 2 and 3. Conclusion: This work demonstrates the key oncogenic role of miR-34a in liver tumors 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 tumor suppressor HNF-4α, which targets cyclin D1, and the induction of a pro-apoptotic program.

International audience; Background: cRel, a subunit of NF‑κB, is implicated in the inflammatory response observed in autoimmune disease. Hence, knocked‑out mice for cRel had a significantly higher mortality, providing new and important functions of cRel in the physiopathology of septic shock. Whether genetic variants in the human REL gene are associated with severity of septic shock is unknown. Methods: We genotyped a population of 1040 ICU patients with septic shock and 855 ICU controls for two known polymorphisms of REL; REL rs842647 and REL rs13031237. Outcome of patients according to the presence of REL variant alleles was compared. Results: The distribution of REL variant alleles was not significantly different between patients and controls. Among the septic shock group, REL rs13031237*T minor allele was not associated with worse outcome. In contrast, REL rs842647*G minor allele was significantly associated with more multi‑organ failure and early death [OR 1.4; 95 % CI (1.02–1.8)]. Conclusion: In a large ICU population, we report a significant clinical association between a variation in the human REL gene and severity and mortality of septic shock, suggesting for the first time a new insight into the role of cRel in response to infection in humans.

Exome sequencing is becoming a standard tool for gene mapping of genetic diseases. Given the vast amount of data generated by Next Generation Sequencing techniques, identification of disease causal variants is like finding a needle in a haystack. The impact assessment and the prioritization of potential pathogenic variants are expected to reduce work in biological validation, which is long and costly. One of the possible approaches to determine the most probable deleterious variants in individual exomes is to use protein function alteration prediction. We propose in this paper to use a machine learning approach, the random forest to build a new meta-score based on five previously described scores (SIFT, Polyphen2, LRT, PhyloP and MutationTaster) and compiled in the dbNSFP database. The functional meta-score was trained on a dataset of 61500 non-synonymous Single Nucleotide Polymorphisms (SNPs). The random forest method (rfPred) appears to be globally better than each of the classifiers separately or in combination in a logistic regression model, and better than a newly described score (CADD) on independent validation sets. RfPred scores have been pre-calculated for all the possible non-synonymous SNPs of human exome and are freely accessible at the web-server http://www.sbim.fr/rfPred/

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.

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.

International audience; Background The RET/GDNF signalling pathway plays a crucial role during development of kidneys and 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 better characterize the involvement of RET and GDNF in kidney development defects, we studied a series of 105 fetuses with bilateral defects including renal agenesis, severe hypodysplasia or multicystic dysplastic kidney. RET and GDNF coding sequences, evolutionary conserved non-coding regions (ECRs) in promoters, 3'UTRs and RET intron 1 were analysed. Copy number variations (CNVs) at these loci were also investigated. Results We identified: (i) a low frequency (< 7%) of potential mutations in the RET coding sequence, with inheritance from the healthy father for four of them; (ii) no GDNF mutation; (iii) similar allele frequencies in patients and controls for most SNP 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); (iv) distribution of the few rare RET variants unidentified in controls into the various 5'-ECRs; (v) absence of CNVs. 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.

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 (Ca 2+ ) from the cytosol into the extracellular space. Multiple lines of evidence support excessive intracellular Ca 2+ 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 Ca 2+ 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.

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.

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 δ 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 δ between the marker and the disease locus was lower than the maximum possible δ, 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.

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.

Background Autism is a complex disorder characterized by deficits involving communication, social interaction, and repetitive and restrictive patterns of behavior. Twin studies have shown that autism is strongly heritable, suggesting a strong genetic component. In other disease states with a complex etiology, such as type 2 diabetes, cancer and cardiovascular disease, combined analysis of multiple genetic variants in a genetic score has helped to identify individuals at high risk of disease. Genetic scores are designed to test for association of genetic markers with disease. Method The accumulation of multiple risk alleles markedly increases the risk of being affected, and compared with studying polymorphisms individually, it improves the identification of subgroups of individuals at greater risk. In the present study, we show that this approach can be applied to autism by specifically looking at a high-risk population of children who have siblings with autism. A two-sample study design and the generation of a genetic score using multiple independent genes were used to assess the risk of autism in a high-risk population. Results In both samples, odds ratios (ORs) increased significantly as a function of the number of risk alleles, with a genetic score of 8 being associated with an OR of 5.54 (95% confidence interval [CI] 2.45 to 12.49). The sensitivities and specificities for each genetic score were similar in both analyses, and the resultant area under the receiver operating characteristic curves were identical (0.59). Conclusions These results suggest that the accumulation of multiple risk alleles in a genetic score is a useful strategy for assessing the risk of autism in siblings of affected individuals, and may be better than studying single polymorphisms for identifying subgroups of individuals with significantly greater risk.

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.