Contributions of common genetic variants to risk of schizophrenia among individuals of African and Latino ancestry.

Schizophrenia is a common, chronic and debilitating neuropsychiatric syndrome affecting tens of millions of individuals worldwide. While rare genetic variants play a role in the etiology of schizophrenia, most of the currently explained liability is within common variation, suggesting that variation predating the human diaspora out of Africa harbors a large fraction of the common variant attributable heritability. However, common variant association studies in schizophrenia have concentrated mainly on cohorts of European descent. We describe genome-wide association studies of 6152 cases and 3918 controls of admixed African ancestry, and of 1234 cases and 3090 controls of Latino ancestry, representing the largest such study in these populations to date. Combining results from the samples with African ancestry with summary statistics from the Psychiatric Genomics Consortium (PGC) study of schizophrenia yielded seven newly genome-wide significant loci, and we identified an additional eight loci by incorporating the results from samples with Latino ancestry. Leveraging population differences in patterns of linkage disequilibrium, we achieve improved fine-mapping resolution at 22 previously reported and 4 newly significant loci. Polygenic risk score profiling revealed improved prediction based on trans-ancestry meta-analysis results for admixed African (Nagelkerke's R2 = 0.032; liability R2 = 0.017; P < 10-52), Latino (Nagelkerke's R2 = 0.089; liability R2 = 0.021; P < 10-58), and European individuals (Nagelkerke's R2 = 0.089; liability R2 = 0.037; P < 10-113), further highlighting the advantages of incorporating data from diverse human populations.

De novo mutations in schizophrenia implicate synaptic networks.

Inherited alleles account for most of the genetic risk for schizophrenia. However, new (de novo) mutations, in the form of large chromosomal copy number changes, occur in a small fraction of cases and disproportionally disrupt genes encoding postsynaptic proteins. Here we show that small de novo mutations, affecting one or a few nucleotides, are overrepresented among glutamatergic postsynaptic proteins comprising activity-regulated cytoskeleton-associated protein (ARC) and N-methyl-d-aspartate receptor (NMDAR) complexes. Mutations are additionally enriched in proteins that interact with these complexes to modulate synaptic strength, namely proteins regulating actin filament dynamics and those whose messenger RNAs are targets of fragile X mental retardation protein (FMRP). Genes affected by mutations in schizophrenia overlap those mutated in autism and intellectual disability, as do mutation-enriched synaptic pathways. Aligning our findings with a parallel case-control study, we demonstrate reproducible insights into aetiological mechanisms for schizophrenia and reveal pathophysiology shared with other neurodevelopmental disorders.

A polygenic burden of rare disruptive mutations in schizophrenia.

Schizophrenia is a common disease with a complex aetiology, probably involving multiple and heterogeneous genetic factors. Here, by analysing the exome sequences of 2,536 schizophrenia cases and 2,543 controls, we demonstrate a polygenic burden primarily arising from rare (less than 1 in 10,000), disruptive mutations distributed across many genes. Particularly enriched gene sets include the voltage-gated calcium ion channel and the signalling complex formed by the activity-regulated cytoskeleton-associated scaffold protein (ARC) of the postsynaptic density, sets previously implicated by genome-wide association and copy-number variation studies. Similar to reports in autism, targets of the fragile X mental retardation protein (FMRP, product of FMR1) are enriched for case mutations. No individual gene-based test achieves significance after correction for multiple testing and we do not detect any alleles of moderately low frequency (approximately 0.5 to 1 per cent) and moderately large effect. Taken together, these data suggest that population-based exome sequencing can discover risk alleles and complements established gene-mapping paradigms in neuropsychiatric disease.

Clonal hematopoiesis and blood-cancer risk inferred from blood DNA sequence.

BACKGROUND: Cancers arise from multiple acquired mutations, which presumably occur over many years. Early stages in cancer development might be present years before cancers become clinically apparent. METHODS: We analyzed data from whole-exome sequencing of DNA in peripheral-blood cells from 12,380 persons, unselected for cancer or hematologic phenotypes. We identified somatic mutations on the basis of unusual allelic fractions. We used data from Swedish national patient registers to follow health outcomes for 2 to 7 years after DNA sampling. RESULTS: Clonal hematopoiesis with somatic mutations was observed in 10% of persons older than 65 years of age but in only 1% of those younger than 50 years of age. Detectable clonal expansions most frequently involved somatic mutations in three genes (DNMT3A, ASXL1, and TET2) that have previously been implicated in hematologic cancers. Clonal hematopoiesis was a strong risk factor for subsequent hematologic cancer (hazard ratio, 12.9; 95% confidence interval, 5.8 to 28.7). Approximately 42% of hematologic cancers in this cohort arose in persons who had clonality at the time of DNA sampling, more than 6 months before a first diagnosis of cancer. Analysis of bone marrow-biopsy specimens obtained from two patients at the time of diagnosis of acute myeloid leukemia revealed that their cancers arose from the earlier clones. CONCLUSIONS: Clonal hematopoiesis with somatic mutations is readily detected by means of DNA sequencing, is increasingly common as people age, and is associated with increased risks of hematologic cancer and death. A subset of the genes that are mutated in patients with myeloid cancers is frequently mutated in apparently healthy persons; these mutations may represent characteristic early events in the development of hematologic cancers. (Funded by the National Human Genome Research Institute and others.).

Genome-wide significant locus for Research Diagnostic Criteria Schizoaffective Disorder Bipolar type.

Studies have suggested that Research Diagnostic Criteria for Schizoaffective Disorder Bipolar type (RDC-SABP) might identify a more genetically homogenous subgroup of bipolar disorder. Aiming to identify loci associated with RDC-SABP, we have performed a replication study using independent RDC-SABP cases (n = 144) and controls (n = 6,559), focusing on the 10 loci that reached a p-value <10-5 for RDC-SABP in the Wellcome Trust Case Control Consortium (WTCCC) bipolar disorder sample. Combining the WTCCC and replication datasets by meta-analysis (combined RDC-SABP, n = 423, controls, n = 9,494), we observed genome-wide significant association at one SNP, rs2352974, located within the intron of the gene TRAIP on chromosome 3p21.31 (p-value, 4.37 × 10-8 ). This locus did not reach genome-wide significance in bipolar disorder or schizophrenia large Psychiatric Genomic Consortium datasets, suggesting that it may represent a relatively specific genetic risk for the bipolar subtype of schizoaffective disorder.

Mutation intolerant genes and targets of FMRP are enriched for nonsynonymous alleles in schizophrenia.

Risk of schizophrenia is conferred by alleles occurring across the full spectrum of frequencies from common SNPs of weak effect through to ultra rare alleles, some of which may be moderately to highly penetrant. Previous studies have suggested that some of the risk of schizophrenia is attributable to uncommon alleles represented on Illumina exome arrays. Here, we present the largest study of exomic variation in schizophrenia to date, using samples from the United Kingdom and Sweden (10,011 schizophrenia cases and 13,791 controls). Single variants, genes, and gene sets were analyzed for association with schizophrenia. No single variant or gene reached genome-wide significance. Among candidate gene sets, we found significant enrichment for rare alleles (minor allele frequency [MAF] < 0.001) in genes intolerant of loss-of-function (LoF) variation and in genes whose messenger RNAs bind to fragile X mental retardation protein (FMRP). We further delineate the genetic architecture of schizophrenia by excluding a role for uncommon exomic variants (0.01 ≤ MAF ≥ 0.001) that confer a relatively large effect (odds ratio [OR] > 4). We also show risk alleles within this frequency range exist, but confer smaller effects and should be identified by larger studies.

De novo CNVs in bipolar affective disorder and schizophrenia.

An increased rate of de novo copy number variants (CNVs) has been found in schizophrenia (SZ), autism and developmental delay. An increased rate has also been reported in bipolar affective disorder (BD). Here, in a larger BD sample, we aimed to replicate these findings and compare de novo CNVs between SZ and BD. We used Illumina microarrays to genotype 368 BD probands, 76 SZ probands and all their parents. Copy number variants were called by PennCNV and filtered for frequency (<1%) and size (>10 kb). Putative de novo CNVs were validated with the z-score algorithm, manual inspection of log R ratios (LRR) and qPCR probes. We found 15 de novo CNVs in BD (4.1% rate) and 6 in SZ (7.9% rate). Combining results with previous studies and using a cut-off of >100 kb, the rate of de novo CNVs in BD was intermediate between controls and SZ: 1.5% in controls, 2.2% in BD and 4.3% in SZ. Only the differences between SZ and BD and SZ and controls were significant. The median size of de novo CNVs in BD (448 kb) was also intermediate between SZ (613 kb) and controls (338 kb), but only the comparison between SZ and controls was significant. Only one de novo CNV in BD was in a confirmed SZ locus (16p11.2). Sporadic or early onset cases were not more likely to have de novo CNVs. We conclude that de novo CNVs play a smaller role in BD compared with SZ. Patients with a positive family history can also harbour de novo mutations.

CNV analysis in a large schizophrenia sample implicates deletions at 16p12.1 and SLC1A1 and duplications at 1p36.33 and CGNL1.

Large and rare copy number variants (CNVs) at several loci have been shown to increase risk for schizophrenia. Aiming to discover novel susceptibility CNV loci, we analyzed 6882 cases and 11 255 controls genotyped on Illumina arrays, most of which have not been used for this purpose before. We identified genes enriched for rare exonic CNVs among cases, and then attempted to replicate the findings in additional 14 568 cases and 15 274 controls. In a combined analysis of all samples, 12 distinct loci were enriched among cases with nominal levels of significance (P < 0.05); however, none would survive correction for multiple testing. These loci include recurrent deletions at 16p12.1, a locus previously associated with neurodevelopmental disorders (P = 0.0084 in the discovery sample and P = 0.023 in the replication sample). Other plausible candidates include non-recurrent deletions at the glutamate transporter gene SLC1A1, a CNV locus recently suggested to be involved in schizophrenia through linkage analysis, and duplications at 1p36.33 and CGNL1. A burden analysis of large (>500 kb), rare CNVs showed a 1.2% excess in cases after excluding known schizophrenia-associated loci, suggesting that additional susceptibility loci exist. However, even larger samples are required for their discovery.

Structural basis of clade-specific HIV-1 neutralization by humanized anti-V3 monoclonal antibody KD-247.

Humanized monoclonal antibody KD-247 targets the Gly(312)-Pro(313)-Gly(314)-Arg(315) arch of the third hypervariable (V3) loop of the HIV-1 surface glycoprotein. It potently neutralizes many HIV-1 clade B isolates, but not of other clades. To understand the molecular basis of this specificity, we solved a high-resolution (1.55 Å) crystal structure of the KD-247 antigen binding fragment and examined the potential interactions with various V3 loop targets. Unlike most antibodies, KD-247 appears to interact with its target primarily through light chain residues. Several of these interactions involve Arg(315) of the V3 loop. To evaluate the role of light chain residues in the recognition of the V3 loop, we generated 20 variants of KD-247 single-chain variable fragments with mutations in the antigen-binding site. Purified proteins were assessed for V3 loop binding using AlphaScreen technology and for HIV-1 neutralization. Our data revealed that recognition of the clade-specificity defining residue Arg(315) of the V3 loop is based on a network of interactions that involve Tyr(L32), Tyr(L92), and Asn(L27d) that directly interact with Arg(315), thus elucidating the molecular interactions of KD-247 with its V3 loop target.

24/7 Neurocritical Care Nurse Practitioner Coverage Reduced Door-to-Needle Time in Stroke Patients Treated with Tissue Plasminogen Activator.

BACKGROUND: Stroke centers with limited on-site neurovascular physician coverage may experience delays in acute stroke treatment. We sought to assess the impact of providing 24/7 neurocritical care acute care nurse practitioner (ACNP) "stroke code" first responder coverage on treatment delays in acute stroke patients who received tissue plasminogen activator (tPA). METHODS: Consecutive acute ischemic stroke patients treated with intravenous tPA at a primary stroke center on Oahu between 2009 and 2014were retrospectively studied. 24/7 ACNP stroke code coverage (intervention) was introduced on July 1, 2011. The tPA utilization, door-to-needle (DTN) time, imaging-to-needle (ITN) time, and independent ambulation at hospital discharge were compared between the preintervention period (24 months) and the postintervention period (33 months). RESULTS: We studied 166 stroke code patients who were treated with intravenous tPA, 44 of whom were treated during the preintervention period and 122 of whom were treated during the postintervention period. After the intervention, the median DTN time was reduced from 53 minutes (interquartile range [IQR] 45-73) to 45 minutes (IQR 35-58) (P = .001), and the median ITN time was reduced from 36 minutes (IQR 28-64) to 21 minutes (IQR 16-31) (P < .0001). Compliance with the 60-minute target DTN improved from 61.4% (27 of 44 patients) in the preintervention period to 81.2% (99 of 122 patients) in the postintervention period (P = .004). The tPA treatment rates were similar between the preintervention and postintervention periods (P = .60). CONCLUSIONS: Addition of 24/7 on-site neurocritical care ACNP first responder coverage for acute stroke code significantly reduced the DTN time among acute stroke patients treated with tPA.

Discovery and statistical genotyping of copy-number variation from whole-exome sequencing depth.

Sequencing of gene-coding regions (the exome) is increasingly used for studying human disease, for which copy-number variants (CNVs) are a critical genetic component. However, detecting copy number from exome sequencing is challenging because of the noncontiguous nature of the captured exons. This is compounded by the complex relationship between read depth and copy number; this results from biases in targeted genomic hybridization, sequence factors such as GC content, and batching of samples during collection and sequencing. We present a statistical tool (exome hidden Markov model [XHMM]) that uses principal-component analysis (PCA) to normalize exome read depth and a hidden Markov model (HMM) to discover exon-resolution CNV and genotype variation across samples. We evaluate performance on 90 schizophrenia trios and 1,017 case-control samples. XHMM detects a median of two rare (<1%) CNVs per individual (one deletion and one duplication) and has 79% sensitivity to similarly rare CNVs overlapping three or more exons discovered with microarrays. With sensitivity similar to state-of-the-art methods, XHMM achieves higher specificity by assigning quality metrics to the CNV calls to filter out bad ones, as well as to statistically genotype the discovered CNV in all individuals, yielding a trio call set with Mendelian-inheritance properties highly consistent with expectation. We also show that XHMM breakpoint quality scores enable researchers to explicitly search for novel classes of structural variation. For example, we apply XHMM to extract those CNVs that are highly likely to disrupt (delete or duplicate) only a portion of a gene.

Analysis of copy number variations at 15 schizophrenia-associated loci.

BACKGROUND: A number of copy number variants (CNVs) have been suggested as susceptibility factors for schizophrenia. For some of these the data remain equivocal, and the frequency in individuals with schizophrenia is uncertain. AIMS: To determine the contribution of CNVs at 15 schizophrenia-associated loci (a) using a large new data-set of patients with schizophrenia (n = 6882) and controls (n = 6316), and (b) combining our results with those from previous studies. METHOD: We used Illumina microarrays to analyse our data. Analyses were restricted to 520 766 probes common to all arrays used in the different data-sets. RESULTS: We found higher rates in participants with schizophrenia than in controls for 13 of the 15 previously implicated CNVs. Six were nominally significantly associated (P<0.05) in this new data-set: deletions at 1q21.1, NRXN1, 15q11.2 and 22q11.2 and duplications at 16p11.2 and the Angelman/Prader-Willi Syndrome (AS/PWS) region. All eight AS/PWS duplications in patients were of maternal origin. When combined with published data, 11 of the 15 loci showed highly significant evidence for association with schizophrenia (P<4.1×10(-4)). CONCLUSIONS: We strengthen the support for the majority of the previously implicated CNVs in schizophrenia. About 2.5% of patients with schizophrenia and 0.9% of controls carry a large, detectable CNV at one of these loci. Routine CNV screening may be clinically appropriate given the high rate of known deleterious mutations in the disorder and the comorbidity associated with these heritable mutations.

Cholesterol metabolism is required for intracellular hedgehog signal transduction in vivo.

We describe the rudolph mouse, a mutant with striking defects in both central nervous system and skeletal development. Rudolph is an allele of the cholesterol biosynthetic enzyme, hydroxysteroid (17-beta) dehydrogenase 7, which is an intriguing finding given the recent implication of oxysterols in mediating intracellular Hedgehog (Hh) signaling. We see an abnormal sterol profile and decreased Hh target gene induction in the rudolph mutant, both in vivo and in vitro. Reduced Hh signaling has been proposed to contribute to the phenotypes of congenital diseases of cholesterol metabolism. Recent in vitro and pharmacological data also indicate a requirement for intracellular cholesterol synthesis for proper regulation of Hh activity via Smoothened. The data presented here are the first in vivo genetic evidence supporting both of these hypotheses, revealing a role for embryonic cholesterol metabolism in both CNS development and normal Hh signaling.

Lethal skeletal dysplasia in mice and humans lacking the golgin GMAP-210.

BACKGROUND: Establishing the genetic basis of phenotypes such as skeletal dysplasia in model organisms can provide insights into biologic processes and their role in human disease. METHODS: We screened mutagenized mice and observed a neonatal lethal skeletal dysplasia with an autosomal recessive pattern of inheritance. Through genetic mapping and positional cloning, we identified the causative mutation. RESULTS: Affected mice had a nonsense mutation in the thyroid hormone receptor interactor 11 gene (Trip11), which encodes the Golgi microtubule-associated protein 210 (GMAP-210); the affected mice lacked this protein. Golgi architecture was disturbed in multiple tissues, including cartilage. Skeletal development was severely impaired, with chondrocytes showing swelling and stress in the endoplasmic reticulum, abnormal cellular differentiation, and increased cell death. Golgi-mediated glycosylation events were altered in fibroblasts and chondrocytes lacking GMAP-210, and these chondrocytes had intracellular accumulation of perlecan, an extracellular matrix protein, but not of type II collagen or aggrecan, two other extracellular matrix proteins. The similarities between the skeletal and cellular phenotypes in these mice and those in patients with achondrogenesis type 1A, a neonatal lethal form of skeletal dysplasia in humans, suggested that achondrogenesis type 1A may be caused by GMAP-210 deficiency. Sequence analysis revealed loss-of-function mutations in the 10 unrelated patients with achondrogenesis type 1A whom we studied. CONCLUSIONS: GMAP-210 is required for the efficient glycosylation and cellular transport of multiple proteins. The identification of a mutation affecting GMAP-210 in mice, and then in humans, as the cause of a lethal skeletal dysplasia underscores the value of screening for abnormal phenotypes in model organisms and identifying the causative mutations.

zCall: a rare variant caller for array-based genotyping: genetics and population analysis.

SUMMARY: zCall is a variant caller specifically designed for calling rare single-nucleotide polymorphisms from array-based technology. This caller is implemented as a post-processing step after a default calling algorithm has been applied. The algorithm uses the intensity profile of the common allele homozygote cluster to define the location of the other two genotype clusters. We demonstrate improved detection of rare alleles when applying zCall to samples that have both Illumina Infinium HumanExome BeadChip and exome sequencing data available. AVAILABILITY: http://atguweb.mgh.harvard.edu/apps/zcall. CONTACT: bneale@broadinstitute.org SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

A genome wide survey supports the involvement of large copy number variants in schizophrenia with and without intellectual disability.

BACKGROUND: Copy number variants (CNVs) have been shown to play a role in schizophrenia and intellectual disability. METHODS: We compared the CNV burden in 66 patients with intellectual disability and no symptoms of psychosis (ID-only) with the burden in 64 patients with intellectual disability and schizophrenia (ID + SCZ). Samples were genotyped on three plates by the Broad Institute using the Affymetrix 6.0 array. RESULTS: For CNVs larger than 100 kb, there was no difference in the CNV burden of ID-only and ID + SCZ. In contrast, the number of duplications larger than 1 Mb was increased in ID + SCZ compared to ID-only. We detected seven large duplications and two large deletions at chromosome 15q11.2 (18.5-20.1 Mb) which were all present in patients with ID + SCZ. The involvement of this region in schizophrenia was confirmed in Scottish samples from the ISC study (N = 2,114; 1,130 cases and 984 controls). Finally, one of the patients with schizophrenia and low IQ carrying a duplication at 15q11.2, is a member of a previously described pedigree with multiple cases of mild intellectual disability, schizophrenia, hearing impairment, retinitis pigmentosa and cataracts. DNA samples were available for 11 members of this family and the duplication was present in all 10 affected individuals and was absent in an unaffected individual. CONCLUSIONS: Duplications at 15q11.2 (18.5-20.1 Mb) are highly prevalent in a severe group of patients characterized by intellectual disability and comorbid schizophrenia. It is also associated with a phenotype that includes schizophrenia, low IQ, hearing and visual impairments resembling the spectrum of symptoms described in "ciliopathies."

ENU mutagenesis in mice identifies candidate genes for hypogonadism.

Genome-wide mutagenesis was performed in mice to identify candidate genes for male infertility, for which the predominant causes remain idiopathic. Mice were mutagenized using N-ethyl-N-nitrosourea (ENU), bred, and screened for phenotypes associated with the male urogenital system. Fifteen heritable lines were isolated and chromosomal loci were assigned using low-density genome-wide SNP arrays. Ten of the 15 lines were pursued further using higher-resolution SNP analysis to narrow the candidate gene regions. Exon sequencing of candidate genes identified mutations in mice with cystic kidneys (Bicc1), cryptorchidism (Rxfp2), restricted germ cell deficiency (Plk4), and severe germ cell deficiency (Prdm9). In two other lines with severe hypogonadism, candidate sequencing failed to identify mutations, suggesting defects in genes with previously undocumented roles in gonadal function. These genomic intervals were sequenced in their entirety and a candidate mutation was identified in SnrpE in one of the two lines. The line harboring the SnrpE variant retains substantial spermatogenesis despite small testis size, an unusual phenotype. In addition to the reproductive defects, heritable phenotypes were observed in mice with ataxia (Myo5a), tremors (Pmp22), growth retardation (unknown gene), and hydrocephalus (unknown gene). These results demonstrate that the ENU screen is an effective tool for identifying potential causes of male infertility.

Agouti C57BL/6N embryonic stem cells for mouse genetic resources.

We report the characterization of a highly germline competent C57BL/6N mouse embryonic stem cell line, JM8. To simplify breeding schemes, the dominant agouti coat color gene was restored in JM8 cells by targeted repair of the C57BL/6 nonagouti mutation. These cells provide a robust foundation for large-scale mouse knockout programs that aim to provide a public resource of targeted mutations in the C57BL/6 genetic background.

Concurrent Lpin1 and Nrcam mouse mutations result in severe peripheral neuropathy with transitory hindlimb paralysis.

Peripheral neuropathy is a broad category of disorders with a diverse etiology, grouped together by their common pathogenic effect on the peripheral nervous system (PNS). Because of the heterogeneity observed to be responsible for these disorders, a forward genetics method of gene discovery was used to identify additional affected pathways. In this report, we describe the mutant mouse line 20884, generated by N-ethyl-N-nitrosourea mutagenesis, which is characterized by adult-onset transitory hindlimb paralysis. Linkage mapping revealed that two point mutations are responsible for the phenotype: a partial loss-of-function mutation in the gene for phosphatidate phosphatase Lpin1 and a truncation mutation in the gene that encodes the neuronal cell adhesion molecule NrCAM. To investigate how the 20884 Lpin1 and Nrcam mutations interact to produce the paralysis phenotype, the double mutant and both single mutants were analyzed by quantitative behavioral, histological, and electrophysiological means. The Lpin1(20884) mutant and the double mutant are characterized by similar levels of demyelination and aberrant myelin structures. Nevertheless, the double mutant exhibits more severe electrophysiological abnormalities than the Lpin1(20884) mutant. The Nrcam(20884) mutant is characterized by normal sciatic nerve morphology and a mild electrophysiological defect. Comparison of the double mutant phenotype with the two single mutants does not point to an additive relationship between the two defects; rather, the Lpin1(20884) and Nrcam(20884) defects appear to act synergistically to produce the 20884 phenotype. It is proposed that the absence of NrCAM in a demyelinating environment has a deleterious effect, possibly by impairing the process of remyelination.

High resolution mapping and positional cloning of ENU-induced mutations in the Rw region of mouse chromosome 5.

BACKGROUND: Forward genetic screens in mice provide an unbiased means to identify genes and other functional genetic elements in the genome. Previously, a large scale ENU mutagenesis screen was conducted to query the functional content of a ~50 Mb region of the mouse genome on proximal Chr 5. The majority of phenotypic mutants recovered were embryonic lethals. RESULTS: We report the high resolution genetic mapping, complementation analyses, and positional cloning of mutations in the target region. The collection of identified alleles include several with known or presumed functions for which no mutant models have been reported (Tbc1d14, Nol14, Tyms, Cad, Fbxl5, Haus3), and mutations in genes we or others previously reported (Tapt1, Rest, Ugdh, Paxip1, Hmx1, Otoe, Nsun7). We also confirmed the causative nature of a homeotic mutation with a targeted allele, mapped a lethal mutation to a large gene desert, and localized a spermiogenesis mutation to a region in which no annotated genes have coding mutations. The mutation in Tbc1d14 provides the first implication of a critical developmental role for RAB-GAP-mediated protein transport in early embryogenesis. CONCLUSION: This collection of alleles contributes to the goal of assigning biological functions to all known genes, as well as identifying novel functional elements that would be missed by reverse genetic approaches.