The Conserved Domain Database (CDD) is the protein classification component of NCBI's Entrez query and retrieval system. CDD is linked to other Entrez databases such as Proteins, Taxonomy and PubMed®, and can be accessed at http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=cdd. CD-Search, which is available at http://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi, is a fast, interactive tool to identify conserved domains in new protein sequences. CD-Search results for protein sequences in Entrez are pre-computed to provide links between proteins and domain models, and computational annotation visible upon request. Protein–protein queries submitted to NCBI's BLAST search service at http://www.ncbi.nlm.nih.gov/BLAST are scanned for the presence of conserved domains by default. While CDD started out as essentially a mirror of publicly available domain alignment collections, such as SMART, Pfam and COG, we have continued an effort to update, and in some cases replace these models with domain hierarchies curated at the NCBI. Here, we report on the progress of the curation effort and associated improvements in the functionality of the CDD information retrieval system.

Ciliary dysfunction leads to a broad range of overlapping phenotypes, termed collectively as ciliopathies. This grouping is underscored by genetic overlap, where causal genes can also contribute modifying alleles to clinically distinct disorders. Here we show that mutations in TTC21B/IFT139, encoding a retrograde intraflagellar transport (IFT) protein, cause both isolated nephronophthisis (NPHP) and syndromic Jeune Asphyxiating Thoracic Dystrophy (JATD). Moreover, although systematic medical resequencing of a large, clinically diverse ciliopathy cohort and matched controls showed a similar frequency of rare changes, in vivo and in vitro evaluations unmasked a significant enrichment of pathogenic alleles in cases, suggesting that TTC21B contributes pathogenic alleles to ∼5% of ciliopathy patients. Our data illustrate how genetic lesions can be both causally associated with diverse ciliopathies, as well as interact in trans with other disease-causing genes, and highlight how saturated resequencing followed by functional analysis of all variants informs the genetic architecture of disorders.

ClinSeq is a pilot project to investigate the use of whole-genome sequencing as a tool for clinical research. By piloting the acquisition of large amounts of DNA sequence data from individual human subjects, we are fostering the development of hypothesis-generating approaches for performing research in genomic medicine, including the exploration of issues related to the genetic architecture of disease, implementation of genomic technology, informed consent, disclosure of genetic information, and archiving, analyzing, and displaying sequence data. In the initial phase of ClinSeq, we are enrolling roughly 1000 participants; the evaluation of each includes obtaining a detailed family and medical history, as well as a clinical evaluation. The participants are being consented broadly for research on many traits and for wholegenome sequencing. Initially, Sanger-based sequencing of 300-400 genes thought to be relevant to atherosclerosis is being performed, with the resulting data analyzed for rare, high-penetrance variants associated with specific clinical traits. The participants are also being consented to allow the contact of family members for additional studies of sequence variants to explore their potential association with specific phenotypes. Here, we present the general considerations in designing ClinSeq, preliminary results based on the generation of an initial 826 Mb of sequence data, the findings for several genes that serve as positive controls for the project, and our views about the potential implications of ClinSeq. The early experiences with ClinSeq illustrate how large-scale medical sequencing can be a practical, productive, and critical component of research in genomic medicine.

Charcot-Marie-Tooth (CMT) disease comprises a genetically and clinically heterogeneous group of peripheral nerve disorders characterized by impaired distal motor and sensory function. Mutations in three genes encoding aminoacyl-tRNA synthetases (ARSs) have been implicated in CMT disease primarily associated with an axonal pathology. ARSs are ubiquitously expressed, essential enzymes responsible for charging tRNA molecules with their cognate amino acids. To further explore the role of ARSs in CMT disease, we performed a large-scale mutation screen of the 37 human ARS genes in a cohort of 355 patients with a phenotype consistent with CMT. Here we describe three variants (p.Leu133His, p.Tyr173SerfsX7, and p.Ile302Met) in the lysyl-tRNA synthetase (KARS) gene in two patients from this cohort. Functional analyses revealed that two of these mutations (p.Leu133His and p.Tyr173SerfsX7) severely affect enzyme activity. Interestingly, both functional variants were found in a single patient with CMT disease and additional neurological and non-neurological sequelae. Based on these data, KARS becomes the fourth ARS gene associated with CMT disease, indicating that this family of enzymes is specifically critical for axon function.

Erythropoiesis is dependent on the activity of transcription factors, including the erythroid-specific erythroid Kruppel-like factor (EKLF). ChIP followed by massively parallel sequencing (ChIP-Seq) is a powerful, unbiased method to map transfactor occupancy. We used ChIP-Seq to study the interactome of EKLF in mouse erythroid progenitor cells and more differentiated erythroblasts. We correlated these results with the nuclear distribution of EKLF, RNA-Seq analysis of the transcriptome, and the occupancy of other erythroid transcription factors. In progenitor cells, EKLF is found predominantly at the periphery of the nucleus, where EKLF primarily occupies the promoter regions of genes and acts as a transcriptional activator. In erythroblasts, EKLF is distributed throughout the nucleus, and erythroblast-specific EKLF occupancy is predominantly in intragenic regions. In progenitor cells, EKLF modulates general cell growth and cell cycle regulatory pathways, whereas in erythroblasts EKLF is associated with repression of these pathways. The EKLF interactome shows very little overlap with the interactomes of GATA1, GATA2, or TAL1, leading to a model in which EKLF directs programs that are independent of those regulated by the GATA factors or TAL1. (Blood. 2011;118(17):e139-e148) IntroductionMore than 2 million red blood cells are released into the circulation every second by a multistep process known as erythropoiesis, 1 which is accompanied by significant changes in the RNA expression profile. 2 The erythroid-specific DNA binding protein erythroid Kruppel-like factor (Klf1; EKLF), the founding member of the mammalian Kruppel/Sp1-like family of C2H2-type zinc finger DNA binding proteins, plays a significant role in this process. 3,4 Klf1 mRNA and EKLF are expressed in both erythroid progenitor cells and terminally differentiating erythrocytes. 5,6 In cell lines and in vitro, EKLF has been demonstrated to physically and functionally interact with DNA at a conserved CCNCNCCCN motif and with additional protein cofactors. 7,8 In EKLF-deficient (Klf1 Ϫ/Ϫ ) mice, definitive fetal liver erythroid progenitor cells fail to progress to erythroblasts, 9 leading to anemia and lethality by embryonic day 15. 10,11 Genome-wide mRNA profiling of WT and Klf1 Ϫ/Ϫ erythroid cells has revealed significant dysregulation of Ͼ 3000 mRNAs. 9,12-15 Two-thirds of the dysregulated transcripts were present at reduced levels in Klf1 Ϫ/Ϫ cells, consistent with the role of EKLF as a transcriptional activator, whereas one-third were present at increased levels, consistent with the role of EKLF as a transcriptional repressor. 16 Chromatin immunoprecipitation (ChIP) demonstrated EKLF occupancy at selected sites in genes encoding erythrocyte proteins AHSP, ankyrin, ␤-spectrin, Band 3, and dematin, all of which are down-regulated in EKLF-deficient cells. 9,12-15 EKLF also participates in chromatin modification and DNase hypersensitive site formation through interactions with CBP/p300 and an SWI/SNFrelated chromatin remodeling complex. 17,18 For example,...

Micrognathia, glossoptosis, and cleft palate comprise one of the most common malformation sequences, Robin sequence. It is a component of the TARP syndrome, talipes equinovarus, atrial septal defect, Robin sequence, and persistent left superior vena cava. This disorder is X-linked and severe, with apparently 100% pre- or postnatal lethality in affected males. Here we characterize a second family with TARP syndrome, confirm linkage to Xp11.23-q13.3, perform massively parallel sequencing of X chromosome exons, filter the results via a number of criteria including the linkage region, use a unique algorithm to characterize sequence changes, and show that TARP syndrome is caused by mutations in the RBM10 gene, which encodes RNA binding motif 10. We further show that this previously uncharacterized gene is expressed in midgestation mouse embryos in the branchial arches and limbs, consistent with the human phenotype. We conclude that massively parallel sequencing is useful to characterize large candidate linkage intervals and that it can be used successfully to allow identification of disease-causing gene mutations.

We report an early onset spastic ataxia-neuropathy syndrome in two brothers of a consanguineous family characterized clinically by lower extremity spasticity, peripheral neuropathy, ptosis, oculomotor apraxia, dystonia, cerebellar atrophy, and progressive myoclonic epilepsy. Whole-exome sequencing identified a homozygous missense mutation (c.1847G>A; p.Y616C) in AFG3L2, encoding a subunit of an m-AAA protease. m-AAA proteases reside in the mitochondrial inner membrane and are responsible for removal of damaged or misfolded proteins and proteolytic activation of essential mitochondrial proteins. AFG3L2 forms either a homo-oligomeric isoenzyme or a hetero-oligomeric complex with paraplegin, a homologous protein mutated in hereditary spastic paraplegia type 7 (SPG7). Heterozygous loss-of-function mutations in AFG3L2 cause autosomal-dominant spinocerebellar ataxia type 28 (SCA28), a disorder whose phenotype is strikingly different from that of our patients. As defined in yeast complementation assays, the AFG3L2Y616C gene product is a hypomorphic variant that exhibited oligomerization defects in yeast as well as in patient fibroblasts. Specifically, the formation of AFG3L2Y616C complexes was impaired, both with itself and to a greater extent with paraplegin. This produced an early-onset clinical syndrome that combines the severe phenotypes of SPG7 and SCA28, in additional to other “mitochondrial” features such as oculomotor apraxia, extrapyramidal dysfunction, and myoclonic epilepsy. These findings expand the phenotype associated with AFG3L2 mutations and suggest that AFG3L2-related disease should be considered in the differential diagnosis of spastic ataxias.

G protein-coupled receptors (GPCRs), the largest human gene family, are important regulators of signaling pathways. However, knowledge of their genetic alterations is limited. In this study, we used exon capture and massively parallel sequencing methods to analyze the mutational status of 734 GPCRs in melanoma. This investigation revealed that one family member, GRM3, was frequently mutated and that one of its mutations clustered within one position. Biochemical analysis of GRM3 alterations revealed that mutant GRM3 selectively regulated the phosphorylation of MEK, leading to increased anchorage-independent growth and migration. Melanoma cells expressing mutant GRM3 had reduced cell growth and cellular migration after short hairpin RNA–mediated knockdown of GRM3 or treatment with a selective MEK inhibitor, AZD-6244, which is currently being used in phase 2 clinical trials. Our study yields the most comprehensive map of genetic alterations in the GPCR gene family.