Faculty Bibliography

The red blood cell related traits are highly heritable but their genetics are poorly defined. Only 5-10% of the total observed variance is explained by the genetic loci found to date, suggesting that additional loci should be searched using approaches alternative to large meta analysis. GWAS (Genome Wide Association Study) for red blood cell traits in a founder population cohort from Northern Italy identified a new locus for mean corpuscular hemoglobin concentration (MCHC) in the TAF3 gene. The association was replicated in two cohorts (rs1887582, P = 4.25E-09). TAF3 encodes a transcription cofactor that participates in core promoter recognition complex, and is involved in zebrafish and mouse erythropoiesis. We show here that TAF3 is required for transcription of the SPTA1 gene, encoding alpha spectrin, one of the proteins that link the plasma membrane to the actin cytoskeleton. Mutations in SPTA1 are responsible for hereditary spherocytosis, a monogenic disorder of MCHC, as well as for the normal MCHC level. Based on our results, we propose that TAF3 is required for normal erythropoiesis in human and that it might have a role in controlling the ratio between hemoglobin (Hb) and cell volume and in the dynamics of RBC maturation in healthy individuals. Finally, TAF3 represents a potential candidate or a modifier gene for disorders of red cell membrane.

Multiple lines of genetic evidence suggest a role for CNTNAP2 in autism. To assess its population impact we studied 2148 common single nucleotide polymorphisms (SNPs) using transmission disequilibrium test (TDT) across the entire ~3.3 Mb CNTNAP2 locus in 186 (408 trios) multiplex and 323 simplex families with autistic spectrum disorder (ASD). This analysis yielded two SNPs with nominal statistical significance (rs17170073, p = 2.0 x 10(-4); rs2215798, p = 1.6 x 10(-4)) that did not survive multiple testing. In a combined analysis of all families, two highly correlated (r (2) = 0.99) SNPs in intron 14 showed significant association with autism (rs2710093, p = 9.0 x 10(-6); rs2253031, p = 2.5 x 10(-5)). To validate these findings and associations at SNPs from previous autism studies (rs7794745, rs2710102 and rs17236239) we genotyped 2051 additional families (572 multiplex and 1479 simplex). None of these variants were significantly associated with ASD after corrections for multiple testing. The analysis of Mendelian errors within each family did not indicate any segregating deletions. Nevertheless, a study of CNTNAP2 gene expression in brains of autistic patients and of normal controls, demonstrated altered expression in a subset of patients (p = 1.9 x10(-5)). Consequently, this study suggests that although CNTNAP2 dysregulation plays a role in some cases, its population contribution to autism susceptibility is limited.

Although variants in many genes have previously been shown to be associated with blood pressure (BP) levels, the molecular mechanism underlying these associations are mostly unknown. We identified a multi-allelic T-rich sequence (TRS) in the 3'UTR of ATP1B1 that varies in length and sequence composition (T22-27 and T12GT 3GT6). The 3'UTR of ATP1B1 contains 2 functional polyadenylation signals and the TRS is downstream of the proximal polyadenylation site (A2). Therefore, we hypothesized that alleles of this TRS might influence ATP1B1 expression by regulating alternative polyadenylation. In vitro, the T12GT 3GT6 allele increases polyadenylation at the A2 polyadenylation site as compared to the T23 allele. Consistent with our hypothesis, the relative abundance of the A2-polyadenylated ATP1B1 mRNA was higher in human kidneys with at least one copy of the T12GT 3GT6 allele than in those lacking this allele. The T12GT 3GT6 allele is also associated with higher systolic BP (beta = 3.3 mmHg, p = 0.014) and diastolic BP (beta = 2.4 mmHg, p = 0.003) in a European-American population. Therefore, we have identified a novel multi-allelic TRS in the 3'UTR of ATP1B1 that is associated with higher BP and may mediate its effect by regulating the polyadenylation of the ATP1B1 mRNA.

Hirschsprung disease (HSCR) genetics is a paradigm for the study and understanding of multigenic disorders. Association between Down syndrome and HSCR suggests that genetic factors that predispose to HSCR map to chromosome 21. To identify these additional factors, we performed a dose-dependent association study on chromosome 21 in Down syndrome patients with HSCR. Assessing 10,895 SNPs in 26 Caucasian cases and their parents led to identify two associated SNPs (rs2837770 and rs8134673) at chromosome-wide level. Those SNPs, which were located in intron 3 of the DSCAM gene within a 19 kb-linkage disequilibrium block region were in complete association and are consistent with DSCAM expression during enteric nervous system development. We replicated the association of HSCR with this region in an independent sample of 220 non-syndromic HSCR Caucasian patients and their parents. At last, we provide the functional rationale to the involvement of DSCAM by network analysis and assessment of SOX10 regulation. Our results reveal the involvement of DSCAM as a HSCR susceptibility locus, both in Down syndrome and HSCR isolated cases. This study further ascertains the chromosome-scan dose-dependent methodology used herein as a mean to map the genetic bases of other sub-phenotypes both in Down syndrome and other aneuploidies.

Microcephaly is a neurodevelopmental disorder causing significantly reduced cerebral cortex size. Many known microcephaly gene products localize to centrosomes, regulating cell fate and proliferation. Here, we identify and characterize a nuclear zinc finger protein, ZNF335/NIF-1, as a causative gene for severe microcephaly, small somatic size, and neonatal death. Znf335 null mice are embryonically lethal, and conditional knockout leads to severely reduced cortical size. RNA-interference and postmortem human studies show that ZNF335 is essential for neural progenitor self-renewal, neurogenesis, and neuronal differentiation. ZNF335 is a component of a vertebrate-specific, trithorax H3K4-methylation complex, directly regulating REST/NRSF, a master regulator of neural gene expression and cell fate, as well as other essential neural-specific genes. Our results reveal ZNF335 as an essential link between H3K4 complexes and REST/NRSF and provide the first direct genetic evidence that this pathway regulates human neurogenesis and neuronal differentiation.

By characterizing the geographic and functional spectrum of human genetic variation, the 1000 Genomes Project aims to build a resource to help to understand the genetic contribution to disease. Here we describe the genomes of 1,092 individuals from 14 populations, constructed using a combination of low-coverage whole-genome and exome sequencing. By developing methods to integrate information across several algorithms and diverse data sources, we provide a validated haplotype map of 38 million single nucleotide polymorphisms, 1.4 million short insertions and deletions, and more than 14,000 larger deletions. We show that individuals from different populations carry different profiles of rare and common variants, and that low-frequency variants show substantial geographic differentiation, which is further increased by the action of purifying selection. We show that evolutionary conservation and coding consequence are key determinants of the strength of purifying selection, that rare-variant load varies substantially across biological pathways, and that each individual contains hundreds of rare non-coding variants at conserved sites, such as motif-disrupting changes in transcription-factor-binding sites. This resource, which captures up to 98% of accessible single nucleotide polymorphisms at a frequency of 1% in related populations, enables analysis of common and low-frequency variants in individuals from diverse, including admixed, populations.

Coffee consumption is a model for addictive behavior. We performed a meta-analysis of genome-wide association studies (GWASs) on coffee intake from 8 Caucasian cohorts (N=18 176) and sought replication of our top findings in a further 7929 individuals. We also performed a gene expression analysis treating different cell lines with caffeine. Genome-wide significant association was observed for two single-nucleotide polymorphisms (SNPs) in the 15q24 region. The two SNPs rs2470893 and rs2472297 (P-values=1.6 × 10(-11) and 2.7 × 10(-11)), which were also in strong linkage disequilibrium (r(2)=0.7) with each other, lie in the 23-kb long commonly shared 5' flanking region between CYP1A1 and CYP1A2 genes. CYP1A1 was found to be downregulated in lymphoblastoid cell lines treated with caffeine. CYP1A1 is known to metabolize polycyclic aromatic hydrocarbons, which are important constituents of coffee, whereas CYP1A2 is involved in the primary metabolism of caffeine. Significant evidence of association was also detected at rs382140 (P-value=3.9 × 10(-09)) near NRCAM-a gene implicated in vulnerability to addiction, and at another independent hit rs6495122 (P-value=7.1 × 10(-09))-an SNP associated with blood pressure-in the 15q24 region near the gene ULK3, in the meta-analysis of discovery and replication cohorts. Our results from GWASs and expression analysis also strongly implicate CAB39L in coffee drinking. Pathway analysis of differentially expressed genes revealed significantly enriched ubiquitin proteasome (P-value=2.2 × 10(-05)) and Parkinson's disease pathways (P-value=3.6 × 10(-05)).

A major unanswered question in neuroscience is whether there exists genomic variability between individual neurons of the brain, contributing to functional diversity or to an unexplained burden of neurological disease. To address this question, we developed a method to amplify genomes of single neurons from human brains. Because recent reports suggest frequent LINE-1 (L1) retrotransposition in human brains, we performed genome-wide L1 insertion profiling of 300 single neurons from cerebral cortex and caudate nucleus of three normal individuals, recovering >80% of germline insertions from single neurons. While we find somatic L1 insertions, we estimate <0.6 unique somatic insertions per neuron, and most neurons lack detectable somatic insertions, suggesting that L1 is not a major generator of neuronal diversity in cortex and caudate. We then genotyped single cortical cells to characterize the mosaicism of a somatic AKT3 mutation identified in a child with hemimegalencephaly. Single-neuron sequencing allows systematic assessment of genomic diversity in the human brain.

Hirschsprung disease (HSCR, aganglionic megacolon) is a complex and heterogeneous disease with an incidence of 1 in 5000 live births. Despite the multifactorial determination of HSCR in the vast majority of cases, there is a monogenic subgroup for which private rare RET coding sequence mutations with high penetrance are found (45% of HSCR familial cases). An asymmetrical parental origin is observed for RET coding sequence mutations with a higher maternal inheritance. A parent-of-origin effect is usually assumed. Here we show that a differential reproductive rate for males and females also leads to an asymmetrical parental origin, which was never considered as a possible explanation till now. In the case of HSCR, we show a positive association between penetrance of the mutation and parental transmission asymmetry: no parental transmission asymmetry is observed in sporadic RET CDS mutation carrier cases for which penetrance of the mutation is low, whereas a parental transmission asymmetry is observed in affected sib-pairs for which penetrance of the mutation is higher. This allows us to conclude that the explanation for this parental asymmetry is that more severe mutations have resulted in a differential reproductive rate between male and female carriers.