Investigation of Type 2 Diabetes Risk Alleles Support CDKN2A/B,CDKAL1, and TCF7L2 As Susceptibility Genes in a Han Chinese Cohort

Background

Recent genome-wide association studies (GWASs) have reported several genetic variants to be reproducibly associated with type 2 diabetes. Additional variants have also been detected from a metaanalysis of three GWASs, performed in populations of European ancestry. In the present study, we evaluated the influence of 17 genetic variants from 15 candidate loci, identified in type 2 diabetes GWASs and the metaanalysis, in a Han Chinese cohort.

Methodology/Principal Findings

Selected type 2 diabetes–associated genetic variants were genotyped in 1,165 type 2 diabetic patients and 1,136 normoglycemic control individuals of Southern Han Chinese ancestry. The OR for risk of developing type 2 diabetes was calculated using a logistic regression model adjusted for age, sex, and BMI. Genotype-phenotype associations were tested using a multivariate linear regression model. Genetic variants in CDKN2A/B, CDKAL1, TCF7L2, TCF2, MC4R, and PPARG showed a nominal association with type 2 diabetes (P≤0.05), of whom the three first would stand correction for multiple testing: CDKN2A/B rs10811661, OR: 1.26 (1.12–1.43) P = 1.8*10⁻⁴; CDKAL1 rs10946398, OR: 1.23 (1.09–1.39); P = 7.1*10⁻⁴, and TCF7L2 rs7903146, OR: 1.61 (1.19–2.18) P = 2.3 * 10⁻³. Only nominal phenotype associations were observed, notably for rs8050136 in FTO and fasting plasma glucose (P = 0.002), postprandial plasma glucose (P = 0.002), and fasting C-peptide levels (P = 0.006) in the diabetic patients, and with BMI in controls (P = 0.033).

Conclusions/Significance

We have identified significant association between variants in CDKN2A/B, CDKAL1 and TCF7L2, and type 2 diabetes in a Han Chinese cohort, indicating these genes as strong candidates conferring susceptibility to type 2 diabetes across different ethnicities.     

Association of Genetic Variants with Isolated Fasting Hyperglycaemia and Isolated Postprandial Hyperglycaemia in a Han Chinese Population

Background

Though multiple single nucleotide polymorphisms (SNPs) associated with type 2 diabetes have been identified, the genetic bases of isolated fasting hyperglycaemia (IFH) and isolated postprandial hyperglycaemia (IPH) were still unclear. In present study, we aimed to investigate the association of genome-wide association study-validated genetic variants and IFH or IPH in Han Chinese.

Methods/Principal Findings

We genotyped 27 validated SNPs in 6,663 unrelated individuals comprising 341 IFH, 865 IPH, 1,203 combined fasting hyperglycaemia and postprandial hyperglycaemia, and 4,254 normal glycaemic subjects of Han ancestry. The distributions of genotype frequencies of FTO, CDKAL1 and GCKR were significant different between individuals with IFH and those with IPH (SNP(p_trend): rs8050136(0.0024), rs9939609(0.0049), rs7756992(0.0122), rs780094(0.0037)). Risk allele of FTO specifically increased the risk of IFH (rs8050136: OR 1.403 [95% CI 1.125–1.750], p = 0.0027; rs9939609: 1.398 [1.120–1.744], p = 0.0030). G allele of CDKAL1 specifically increased the risk of IPH (1.217 [1.092–1.355], p = 0.0004). G allele of GCKR increased the risk of IFH (1.167 [0.999–1.362], p = 0.0513), but decreased the risk of IPH (0.891 [0.801–0.991], p = 0.0331). In addition, TCF7L2 and KCNQ1 increased the risk of both IFH and IPH. When combined, each additional risk allele associated with IFH increased the risk for IFH by 1.246-fold (p<0.0001), while each additional risk allele associated with IPH increased the risk for IPH by 1.190-fold (p<0.0001).

Conclusion/Significance

Our results indicate that genotype distributions of variants from FTO, GCKR, CDKAL1 were different between IPH and IFH in Han Chinese. Variants of genes modulating insulin sensitivity (FTO, GCKR) contributed to the risk of IFH, while variants of genes related to beta cell function (CDKAL1) increase the risk of IPH.     

Whole-Genome Resequencing of Worldwide Wild and Domestic Sheep Elucidates Genetic Diversity,Introgression, and Agronomically Important Loci

Domestic sheep and their wild relatives harbor substantial genetic variants that can form the backbone of molecular breeding, but their genome landscapes remain understudied. Here, we present a comprehensive genome resource for wild ovine species, landraces and improved breeds of domestic sheep, comprising high-coverage (∼16.10×) whole genomes of 810 samples from 7 wild species and 158 diverse domestic populations. We detected, in total, ∼121.2 million single nucleotide polymorphisms, ∼61 million of which are novel. Some display significant (P < 0.001) differences in frequency between wild and domestic species, or are private to continent-wide or individual sheep populations. Retained or introgressed wild gene variants in domestic populations have contributed to local adaptation, such as the variation in the HBB associated with plateau adaptation. We identified novel and previously reported targets of selection on morphological and agronomic traits such as stature, horn, tail configuration, and wool fineness. We explored the genetic basis of wool fineness and unveiled a novel mutation (chr25: T7,068,586C) in the 3′-UTR of IRF2BP2 as plausible causal variant for fleece fiber diameter. We reconstructed prehistorical migrations from the Near Eastern domestication center to South-and-Southeast Asia and found two main waves of migrations across the Eurasian Steppe and the Iranian Plateau in the Early and Late Bronze Ages. Our findings refine our understanding of genome variation as shaped by continental migrations, introgression, adaptation, and selection of sheep.     

Cofactor Activity in Factor VIIIa of the Blood Clotting Pathway Is Stabilized by an Interdomain Bond between His281 and Ser524 Formed in Factor VIII

The factor VIII (FVIII) crystal structure suggests a possible bonding interaction of His²⁸¹ (A1 domain) with Ser⁵²⁴ (A2 domain), although the resolution of the structure (∼4 Å) does not firmly establish this bonding. To establish that side chains of these residues participate in an interdomain bond, we prepared and examined the functional properties of a residue swap variant (H281S/S524H) where His²⁸¹ and Ser⁵²⁴ residues were exchanged with one another and a disulfide-bridged variant (H281C/S524C) where the two residues were replaced with Cys. The latter variant showed efficient disulfide bonding of the A1 and A2 domains. The swap variant showed WT-like FVIII and FVIIIa stability, which were markedly reduced for H281A and S524A variants in an earlier study. The disulfide-bridged variant showed ∼20% increased FVIII stability, and FVIIIa did not decay during the time course measured. This variant also yielded 35% increased thrombin peak values compared with WT in a plasma-based thrombin generation assay. Binding analyses of H281S-A1/A3C1C2 dimer with S524H-A2 subunit yielded a near WT-like affinity value, whereas combining the variant dimer or A2 subunit with the WT complement yielded ∼5- and ∼10-fold reductions, respectively, in affinity. Other functional properties including thrombin generation potential, FIXa binding affinity, K_m for FX of FXase complexes, thrombin activation efficiency, and down-regulation by activated protein C showed similar results for the two variants compared with WT FVIII. These results indicate that the side chains of His²⁸¹ and Ser⁵²⁴ are in close proximity and contribute to a bonding interaction in FVIII that is retained in FVIIIa.     

Deletion of CDKAL1 affects mitochondrial ATP generation and first-phase insulin exocytosis

Background

A variant of the CDKAL1 gene was reported to be associated with type 2 diabetes and reduced insulin release in humans; however, the role of CDKAL1 in β cells is largely unknown. Therefore, to determine the role of CDKAL1 in insulin release from β cells, we studied insulin release profiles in CDKAL1 gene knockout (CDKAL1 KO) mice.

Principal Findings

Total internal reflection fluorescence imaging of CDKAL1 KO β cells showed that the number of fusion events during first-phase insulin release was reduced. However, there was no significant difference in the number of fusion events during second-phase release or high K⁺-induced release between WT and KO cells. CDKAL1 deletion resulted in a delayed and slow increase in cytosolic free Ca²⁺ concentration during high glucose stimulation. Patch-clamp experiments revealed that the responsiveness of ATP-sensitive K⁺ (K_ATP) channels to glucose was blunted in KO cells. In addition, glucose-induced ATP generation was impaired. Although CDKAL1 is homologous to cyclin-dependent kinase 5 (CDK5) regulatory subunit-associated protein 1, there was no difference in the kinase activity of CDK5 between WT and CDKAL1 KO islets.

Conclusions/Significance

We provide the first report describing the function of CDKAL1 in β cells. Our results indicate that CDKAL1 controls first-phase insulin exocytosis in β cells by facilitating ATP generation, K_ATP channel responsiveness and the subsequent activity of Ca²⁺ channels through pathways other than CDK5-mediated regulation.     

Deletion of CDKAL1 Affects High-Fat Diet–Induced Fat Accumulation and Glucose-Stimulated Insulin Secretion in Mice,Indicating Relevance to Diabetes

Background/Objective

The CDKAL1 gene is among the best-replicated susceptibility loci for type 2 diabetes, originally identified by genome-wide association studies in humans. To clarify a physiological importance of CDKAL1, we examined effects of a global Cdkal1-null mutation in mice and also evaluated the influence of a CDKAL1 risk allele on body mass index (BMI) in Japanese subjects.

Methods

In Cdkal1-deficient (Cdkal1 ^−/−) mice, we performed oral glucose tolerance test, insulin tolerance test, and perfusion experiments with and without high-fat feeding. Based on the findings in mice, we tested genetic association of CDKAL1 variants with BMI, as a measure of adiposity, and type 2 diabetes in Japanese.

Principal Findings

On a standard diet, Cdkal1 ^−/− mice were modestly lighter in weight than wild-type littermates without major alterations in glucose metabolism. On a high fat diet, Cdkal1 ^−/− mice showed significant reduction in fat accumulation (17% reduction in %intraabdominal fat, P = 0.023 vs. wild-type littermates) with less impaired insulin sensitivity at an early stage. High fat feeding did not potentiate insulin secretion in Cdkal1 ^−/− mice (1.0-fold), contrary to the results in wild-type littermates (1.6-fold, P<0.01). Inversely, at a later stage, Cdkal1 ^−/− mice showed more prominent impairment of insulin sensitivity and glucose tolerance. mRNA expression analysis indicated that Scd1 might function as a critical mediator of the altered metabolism in Cdkal1 ^−/− mice. In accordance with the findings in mice, a nominally significant (P<0.05) association between CDKAL1 rs4712523 and BMI was replicated in 2 Japanese general populations comprising 5,695 and 12,569 samples; the risk allele for type 2 diabetes was also associated with decreased BMI.

Conclusions

Cdkal1 gene deletion is accompanied by modestly impaired insulin secretion and longitudinal fluctuations in insulin sensitivity during high-fat feeding in mice. CDKAL1 may affect such compensatory mechanisms regulating glucose homeostasis through interaction with diet.     

Pooled Sequencing of 531 Genes in Inflammatory Bowel Disease Identifies an Associated Rare Variant in BTNL2 and Implicates Other Immune Related Genes

The contribution of rare coding sequence variants to genetic susceptibility in complex disorders is an important but unresolved question. Most studies thus far have investigated a limited number of genes from regions which contain common disease associated variants. Here we investigate this in inflammatory bowel disease by sequencing the exons and proximal promoters of 531 genes selected from both genome-wide association studies and pathway analysis in pooled DNA panels from 474 cases of Crohn’s disease and 480 controls. 80 variants with evidence of association in the sequencing experiment or with potential functional significance were selected for follow up genotyping in 6,507 IBD cases and 3,064 population controls. The top 5 disease associated variants were genotyped in an extension panel of 3,662 IBD cases and 3,639 controls, and tested for association in a combined analysis of 10,147 IBD cases and 7,008 controls. A rare coding variant p.G454C in the BTNL2 gene within the major histocompatibility complex was significantly associated with increased risk for IBD (p = 9.65x10⁻¹⁰, OR = 2.3[95% CI = 1.75–3.04]), but was independent of the known common associated CD and UC variants at this locus. Rare (<1%) and low frequency (1–5%) variants in 3 additional genes showed suggestive association (p<0.005) with either an increased risk (ARIH2 c.338-6C>T) or decreased risk (IL12B p.V298F, and NICN p.H191R) of IBD. These results provide additional insights into the involvement of the inhibition of T cell activation in the development of both sub-phenotypes of inflammatory bowel disease. We suggest that although rare coding variants may make a modest overall contribution to complex disease susceptibility, they can inform our understanding of the molecular pathways that contribute to pathogenesis.     

A Genetic Variant in 12q13, a Possible Risk Factor for Bipolar Disorder,Is Associated with Depressive State,Accounting for Stressful Life Events

Genome-wide association studies (GWASs) have identified a number of susceptibility genes for schizophrenia (SCZ) and bipolar disorder (BD). However, the identification of risk genes for major depressive disorder (MDD) has been unsuccessful because the etiology of MDD is more influenced by environmental factors; thus, gene–environment (G×E) interactions are important, such as interplay with stressful life events (SLEs). We assessed the G×E interactions and main effects of genes targeting depressive symptoms. Using a case–control design, 922 hospital staff members were evaluated for depressive symptoms according to Beck Depressive Inventory (BDI; “depression” and “control” groups were classified by scores of 10 in the BDI test), SLEs, and personality. A total of sixty-three genetic variants were selected on the basis of previous GWASs of MDD, SCZ, and BD as well as candidate-gene (SLC6A4, BDNF, DBH, and FKBP5) studies. Logistic regression analysis revealed a marginally significant interaction (genetic variant × SLE) at rs4523957 (P_uncorrected = 0.0034) with depression and a significant association of single nucleotide polymorphism identified from evidence of BD GWAS (rs7296288, downstream of DHH at 12q13.1) with depression as the main effect (P_uncorrected = 9.4×10⁻⁴, P_corrected = 0.0424). We also found that SLEs had a larger impact on depression (odds ratio∼3), as reported previously. These results suggest that DHH plays a possible role in depression etiology; however, variants from MDD or SCZ GWAS evidence or candidate genes showed no significant associations or minimal effects of interactions with SLEs on depression.     

Naturally Occurring Variants of the Dysglycemic Peptide Pancreastatin: DIFFERENTIAL POTENCIES FOR MULTIPLE CELLULAR FUNCTIONS AND STRUCTURE-FUNCTION CORRELATION*

Pancreastatin (PST), a chromogranin A-derived peptide, is a potent physiological inhibitor of glucose-induced insulin secretion. PST also triggers glycogenolysis in liver and reduces glucose uptake in adipocytes and hepatocytes. Here, we probed for genetic variations in PST sequence and identified two variants within its functionally important carboxyl terminus domain: E287K and G297S. To understand functional implications of these amino acid substitutions, we tested the effects of wild-type (PST-WT), PST-287K, and PST-297S peptides on various cellular processes/events. The rank order of efficacy to inhibit insulin-stimulated glucose uptake was: PST-297S > PST-287K > PST-WT. The PST peptides also displayed the same order of efficacy for enhancing intracellular nitric oxide and Ca²⁺ levels in various cell types. In addition, PST peptides activated gluconeogenic genes in the following order: PST-297S ≈ PST-287K > PST-WT. Consistent with these in vitro results, the common PST variant allele Ser-297 was associated with significantly higher (by ∼17 mg/dl, as compared with the wild-type Gly-297 allele) plasma glucose level in our study population (n = 410). Molecular modeling and molecular dynamics simulations predicted the following rank order of α-helical content: PST-297S > PST-287K > PST-WT. Corroboratively, circular dichroism analysis of PST peptides revealed significant differences in global structures (e.g. the order of propensity to form α-helix was: PST-297S ≈ PST-287K > PST-WT). This study provides a molecular basis for enhanced potencies/efficacies of human PST variants (likely to occur in ∼300 million people worldwide) and has quantitative implications for inter-individual variations in glucose/insulin homeostasis.     

Germline Heterozygous Variants in SEC23B Are Associated with Cowden Syndrome and Enriched in Apparently Sporadic Thyroid Cancer

Cancer-predisposing genes associated with inherited cancer syndromes help explain mechanisms of sporadic carcinogenesis and often inform normal development. Cowden syndrome (CS) is an autosomal-dominant disorder characterized by high lifetime risks of epithelial cancers, such that ∼50% of affected individuals are wild-type for known cancer-predisposing genes. Using whole-exome and Sanger sequencing of a multi-generation CS family affected by thyroid and other cancers, we identified a pathogenic missense heterozygous SEC23B variant (c.1781T>G [p.Val594Gly]) that segregates with the phenotype. We also found germline heterozygous SEC23B variants in 3/96 (3%) unrelated mutation-negative CS probands with thyroid cancer and in The Cancer Genome Atlas (TCGA), representing apparently sporadic cancers. We note that the TCGA thyroid cancer dataset is enriched with unique germline deleterious SEC23B variants associated with a significantly younger age of onset. SEC23B encodes Sec23 homolog B (S. cerevisiae), a component of coat protein complex II (COPII), which transports proteins from the endoplasmic reticulum (ER) to the Golgi apparatus. Interestingly, germline homozygous or compound-heterozygous SEC23B mutations cause an unrelated disorder, congenital dyserythropoietic anemia type II, and SEC23B-deficient mice suffer from secretory organ degeneration due to ER-stress-associated apoptosis. By characterizing the p.Val594Gly variant in a normal thyroid cell line, we show that it is a functional alteration that results in ER-stress-mediated cell-colony formation and survival, growth, and invasion, which reflect aspects of a cancer phenotype. Our findings suggest a different role for SEC23B, whereby germline heterozygous variants associate with cancer predisposition potentially mediated by ER stress “addiction.”     

Evolutionary Genomics of Structural Variation in Asian Rice (Oryza sativa) Domestication

Structural variants (SVs) are a largely unstudied feature of plant genome evolution, despite the fact that SVs contribute substantially to phenotypes. In this study, we discovered SVs across a population sample of 347 high-coverage, resequenced genomes of Asian rice (Oryza sativa) and its wild ancestor (O. rufipogon). In addition to this short-read data set, we also inferred SVs from whole-genome assemblies and long-read data. Comparisons among data sets revealed different features of genome variability. For example, genome alignment identified a large (∼4.3 Mb) inversion in indica rice varieties relative to japonica varieties, and long-read analyses suggest that ∼9% of genes from the outgroup (O. longistaminata) are hemizygous. We focused, however, on the resequencing sample to investigate the population genomics of SVs. Clustering analyses with SVs recapitulated the rice cultivar groups that were also inferred from SNPs. However, the site-frequency spectrum of each SV type—which included inversions, duplications, deletions, translocations, and mobile element insertions—was skewed toward lower frequency variants than synonymous SNPs, suggesting that SVs may be predominantly deleterious. Among transposable elements, SINE and mariner insertions were found at especially low frequency. We also used SVs to study domestication by contrasting between rice and O. rufipogon. Cultivated genomes contained ∼25% more derived SVs and mobile element insertions than O. rufipogon, indicating that SVs contribute to the cost of domestication in rice. Peaks of SV divergence were enriched for known domestication genes, but we also detected hundreds of genes gained and lost during domestication, some of which were enriched for traits of agronomic interest.     

Complete telomere-to-telomere de novo assembly of the Plasmodium falciparum genome through long-read (>11?kb), single molecule,real-time sequencing

The application of next-generation sequencing to estimate genetic diversity of Plasmodium falciparum, the most lethal malaria parasite, has proved challenging due to the skewed AT-richness [∼80.6% (A + T)] of its genome and the lack of technology to assemble highly polymorphic subtelomeric regions that contain clonally variant, multigene virulence families (Ex: var and rifin). To address this, we performed amplification-free, single molecule, real-time sequencing of P. falciparum genomic DNA and generated reads of average length 12 kb, with 50% of the reads between 15.5 and 50 kb in length. Next, using the Hierarchical Genome Assembly Process, we assembled the P. falciparum genome de novo and successfully compiled all 14 nuclear chromosomes telomere-to-telomere. We also accurately resolved centromeres [∼90–99% (A + T)] and subtelomeric regions and identified large insertions and duplications that add extra var and rifin genes to the genome, along with smaller structural variants such as homopolymer tract expansions. Overall, we show that amplification-free, long-read sequencing combined with de novo assembly overcomes major challenges inherent to studying the P. falciparum genome. Indeed, this technology may not only identify the polymorphic and repetitive subtelomeric sequences of parasite populations from endemic areas but may also evaluate structural variation linked to virulence, drug resistance and disease transmission.     

Diversity of Extended HLA-DRB1 Haplotypes in the Finnish Population

The Major Histocompatibility Complex (MHC, 6p21) codes for traditional HLA and other host response related genes. The polymorphic HLA-DRB1 gene in MHC Class II has been associated with several complex diseases. In this study we focus on MHC haplotype structures in the Finnish population. We explore the variability of extended HLA-DRB1 haplotypes in relation to the other traditional HLA genes and a selected group of MHC class III genes. A total of 150 healthy Finnish individuals were included in the study. Subjects were genotyped for HLA alleles (HLA-A, -B, -DRB1, -DQB1, and -DPB1). The polymorphism of TNF, LTA, C4, BTNL2 and HLA-DRA genes was studied with 74 SNPs (single nucleotide polymorphism). The C4A and C4B gene copy numbers and a 2-bp silencing insertion at exon 29 in C4A gene were analysed with quantitative genomic realtime-PCR. The allele frequencies for each locus were calculated and haplotypes were constructed using both the traditional HLA alleles and SNP blocks. The most frequent Finnish A∼B∼DR -haplotype, uncommon in elsewhere in Europe, was A*03∼B*35∼DRB1*01∶01. The second most common haplotype was a common European ancestral haplotype AH 8.1 (A*01∼B*08∼DRB1*03∶01). Extended haplotypes containing HLA-B, TNF block, C4 and HLA-DPB1 strongly increased the number of HLA-DRB1 haplotypes showing variability in the extended HLA-DRB1 haplotype structures. On the contrary, BTNL2 block and HLA-DQB1 were more conserved showing linkage with the HLA-DRB1 alleles. We show that the use of HLA-DRB1 haplotypes rather than single HLA-DRB1 alleles is advantageous when studying the polymorphisms and LD patters of the MHC region. For disease association studies the HLA-DRB1 haplotypes with various MHC markers allows us to cluster haplotypes with functionally important gene variants such as C4 deficiency and cytokines TNF and LTA, and provides hypotheses for further assessment. Our study corroborates the importance of studying population-specific MHC haplotypes.     

Hypoxia Adaptations in the Grey Wolf (Canis lupus chanco) from Qinghai-Tibet Plateau

The Tibetan grey wolf (Canis lupus chanco) occupies habitats on the Qinghai-Tibet Plateau, a high altitude (>3000 m) environment where low oxygen tension exerts unique selection pressure on individuals to adapt to hypoxic conditions. To identify genes involved in hypoxia adaptation, we generated complete genome sequences of nine Chinese wolves from high and low altitude populations at an average coverage of 25× coverage. We found that, beginning about 55,000 years ago, the highland Tibetan grey wolf suffered a more substantial population decline than lowland wolves. Positively selected hypoxia-related genes in highland wolves are enriched in the HIF signaling pathway (P = 1.57E-6), ATP binding (P = 5.62E-5), and response to an oxygen-containing compound (P≤5.30E-4). Of these positively selected hypoxia-related genes, three genes (EPAS1, ANGPT1, and RYR2) had at least one specific fixed non-synonymous SNP in highland wolves based on the nine genome data. Our re-sequencing studies on a large panel of individuals showed a frequency difference greater than 58% between highland and lowland wolves for these specific fixed non-synonymous SNPs and a high degree of LD surrounding the three genes, which imply strong selection. Past studies have shown that EPAS1 and ANGPT1 are important in the response to hypoxic stress, and RYR2 is involved in heart function. These three genes also exhibited significant signals of natural selection in high altitude human populations, which suggest similar evolutionary constraints on natural selection in wolves and humans of the Qinghai-Tibet Plateau.     

Association and Mutation Analyses of 16p11.2 Autism Candidate Genes

Background

Autism is a complex childhood neurodevelopmental disorder with a strong genetic basis. Microdeletion or duplication of a ∼500–700-kb genomic rearrangement on 16p11.2 that contains 24 genes represents the second most frequent chromosomal disorder associated with autism. The role of common and rare 16p11.2 sequence variants in autism etiology is unknown.

Methodology/Principal Findings

To identify common 16p11.2 variants with a potential role in autism, we performed association studies using existing data generated from three microarray platforms: Affymetrix 5.0 (777 families), Illumina 550 K (943 families), and Affymetrix 500 K (60 families). No common variants were identified that were significantly associated with autism. To look for rare variants, we performed resequencing of coding and promoter regions for eight candidate genes selected based on their known expression patterns and functions. In total, we identified 26 novel variants in autism: 13 exonic (nine non-synonymous, three synonymous, and one untranslated region) and 13 promoter variants. We found a significant association between autism and a coding variant in the seizure-related gene SEZ6L2 (12/1106 autism vs. 3/1161 controls; p = 0.018). Sez6l2 expression in mouse embryos was restricted to the spinal cord and brain. SEZ6L2 expression in human fetal brain was highest in post-mitotic cortical layers, hippocampus, amygdala, and thalamus. Association analysis of SEZ6L2 in an independent sample set failed to replicate our initial findings.

Conclusions/Significance

We have identified sequence variation in at least one candidate gene in 16p11.2 that may represent a novel genetic risk factor for autism. However, further studies are required to substantiate these preliminary findings.