ABSTRACT
Immunity that controls parasitemia and inflammation during Plasmodium falciparum (Pf) malaria can be acquired with repeated infections. A limited understanding of this complex immune response impedes the development of vaccines and adjunctive therapies. We conducted a prospective systems biology study of children who differed in their ability to control parasitemia and fever following Pf infection. By integrating whole-blood transcriptomics, flow-cytometric analysis, and plasma cytokine and antibody profiles, we demonstrate that a pre-infection signature of B cell enrichment, upregulation of T helper type 1 (Th1) and Th2 cell-associated pathways, including interferon responses, and p53 activation associated with control of malarial fever and coordinated with Pf-specific immunoglobulin G (IgG) and Fc receptor activation to control parasitemia. Our hypothesis-generating approach identified host molecules that may contribute to differential clinical outcomes during Pf infection. As a proof of concept, we have shown that enhanced p53 expression in monocytes attenuated Plasmodium-induced inflammation and predicted protection from fever.
Subject(s)
B-Lymphocytes/immunology , Blood Proteins/metabolism , Inflammation/metabolism , Malaria, Falciparum/metabolism , Plasmodium falciparum/physiology , Th1 Cells/immunology , Th2 Cells/immunology , Tumor Suppressor Protein p53/metabolism , Adolescent , Adult , Animals , Antibodies, Protozoan/metabolism , Child , Child, Preschool , Disease Resistance , Female , Gene Expression Profiling , Humans , Infant , Interferons/metabolism , Male , Mice , Mice, Inbred C57BL , Prospective Studies , Receptors, Fc/metabolism , Signal Transduction , Tumor Suppressor Protein p53/genetics , Young AdultABSTRACT
Genome sequencing has established clinical utility for rare disease diagnosis. While increasing numbers of individuals have undergone elective genome sequencing, a comprehensive study surveying genome-wide disease-associated genes in adults with deep phenotyping has not been reported. Here we report the results of a 3-y precision medicine study with a goal to integrate whole-genome sequencing with deep phenotyping. A cohort of 1,190 adult participants (402 female [33.8%]; mean age, 54 y [range 20 to 89+]; 70.6% European) had whole-genome sequencing, and were deeply phenotyped using metabolomics, advanced imaging, and clinical laboratory tests in addition to family/medical history. Of 1,190 adults, 206 (17.3%) had at least 1 genetic variant with pathogenic (P) or likely pathogenic (LP) assessment that suggests a predisposition of genetic risk. A multidisciplinary clinical team reviewed all reportable findings for the assessment of genotype and phenotype associations, and 137 (11.5%) had genotype and phenotype associations. A high percentage of genotype and phenotype associations (>75%) was observed for dyslipidemia (n = 24), cardiomyopathy, arrhythmia, and other cardiac diseases (n = 42), and diabetes and endocrine diseases (n = 17). A lack of genotype and phenotype associations, a potential burden for patient care, was observed in 69 (5.8%) individuals with P/LP variants. Genomics and metabolomics associations identified 61 (5.1%) heterozygotes with phenotype manifestations affecting serum metabolite levels in amino acid, lipid and cofactor, and vitamin pathways. Our descriptive analysis provides results on the integration of whole-genome sequencing and deep phenotyping for clinical assessments in adults.
Subject(s)
Diagnostic Imaging , Metabolomics , Precision Medicine/methods , Whole Genome Sequencing , Adult , Aged , Aged, 80 and over , Cohort Studies , Female , Genetic Predisposition to Disease/genetics , Genotype , Heart Diseases/genetics , Humans , Male , Middle Aged , Phenotype , Young AdultABSTRACT
Short tandem repeats (STRs) are hyper-mutable sequences in the human genome. They are often used in forensics and population genetics and are also the underlying cause of many genetic diseases. There are challenges associated with accurately determining the length polymorphism of STR loci in the genome by next-generation sequencing (NGS). In particular, accurate detection of pathological STR expansion is limited by the sequence read length during whole-genome analysis. We developed TREDPARSE, a software package that incorporates various cues from read alignment and paired-end distance distribution, as well as a sequence stutter model, in a probabilistic framework to infer repeat sizes for genetic loci, and we used this software to infer repeat sizes for 30 known disease loci. Using simulated data, we show that TREDPARSE outperforms other available software. We sampled the full genome sequences of 12,632 individuals to an average read depth of approximately 30× to 40× with Illumina HiSeq X. We identified 138 individuals with risk alleles at 15 STR disease loci. We validated a representative subset of the samples (n = 19) by Sanger and by Oxford Nanopore sequencing. Additionally, we validated the STR calls against known allele sizes in a set of GeT-RM reference cell-line materials (n = 6). Several STR loci that are entirely guanine or cytosines (G or C) have insufficient read evidence for inference and therefore could not be assayed precisely by TREDPARSE. TREDPARSE extends the limit of STR size detection beyond the physical sequence read length. This extension is critical because many of the disease risk cutoffs are close to or beyond the short sequence read length of 100 to 150 bases.
Subject(s)
Genome, Human/genetics , Microsatellite Repeats/genetics , Adolescent , Adult , Alleles , Child , Female , Genetics, Population/methods , High-Throughput Nucleotide Sequencing/methods , Humans , Male , Middle Aged , Polymorphism, Genetic/genetics , Sequence Analysis, DNA/methods , SoftwareABSTRACT
The HLA gene complex on human chromosome 6 is one of the most polymorphic regions in the human genome and contributes in large part to the diversity of the immune system. Accurate typing of HLA genes with short-read sequencing data has historically been difficult due to the sequence similarity between the polymorphic alleles. Here, we introduce an algorithm, xHLA, that iteratively refines the mapping results at the amino acid level to achieve 99-100% four-digit typing accuracy for both class I and II HLA genes, taking only [Formula: see text]3 min to process a 30× whole-genome BAM file on a desktop computer.
Subject(s)
Histocompatibility Testing/methods , Algorithms , Benchmarking , HumansABSTRACT
We report on the sequencing of 10,545 human genomes at 30×-40× coverage with an emphasis on quality metrics and novel variant and sequence discovery. We find that 84% of an individual human genome can be sequenced confidently. This high-confidence region includes 91.5% of exon sequence and 95.2% of known pathogenic variant positions. We present the distribution of over 150 million single-nucleotide variants in the coding and noncoding genome. Each newly sequenced genome contributes an average of 8,579 novel variants. In addition, each genome carries on average 0.7 Mb of sequence that is not found in the main build of the hg38 reference genome. The density of this catalog of variation allowed us to construct high-resolution profiles that define genomic sites that are highly intolerant of genetic variation. These results indicate that the data generated by deep genome sequencing is of the quality necessary for clinical use.
Subject(s)
Genome, Human , Genomics , Whole Genome Sequencing , Chromosome Mapping , Computational Biology/methods , Databases, Nucleic Acid , Genetic Predisposition to Disease , Genetic Variation , Genomics/methods , Humans , Open Reading Frames , Polymorphism, Single Nucleotide , Reproducibility of Results , Untranslated RegionsABSTRACT
Human high-altitude (HA) adaptation or mal-adaptation is explored to understand the physiology, pathophysiology, and molecular mechanisms that underlie long-term exposure to hypoxia. Here, we report the results of an analysis of the largest whole-genome-sequencing of Chronic Mountain Sickness (CMS) and nonCMS individuals, identified candidate genes and functionally validated these candidates in a genetic model system (Drosophila). We used PreCIOSS algorithm that uses Haplotype Allele Frequency score to separate haplotypes carrying the favored allele from the noncarriers and accordingly, prioritize genes associated with the CMS or nonCMS phenotype. Haplotypes in eleven candidate regions, with SNPs mostly in nonexonic regions, were significantly different between CMS and nonCMS subjects. Closer examination of individual genes in these regions revealed the involvement of previously identified candidates (e.g., SENP1) and also unreported ones SGK3, COPS5, PRDM1, and IFT122 in CMS. Remarkably, in addition to genes like SENP1, SGK3, and COPS5 which are HIF-dependent, our study reveals for the first time HIF-independent gene PRDM1, indicating an involvement of wider, nonHIF pathways in HA adaptation. Finally, we observed that down-regulating orthologs of these genes in Drosophila significantly enhanced their hypoxia tolerance. Taken together, the PreCIOSS algorithm, applied on a large number of genomes, identifies the involvement of both new and previously reported genes in selection sweeps, highlighting the involvement of multiple hypoxia response systems. Since the overwhelming majority of SNPs are in nonexonic (and possibly regulatory) regions, we speculate that adaptation to HA necessitates greater genetic flexibility allowing for transcript variability in response to graded levels of hypoxia.
Subject(s)
Acclimatization/genetics , Altitude Sickness/genetics , Adaptation, Physiological/genetics , Adult , Alleles , Altitude , Altitude Sickness/metabolism , Altitude Sickness/physiopathology , Animals , Chronic Disease , Drosophila/genetics , Evolution, Molecular , Gene Frequency/genetics , Haplotypes/genetics , Humans , Hypoxia/genetics , Hypoxia/physiopathology , Male , Peru , Polymorphism, Single Nucleotide/genetics , Positive Regulatory Domain I-Binding Factor 1/genetics , Positive Regulatory Domain I-Binding Factor 1/metabolism , Whole Genome Sequencing/methodsABSTRACT
Defined transcription factors can induce epigenetic reprogramming of adult mammalian cells into induced pluripotent stem cells. Although DNA factors are integrated during some reprogramming methods, it is unknown whether the genome remains unchanged at the single nucleotide level. Here we show that 22 human induced pluripotent stem (hiPS) cell lines reprogrammed using five different methods each contained an average of five protein-coding point mutations in the regions sampled (an estimated six protein-coding point mutations per exome). The majority of these mutations were non-synonymous, nonsense or splice variants, and were enriched in genes mutated or having causative effects in cancers. At least half of these reprogramming-associated mutations pre-existed in fibroblast progenitors at low frequencies, whereas the rest occurred during or after reprogramming. Thus, hiPS cells acquire genetic modifications in addition to epigenetic modifications. Extensive genetic screening should become a standard procedure to ensure hiPS cell safety before clinical use.
Subject(s)
Cellular Reprogramming/genetics , Induced Pluripotent Stem Cells/metabolism , Mutagenesis/genetics , Point Mutation/genetics , Cells, Cultured , DNA Mutational Analysis , Epistasis, Genetic/genetics , Fibroblasts/cytology , Fibroblasts/metabolism , Humans , Induced Pluripotent Stem Cells/cytology , Male , Middle Aged , Models, Genetic , Open Reading Frames/geneticsABSTRACT
There is increasing evidence that the phenotypic effects of genomic sequence variants are best understood in terms of variant haplotypes rather than as isolated polymorphisms. Haplotype analysis is also critically important for uncovering population histories and for the study of evolutionary genetics. Although the sequencing of individual human genomes to reveal personal collections of sequence variants is now well established, there has been slower progress in the phasing of these variants into pairs of haplotypes along each pair of chromosomes. Here, we have developed a distinct approach to haplotyping that can yield chromosome-length haplotypes, including the vast majority of heterozygous single-nucleotide polymorphisms (SNPs) in an individual human genome. This approach exploits the haploid nature of sperm cells and employs a combination of genotyping and low-coverage sequencing on a short-read platform. In addition to generating chromosome-length haplotypes, the approach can directly identify recombination events (averaging 1.1 per chromosome) with a median resolution of <100 kb.
Subject(s)
Genome, Human , Haplotypes/genetics , Spermatozoa , Chromosome Mapping , Genotype , Humans , Male , Polymorphism, Single Nucleotide , Sequence Analysis, DNAABSTRACT
The freshwater cnidarian Hydra was first described in 1702 and has been the object of study for 300 years. Experimental studies of Hydra between 1736 and 1744 culminated in the discovery of asexual reproduction of an animal by budding, the first description of regeneration in an animal, and successful transplantation of tissue between animals. Today, Hydra is an important model for studies of axial patterning, stem cell biology and regeneration. Here we report the genome of Hydra magnipapillata and compare it to the genomes of the anthozoan Nematostella vectensis and other animals. The Hydra genome has been shaped by bursts of transposable element expansion, horizontal gene transfer, trans-splicing, and simplification of gene structure and gene content that parallel simplification of the Hydra life cycle. We also report the sequence of the genome of a novel bacterium stably associated with H. magnipapillata. Comparisons of the Hydra genome to the genomes of other animals shed light on the evolution of epithelia, contractile tissues, developmentally regulated transcription factors, the Spemann-Mangold organizer, pluripotency genes and the neuromuscular junction.
Subject(s)
Genome/genetics , Hydra/genetics , Animals , Anthozoa/genetics , Comamonadaceae/genetics , DNA Transposable Elements/genetics , Gene Transfer, Horizontal/genetics , Genome, Bacterial/genetics , Hydra/microbiology , Hydra/ultrastructure , Molecular Sequence Data , Neuromuscular Junction/ultrastructureABSTRACT
BACKGROUND: In humans it is unknown if the composition of the gut microbiota alters the risk of Plasmodium falciparum infection or the risk of developing febrile malaria once P. falciparum infection is established. Here we collected stool samples from a cohort composed of 195 Malian children and adults just prior to an intense P. falciparum transmission season. We assayed these samples using massively parallel sequencing of the 16S ribosomal RNA gene to identify the composition of the gut bacterial communities in these individuals. During the ensuing 6-month P. falciparum transmission season we examined the relationship between the stool microbiota composition of individuals in this cohort and their prospective risk of both P. falciparum infection and febrile malaria. RESULTS: Consistent with prior studies, stool microbial diversity in the present cohort increased with age, although the overall microbiota profile was distinct from cohorts in other regions of Africa, Asia and North America. Age-adjusted Cox regression analysis revealed a significant association between microbiota composition and the prospective risk of P. falciparum infection; however, no relationship was observed between microbiota composition and the risk of developing febrile malaria once P. falciparum infection was established. CONCLUSIONS: These findings underscore the diversity of gut microbiota across geographic regions, and suggest that strategic modulation of gut microbiota composition could decrease the risk of P. falciparum infection in malaria-endemic areas, potentially as an adjunct to partially effective malaria vaccines.
Subject(s)
Bacteria/classification , Feces/microbiology , High-Throughput Nucleotide Sequencing/methods , Malaria, Falciparum/parasitology , Sequence Analysis, RNA/methods , Adolescent , Bacteria/isolation & purification , Child , Child, Preschool , Female , Humans , Infant , Malaria, Falciparum/blood , Malaria, Falciparum/transmission , Male , Mali/epidemiology , Microbiota , Prospective Studies , RNA, Bacterial/analysis , RNA, Ribosomal, 16S/analysis , Risk Factors , Young AdultABSTRACT
UNLABELLED: Rotaviruses (RVs) are leading causes of severe diarrhea and vomiting in infants and young children. RVs with G10P[11] genotype specificity have been associated with symptomatic and asymptomatic neonatal infections in Vellore, India. To identify possible viral genetic determinants responsible for differences in symptomology, the genome sequences of G10P[11] RVs in stool samples of 19 neonates with symptomatic infections and 20 neonates with asymptomatic infections were determined by Sanger and next-generation sequencing. The data showed that all 39 viruses had identical genotype constellations (G10-P[11]-I2-R2-C2-M2-A1-N1-T1-E2-H3), the same as those of the previously characterized symptomatic N155 Vellore isolate. The data also showed that the RNA and deduced protein sequences of all the Vellore G10P[11] viruses were nearly identical; no nucleotide or amino acid differences were found that correlated with symptomatic versus asymptomatic infection. Next-generation sequencing data revealed that some stool samples, both from neonates with symptomatic infections and from neonates with asymptomatic infections, also contained one or more positive-strand RNA viruses (Aichi virus, astrovirus, or salivirus/klassevirus) suspected of being potential causes of pediatric gastroenteritis. However, none of the positive-strand RNA viruses could be causally associated with the development of symptoms. These results indicate that the diversity of clinical symptoms in Vellore neonates does not result from genetic differences among G10P[11] RVs; instead, other undefined factors appear to influence whether neonates develop gastrointestinal disease symptoms. IMPORTANCE: Rotavirus (RV) strains have been identified that preferentially replicate in neonates, in some cases, without causing gastrointestinal disease. Surveillance studies have established that G10P[11] RVs are a major cause of neonatal infection in Vellore, India, with half of infected neonates exhibiting symptoms. We used Sanger and next-generation sequencing technologies to contrast G10P[11] RVs recovered from symptomatic and asymptomatic neonates. Remarkably, the data showed that the RNA genomes of the viruses were virtually indistinguishable and lacked any differences that could explain the diversity of clinical outcomes among infected Vellore neonates. The sequencing results also indicated that some symptomatic and some asymptomatic Vellore neonates were infected with other enteric viruses (Aichi virus, astrovirus, salvirus/klassevirus); however, none could be correlated with the presence of symptoms in neonates. Together, our findings suggest that other poorly defined factors, not connected to the genetic makeup of the Vellore G10P[11] viruses, influence whether neonates develop gastrointestinal disease symptoms.
Subject(s)
Diarrhea/virology , Rotavirus Infections/virology , Rotavirus/genetics , Feces/virology , Gastroenteritis/virology , Genotype , Humans , India , Infant, Newborn , Kobuvirus/geneticsABSTRACT
UNLABELLED: Rotaviruses (RVs) are 11-segmented, double-stranded RNA viruses that cause severe gastroenteritis in children. In addition to an error-prone genome replication mechanism, RVs can increase their genetic diversity by reassorting genes during host coinfection. Such exchanges allow RVs to acquire advantageous genes and adapt in the face of selective pressures. However, reassortment may also impose fitness costs if it unlinks genes/proteins that have accumulated compensatory, coadaptive mutations and that operate best when kept together. To better understand human RV evolutionary dynamics, we analyzed the genome sequences of 135 strains (genotype G1/G3/G4-P[8]-I1-C1-R1-A1-N1-T1-E1-H1) that were collected at a single location in Washington, DC, during the years 1974 to 1991. Intragenotypic phylogenetic trees were constructed for each viral gene using the nucleotide sequences, thereby defining novel allele level gene constellations (GCs) and illuminating putative reassortment events. The results showed that RVs with distinct GCs cocirculated during the vast majority of the collection years and that some of these GCs persisted in the community unchanged by reassortment. To investigate the influence of protein coadaptation on GC maintenance, we performed a mutual information-based analysis of the concatenated amino acid sequences and identified an extensive covariance network. Unexpectedly, amino acid covariation was highest between VP4 and VP2, which are structural components of the RV virion that are not thought to directly interact. These results suggest that GCs may be influenced by the selective constraints placed on functionally coadapted, albeit noninteracting, viral proteins. This work raises important questions about mutation-reassortment interplay and its impact on human RV evolution. IMPORTANCE: Rotaviruses are devastating human pathogens that cause severe diarrhea and kill >450,000 children each year. The virus can evolve by accumulating mutations and by acquiring new genes from other strains via a process called reassortment. However, little is known about the relationship between mutation accumulation and gene reassortment for rotaviruses and how it impacts viral evolution. In this study, we analyzed the genome sequences of human strains found in clinical fecal specimens that were collected at a single hospital over an 18-year time span. We found that many rotaviruses did not reassort their genes but instead maintained them as specific sets (i.e., constellations). By analyzing the encoded proteins, we discovered concurrent amino acid changes among them, which suggests that they are functionally coadapted to operate best when kept together. This study increases our understanding of how rotaviruses evolve over time in the human population.
Subject(s)
Evolution, Molecular , Rotavirus/genetics , Rotavirus/isolation & purification , Viral Proteins/genetics , Adaptation, Biological , Child, Preschool , Cluster Analysis , District of Columbia , Genome, Viral , Humans , Infant , Molecular Sequence Data , Phylogeny , Rotavirus/classification , Sequence Analysis, DNAABSTRACT
Mutations in isocitrate dehydrogenase 1 (IDH1) have been found in the vast majority of low grade and progressive infiltrating gliomas and are characterized by the production of 2-hydroxyglutarate from α-ketoglutarate. Recent investigations of malignant gliomas have identified additional genetic and chromosomal abnormalities which cluster with IDH1 mutations into two distinct subgroups. The astrocytic subgroup was found to have frequent mutations in ATRX, TP53 and displays alternative lengthening of telomeres. The second subgroup with oligodendrocytic morphology has frequent mutations in CIC or FUBP1, and is linked to co-deletion of the 1p/19q arms. These mutations reflect the development of two distinct molecular pathways representing the majority of IDH1 mutant gliomas. Unfortunately, due to the scarcity of endogenously derived IDH1 mutant models, there is a lack of accurate models to study mechanism and develop new therapy. Here we report the generation of an endogenous IDH1 anaplastic astrocytoma in vivo model with concurrent mutations in TP53, CDKN2A and ATRX. The model has a similar phenotype and histopathology as the original patient tumor, expresses the IDH1 (R132H) mutant protein and exhibits an alternative lengthening of telomeres phenotype. The JHH-273 model is characteristic of anaplastic astrocytoma and represents a valuable tool for investigating the pathogenesis of this distinct molecular subset of gliomas and for preclinical testing of compounds targeting IDH1 mutations or alternative lengthening of telomeres.
Subject(s)
Astrocytoma/genetics , Astrocytoma/pathology , Isocitrate Dehydrogenase/genetics , Mutation , Telomere/pathology , Adult , Animals , DNA Helicases/genetics , Disease Models, Animal , Genes, p16 , Heterografts , Humans , Immunohistochemistry , In Situ Hybridization, Fluorescence , Male , Mice , Neoplasm Transplantation/methods , Nuclear Proteins/genetics , Phenotype , Reverse Transcriptase Polymerase Chain Reaction , Tumor Suppressor Protein p53/genetics , X-linked Nuclear ProteinABSTRACT
Group A rotaviruses (RVs) are 11-segmented, double-stranded RNA viruses and are primary causes of gastroenteritis in young children. Despite their medical relevance, the genetic diversity of modern human RVs is poorly understood, and the impact of vaccine use on circulating strains remains unknown. In this study, we report the complete genome sequence analysis of 58 RVs isolated from children with severe diarrhea and/or vomiting at Vanderbilt University Medical Center (VUMC) in Nashville, TN, during the years spanning community vaccine implementation (2005 to 2009). The RVs analyzed include 36 G1P[8], 18 G3P[8], and 4 G12P[8] Wa-like genogroup 1 strains with VP6-VP1-VP2-VP3-NSP1-NSP2-NSP3-NSP4-NSP5/6 genotype constellations of I1-R1-C1-M1-A1-N1-T1-E1-H1. By constructing phylogenetic trees, we identified 2 to 5 subgenotype alleles for each gene. The results show evidence of intragenogroup gene reassortment among the cocirculating strains. However, several isolates from different seasons maintained identical allele constellations, consistent with the notion that certain RV clades persisted in the community. By comparing the genes of VUMC RVs to those of other archival and contemporary RV strains for which sequences are available, we defined phylogenetic lineages and verified that the diversity of the strains analyzed in this study reflects that seen in other regions of the world. Importantly, the VP4 and VP7 proteins encoded by VUMC RVs and other contemporary strains show amino acid changes in or near neutralization domains, which might reflect antigenic drift of the virus. Thus, this large-scale, comparative genomic study of modern human RVs provides significant insight into how this pathogen evolves during its spread in the community.
Subject(s)
Diarrhea/virology , Gastroenteritis/virology , Genetic Variation , Genome, Viral , Rotavirus/genetics , Rotavirus/isolation & purification , Child , Child, Preschool , Female , Genomics , Humans , Male , Molecular Sequence Data , Phylogeny , Rotavirus/classification , Viral Proteins/geneticsABSTRACT
Although patterns of somatic alterations have been reported for tumor genomes, little is known on how they compare with alterations present in non-tumor genomes. A comparison of the two would be crucial to better characterize the genetic alterations driving tumorigenesis. We sequenced the genomes of a lymphoblastoid (HCC1954BL) and a breast tumor (HCC1954) cell line derived from the same patient and compared the somatic alterations present in both. The lymphoblastoid genome presents a comparable number and similar spectrum of nucleotide substitutions to that found in the tumor genome. However, a significant difference in the ratio of non-synonymous to synonymous substitutions was observed between both genomes (P = 0.031). Protein-protein interaction analysis revealed that mutations in the tumor genome preferentially affect hub-genes (P = 0.0017) and are co-selected to present synergistic functions (P < 0.0001). KEGG analysis showed that in the tumor genome most mutated genes were organized into signaling pathways related to tumorigenesis. No such organization or synergy was observed in the lymphoblastoid genome. Our results indicate that endogenous mutagens and replication errors can generate the overall number of mutations required to drive tumorigenesis and that it is the combination rather than the frequency of mutations that is crucial to complete tumorigenic transformation.
Subject(s)
Breast Neoplasms/genetics , Genetic Variation , Genome, Human , Cell Line, Transformed , Cell Line, Tumor , Chromosome Aberrations , Female , Humans , Lymphocytes , Middle Aged , Mutation , Point Mutation , Protein Interaction Mapping , Sequence Analysis, DNAABSTRACT
As an obligatory parasite of humans, the body louse (Pediculus humanus humanus) is an important vector for human diseases, including epidemic typhus, relapsing fever, and trench fever. Here, we present genome sequences of the body louse and its primary bacterial endosymbiont Candidatus Riesia pediculicola. The body louse has the smallest known insect genome, spanning 108 Mb. Despite its status as an obligate parasite, it retains a remarkably complete basal insect repertoire of 10,773 protein-coding genes and 57 microRNAs. Representing hemimetabolous insects, the genome of the body louse thus provides a reference for studies of holometabolous insects. Compared with other insect genomes, the body louse genome contains significantly fewer genes associated with environmental sensing and response, including odorant and gustatory receptors and detoxifying enzymes. The unique architecture of the 18 minicircular mitochondrial chromosomes of the body louse may be linked to the loss of the gene encoding the mitochondrial single-stranded DNA binding protein. The genome of the obligatory louse endosymbiont Candidatus Riesia pediculicola encodes less than 600 genes on a short, linear chromosome and a circular plasmid. The plasmid harbors a unique arrangement of genes required for the synthesis of pantothenate, an essential vitamin deficient in the louse diet. The human body louse, its primary endosymbiont, and the bacterial pathogens that it vectors all possess genomes reduced in size compared with their free-living close relatives. Thus, the body louse genome project offers unique information and tools to use in advancing understanding of coevolution among vectors, symbionts, and pathogens.
Subject(s)
Genome, Bacterial/genetics , Genome, Insect/genetics , Pediculus/genetics , Pediculus/microbiology , Animals , Enterobacteriaceae/genetics , Genes, Bacterial/genetics , Genes, Insect/genetics , Genomics/methods , Humans , Lice Infestations/parasitology , Molecular Sequence Data , Sequence Analysis, DNA , SymbiosisABSTRACT
Structural variants (SVs) are common in the human genome. Because approximately half of the human genome consists of repetitive, transposable DNA sequences, it is plausible that these elements play an important role in generating SVs in humans. Sequencing of the diploid genome of one individual human (HuRef) affords us the opportunity to assess, for the first time, the impact of mobile elements on SVs in an individual in a thorough and unbiased fashion. In this study, we systematically evaluated more than 8000 SVs to identify mobile element-associated SVs as small as 100 bp and specific to the HuRef genome. Combining computational and experimental analyses, we identified and validated 706 mobile element insertion events (including Alu, L1, SVA elements, and nonclassical insertions), which added more than 305 kb of new DNA sequence to the HuRef genome compared with the Human Genome Project (HGP) reference sequence (hg18). We also identified 140 mobile element-associated deletions, which removed approximately 126 kb of sequence from the HuRef genome. Overall, approximately 10% of the HuRef-specific indels larger than 100 bp are caused by mobile element-associated events. More than one-third of the insertion/deletion events occurred in genic regions, and new Alu insertions occurred in exons of three human genes. Based on the number of insertions and the estimated time to the most recent common ancestor of HuRef and the HGP reference genome, we estimated the Alu, L1, and SVA retrotransposition rates to be one in 21 births, 212 births, and 916 births, respectively. This study presents the first comprehensive analysis of mobile element-related structural variants in the complete DNA sequence of an individual and demonstrates that mobile elements play an important role in generating inter-individual structural variation.
Subject(s)
Computational Biology/methods , DNA Transposable Elements/genetics , Genetic Variation , Genome, Human , Alu Elements , Gene Deletion , Humans , Molecular Sequence Data , Polymorphism, Single Nucleotide , Sequence Analysis, DNAABSTRACT
Owing to their phylogenetic position, cartilaginous fishes (sharks, rays, skates, and chimaeras) provide a critical reference for our understanding of vertebrate genome evolution. The relatively small genome of the elephant shark, Callorhinchus milii, a chimaera, makes it an attractive model cartilaginous fish genome for whole-genome sequencing and comparative analysis. Here, the authors describe survey sequencing (1.4x coverage) and comparative analysis of the elephant shark genome, one of the first cartilaginous fish genomes to be sequenced to this depth. Repetitive sequences, represented mainly by a novel family of short interspersed element-like and long interspersed element-like sequences, account for about 28% of the elephant shark genome. Fragments of approximately 15,000 elephant shark genes reveal specific examples of genes that have been lost differentially during the evolution of tetrapod and teleost fish lineages. Interestingly, the degree of conserved synteny and conserved sequences between the human and elephant shark genomes are higher than that between human and teleost fish genomes. Elephant shark contains putative four Hox clusters indicating that, unlike teleost fish genomes, the elephant shark genome has not experienced an additional whole-genome duplication. These findings underscore the importance of the elephant shark as a critical reference vertebrate genome for comparative analysis of the human and other vertebrate genomes. This study also demonstrates that a survey-sequencing approach can be applied productively for comparative analysis of distantly related vertebrate genomes.
Subject(s)
Genome , Sharks/genetics , Amino Acid Sequence , Animals , Base Sequence , DNA , Humans , Molecular Sequence Data , Phylogeny , Repetitive Sequences, Nucleic AcidABSTRACT
Presented here is a genome sequence of an individual human. It was produced from approximately 32 million random DNA fragments, sequenced by Sanger dideoxy technology and assembled into 4,528 scaffolds, comprising 2,810 million bases (Mb) of contiguous sequence with approximately 7.5-fold coverage for any given region. We developed a modified version of the Celera assembler to facilitate the identification and comparison of alternate alleles within this individual diploid genome. Comparison of this genome and the National Center for Biotechnology Information human reference assembly revealed more than 4.1 million DNA variants, encompassing 12.3 Mb. These variants (of which 1,288,319 were novel) included 3,213,401 single nucleotide polymorphisms (SNPs), 53,823 block substitutions (2-206 bp), 292,102 heterozygous insertion/deletion events (indels)(1-571 bp), 559,473 homozygous indels (1-82,711 bp), 90 inversions, as well as numerous segmental duplications and copy number variation regions. Non-SNP DNA variation accounts for 22% of all events identified in the donor, however they involve 74% of all variant bases. This suggests an important role for non-SNP genetic alterations in defining the diploid genome structure. Moreover, 44% of genes were heterozygous for one or more variants. Using a novel haplotype assembly strategy, we were able to span 1.5 Gb of genome sequence in segments >200 kb, providing further precision to the diploid nature of the genome. These data depict a definitive molecular portrait of a diploid human genome that provides a starting point for future genome comparisons and enables an era of individualized genomic information.
Subject(s)
Chromosome Mapping , Diploidy , Genome, Human , Sequence Analysis, DNA , Base Sequence , Chromosome Mapping/instrumentation , Chromosome Mapping/methods , Chromosomes, Human , Chromosomes, Human, Y/genetics , Gene Dosage , Genotype , Haplotypes , Human Genome Project , Humans , INDEL Mutation , In Situ Hybridization, Fluorescence , Male , Microarray Analysis , Middle Aged , Molecular Sequence Data , Pedigree , Phenotype , Polymorphism, Single Nucleotide , Reproducibility of Results , Sequence Analysis, DNA/instrumentation , Sequence Analysis, DNA/methodsABSTRACT
5-hydroxytryptamine type 3 (5-HT3) receptors are cation-selective transmitter-gated ion channels of the Cys-loop superfamily. The single-channel conductance of human recombinant 5-HT3 receptors assembled as homomers of 5-HT3A subunits, or heteromers of 5-HT3A and 5-HT3B subunits, are markedly different, being 0.4 pS (refs 6, 9) and 16 pS (ref. 7), respectively. Paradoxically, the channel-lining M2 domain of the 5-HT3A subunit would be predicted to promote cation conduction, whereas that of the 5-HT3B subunit would not. Here we describe a determinant of single-channel conductance that can explain these observations. By constructing chimaeric 5-HT3A and 5-HT3B subunits we identified a region (the 'HA-stretch') within the large cytoplasmic loop of the receptor that markedly influences channel conductance. Replacement of three arginine residues unique to the HA-stretch of the 5-HT3A subunit by their 5-HT3B subunit counterparts increased single-channel conductance 28-fold. Significantly, ultrastructural studies of the Torpedo nicotinic acetylcholine receptor indicate that the key residues might frame narrow openings that contribute to the permeation pathway. Our findings solve the conundrum of the anomalously low conductance of homomeric 5-HT3A receptors and indicate an important function for the HA-stretch in Cys-loop transmitter-gated ion channels.