Generation of functional insulin-producing cells in the gut by Foxo1 ablation

Restoration of regulated insulin secretion is the ultimate goal of therapy for type 1 diabetes. Here, we show that, unexpectedly, somatic ablation of Foxo1 in Neurog3(+) enteroendocrine progenitor cells gives rise to gut insulin-positive (Ins(+)) cells that express markers of mature β cells and secrete bioactive insulin as well as C-peptide in response to glucose and sulfonylureas. Lineage tracing experiments showed that gut Ins(+) cells arise cell autonomously from Foxo1-deficient cells. Inducible Foxo1 ablation in adult mice also resulted in the generation of gut Ins(+) cells. Following ablation by the β-cell toxin streptozotocin, gut Ins(+) cells regenerate and produce insulin, reversing hyperglycemia in mice. The data indicate that Neurog3(+) enteroendocrine progenitors require active Foxo1 to prevent differentiation into Ins(+) cells. Foxo1 ablation in gut epithelium may provide an approach to restore insulin production in type 1 diabetes.     

Genome-wide copy number variation study associates metabotropic glutamate receptor gene networks with attention deficit hyperactivity disorder

Attention deficit hyperactivity disorder (ADHD) is a common, heritable neuropsychiatric disorder of unknown etiology. We performed a whole-genome copy number variation (CNV) study on 1,013 cases with ADHD and 4,105 healthy children of European ancestry using 550,000 SNPs. We evaluated statistically significant findings in multiple independent cohorts, with a total of 2,493 cases with ADHD and 9,222 controls of European ancestry, using matched platforms. CNVs affecting metabotropic glutamate receptor genes were enriched across all cohorts (P = 2.1 × 10(-9)). We saw GRM5 (encoding glutamate receptor, metabotropic 5) deletions in ten cases and one control (P = 1.36 × 10(-6)). We saw GRM7 deletions in six cases, and we saw GRM8 deletions in eight cases and no controls. GRM1 was duplicated in eight cases. We experimentally validated the observed variants using quantitative RT-PCR. A gene network analysis showed that genes interacting with the genes in the GRM family are enriched for CNVs in ～10% of the cases (P = 4.38 × 10(-10)) after correction for occurrence in the controls. We identified rare recurrent CNVs affecting glutamatergic neurotransmission genes that were overrepresented in multiple ADHD cohorts.     

Genome-wide association study of flowering time and grain yield traits in a worldwide collection of rice germplasm

A high-density haplotype map recently enabled a genome-wide association study (GWAS) in a population of indica subspecies of Chinese rice landraces. Here we extend this methodology to a larger and more diverse sample of 950 worldwide rice varieties, including the Oryza sativa indica and Oryza sativa japonica subspecies, to perform an additional GWAS. We identified a total of 32 new loci associated with flowering time and with ten grain-related traits, indicating that the larger sample increased the power to detect trait-associated variants using GWAS. To characterize various alleles and complex genetic variation, we developed an analytical framework for haplotype-based de novo assembly of the low-coverage sequencing data in rice. We identified candidate genes for 18 associated loci through detailed annotation. This study shows that the integrated approach of sequence-based GWAS and functional genome annotation has the potential to match complex traits to their causal polymorphisms in rice.     

Estimating the proportion of variation in susceptibility to schizophrenia captured by common SNPs

Schizophrenia is a complex disorder caused by both genetic and environmental factors. Using 9,087 affected individuals, 12,171 controls and 915,354 imputed SNPs from the Schizophrenia Psychiatric Genome-Wide Association Study (GWAS) Consortium (PGC-SCZ), we estimate that 23% (s.e. = 1%) of variation in liability to schizophrenia is captured by SNPs. We show that a substantial proportion of this variation must be the result of common causal variants, that the variance explained by each chromosome is linearly related to its length (r = 0.89, P = 2.6 × 10(-8)), that the genetic basis of schizophrenia is the same in males and females, and that a disproportionate proportion of variation is attributable to a set of 2,725 genes expressed in the central nervous system (CNS; P = 7.6 × 10(-8)). These results are consistent with a polygenic genetic architecture and imply more individual SNP associations will be detected for this disease as sample size increases.     

Common variants near TARDBP and EGR2 are associated with susceptibility to Ewing sarcoma

Ewing sarcoma, a pediatric tumor characterized by EWSR1-ETS fusions, is predominantly observed in populations of European ancestry. We performed a genome-wide association study (GWAS) of 401 French individuals with Ewing sarcoma, 684 unaffected French individuals and 3,668 unaffected individuals of European descent and living in the United States. We identified candidate risk loci at 1p36.22, 10q21 and 15q15. We replicated these loci in two independent sets of cases and controls. Joint analysis identified associations with rs9430161 (P = 1.4 × 10(-20); odds ratio (OR) = 2.2) located 25 kb upstream of TARDBP, rs224278 (P = 4.0 × 10(-17); OR = 1.7) located 5 kb upstream of EGR2 and, to a lesser extent, rs4924410 at 15q15 (P = 6.6 × 10(-9); OR = 1.5). The major risk haplotypes were less prevalent in Africans, suggesting that these loci could contribute to geographical differences in Ewing sarcoma incidence. TARDBP shares structural similarities with EWSR1 and FUS, which encode RNA binding proteins, and EGR2 is a target gene of EWSR1-ETS. Variants at these loci were associated with expression levels of TARDBP, ADO (encoding cysteamine dioxygenase) and EGR2.     

A genome-wide association study in Han Chinese identifies new susceptibility loci for ankylosing spondylitis

To identify susceptibility loci for ankylosing spondylitis, we performed a two-stage genome-wide association study in Han Chinese. In the discovery stage, we analyzed 1,356,350 autosomal SNPs in 1,837 individuals with ankylosing spondylitis and 4,231 controls; in the validation stage, we analyzed 30 suggestive SNPs in an additional 2,100 affected individuals and 3,496 controls. We identified two new susceptibility loci between EDIL3 and HAPLN1 at 5q14.3 (rs4552569; P = 8.77 × 10(-10)) and within ANO6 at 12q12 (rs17095830; P = 1.63 × 10(-8)). We also confirmed previously reported associations in Europeans within the major histocompatibility complex (MHC) region (top SNP, rs13202464; P < 5 × 10(-324)) and at 2p15 (rs10865331; P = 1.98 × 10(-8)). We show that rs13202464 within the MHC region mainly represents the risk effect of HLA-B*27 variants (including HLA-B*2704, HLA-B*2705 and HLA-B*2715) in Chinese. The two newly discovered loci implicate genes related to bone formation and cartilage development, suggesting their potential involvement in the etiology of ankylosing spondylitis.     

Extremely low-coverage sequencing and imputation increases power for genome-wide association studies

Genome-wide association studies (GWAS) have proven to be a powerful method to identify common genetic variants contributing to susceptibility to common diseases. Here, we show that extremely low-coverage sequencing (0.1-0.5×) captures almost as much of the common (>5%) and low-frequency (1-5%) variation across the genome as SNP arrays. As an empirical demonstration, we show that genome-wide SNP genotypes can be inferred at a mean r(2) of 0.71 using off-target data (0.24× average coverage) in a whole-exome study of 909 samples. Using both simulated and real exome-sequencing data sets, we show that association statistics obtained using extremely low-coverage sequencing data attain similar P values at known associated variants as data from genotyping arrays, without an excess of false positives. Within the context of reductions in sample preparation and sequencing costs, funds invested in extremely low-coverage sequencing can yield several times the effective sample size of GWAS based on SNP array data and a commensurate increase in statistical power.     

Dominant missense mutations in ABCC9 cause Cantú syndrome

Cantú syndrome is characterized by congenital hypertrichosis, distinctive facial features, osteochondrodysplasia and cardiac defects. By using family-based exome sequencing, we identified a de novo mutation in ABCC9. Subsequently, we discovered novel dominant missense mutations in ABCC9 in 14 of the 16 individuals with Cantú syndrome examined. The ABCC9 protein is part of an ATP-dependent potassium (K(ATP)) channel that couples the metabolic state of a cell with its electrical activity. All mutations altered amino acids in or close to the transmembrane domains of ABCC9. Using electrophysiological measurements, we show that mutations in ABCC9 reduce the ATP-mediated potassium channel inhibition, resulting in channel opening. Moreover, similarities between the phenotype of individuals with Cantú syndrome and side effects from the K(ATP) channel agonist minoxidil indicate that the mutations in ABCC9 result in channel opening. Given the availability of ABCC9 antagonists, our findings may have direct implications for the treatment of individuals with Cantú syndrome.     

CCDC103 mutations cause primary ciliary dyskinesia by disrupting assembly of ciliary dynein arms

Cilia are essential for fertilization, respiratory clearance, cerebrospinal fluid circulation and establishing laterality. Cilia motility defects cause primary ciliary dyskinesia (PCD, MIM244400), a disorder affecting 1:15,000-30,000 births. Cilia motility requires the assembly of multisubunit dynein arms that drive ciliary bending. Despite progress in understanding the genetic basis of PCD, mutations remain to be identified for several PCD-linked loci. Here we show that the zebrafish cilia paralysis mutant schmalhans (smh(tn222)) encodes the coiled-coil domain containing 103 protein (Ccdc103), a foxj1a-regulated gene product. Screening 146 unrelated PCD families identified individuals in six families with reduced outer dynein arms who carried mutations in CCDC103. Dynein arm assembly in smh mutant zebrafish was rescued by wild-type but not mutant human CCDC103. Chlamydomonas Ccdc103/Pr46b functions as a tightly bound, axoneme-associated protein. These results identify Ccdc103 as a dynein arm attachment factor that causes primary ciliary dyskinesia when mutated.     

Exome sequencing of liver fluke-associated cholangiocarcinoma

Opisthorchis viverrini-related cholangiocarcinoma (CCA), a fatal bile duct cancer, is a major public health concern in areas endemic for this parasite. We report here whole-exome sequencing of eight O. viverrini-related tumors and matched normal tissue. We identified and validated 206 somatic mutations in 187 genes using Sanger sequencing and selected 15 genes for mutation prevalence screening in an additional 46 individuals with CCA (cases). In addition to the known cancer-related genes TP53 (mutated in 44.4% of cases), KRAS (16.7%) and SMAD4 (16.7%), we identified somatic mutations in 10 newly implicated genes in 14.8-3.7% of cases. These included inactivating mutations in MLL3 (in 14.8% of cases), ROBO2 (9.3%), RNF43 (9.3%) and PEG3 (5.6%), and activating mutations in the GNAS oncogene (9.3%). These genes have functions that can be broadly grouped into three biological classes: (i) deactivation of histone modifiers, (ii) activation of G protein signaling and (iii) loss of genome stability. This study provides insight into the mutational landscape contributing to O. viverrini-related CCA.