Search

Antibodies with ontogenies from VH1-2 or VH1-46-germline genes dominate the broadly neutralizing response against the CD4-binding site (CD4bs) on HIV-1. Here, we define with longitudinal sampling from time-of-infection the development of a VH1-46-derived antibody lineage that matured to neutralize 90% of HIV-1 isolates. Structures of lineage antibodies CH235 (week 41 from time-of-infection, 18% breadth), CH235.9 (week 152, 77%), and CH235.12 (week 323, 90%) demonstrated the maturing epitope to focus on the conformationally invariant portion of the CD4bs. Similarities between CH235 lineage and five unrelated CD4bs lineages in epitope focusing, length-of-time to develop breadth, and extraordinary level of somatic hypermutation suggested commonalities in maturation among all CD4bs antibodies. Fortunately, the required CH235-lineage hypermutation appeared substantially guided by the intrinsic mutability of the VH1-46 gene, which closely resembled VH1-2. We integrated our CH235-lineage findings with a second broadly neutralizing lineage and HIV-1 co-evolution to suggest a vaccination strategy for inducing both lineages.

The aim of the study was to assess exome data for preemptive pharmacogenetic screening for 203 clinically relevant pharmacogenetic variant positions from the Pharmacogenomics Knowledgebase and Clinical Pharmacogenetics Implementation Consortium and identify copy-number variants (CNVs) in CYP2D6. We examined the coverage and genotype quality of 203 pharmacogenetic variant positions in 973 exomes compared with five genomes and with five genotyping chip data sets. Then, we determined the agreement of exome and chip genotypes by evaluating concordance in a three-way comparison of exome, genome, and chip-based genotyping at 1,929 variant positions in five individuals. Finally, we evaluated the utility of exomes for detecting CYP2D6 CNVs. For 5 individuals examined for 203 pharmacogenetic variants (5 × 203 = 1,015), 998/1,015 were identified by genome, 849/1,015 were identified by exome, and 295/1,015 by genotyping chip. Thirty-six pharmacogenetic star allele variants with moderate to strong Clinical Pharmacogenetics Implementation Consortium (CPIC) therapeutic recommendations were identified in 973 exomes. Exomes had high (98%) genotype concordance with chip-based genotyping. CYP2D6 CNVs were identified in 57/973 exomes. Exomes outperformed the current chip-based assay in detecting more important pharmacogenetic variant positions and CYP2D6 CNVs for preemptive pharmacogenetic screening. Tools should be developed to derive pharmacogenetic variants from exomes. Genet Med 19 3, 357–361.

Endometrial cancer is the 6th most commonly diagnosed cancer among women worldwide, causing ~74,000 deaths annually 1. Serous endometrial cancers are a clinically aggressive subtype with a poorly defined genetic etiology 2-4. We used whole exome sequencing (WES) to comprehensively search for somatic mutations within ~22,000 protein-encoding genes among 13 primary serous endometrial tumors. We subsequently resequenced 18 genes that were mutated in more than one tumor, and/or were genes that formed an enriched functional grouping, from 40 additional serous tumors. We identified high frequencies of somatic mutations in CHD4 (17%), EP300 (8%), ARID1A (6%), TSPYL2 (6%), FBXW7 (29%), SPOP (8%), MAP3K4 (6%) and ABCC9 (6%). Overall, 36.5% of serous tumors had mutated a chromatin-remodeling gene and 35% had mutated a ubiquitin ligase complex gene, implicating the frequent mutational disruption of these processes in the molecular pathogenesis of one of the deadliest forms of endometrial cancer.

What Kitty Shares with Kings Although long-studied, the underlying basis of mammalian coat patterns remains unclear. By studying a large number of cat species and varieties, Kaelin et al. (p. 1536 ) identified two genes, Taqpep and Edn3 , as critical factors in the development of feline pigment patterns. Mutations in Taqpep are responsible for the blotched tabby pattern in domestic cats and the unusual coat of wild king cheetahs. Gene expression patterns in cat and cheetah skin suggest that Edn3 is a likely regulator of felid hair color. The findings support a common model for coat and pigment pattern formation in domestic and wild cats.

Long noncoding RNAs (lncRNAs) play a broad range of biological roles, including regulation of expression of genes and chromosomes. Here, we present evidence that lncRNAs are involved in vertebrate circadian biology. Differential night/day expression of 112 lncRNAs (0.3 to >50 kb) occurs in the rat pineal gland, which is the source of melatonin, the hormone of the night. Approximately one-half of these changes reflect nocturnal increases. Studies of eight lncRNAs with 2- to >100-fold daily rhythms indicate that, in most cases, the change results from neural stimulation from the central circadian oscillator in the suprachiasmatic nucleus (doubling time = 0.5–1.3 h). Light exposure at night rapidly reverses (halving time = 9–32 min) levels of some of these lncRNAs. Organ culture studies indicate that expression of these lncRNAs is regulated by norepinephrine acting through cAMP. These findings point to a dynamic role of lncRNAs in the circadian system.

Genome- and exome-sequencing costs are continuing to fall, and many individuals are undergoing these assessments as research participants and patients. The issue of secondary (so-called incidental) findings in exome analysis is controversial, and data are needed on methods of detection and their frequency. We piloted secondary variant detection by analyzing exomes for mutations in cancer-susceptibility syndromes in subjects ascertained for atherosclerosis phenotypes. We performed exome sequencing on 572 ClinSeq participants, and in 37 genes, we interpreted variants that cause high-penetrance cancer syndromes by using an algorithm that filtered results on the basis of mutation type, quality, and frequency and that filtered mutation-database entries on the basis of defined categories of causation. We identified 454 sequence variants that differed from the human reference. Exclusions were made on the basis of sequence quality (26 variants) and high frequency in the cohort (77 variants) or dbSNP (17 variants), leaving 334 variants of potential clinical importance. These were further filtered on the basis of curation of literature reports. Seven participants, four of whom were of Ashkenazi Jewish descent and three of whom did not meet family-history-based referral criteria, had deleterious BRCA1 or BRCA2 mutations. One participant had a deleterious SDHC mutation, which causes paragangliomas. Exome sequencing, coupled with multidisciplinary interpretation, detected clinically important mutations in cancer-susceptibility genes; four of such mutations were in individuals without a significant family history of disease. We conclude that secondary variants of high clinical importance will be detected at an appreciable frequency in exomes, and we suggest that priority be given to the development of more efficient modes of interpretation with trials in larger patient groups.

Objective: To utilize high-throughput sequencing to determine the etiology of juvenile-onset neurodegeneration in a 19-year-old woman with progressive motor and cognitive decline. Methods: Exome sequencing identified an initial list of 133,555 variants in the proband9s family, which were filtered using segregation analysis, presence in dbSNP, and an empirically derived gene exclusion list. The filtered list comprised 52 genes: 21 homozygous variants and 31 compound heterozygous variants. These variants were subsequently scrutinized with predicted pathogenicity programs and for association with appropriate clinical syndromes. Results: Exome sequencing data identified 2 GLB1 variants (c.602G>A, p.R201H; c.785G>T, p.G262V). β-Galactosidase enzyme analysis prior to our evaluation was reported as normal; however, subsequent testing was consistent with juvenile-onset GM 1 -gangliosidosis. Urine oligosaccharide analysis was positive for multiple oligosaccharides with terminal galactose residues. Conclusions: We describe a patient with juvenile-onset neurodegeneration that had eluded diagnosis for over a decade. GM 1 -gangliosidosis had previously been excluded from consideration, but was subsequently identified as the correct diagnosis using exome sequencing. Exome sequencing can evaluate genes not previously associated with neurodegeneration, as well as most known neurodegeneration-associated genes. Our results demonstrate the utility of “agnostic” exome sequencing to evaluate patients with undiagnosed disorders, without prejudice from prior testing results.

Genomic technologies, such as whole-exome sequencing, are a powerful tool in genetic research. Such testing yields a great deal of incidental medical information, or medical information not related to the primary research target. We describe the management of incidental medical information derived from whole-exome sequencing in the research context. We performed whole-exome sequencing on a monozygotic twin pair in which only 1 child was affected with congenital anomalies and applied an institutional review board–approved algorithm to determine what genetic information would be returned. Whole-exome sequencing identified 79 525 genetic variants in the twins. Here, we focus on novel variants. After filtering artifacts and excluding known single nucleotide polymorphisms and variants not predicted to be pathogenic, the twins had 32 novel variants in 32 genes that were felt to be likely to be associated with human disease. Eighteen of these novel variants were associated with recessive disease and 18 were associated with dominantly manifesting conditions (variants in some genes were potentially associated with both recessive and dominant conditions), but only 1 variant ultimately met our institutional review board–approved criteria for return of information to the research participants.

The current genetic and recombination maps of the cat have fewer than 3,000 markers and a resolution limit greater than 1 Mb. To complement the first-generation domestic cat maps, support higher resolution mapping studies, and aid genome assembly in specific areas as well as in the whole genome, a 15,000Rad radiation hybrid (RH) panel for the domestic cat was generated. Fibroblasts from the female Abyssinian cat that was used to generate the cat genomic sequence were fused to a Chinese hamster cell line (A23), producing 150 hybrid lines. The clones were initially characterized using 39 short tandem repeats (STRs) and 1,536 SNP markers. The utility of whole-genome amplification in preserving and extending RH panel DNA was also tested using 10 STR markers; no significant difference in retention was observed. The resolution of the 15,000Rad RH panel was established by constructing framework maps across 10 different 1-Mb regions on different feline chromosomes. In these regions, 2-point analysis was used to estimate RH distances, which compared favorably with the estimation of physical distances. The study demonstrates that the 15,000Rad RH panel constitutes a powerful tool for constructing high-resolution maps, having an average resolution of 40.1 kb per marker across the ten 1-Mb regions. In addition, the RH panel will complement existing genomic resources for the domestic cat, aid in the accurate re-assemblies of the forthcoming cat genomic sequence, and support cross-species genomic comparisons.

Large-scale genomic analysis such as whole-exome and whole-genome sequencing is becoming increasingly prevalent in the research arena. Clinically, many potential uses of this technology have been proposed. One such application is the extension or augmentation of newborn screening. In order to explore this application, we examined data from 3 children with normal newborn screens who underwent whole-exome sequencing as part of research participation. We analyzed sequence information for 151 selected genes associated with conditions ascertained by newborn screening. We compared findings with publicly available databases and results from over 500 individuals who underwent whole-exome sequencing at the same facility. Novel variants were confirmed through bidirectional dideoxynucleotide sequencing. High-density microarrays (Illumina Omni1-Quad) were also performed to detect potential copy number variations affecting these genes. We detected an average of 87 genetic variants per individual. After excluding artifacts, 96% of the variants were found to be reported in public databases and have no evidence of pathogenicity. No variants were identified that would predict disease in the tested individuals, which is in accordance with their normal newborn screens. However, we identified 6 previously reported variants and 2 novel variants that, according to published literature, could result in affected offspring if the reproductive partner were also a mutation carrier; other specific molecular findings highlight additional means by which genomic testing could augment newborn screening.

Ciliopathies are caused by alterations in the development and function of cilia. Now Jeffrey Martens and his colleagues demonstrate anatomic and functional rescue of cilia development in mature, differentiated neurons by adenovirus-mediated restoration of expression of the wild-type protein intraflagellar transport protein 88 (Ift88) and show restoration of olfactory function in a mouse model of ciliopathy. A loss-of-function mutation in IFT88 is also identified in individuals with ciliopathies. Cilia are evolutionarily conserved microtubule-based organelles that are crucial for diverse biological functions, including motility, cell signaling and sensory perception1. In humans, alterations in the formation and function of cilia manifest clinically as ciliopathies, a growing class of pleiotropic genetic disorders2,3,4. Despite the substantial progress that has been made in identifying genes that cause ciliopathies, therapies for these disorders are not yet available to patients. Although mice with a hypomorphic mutation in the intraflagellar transport protein IFT88 (Ift88Tg737Rpw mice, also known as ORPK mice)5 have been well studied, the relevance of IFT88 mutations to human pathology is unknown. We show that a mutation in IFT88 causes a hitherto unknown human ciliopathy. In vivo complementation assays in zebrafish and mIMCD3 cells show the pathogenicity of this newly discovered allele. We further show that ORPK mice are functionally anosmic as a result of the loss of cilia on their olfactory sensory neurons (OSNs). Notably, adenoviral-mediated expression of IFT88 in mature, fully differentiated OSNs of ORPK mice is sufficient to restore ciliary structures and rescue olfactory function. These studies are the first to use in vivo therapeutic treatment to reestablish cilia in a mammalian ciliopathy. More broadly, our studies indicate that gene therapy is a viable option for cellular and functional rescue of the complex ciliary organelle in established differentiated cells.

A b s t r ac t Background The Proteus syndrome is characterized by the overgrowth of skin, connective tissue, brain, and other tissues. It has been hypothesized that the syndrome is caused by somatic mosaicism for a mutation that is lethal in the nonmosaic state. Methods We performed exome sequencing of DNA from biopsy samples obtained from patients with the Proteus syndrome and compared the resultant DNA sequences with those of unaffected tissues obtained from the same patients. We confirmed and extended an observed association, using a custom restriction-enzyme assay to analyze the DNA in 158 samples from 29 patients with the Proteus syndrome. We then assayed activation of the AKT protein in affected tissues, using phosphorylationspecific antibodies on Western blots. Results Of 29 patients with the Proteus syndrome, 26 had a somatic activating mutation (c.49G→A, p.Glu17Lys) in the oncogene AKT1, encoding the AKT1 kinase, an enzyme known to mediate processes such as cell proliferation and apoptosis. Tissues and cell lines from patients with the Proteus syndrome harbored admixtures of mutant alleles that ranged from 1% to approximately 50%. Mutant cell lines showed greater AKT phosphorylation than did control cell lines. A pair of single-cell clones that were established from the same starting culture and differed with respect to their mutation status had different levels of AKT phosphorylation. Conclusions The Proteus syndrome is caused by a somatic activating mutation in AKT1, proving the hypothesis of somatic mosaicism and implicating activation of the PI3K–AKT pathway in the characteristic clinical findings of overgrowth and tumor susceptibility in this disorder. (Funded by the Intramural Research Program of the National Human Genome Research Institute.)

Recent advances in sequencing technology have led to a rapid accumulation of mitochondrial DNA (mtDNA) sequences, which now represent the wide spectrum of animal diversity. However, one animal phylum--Ctenophora--has, to date, remained completely unsampled. Ctenophores, a small group of marine animals, are of interest due to their unusual biology, controversial phylogenetic position, and devastating impact as invasive species. Using data from the Mnemiopsis leidyi genome sequencing project, we Polymerase Chain Reaction (PCR) amplified and analyzed its complete mitochondrial (mt-) genome. At just over 10 kb, the mt-genome of M. leidyi is the smallest animal mtDNA ever reported and is among the most derived. It has lost at least 25 genes, including atp6 and all tRNA genes. We show that atp6 has been relocated to the nuclear genome and has acquired introns and a mitochondrial targeting presequence, while tRNA genes have been genuinely lost, along with nuclear-encoded mt-aminoacyl tRNA synthetases. The mt-genome of M. leidyi also displays extremely high rates of sequence evolution, which likely led to the degeneration of both protein and rRNA genes. In particular, encoded rRNA molecules possess little similarity with their homologs in other organisms and have highly reduced secondary structures. At the same time, nuclear encoded mt-ribosomal proteins have undergone expansions, likely to compensate for the reductions in mt-rRNA. The unusual features identified in M. leidyi mtDNA make this organism an interesting system for the study of various aspects of mitochondrial biology, particularly protein and tRNA import and mt-ribosome structures, and add to its value as an emerging model species. Furthermore, the fast-evolving M. leidyi mtDNA should be a convenient molecular marker for species- and population-level studies.

Massively parallel DNA sequencing technologies have greatly increased our ability to generate large amounts of sequencing data at a rapid pace. Several methods have been developed to enrich for genomic regions of interest for targeted sequencing. We have compared three of these methods: Molecular Inversion Probes (MIP), Solution Hybrid Selection (SHS), and Microarray-based Genomic Selection (MGS). Using HapMap DNA samples, we compared each of these methods with respect to their ability to capture an identical set of exons and evolutionarily conserved regions associated with 528 genes (2.61 Mb). For sequence analysis, we developed and used a novel Bayesian genotype-assigning algorithm, Most Probable Genotype (MPG). All three capture methods were effective, but sensitivities (percentage of targeted bases associated with high-quality genotypes) varied for an equivalent amount of pass-filtered sequence: for example, 70% (MIP), 84% (SHS), and 91% (MGS) for 400 Mb. In contrast, all methods yielded similar accuracies of >99.84% when compared to Infinium 1M SNP BeadChip-derived genotypes and >99.998% when compared to 30-fold coverage whole-genome shotgun sequencing data. We also observed a low false-positive rate with all three methods; of the heterozygous positions identified by each of the capture methods, >99.57% agreed with 1M SNP BeadChip, and >98.840% agreed with the whole-genome shotgun data. In addition, we successfully piloted the genomic enrichment of a set of 12 pooled samples via the MGS method using molecular bar codes. We find that these three genomic enrichment methods are highly accurate and practical, with sensitivities comparable to that of 30-fold coverage whole-genome shotgun data.

The development of massively parallel sequencing technologies, coupled with new massively parallel DNA enrichment technologies (genomic capture), has allowed the sequencing of targeted regions of the human genome in rapidly increasing numbers of samples. Genomic capture can target specific areas in the genome, including genes of interest and linkage regions, but this limits the study to what is already known. Exome capture allows an unbiased investigation of the complete protein-coding regions in the genome. Researchers can use exome capture to focus on a critical part of the human genome, allowing larger numbers of samples than are currently practical with whole-genome sequencing. In this review, we briefly describe some of the methodologies currently used for genomic and exome capture and highlight recent applications of this technology.

Summary Virtually the entire surface of the HIV-1-envelope trimer is recognized by neutralizing antibodies, except for a highly glycosylated region at the center of the "silent face" on the gp120 subunit. From an HIV-1-infected donor, #74, we identified antibody VRC-PG05, which neutralized 27% of HIV-1 strains. The crystal structure of the antigen-binding fragment of VRC-PG05 in complex with gp120 revealed an epitope comprised primarily of N-linked glycans from N262, N295, and N448 at the silent face center. Somatic hypermutation occurred preferentially at antibody residues that interacted with these glycans, suggesting somatic development of glycan recognition. Resistance to VRC-PG05 in donor #74 involved shifting of glycan-N448 to N446 or mutation of glycan-proximal residue E293. HIV-1 neutralization can thus be achieved at the silent face center by glycan-recognizing antibody; along with other known epitopes, the VRC-PG05 epitope completes coverage by neutralizing antibody of all major exposed regions of the prefusion closed trimer.

Human adenylate kinase 2 deficiency causes the most profound human haematopoietic defect associated with sensorineural deafness. Reticular dysgenesis (RD) is an autosomal recessive form of human severe combined immunodeficiency (SCID) characterized by an early differentiation block in the myeloid lineage and a profound impairment in lymphoid maturation associated with bilateral sensorineural deafness. The lack of polymorphonuclear neutrophils in affected newborns is responsible for the occurrence of severe infections earlier than usually observed in the other forms of SCID. Furthermore, RD associated neutropenia is characterized by the lack of responsiveness to G-CSF. We have identified bi-allelic mutations in the adenylate kinase 2 (AK2) gene is seven patients affected with RD. These mutations resulted in the absence or a strong decrease in protein expression. Restoration of AK2 expression in the bone marrow cells of RD patients overcomes the neutrophil differentiation arrest, underlying its specific requirement in the development of a restricted set of haematopoietic lineages. Lastly, we established that AK2 is specifically expressed in the stria vascularis region of the inner ear. The function of this gene in the differentiation of a given set of cell lineages is rapidly in progress, showing that studies of primary immunodeficiencies continue to provide key information on human lympho-haematopoietic development. Moreover the AK2 enzyme seems a key molecule in different biological systems such as the lympho-haematopoiesis and the brain development.

Genetic variation among individual humans occurs on many different scales, ranging from gross alterations in the human karyotype to single nucleotide changes. Here we explore variation on an intermediate scale--particularly insertions, deletions and inversions affecting from a few thousand to a few million base pairs. We employed a clone-based method to interrogate this intermediate structural variation in eight individuals of diverse geographic ancestry. Our analysis provides a comprehensive overview of the normal pattern of structural variation present in these genomes, refining the location of 1,695 structural variants. We find that 50% were seen in more than one individual and that nearly half lay outside regions of the genome previously described as structurally variant. We discover 525 new insertion sequences that are not present in the human reference genome and show that many of these are variable in copy number between individuals. Complete sequencing of 261 structural variants reveals considerable locus complexity and provides insights into the different mutational processes that have shaped the human genome. These data provide the first high-resolution sequence map of human structural variation--a standard for genotyping platforms and a prelude to future individual genome sequencing projects.

Systemic autoinflammatory diseases are driven by abnormal activation of innate immunity. Herein we describe a new disease caused by high-penetrance heterozygous germline mutations in TNFAIP3, which encodes the NF-κB regulatory protein A20, in six unrelated families with early-onset systemic inflammation. The disorder resembles Behçet's disease, which is typically considered a polygenic disorder with onset in early adulthood. A20 is a potent inhibitor of the NF-κB signaling pathway. Mutant, truncated A20 proteins are likely to act through haploinsufficiency because they do not exert a dominant-negative effect in overexpression experiments. Patient-derived cells show increased degradation of IκBα and nuclear translocation of the NF-κB p65 subunit together with increased expression of NF-κB-mediated proinflammatory cytokines. A20 restricts NF-κB signals via its deubiquitinase activity. In cells expressing mutant A20 protein, there is defective removal of Lys63-linked ubiquitin from TRAF6, NEMO and RIP1 after stimulation with tumor necrosis factor (TNF). NF-κB-dependent proinflammatory cytokines are potential therapeutic targets for the patients with this disease.

Background In chicken, loss of TALPID3 results in non-functional cilia and short-rib polydactyly syndrome. This phenotype is caused by a frameshift mutation in the chicken ortholog of the human KIAA0586 gene, which encodes a novel coiled-coil domain protein essential for primary ciliogenesis, suggesting that KIAA0586 can be associated with ciliopathy in human beings. Methods In our patients with ciliopathy (http://www.clinicaltrials.gov: NCT00068224), we have collected extensive clinical and neuroimaging data from affected individuals, and performed whole exome sequencing on DNA from affected individuals and their parents. We analysed gene expression on fibroblast cell line, and determined the effect of gene mutation on ciliogenesis in cells derived from patients. Results We identified biallelic mutations in the human TALPID3 ortholog, KIAA0586 , in six children with findings of overlapping Jeune and Joubert syndromes. Fibroblasts cultured from one of the patients with Jeune–Joubert syndrome exhibited more severe cilia defects than fibroblasts from patients with only Joubert syndrome; this difference was reflected in KIAA0586 RNA expression levels. Rescue of the cilia defect with full-length wild type KIAA0586 indicated a causal link between cilia formation and KIAA0586 function. Conclusions Our results show that biallelic deleterious mutations in KIAA0586 lead to Joubert syndrome with or without Jeune asphyxiating thoracic dystrophy. Furthermore, our results confirm that KIAA0586/TALPID3 is essential in cilia formation in human beings, expand the KIAA0586 phenotype to include features of Jeune syndrome and provide a pathogenetic connection between Joubert and Jeune syndromes, based on aberrant ciliogenesis.

The term neurotranscriptomics is used here to describe genome-wide analysis of neural control of transcriptomes. In this report, next-generation RNA sequencing was using to analyze the effects of neonatal (5-days-of-age) surgical stimulus deprivation on the adult rat pineal transcriptome. In intact animals, more than 3000 coding genes were found to exhibit differential expression (adjusted-p < 0.001) on a night/day basis in the pineal gland (70% of these increased at night, 376 genes changed more than 4-fold in either direction). Of these, more than two thousand genes were not previously known to be differentially expressed on a night/day basis. The night/day changes in expression were almost completely eliminated by neonatal removal (SCGX) or decentralization (DCN) of the superior cervical ganglia (SCG), which innervate the pineal gland. Other than the loss of rhythmic variation, surgical stimulus deprivation had little impact on the abundance of most genes; of particular interest, expression levels of the melatonin-synthesis-related genes Tph1, Gch1, and Asmt displayed little change (less than 35%) following DCN or SCGX. However, strong and consistent changes were observed in the expression of a small number of genes including the gene encoding Serpina1, a secreted protease inhibitor that might influence extracellular architecture. Many of the genes that exhibited night/day differential expression in intact animals also exhibited similar changes following in vitro treatment with norepinephrine, a superior cervical ganglia transmitter, or with an analog of cyclic AMP, a norepinephrine second messenger in this tissue. These findings are of significance in that they establish that the pineal-defining transcriptome is established prior to the neonatal period. Further, this work expands our knowledge of the biological process under neural control in this tissue and underlines the value of RNA sequencing in revealing how neurotransmission influences cell biology.

Background Laminins are heterotrimeric complexes, consisting of α, β and γ subunits that form a major component of basement membranes and extracellular matrix. Laminin complexes have different, but often overlapping, distributions and functions. Methods Under our clinical protocol, NCT00068224, we have performed extensive clinical and neuropsychiatric phenotyping, neuroimaging and molecular analysis in patients with laminin α1 ( LAMA1 )-associated lamininopathy. We investigated the consequence of mutations in LAMA1 using patient-derived fibroblasts and neuronal cells derived from neuronal stem cells. Results In this paper we describe individuals with biallelic mutations in LAMA1 , all of whom had the cerebellar dysplasia, myopia and retinal dystrophy, in addition to obsessive compulsive traits, tics and anxiety. Patient-derived fibroblasts have impaired adhesion, reduced migration, abnormal morphology and increased apoptosis due to impaired activation of Cdc42, a member of the Rho family of GTPases that is involved in cytoskeletal dynamics. LAMA1 knockdown in human neuronal cells also showed abnormal morphology and filopodia formation, supporting the importance of LAMA1 in neuronal migration, and marking these cells potentially useful tools for disease modelling and therapeutic target discovery. Conclusion This paper broadens the phenotypes associated with LAMA1 mutations. We demonstrate that LAMA1 deficiency can lead to alteration in cytoskeletal dynamics, which may invariably lead to alteration in dendrite growth and axonal formation. Estimation of disease prevalence based on population studies in LAMA1 reveals a prevalence of 1–20 in 1 000 000. Trial registration number: NCT00068224

SummaryThe site on the HIV-1 gp120 glycoprotein that binds the CD4 receptor is recognized by broadly reactive antibodies, several of which neutralize over 90% of HIV-1 strains. To understand how antibodies achieve such neutralization, we isolated CD4-binding-site (CD4bs) antibodies and analyzed 16 co-crystal structures –8 determined here– of CD4bs antibodies from 14 donors. The 16 antibodies segregated by recognition mode and developmental ontogeny into two types: CDR H3-dominated and VH-gene-restricted. Both could achieve greater than 80% neutralization breadth, and both could develop in the same donor. Although paratope chemistries differed, all 16 gp120-CD4bs antibody complexes showed geometric similarity, with antibody-neutralization breadth correlating with antibody-angle of approach relative to the most effective antibody of each type. The repertoire for effective recognition of the CD4 supersite thus comprises antibodies with distinct paratopes arrayed about two optimal geometric orientations, one achieved by CDR H3 ontogenies and the other achieved by VH-gene-restricted ontogenies.

Limb body wall complex (LBWC) and amniotic band sequence (ABS) are multiple congenital anomaly conditions with craniofacial, limb, and ventral wall defects. LBWC and ABS are considered separate entities by some, and a continuum of severity of the same condition by others. The etiology of LBWC/ABS remains unknown and multiple hypotheses have been proposed. One individual with features of LBWC and his unaffected parents were whole exome sequenced and Sanger sequenced as confirmation of the mutation. Functional studies were conducted using morpholino knockdown studies followed by human mRNA rescue experiments. Using whole exome sequencing, a de novo heterozygous mutation was found in the gene IQCK: c.667C>G; p.Q223E and confirmed by Sanger sequencing in an individual with LBWC. Morpholino knockdown of iqck mRNA in the zebrafish showed ventral defects including failure of ventral fin to develop and cardiac edema. Human wild-type IQCK mRNA rescued the zebrafish phenotype, whereas human p.Q223E IQCK mRNA did not, but worsened the phenotype of the morpholino knockdown zebrafish. This study supports a genetic etiology for LBWC/ABS, or potentially a new syndrome.

Inherited genetic variation has a critical but as yet largely uncharacterized role in human disease. Here we report a public database of common variation in the human genome: more than one million single nucleotide polymorphisms (SNPs) for which accurate and complete genotypes have been obtained in 269 DNA samples from four populations, including ten 500-kilobase regions in which essentially all information about common DNA variation has been extracted. These data document the generality of recombination hotspots, a block-like structure of linkage disequilibrium and low haplotype diversity, leading to substantial correlations of SNPs with many of their neighbours. We show how the HapMap resource can guide the design and analysis of genetic association studies, shed light on structural variation and recombination, and identify loci that may have been subject to natural selection during human evolution. © 2005 Nature Publishing Group., link_to_OA_fulltext

Inbred mouse strains provide the foundation for mouse genetics. By selecting for phenotypic features of interest, inbreeding drives genomic evolution and eliminates individual variation, while fixing certain sets of alleles that are responsible for the trait characteristics of the strain. Mouse strains 129Sv (129S5) and C57BL/6J, two of the most widely used inbred lines, diverged from common ancestors within the last century, yet very little is known about the genomic differences between them. By comparative genomic hybridization and sequence analysis of 129S5 short insert libraries, we identified substantial structural variation, a complex fine-scale haplotype pattern with a continuous distribution of diversity blocks, and extensive nucleotide variation, including nonsynonymous coding SNPs and stop codons. Collectively, these genomic changes denote the level and direction of allele fixation that has occurred during inbreeding and provide a basis for defining what makes these mouse strains unique.

The human X chromosome has a unique biology that was shaped by its evolution as the sex chromosome shared by males and females. We have determined 99.3% of the euchromatic sequence of the X chromosome. Our analysis illustrates the autosomal origin of the mammalian sex chromosomes, the stepwise process that led to the progressive loss of recombination between X and Y, and the extent of subsequent degradation of the Y chromosome. LINE1 repeat elements cover one-third of the X chromosome, with a distribution that is consistent with their proposed role as way stations in the process of X-chromosome inactivation. We found 1,098 genes in the sequence, of which 99 encode proteins expressed in testis and in various tumour types. A disproportionately high number of mendelian diseases are documented for the X chromosome. Of this number, 168 have been explained by mutations in 113 X-linked genes, which in many cases were characterized with the aid of the DNA sequence.

Chromosome 9 is highly structurally polymorphic. It contains the largest autosomal block of heterochromatin, which is heteromorphic in 6–8% of humans, whereas pericentric inversions occur in more than 1% of the population. The finished euchromatic sequence of chromosome 9 comprises 109,044,351 base pairs and represents >99.6% of the region. Analysis of the sequence reveals many intra- and interchromosomal duplications, including segmental duplications adjacent to both the centromere and the large heterochromatic block. We have annotated 1,149 genes, including genes implicated in male-to-female sex reversal, cancer and neurodegenerative disease, and 426 pseudogenes. The chromosome contains the largest interferon gene cluster in the human genome. There is also a region of exceptionally high gene and G + C content including genes paralogous to those in the major histocompatibility complex. We have also detected recently duplicated genes that exhibit different rates of sequence divergence, presumably reflecting natural selection.

Analysis of the human genome sequence has identified ≈25,000–30,000 protein-coding genes, but little is known about how most of these are regulated. Mapping DNase I hypersensitive (HS) sites has traditionally represented the gold-standard experimental method for identifying regulatory elements, but the labor-intensive nature of this technique has limited its application to only a small number of human genes. We have developed a protocol to generate a genome-wide library of gene regulatory sequences by cloning DNase HS sites. We generated a library of DNase HS sites from quiescent primary human CD4 + T cells and analyzed ≈5,600 of the resulting clones. Compared to sequences from randomly generated in silico libraries, sequences from these clones were found to map more frequently to regions of the genome known to contain regulatory elements, such as regions upstream of genes, within CpG islands, and in sequences that align between mouse and human. These cloned sites also tend to map near genes that have detectable transcripts in CD4 + T cells, demonstrating that transcriptionally active regions of the genome are being selected. Validation of putative regulatory elements was achieved by repeated recovery of the same sequence and real-time PCR. This cloning strategy, which can be scaled up and applied to any cell line or tissue, will be useful in identifying regulatory elements controlling global expression differences that delineate tissue types, stages of development, and disease susceptibility.

The existence of reliable mtDNA reference sequences for each species is of great relevance in a variety of fields, from phylogenetic and population genetics studies to pathogenetic determination of mtDNA variants in humans or in animal models of mtDNA-linked diseases. We present compelling evidence for the existence of sequencing errors on the current mouse mtDNA reference sequence. This includes the deletion of a full codon in two genes, the substitution of one amino acid on five occasions and also the involvement of tRNA and rRNA genes. The conclusions are supported by: (i) the re-sequencing of the original cell line used by Bibb and Clayton, the LA9 cell line, (ii) the sequencing of a second L-derivative clone (L929), and (iii) the comparison with 12 other mtDNA sequences from live mice, 10 of them maternally related with the mouse from which the L cells were generated. Two of the latest sequences are reported for the first time in this study (Balb/cJ and C57BL/6J). In addition, we found that both the LA9 and L929 mtDNAs also contain private clone polymorphic variants that, at least in the case of L929, promote functional impairment of the oxidative phosphorylation system. Conse quently, the mtDNA of the strain used for the mouse genome project (C57BL/6J) is proposed as the new standard for the mouse mtDNA sequence.

Most inbred laboratory mouse strains are known to have originated from a mixed but limited founder population in a few laboratories. However, the effect of this breeding history on patterns of genetic variation among these strains and the implications for their use are not well understood. Here we present an analysis of the fine structure of variation in the mouse genome, using single nucleotide polymorphisms (SNPs). When the recently assembled genome sequence from the C57BL/6J strain is aligned with sample sequence from other strains, we observe long segments of either extremely high (approximately 40 SNPs per 10 kb) or extremely low (approximately 0.5 SNPs per 10 kb) polymorphism rates. In all strain-to-strain comparisons examined, only one-third of the genome falls into long regions (averaging >1 Mb) of a high SNP rate, consistent with estimated divergence rates between Mus musculus domesticus and either M. m. musculus or M. m. castaneus. These data suggest that the genomes of these inbred strains are mosaics with the vast majority of segments derived from domesticus and musculus sources. These observations have important implications for the design and interpretation of positional cloning experiments.

Variation in the human genome sequence is key to understanding susceptibility to disease in modern populations and the history of ancestral populations. Unlocking this information requires knowledge of the patterns and underlying causes of human sequence diversity. By applying a new population-genetic framework to two genome-wide polymorphism surveys, we find that the human genome contains sizeable regions (stretching over tens of thousands of base pairs) that have intrinsically high and low rates of sequence variation. We show that the primary determinant of these patterns is shared genealogical history. Only a fraction of the variation (at most 25%) is due to the local mutation rate. By measuring the average distance over which genealogical histories are typically preserved, these data provide the first genome-wide estimate of the average extent of correlation among variants (linkage disequilibrium). The results are best explained by extreme variability in the recombination rate at a fine scale, and provide the first empirical evidence that such recombination 'hot spots' are a general feature of the human genome and have a principal role in shaping genetic variation in the human population.

BACKGROUND: Neurofibromatosis type 1 (NF1) is a genetic tumor predisposition disorder caused by germline mutations in tumor suppressor NF1. Plexiform neurofibromas (PN) are benign tumors that arise prenatally or early in childhood and affect 30-50% of NF1 population. Somatic inactivation of second copy of NF1 is believed to be primary genetic event leading to PN initiation. NF1 patients have 8-12% lifetime risk of developing malignant peripheral nerve sheath tumor (MPNST), a highly aggressive soft tissue sarcoma, often arising from pre-existing PN and atypical NF (ANF). ANF are pre-malignant tumors that often arise within PN and can transform into MPNST. They are distinct from both PN and MPNST clinically and histologically, thus representing an intermediate step in malignant transformation. Several studies identified deletion of the CDKN2A/2B locus as the most frequent genetic event in ANF, however it is not clear whether other genes or pathways play role in PN transformation into ANF and further into MPNST. In this study, we performed genomic analysis of 16 ANF and 4 MPNST matched with normal DNA obtained from 14 and 4 patients, respectively. METHODS: We performed whole exome sequencing and whole transcriptome RNASeq analyses on Illumina Hi-Seq 2500 platform and copy-number variant (CNV) analysis on Illumina HumanOmniExpressExome-8 SNP-arrays. In addition, we performed deep sequencing of NF1 and validation of select mutations on IonTorrent platform. For select tumors we estimated growth rate and metabolic activity by using volumetric MRI and FDG-PET. RESULTS: We identified inactivation of NF1 in the majority of ANF and all MPNST. We also detected CDKN2A/B locus deletion in the majority of ANF and MPNST (heterozygous in ANF and mostly homozygous in MPNST). We determined that PRC2 genes (EED and SUZ12) were mutated in multiple MPNST but never in ANF. We identified a low number of point mutations and small indels in the genomes of ANF (median 1, range 0-4) and somewhat elevated mutation burden in MPNST (median 23, range 18-31), however none of these mutations were recurrent and none of the mutant genes (other than NF1 and CDKN2A) were present in multiple samples. We found 93 CNV per tumor (median) in ANF that constituted ~2% of their genomes. In comparison, we observed 2,249 CNV (median) in MPNST that comprised ~75% of their genomes. We didn’t detect significant correlation between growth rate or metabolic activity and the degree of genomic instability or mutation burden in the tumors, however the size of the sample set was modest. RNAseq data analysis is pending. CONCLUSIONS: It appears that PN-ANF transition is predominantly if not exclusively driven by heterozygous deletion of the CDKN2A/2B locus. Further progression to MPNST likely involves homozygous loss of CDKN2A/B and complete inactivation of the PRC2 complex. Widespread LOH in MPNST may accelerate inactivation of key gatekeepers. Citation Format: Alexander Pemov, Nancy F. Hansen, Rajesh Patidar, Christine Higham, Eva Dombi, Joseph F. Boland, Settara C. Chandrasekharappa, NIH Intramural Sequencing Center, James C. Mullikin, Margaret Wallace, Javed Khan, Eric Legius, Brigitte Widemann, Douglas R. Stewart. Identifying early genetic steps in malignant transformation of neurofibromatosis type 1- associated plexiform neurofibromas [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 1456. doi:10.1158/1538-7445.AM2017-1456

While genomic sequencing methods are powerful tools in the discovery of the genetic underpinnings of human disease, incidentally-revealed novel genomic risk factors may be equally important, both scientifically, and as relates to direct patient care. We performed whole-exome sequencing on a child with VACTERL association who suffered severe post-surgical neonatal pulmonary hypertension, and identified a potential novel genetic risk factor for this complication: a heterozygous mutation in CPSI. Newborn screening results from this patient’s monozygotic twin provided evidence that this mutation, in combination with an environmental trigger (in this case, surgery), may have resulted in pulmonary artery hypertension due to inadequate nitric oxide production. Identification of this genetic risk factor allows for targeted medical preventative measures in this patient as well as relatives with the same mutation, and illustrates the power of incidental medical information unearthed by whole-exome sequencing.

Background The Hawaiian honeycreepers are an avian adaptive radiation containing many endangered and extinct species. They display a dramatic range of phenotypic variation and are a model system for studies of evolution, conservation, disease dynamics and population genetics. Development of a genome-scale resources for this group would augment the quality of research focusing on Hawaiian honeycreepers and facilitate comparative avian genomic research. Results We assembled the genome sequence of a Hawaii amakihi (Hemignathus virens),and identified ~3.9 million single nucleotide polymorphisms (SNPs) in the genome. Using the amakihi genome as a reference, we also identified ~156,000 SNPs in RAD tag (restriction site associated DNA) sequencing of five honeycreeper species (palila [Loxioides bailleui], Nihoa finch [Telespiza ultima], iiwi [Vestiaria coccinea], apapane [Himatione sanguinea], and amakihi). SNPs are distributed throughout the amakihi genome, and the individual sequenced shows several large regions of low heterozygosity on chromosomes 1, 5, 6, 8 and 11. SNPs from RAD tag sequencing were also found throughout the genome but were found to be more densely located on microchromosomes, apparently a result of differential distribution of the particular site recognized by restriction enzyme BseXI. Conclusions The amakihi genome sequence will be useful for comparative avian genomics research and provides a significant resource for studies in such areas as disease ecology, evolution, and conservation genetics. The genome sequences will enable mapping of transcriptome data for honeycreepers and comparison of gene sequences between avian taxa. Researchers will be able to use the large number of SNP markers to genotype honeycreepers in regions of interest or across the whole genome. There are enough markers to enable use of methods such as genome-wide association studies (GWAS) that will allow researchers to make connections between phenotypic diversity of honeycreepers and specific genetic variants. Genome-wide markers will also help resolve phylogenetic and population genetic questions in honeycreepers. Electronic supplementary material The online version of this article (doi:10.1186/1471-2164-15-1098) contains supplementary material, which is available to authorized users.

Neurofibromatosis type 1 (NF1) is an autosomal dominant, monogenic disorder of dysregulated neurocutaneous tissue growth. Pleiotropy, variable expressivity and few NF1 genotype-phenotype correlates limit clinical prognostication in NF1. Phenotype complexity in NF1 is hypothesized to derive in part from genetic modifiers unlinked to the NF1 locus. In this study, we hypothesized that normal variation in germline gene expression confers risk for certain phenotypes in NF1. In a set of 79 individuals with NF1, we examined the association between gene expression in lymphoblastoid cell lines with NF1-associated phenotypes and sequenced select genes with significant phenotype/expression correlations. In a discovery cohort of 89 self-reported European-Americans with NF1 we examined the association between germline sequence variants of these genes with café-au-lait macule (CALM) count, a tractable, tumor-like phenotype in NF1. Two correlated, common SNPs (rs4660761 and rs7161) between DPH2 and ATP6V0B were significantly associated with the CALM count. Analysis with tiled regression also identified SNP rs4660761 as significantly associated with CALM count. SNP rs1800934 and 12 rare variants in the mismatch repair gene MSH6 were also associated with CALM count. Both SNPs rs7161 and rs4660761 (DPH2 and ATP6V0B) were highly significant in a mega-analysis in a combined cohort of 180 self-reported European-Americans; SNP rs1800934 (MSH6) was near-significant in a meta-analysis assuming dominant effect of the minor allele. SNP rs4660761 is predicted to regulate ATP6V0B, a gene associated with melanosome biology. Individuals with homozygous mutations in MSH6 can develop an NF1-like phenotype, including multiple CALMs. Through a multi-platform approach, we identified variants that influence NF1 CALM count., Author Summary Neurofibromatosis type 1 (NF1) is a relatively common genetic disease that increases the chance to develop a variety of benign and malignant tumors. People with NF1 also typically feature a large number of birthmarks called café-au-lait macules. It is difficult to predict severity or specific problems in NF1. We sought to identify genes (other than NF1, the gene that causes the disease) that influence severity in NF1. We determined the number of café-au-lait macules in two groups of people with NF1. We measured the gene expression of about 10,000 genes in the cultured white blood cells from one group of people. We then sequenced a group of genes whose expression level was increased in people with higher numbers of café-au-lait macules. In the first group, we found common variants in genes MSH6 and near DPH2 and ATP6V0B that were significantly associated with the number of café-au-lait macules. Some of these variants were close to significant in the second group of people. The two variants near DPH2 and ATP6V0B were very significant when analysed in both groups combined. Our work is among the first to identify genetic variants that influence the severity of NF1.

Public health officials have raised concerns that plasmid transfer between Enterobacteriaceae species may spread resistance to carbapenems, an antibiotic class of last resort, thereby rendering common health care–associated infections nearly impossible to treat. To determine the diversity of carbapenemase-encoding plasmids and assess their mobility among bacterial species, we performed comprehensive surveillance and genomic sequencing of carbapenem-resistant Enterobacteriaceae in the National Institutes of Health (NIH) Clinical Center patient population and hospital environment. We isolated a repertoire of carbapenemase-encoding Enterobacteriaceae, including multiple strains of Klebsiella pneumoniae , Klebsiella oxytoca , Escherichia coli , Enterobacter cloacae , Citrobacter freundii , and Pantoea species. Long-read genome sequencing with full end-to-end assembly revealed that these organisms carry the carbapenem resistance genes on a wide array of plasmids. K. pneumoniae and E. cloacae isolated simultaneously from a single patient harbored two different carbapenemase-encoding plasmids, indicating that plasmid transfer between organisms was unlikely within this patient. We did, however, find evidence of horizontal transfer of carbapenemase-encoding plasmids between K. pneumoniae , E. cloacae , and C. freundii in the hospital environment. Our data, including full plasmid identification, challenge assumptions about horizontal gene transfer events within patients and identify possible connections between patients and the hospital environment. In addition, we identified a new carbapenemase-encoding plasmid of potentially high clinical impact carried by K. pneumoniae , E. coli , E. cloacae , and Pantoea species, in unrelated patients and in the hospital environment.

Background Endometrial cancer (EC) is the 8th leading cause of cancer death amongst American women. Most ECs are endometrioid, serous, or clear cell carcinomas, or an admixture of histologies. Serous and clear ECs are clinically aggressive tumors for which alternative therapeutic approaches are needed. The purpose of this study was to search for somatic mutations in the tyrosine kinome of serous and clear cell ECs, because mutated kinases can point to potential therapeutic targets. Methods In a mutation discovery screen, we PCR amplified and Sanger sequenced the exons encoding the catalytic domains of 86 tyrosine kinases from 24 serous, 11 clear cell, and 5 mixed histology ECs. For somatically mutated genes, we next sequenced the remaining coding exons from the 40 discovery screen tumors and sequenced all coding exons from another 72 ECs (10 clear cell, 21 serous, 41 endometrioid). We assessed the copy number of mutated kinases in this cohort of 112 tumors using quantitative real time PCR, and we used immunoblotting to measure expression of these kinases in endometrial cancer cell lines. Results Overall, we identified somatic mutations in TNK2 (tyrosine kinase non-receptor, 2) and DDR1 (discoidin domain receptor tyrosine kinase 1) in 5.3% (6 of 112) and 2.7% (3 of 112) of ECs. Copy number gains of TNK2 and DDR1 were identified in another 4.5% and 0.9% of 112 cases respectively. Immunoblotting confirmed TNK2 and DDR1 expression in endometrial cancer cell lines. Three of five missense mutations in TNK2 and one of two missense mutations in DDR1 are predicted to impact protein function by two or more in silico algorithms. The TNK2P761Rfs*72 frameshift mutation was recurrent in EC, and the DDR1R570Q missense mutation was recurrent across tumor types. Conclusions This is the first study to systematically search for mutations in the tyrosine kinome in clear cell endometrial tumors. Our findings indicate that high-frequency somatic mutations in the catalytic domains of the tyrosine kinome are rare in clear cell ECs. We uncovered ten new mutations in TNK2 and DDR1 within serous and endometrioid ECs, thus providing novel insights into the mutation spectrum of each gene in EC. Electronic supplementary material The online version of this article (doi:10.1186/1471-2407-14-884) contains supplementary material, which is available to authorized users.

We describe an algorithm, SSAHA (SequenceSearch and Alignment by HashingAlgorithm), for performing fast searches on databases containing multiple gigabases of DNA. Sequences in the database are preprocessed by breaking them into consecutive k-tuples ofk contiguous bases and then using a hash table to store the position of each occurrence of each k-tuple. Searching for a query sequence in the database is done by obtaining from the hash table the “hits” for each k-tuple in the query sequence and then performing a sort on the results. We discuss the effect of the tuple length k on the search speed, memory usage, and sensitivity of the algorithm and present the results of computational experiments which show that SSAHA can be three to four orders of magnitude faster than BLAST or FASTA, while requiring less memory than suffix tree methods. The SSAHAalgorithm is used for high-throughput single nucleotide polymorphism (SNP) detection and very large scale sequence assembly. Also, it provides Web-based sequence search facilities for Ensembl projects.

An understanding of ctenophore biology is critical for reconstructing events that occurred early in animal evolution. Toward this goal, we have sequenced, assembled, and annotated the genome of the ctenophore Mnemiopsis leidyi. Our phylogenomic analyses of both amino acid positions and gene content suggest that ctenophores rather than sponges are the sister lineage to all other animals. Mnemiopsis lacks many of the genes found in bilaterian mesodermal cell types, suggesting that these cell types evolved independently. The set of neural genes in Mnemiopsis is similar to that of sponges, indicating that sponges may have lost a nervous system. These results present a newly supported view of early animal evolution that accounts for major losses and/or gains of sophisticated cell types, including nerve and muscle cells.

Next-generation sequencing of antibody transcripts provides a wealth of data, but the ability to identify function-specific antibodies solely on the basis of sequence has remained elusive. We previously characterized the VRC01 class of antibodies, which target the CD4-binding site on gp120, appear in multiple donors, and broadly neutralize HIV-1. Antibodies of this class have developmental commonalities, but typically share only ∼50% amino acid sequence identity among different donors. Here we apply next-generation sequencing to identify VRC01 class antibodies in a new donor, C38, directly from B cell transcript sequences. We first tested a lineage rank approach, but this was unsuccessful, likely because VRC01 class antibody sequences were not highly prevalent in this donor. We next identified VRC01 class heavy chains through a phylogenetic analysis that included thousands of sequences from C38 and a few known VRC01 class sequences from other donors. This "cross-donor analysis" yielded heavy chains with little sequence homology to previously identified VRC01 class heavy chains. Nonetheless, when reconstituted with the light chain from VRC01, half of the heavy chain chimeric antibodies showed substantial neutralization potency and breadth. We then identified VRC01 class light chains through a five-amino-acid sequence motif necessary for VRC01 light chain recognition. From over a million light chain sequences, we identified 13 candidate VRC01 class members. Pairing of these light chains with the phylogenetically identified C38 heavy chains yielded functional antibodies that effectively neutralized HIV-1. Bioinformatics analysis can thus directly identify functional HIV-1-neutralizing antibodies of the VRC01 class from a sequenced antibody repertoire.

We present a high-quality genome sequence of a Neanderthal woman from Siberia. We show that her parents were related at the level of half-siblings and that mating among close relatives was common among her recent ancestors. We also sequenced the genome of a Neanderthal from the Caucasus to low coverage. An analysis of the relationships and population history of available archaic genomes and 25 present-day human genomes shows that several gene flow events occurred among Neanderthals, Denisovans and early modern humans, possibly including gene flow into Denisovans from an unknown archaic group. Thus, interbreeding, albeit of low magnitude, occurred among many hominin groups in the Late Pleistocene. In addition, the high-quality Neanderthal genome allows us to establish a definitive list of substitutions that became fixed in modern humans after their separation from the ancestors of Neanderthals and Denisovans.

Although a considerable proportion of serum lipids loci identified in European ancestry individuals (EA) replicate in African Americans (AA), interethnic differences in the distribution of serum lipids suggest that some genetic determinants differ by ethnicity. We conducted a comprehensive evaluation of five lipid candidate genes to identify variants with ethnicity-specific effects. We sequenced ABCA1, LCAT, LPL, PON1, and SERPINE1 in 48 AA individuals with extreme serum lipid concentrations (high HDLC/low TG or low HDLC/high TG). Identified variants were genotyped in the full population-based sample of AA (n = 1694) and tested for an association with serum lipids. rs328 (LPL) and correlated variants were associated with higher HDLC and lower TG. Interestingly, a stronger effect was observed on a “European” vs. “African” genetic background at this locus. To investigate this effect, we evaluated the region among West Africans (WA). For TG, the effect size among WA was the same in AA with only African local ancestry (2–3% lower TG), while the larger association among AA with local European ancestry matched previous reports in EA (10%). For HDLC, there was no association with rs328 in AA with only African local ancestry or in WA, while the association among AA with European local ancestry was much greater than what has been observed for EA (15 vs. ∼5 mg/dl), suggesting an interaction with an environmental or genetic factor that differs by ethnicity. Beyond this ancestry effect, the importance of African ancestry-focused, sequence-based work was also highlighted by serum lipid associations of variants that were in higher frequency (or present only) among those of African ancestry. By beginning our study with the sequence variation present in AA individuals, investigating local ancestry effects, and seeking replication in WA, we were able to comprehensively evaluate the role of a set of candidate genes in serum lipids in AA., Author Summary Most of the work on the genetic epidemiology of serum lipids in African Americans (AA) has focused on replicating findings that were identified in European ancestry individuals. While this can be very informative about the generalizability of lipids loci across populations, African ancestry-specific variation will be missed using this approach. Our aim was to comprehensively evaluate five lipid candidate genes in an AA population, from the identification of variants of interest to population-level analysis of high-density lipoprotein cholesterol (HDLC) and triglycerides (TG). We sequenced five genes in individuals with extreme lipids (n = 48) drawn from a population-based study of AA. The variants identified were genotyped in 1,694 AA and analyzed. Notable among the findings were the observation of ancestry specific effect for several variants in the LPL gene among these admixed individuals, with a greater effect observed among those with European ancestry in this region. These associations were further elucidated by replication in West Africans. By beginning with the sequence variation present among AA, investigating ancestry effects, and seeking replication in West Africans, we were able to comprehensively evaluate these candidate genes with a focus on African ancestry individuals.

Genome-wide association studies (GWAS) are a powerful means of identifying genes with disease-associated common variants, but they are not well-suited to detecting genes with disease-associated rare and low-frequency variants. In the current study of Behçet disease (BD), nonsynonymous variants (NSVs) identified by deep exonic resequencing of 10 genes found by GWAS ( IL10, IL23R , CCR1 , STAT4 , KLRK1 , KLRC1, KLRC2, KLRC3, KLRC4 , and ERAP1 ) and 11 genes selected for their role in innate immunity (IL1B, IL1R1, IL1RN, NLRP3, MEFV, TNFRSF1A, PSTPIP1, CASP1, PYCARD, NOD2, and TLR4 ) were evaluated for BD association. A differential distribution of the rare and low-frequency NSVs of a gene in 2,461 BD cases compared with 2,458 controls indicated their collective association with disease. By stringent criteria requiring at least a single burden test with study-wide significance and a corroborating test with at least nominal significance, rare and low-frequency NSVs in one GWAS-identified gene, IL23R ( P = 6.9 × 10 −5 ), and one gene involved in innate immunity, TLR4 ( P = 8.0 × 10 −4 ), were associated with BD. In addition, damaging or rare damaging NOD2 variants were nominally significant across all three burden tests applied ( P = 0.0063–0.045). Furthermore, carriage of the familial Mediterranean fever gene ( MEFV ) mutation Met694Val, which is known to cause recessively inherited familial Mediterranean fever, conferred BD risk in the Turkish population (OR, 2.65; P = 1.8 × 10 −12 ). The disease-associated NSVs in MEFV and TLR4 implicate innate immune and bacterial sensing mechanisms in BD pathogenesis.

Core duplicons in the human genome represent ancestral duplication modules shared by the majority of intrachromosomal duplication blocks within a given chromosome. These cores are associated with the emergence of novel gene families in the hominoid lineage, but their genomic organization and gene characterization among other primates are largely unknown. Here, we investigate the genomic organization and expression of the core duplicon on chromosome 17 that led to the expansion of LRRC37 during primate evolution. A comparison of the LRRC37 gene family organization in human, orangutan, macaque, marmoset, and lemur genomes shows the presence of both orthologous and species-specific gene copies in all primate lineages. Expression profiling in mouse, macaque, and human tissues reveals that the ancestral expression of LRRC37 was restricted to the testis. In the hominid lineage, the pattern of LRRC37 became increasingly ubiquitous, with significantly higher levels of expression in the cerebellum and thymus, and showed a remarkable diversity of alternative splice forms. Transfection studies in HeLa cells indicate that the human FLAG-tagged recombinant LRRC37 protein is secreted after cleavage of a transmembrane precursor and its overexpression can induce filipodia formation., This work was supported by National Institutes of Health (NIH) grants HG002385 and GM058815 to E.E.E. and in part by the Intramural Research Program of the National Human Genome Research Institute, NIH. E.E.E. is an investigator of the Howard Hughes Medical Institute.

Antibody 10E8 targets the membrane-proximal external region (MPER) of HIV-1 gp41, neutralizes >97% of HIV-1 isolates, and lacks the auto-reactivity often associated with MPER-directed antibodies. The developmental pathway of 10E8 might therefore serve as a promising template for vaccine design, but samples from time-of-infection-often used to infer the B cell record-are unavailable. In this study, we used crystallography, next-generation sequencing (NGS), and functional assessments to infer the 10E8 developmental pathway from a single time point. Mutational analysis indicated somatic hypermutation of the 2nd-heavy chain-complementarity determining region (CDR H2) to be critical for neutralization, and structures of 10E8 variants with V-gene regions reverted to genomic origin for heavy-and-light chains or heavy chain-only showed structural differences >2 Å relative to mature 10E8 in the CDR H2 and H3. To understand these developmental changes, we used bioinformatic sieving, maximum likelihood, and parsimony analyses of immunoglobulin transcripts to identify 10E8-lineage members, to infer the 10E8-unmutated common ancestor (UCA), and to calculate 10E8-developmental intermediates. We were assisted in this analysis by the preservation of a critical D-gene segment, which was unmutated in most 10E8-lineage sequences. UCA and early intermediates weakly bound a 26-residue-MPER peptide, whereas HIV-1 neutralization and epitope recognition in liposomes were only observed with late intermediates. Antibody 10E8 thus develops from a UCA with weak MPER affinity and substantial differences in CDR H2 and H3 from the mature 10E8; only after extensive somatic hypermutation do 10E8-lineage members gain recognition in the context of membrane and HIV-1 neutralization.

Although RNA-Seq data provide unprecedented isoform-level expression information, detection of alternative isoform regulation (AIR) remains difficult, particularly when working with an incomplete transcript annotation. We introduce JunctionSeq, a new method that builds on the statistical techniques used by the well-established DEXSeq package to detect differential usage of both exonic regions and splice junctions. In particular, JunctionSeq is capable of detecting differential usage of novel splice junctions without the need for an additional isoform assembly step, greatly improving performance when the available transcript annotation is flawed or incomplete. JunctionSeq also provides a powerful and streamlined visualization toolset that allows bioinformaticians to quickly and intuitively interpret their results. We tested our method on publicly available data from several experiments performed on the rat pineal gland and Toxoplasma gondii, successfully detecting known and previously validated AIR genes in 19 out of 19 gene-level hypothesis tests. Due to its ability to query novel splice sites, JunctionSeq is still able to detect these differences even when all alternative isoforms for these genes were not included in the transcript annotation. JunctionSeq thus provides a powerful method for detecting alternative isoform regulation even with low-quality annotations. An implementation of JunctionSeq is available as an R/Bioconductor package.

Early-onset myopathy, areflexia, respiratory distress and dysphagia (EMARDD) is a myopathic disorder associated with mutations in MEGF10. By novel analysis of SNP array hybridization and exome sequence coverage, we diagnosed a 10-years old girl with EMARDD following identification of a novel homozygous deletion of exon 7 in MEGF10. In contrast to previously reported EMARDD patients, her weakness was more prominent proximally than distally, and involved her legs more than her arms. MRI of her pelvis and thighs showed muscle atrophy and fatty replacement. Ultrasound of several muscle groups revealed dense homogenous increases in echogenicity. Cloning and sequencing of the deletion breakpoint identified features suggesting the mutation arose by fork stalling and template switching. These findings constitute the first genomic deletion causing EMARDD, expand the clinical phenotype, and provide new insight into the pattern and histology of its muscular pathology.