Comprehensive Characterization of Human Genome Variation by High Coverage Whole-Genome Sequencing of Forty Four Caucasians

Whole genome sequencing studies are essential to obtain a comprehensive understanding of the vast pattern of human genomic variations. Here we report the results of a high-coverage whole genome sequencing study for 44 unrelated healthy Caucasian adults, each sequenced to over 50-fold coverage (averaging 65.8×). We identified approximately 11 million single nucleotide polymorphisms (SNPs), 2.8 million short insertions and deletions, and over 500,000 block substitutions. We showed that, although previous studies, including the 1000 Genomes Project Phase 1 study, have catalogued the vast majority of common SNPs, many of the low-frequency and rare variants remain undiscovered. For instance, approximately 1.4 million SNPs and 1.3 million short indels that we found were novel to both the dbSNP and the 1000 Genomes Project Phase 1 data sets, and the majority of which (∼96%) have a minor allele frequency less than 5%. On average, each individual genome carried ∼3.3 million SNPs and ∼492,000 indels/block substitutions, including approximately 179 variants that were predicted to cause loss of function of the gene products. Moreover, each individual genome carried an average of 44 such loss-of-function variants in a homozygous state, which would completely “knock out” the corresponding genes. Across all the 44 genomes, a total of 182 genes were “knocked-out” in at least one individual genome, among which 46 genes were “knocked out” in over 30% of our samples, suggesting that a number of genes are commonly “knocked-out” in general populations. Gene ontology analysis suggested that these commonly “knocked-out” genes are enriched in biological process related to antigen processing and immune response. Our results contribute towards a comprehensive characterization of human genomic variation, especially for less-common and rare variants, and provide an invaluable resource for future genetic studies of human variation and diseases.     

Diversity of “Candidatus Liberibacter asiaticus,” Based on the omp Gene Sequence

Huanglongbing (yellow dragon disease) is a destructive disease of citrus. The etiological agent is a noncultured, phloem-restricted alpha-proteobacterium, “Candidatus Liberibacter africanus” in Africa and “Candidatus Liberibacter asiaticus” in Asia. In this study, we used an omp-based PCR-restriction fragment length polymorphism (RFLP) approach to analyze the genetic variability of “Ca. Liberibacter asiaticus” isolates. By using five different enzymes, each the 10 isolates tested could be associated with a specific combination of restriction profiles. The results indicate that the species “Ca. Liberibacter asiaticus,” even within a given region, may comprise several different variants. Thus, omp-based PCR-RFLP analysis is a simple method for detecting and differentiating “Ca. Liberibacter asiaticus” isolates.     

Bayesian Analysis of Congruence of Core Genes in Prochlorococcus and Synechococcus and Implications on Horizontal Gene Transfer

It is often suggested that horizontal gene transfer is so ubiquitous in microbes that the concept of a phylogenetic tree representing the pattern of vertical inheritance is oversimplified or even positively misleading. “Universal proteins” have been used to infer the organismal phylogeny, but have been criticized as being only the “tree of one percent.” Currently, few options exist for those wishing to rigorously assess how well a universal protein phylogeny, based on a relative handful of well-conserved genes, represents the phylogenetic histories of hundreds of genes. Here, we address this problem by proposing a visualization method and a statistical test within a Bayesian framework. We use the genomes of marine cyanobacteria, a group thought to exhibit substantial amounts of HGT, as a test case. We take 379 orthologous gene families from 28 cyanobacteria genomes and estimate the Bayesian posterior distributions of trees – a “treecloud” – for each, as well as for a concatenated dataset based on putative “universal proteins.” We then calculate the average distance between trees within and between all treeclouds on various metrics and visualize this high-dimensional space with non-metric multidimensional scaling (NMMDS). We show that the tree space is strongly clustered and that the universal protein treecloud is statistically significantly closer to the center of this tree space than any individual gene treecloud. We apply several commonly-used tests for incongruence/HGT and show that they agree HGT is rare in this dataset, but make different choices about which genes were subject to HGT. Our results show that the question of the representativeness of the “tree of one percent” is a quantitative empirical question, and that the phylogenetic central tendency is a meaningful observation even if many individual genes disagree due to the various sources of incongruence.     

Whole Exome Re-Sequencing Implicates CCDC38 and Cilia Structure and Function in Resistance to Smoking Related Airflow Obstruction

Chronic obstructive pulmonary disease (COPD) is a leading cause of global morbidity and mortality and, whilst smoking remains the single most important risk factor, COPD risk is heritable. Of 26 independent genomic regions showing association with lung function in genome-wide association studies, eleven have been reported to show association with airflow obstruction. Although the main risk factor for COPD is smoking, some individuals are observed to have a high forced expired volume in 1 second (FEV₁) despite many years of heavy smoking. We hypothesised that these “resistant smokers” may harbour variants which protect against lung function decline caused by smoking and provide insight into the genetic determinants of lung health. We undertook whole exome re-sequencing of 100 heavy smokers who had healthy lung function given their age, sex, height and smoking history and applied three complementary approaches to explore the genetic architecture of smoking resistance. Firstly, we identified novel functional variants in the “resistant smokers” and looked for enrichment of these novel variants within biological pathways. Secondly, we undertook association testing of all exonic variants individually with two independent control sets. Thirdly, we undertook gene-based association testing of all exonic variants. Our strongest signal of association with smoking resistance for a non-synonymous SNP was for rs10859974 (P = 2.34×10⁻⁴) in CCDC38, a gene which has previously been reported to show association with FEV₁/FVC, and we demonstrate moderate expression of CCDC38 in bronchial epithelial cells. We identified an enrichment of novel putatively functional variants in genes related to cilia structure and function in resistant smokers. Ciliary function abnormalities are known to be associated with both smoking and reduced mucociliary clearance in patients with COPD. We suggest that genetic influences on the development or function of cilia in the bronchial epithelium may affect growth of cilia or the extent of damage caused by tobacco smoke.     

Modeling the Evolutionary Architectures of Transcribed Human Enhancer Sequences Reveals Distinct Origins,Functions, and Associations with Human Trait Variation

Despite the importance of gene regulatory enhancers in human biology and evolution, we lack a comprehensive model of enhancer evolution and function. This substantially limits our understanding of the genetic basis of species divergence and our ability to interpret the effects of noncoding variants on human traits.To explore enhancer sequence evolution and its relationship to regulatory function, we traced the evolutionary origins of transcribed human enhancer sequences with activity across diverse tissues and cellular contexts from the FANTOM5 consortium. The transcribed enhancers are enriched for sequences of a single evolutionary age (“simple” evolutionary architectures) compared with enhancers that are composites of sequences of multiple evolutionary ages (“complex” evolutionary architectures), likely indicating constraint against genomic rearrangements. Complex enhancers are older, more pleiotropic, and more active across species than simple enhancers. Genetic variants within complex enhancers are also less likely to associate with human traits and biochemical activity. Transposable-element-derived sequences (TEDS) have made diverse contributions to enhancers of both architectures; the majority of TEDS are found in enhancers with simple architectures, while a minority have remodeled older sequences to create complex architectures. Finally, we compare the evolutionary architectures of transcribed enhancers with histone-mark-defined enhancers.Our results reveal that most human transcribed enhancers are ancient sequences of a single age, and thus the evolution of most human enhancers was not driven by increases in evolutionary complexity over time. Our analyses further suggest that considering enhancer evolutionary histories provides context that can aid interpretation of the effects of variants on enhancer function. Based on these results, we propose a framework for analyzing enhancer evolutionary architecture.     

Selecting Essential Information for Biosurveillance—A Multi-Criteria Decision Analysis

The National Strategy for Biosurveillancedefines biosurveillance as “the process of gathering, integrating, interpreting, and communicating essential information related to all-hazards threats or disease activity affecting human, animal, or plant health to achieve early detection and warning, contribute to overall situational awareness of the health aspects of an incident, and to enable better decision-making at all levels.” However, the strategy does not specify how “essential information” is to be identified and integrated into the current biosurveillance enterprise, or what the metrics qualify information as being “essential”. Thequestion of data stream identification and selection requires a structured methodology that can systematically evaluate the tradeoffs between the many criteria that need to be taken in account. Multi-Attribute Utility Theory, a type of multi-criteria decision analysis, can provide a well-defined, structured approach that can offer solutions to this problem. While the use of Multi-Attribute Utility Theoryas a practical method to apply formal scientific decision theoretical approaches to complex, multi-criteria problems has been demonstrated in a variety of fields, this method has never been applied to decision support in biosurveillance.We have developed a formalized decision support analytic framework that can facilitate identification of “essential information” for use in biosurveillance systems or processes and we offer this framework to the global BSV community as a tool for optimizing the BSV enterprise. To demonstrate utility, we applied the framework to the problem of evaluating data streams for use in an integrated global infectious disease surveillance system.     

Evaluating the Effect of 3′-UTR Variants in DICER1 and DROSHA on Their Tissue-Specific Expression by miRNA Target Prediction

Untranslated gene regions (UTRs) play an important role in controlling gene expression. 3′-UTRs are primarily targeted by microRNA (miRNA) molecules that form complex gene regulatory networks. Cancer genomes are replete with non-coding mutations, many of which are connected to changes in tumor gene expression that accompany the development of cancer and are associated with resistance to therapy. Therefore, variants that occurred in 3′-UTR under cancer progression should be analysed to predict their phenotypic effect on gene expression, e.g., by evaluating their impact on miRNA target sites. Here, we analyze 3′-UTR variants in DICER1 and DROSHA genes in the context of myelodysplastic syndrome (MDS) development. The key features of this analysis include an assessment of both “canonical” and “non-canonical” types of mRNA-miRNA binding and tissue-specific profiling of miRNA interactions with wild-type and mutated genes. As a result, we obtained a list of DICER1 and DROSHA variants likely altering the miRNA sites and, therefore, potentially leading to the observed tissue-specific gene downregulation. All identified variants have low population frequency consistent with their potential association with pathology progression.     

The Characterization of Twenty Sequenced Human Genomes

We present the analysis of twenty human genomes to evaluate the prospects for identifying rare functional variants that contribute to a phenotype of interest. We sequenced at high coverage ten “case” genomes from individuals with severe hemophilia A and ten “control” genomes. We summarize the number of genetic variants emerging from a study of this magnitude, and provide a proof of concept for the identification of rare and highly-penetrant functional variants by confirming that the cause of hemophilia A is easily recognizable in this data set. We also show that the number of novel single nucleotide variants (SNVs) discovered per genome seems to stabilize at about 144,000 new variants per genome, after the first 15 individuals have been sequenced. Finally, we find that, on average, each genome carries 165 homozygous protein-truncating or stop loss variants in genes representing a diverse set of pathways.     

Mitochondrial DNA Variant Discovery and Evaluation in Human Cardiomyopathies through Next-Generation Sequencing

Mutations in mitochondrial DNA (mtDNA) may cause maternally-inherited cardiomyopathy and heart failure. In homoplasmy all mtDNA copies contain the mutation. In heteroplasmy there is a mixture of normal and mutant copies of mtDNA. The clinical phenotype of an affected individual depends on the type of genetic defect and the ratios of mutant and normal mtDNA in affected tissues. We aimed at determining the sensitivity of next-generation sequencing compared to Sanger sequencing for mutation detection in patients with mitochondrial cardiomyopathy. We studied 18 patients with mitochondrial cardiomyopathy and two with suspected mitochondrial disease. We “shotgun” sequenced PCR-amplified mtDNA and multiplexed using a single run on Roche''s 454 Genome Sequencer. By mapping to the reference sequence, we obtained 1,300× average coverage per case and identified high-confidence variants. By comparing these to >400 mtDNA substitution variants detected by Sanger, we found 98% concordance in variant detection. Simulation studies showed that >95% of the homoplasmic variants were detected at a minimum sequence coverage of 20× while heteroplasmic variants required >200× coverage. Several Sanger “misses” were detected by 454 sequencing. These included the novel heteroplasmic 7501T>C in tRNA serine 1 in a patient with sudden cardiac death. These results support a potential role of next-generation sequencing in the discovery of novel mtDNA variants with heteroplasmy below the level reliably detected with Sanger sequencing. We hope that this will assist in the identification of mtDNA mutations and key genetic determinants for cardiomyopathy and mitochondrial disease.     

The Epidemiologist's Dilemma: A Position Paper on the Role of Epidemiology in Canadian Research

Canadian epidemiology is currently in transition from being primarily interested in infectious disease to becoming active in investigating causes of non-infectious disease. Generally, epidemiologists limit their work to “field”, “basic”, or “theoretical” epidemiology. In all three fields there appears to be a shortage of qualified personnel which is likely to become even more acute because of new roles that epidemiologists will probably play in the future. A minimum of 40 full-time epidemiologists is currently required in Canada.Departments other than departments of preventive medicine at medical schools and teaching hospitals are currently spending 79.5% of all funds allocated for non-microbiological epidemiological research in Canada. Since epidemiology is by its very nature population orientated, rather than clinically orientated, clinicians require consultative advice from epidemiologists at many stages of their research. Epidemiological facilities in departments of preventive medicine should therefore be strengthened in order to provide research training and didactic courses needed as the numbers of physicians, nurses and paramedical persons are increased.     

Evolutionary Relationships of “Candidatus Riesia spp.,” Endosymbiotic Enterobacteriaceae Living within Hematophagous Primate Lice

The primary endosymbiotic bacteria from three species of parasitic primate lice were characterized molecularly. We have confirmed the characterization of the primary endosymbiont (P-endosymbiont) of the human head/body louse Pediculus humanus and provide new characterizations of the P-endosymbionts from Pediculus schaeffi from chimpanzees and Pthirus pubis, the pubic louse of humans. The endosymbionts show an average percent sequence divergence of 11 to 15% from the most closely related known bacterium “Candidatus Arsenophonus insecticola.” We propose that two additional species be added to the genus “Candidatus Riesia.” The new species proposed within “Candidatus Riesia” have sequence divergences of 3.4% and 10 to 12% based on uncorrected pairwise differences. Our Bayesian analysis shows that the branching pattern for the primary endosymbionts was the same as that for their louse hosts, suggesting a long coevolutionary history between primate lice and their primary endosymbionts. We used a calibration of 5.6 million years to date the divergence between endosymbionts from human and chimpanzee lice and estimated an evolutionary rate of nucleotide substitution of 0.67% per million years, which is 15 to 30 times faster than previous estimates calculated for Buchnera, the primary endosymbiont in aphids. Given the evidence for cospeciation with primate lice and the evidence for fast evolutionary rates, this lineage of endosymbiotic bacteria can be evaluated as a fast-evolving marker of both louse and primate evolutionary histories.     

Monitoring DNA Contamination in Handled vs. Directly Excavated Ancient Human Skeletal Remains

Bones, teeth and hair are often the only physical evidence of human or animal presence at an archaeological site; they are also the most widely used sources of samples for ancient DNA (aDNA) analysis. Unfortunately, the DNA extracted from ancient samples, already scarce and highly degraded, is widely susceptible to exogenous contaminations that can affect the reliability of aDNA studies. We evaluated the molecular effects of sample handling on five human skeletons freshly excavated from a cemetery dated between the 11 to the 14^th century. We collected specimens from several skeletal areas (teeth, ribs, femurs and ulnas) from each individual burial. We then divided the samples into two different sets: one labeled as “virgin samples” (i.e. samples that were taken by archaeologists under contamination-controlled conditions and then immediately sent to the laboratory for genetic analyses), and the second called “lab samples”(i.e. samples that were handled without any particular precautions and subject to normal washing, handling and measuring procedures in the osteological lab). Our results show that genetic profiles from “lab samples” are incomplete or ambiguous in the different skeletal areas while a different outcome is observed in the “virgin samples” set. Generally, all specimens from different skeletal areas in the exception of teeth present incongruent results between “lab” and “virgin” samples. Therefore teeth are less prone to contamination than the other skeletal areas we analyzed and may be considered a material of choice for classical aDNA studies. In addition, we showed that bones can also be a good candidate for human aDNA analysis if they come directly from the excavation site and are accompanied by a clear taphonomic history.     

Categorizing Ideas about Trees: A Tree of Trees

The aim of this study is to explore whether matrices and MP trees used to produce systematic categories of organisms could be useful to produce categories of ideas in history of science. We study the history of the use of trees in systematics to represent the diversity of life from 1766 to 1991. We apply to those ideas a method inspired from coding homologous parts of organisms. We discretize conceptual parts of ideas, writings and drawings about trees contained in 41 main writings; we detect shared parts among authors and code them into a 91-characters matrix and use a tree representation to show who shares what with whom. In other words, we propose a hierarchical representation of the shared ideas about trees among authors: this produces a “tree of trees.” Then, we categorize schools of tree-representations. Classical schools like “cladists” and “pheneticists” are recovered but others are not: “gradists” are separated into two blocks, one of them being called here “grade theoreticians.” We propose new interesting categories like the “buffonian school,” the “metaphoricians,” and those using “strictly genealogical classifications.” We consider that networks are not useful to represent shared ideas at the present step of the study. A cladogram is made for showing who is sharing what with whom, but also heterobathmy and homoplasy of characters. The present cladogram is not modelling processes of transmission of ideas about trees, and here it is mostly used to test for proximity of ideas of the same age and for categorization.     

Dissecting the Dynamics of HIV-1 Protein Sequence Diversity

The rapid mutation of human immunodeficiency virus-type 1 (HIV-1) and the limited characterization of the composition and incidence of the variant population are major obstacles to the development of an effective HIV-1 vaccine. This issue was addressed by a comprehensive analysis of over 58,000 clade B HIV-1 protein sequences reported over at least 26 years. The sequences were aligned and the 2,874 overlapping nonamer amino acid positions of the viral proteome, each a possible core binding domain for human leukocyte antigen molecules and T-cell receptors, were quantitatively analyzed for four patterns of sequence motifs: (1) “index”, the most prevalent sequence; (2) “major” variant, the most common variant sequence; (3) “minor” variants, multiple different sequences, each with an incidence less than that of the major variant; and (4) “unique” variants, each observed only once in the alignment. The collective incidence of the major, minor, and unique variants at each nonamer position represented the total variant population for the position. Positions with more than 50% total variants contained correspondingly reduced incidences of index and major variant sequences and increased minor and unique variants. Highly diverse positions, with 80 to 98% variant nonamer sequences, were present in each protein, including 5% of Gag, and 27% of Env and Nef, each. The multitude of different variant nonamer sequences (i.e. nonatypes; up to 68%) at the highly diverse positions, represented by the major, multiple minor, and multiple unique variants likely supported variants function both in immune escape and as altered peptide ligands with deleterious T-cell responses. The patterns of mutational change were consistent with the sequences of individual HXB2 and C1P viruses and can be considered applicable to all HIV-1 viruses. This characterization of HIV-1 protein mutation provides a foundation for the design of peptide-based vaccines and therapeutics.     

What can experience add to early medical education? Consensus survey

Objective To provide a rationale for integrating experience into early medical education (“early experience”).Design Small group discussions to obtain stakeholders'' views. Grounded theory analysis with respondent, internal, and external validation.Setting Problem based, undergraduate medical curriculum that is not vertically integrated.Participants A purposive sample of 64 students, staff, and curriculum leaders from three university medical schools in the United Kingdom.Results Without early experience, the curriculum was socially isolating and divorced from clinical practice. The abruptness of students'' transition to the clinical environment in year 3 generated positive and negative emotions. The rationale for early experience would be to ease the transition; orientate the curriculum towards the social context of practice; make students more confident to approach patients; motivate them; increase their awareness of themselves and others; strengthen, deepen, and contextualise their theoretical knowledge; teach intellectual skills; strengthen learning of behavioural and social sciences; and teach them about the role of health professionals.Conclusion A rationale for early experience would be to strengthen and deepen cognitively, broaden affectively, contextualise, and integrate medical education. This is partly a process of professional socialisation that should start earlier to avoid an abrupt transition. “Experience” can be defined as “authentic human contact in a social or clinical context that enhances learning of health, illness or disease, and the role of the health professional.”     

Wet-dry-wet drug screen leads to the synthesis of TS1, a novel compound reversing lung fibrosis through inhibition of myofibroblast differentiation

Therapies halting the progression of fibrosis are ineffective and limited. Activated myofibroblasts are emerging as important targets in the progression of fibrotic diseases. Previously, we performed a high-throughput screen on lung fibroblasts and subsequently demonstrated that the inhibition of myofibroblast activation is able to prevent lung fibrosis in bleomycin-treated mice. High-throughput screens are an ideal method of repurposing drugs, yet they contain an intrinsic limitation, which is the size of the library itself. Here, we exploited the data from our “wet” screen and used “dry” machine learning analysis to virtually screen millions of compounds, identifying novel anti-fibrotic hits which target myofibroblast differentiation, many of which were structurally related to dopamine. We synthesized and validated several compounds ex vivo (“wet”) and confirmed that both dopamine and its derivative TS1 are powerful inhibitors of myofibroblast activation. We further used RNAi-mediated knock-down and demonstrated that both molecules act through the dopamine receptor 3 and exert their anti-fibrotic effect by inhibiting the canonical transforming growth factor β pathway. Furthermore, molecular modelling confirmed the capability of TS1 to bind both human and mouse dopamine receptor 3. The anti-fibrotic effect on human cells was confirmed using primary fibroblasts from idiopathic pulmonary fibrosis patients. Finally, TS1 prevented and reversed disease progression in a murine model of lung fibrosis. Both our interdisciplinary approach and our novel compound TS1 are promising tools for understanding and combating lung fibrosis.Subject terms: Drug discovery, Respiratory tract diseases     

Dogs (Canis familiaris) Evaluate Humans on the Basis of Direct Experiences Only

Reputation formation is a key component in the social interactions of many animal species. An evaluation of reputation is drawn from two principal sources: direct experience of an individual and indirect experience from observing that individual interacting with a third party. In the current study we investigated whether dogs use direct and/or indirect experience to choose between two human interactants. In the first experiment, subjects had direct interaction either with a “nice” human (who played with, talked to and stroked the dog) or with an “ignoring” experimenter who ignored the dog completely. Results showed that the dogs stayed longer close to the “nice” human. In a second experiment the dogs observed a “nice” or “ignoring” human interacting with another dog. This indirect experience, however, did not lead to a preference between the two humans. These results suggest that the dogs in our study evaluated humans solely on the basis of direct experience.