Screening the human exome: a comparison of whole genome and whole transcriptome sequencing.

BACKGROUND: There is considerable interest in the development of methods to efficiently identify all coding variants present in large sample sets of humans. There are three approaches possible: whole-genome sequencing, whole-exome sequencing using exon capture methods, and RNA-Seq. While whole-genome sequencing is the most complete, it remains sufficiently expensive that cost effective alternatives are important. RESULTS: Here we provide a systematic exploration of how well RNA-Seq can identify human coding variants by comparing variants identified through high coverage whole-genome sequencing to those identified by high coverage RNA-Seq in the same individual. This comparison allowed us to directly evaluate the sensitivity and specificity of RNA-Seq in identifying coding variants, and to evaluate how key parameters such as the degree of coverage and the expression levels of genes interact to influence performance. We find that although only 40% of exonic variants identified by whole genome sequencing were captured using RNA-Seq; this number rose to 81% when concentrating on genes known to be well-expressed in the source tissue. We also find that a high false positive rate can be problematic when working with RNA-Seq data, especially at higher levels of coverage. CONCLUSIONS: We conclude that as long as a tissue relevant to the trait under study is available and suitable quality control screens are implemented, RNA-Seq is a fast and inexpensive alternative approach for finding coding variants in genes with sufficiently high expression levels.

Display of cell surface sites for fibronectin assembly is modulated by cell adherence to (1)F3 and C-terminal modules of fibronectin.

BACKGROUND: Fibronectin-null cells assemble soluble fibronectin shortly after adherence to a substrate coated with intact fibronectin but not when adherent to the cell-binding domain of fibronectin (modules (7)F3-(10)F3). Interactions of adherent cells with regions of adsorbed fibronectin other than modules (7)F3-(10)F3, therefore, are required for early display of the cell surface sites that initiate and direct fibronectin assembly. METHODOLOGY/PRINCIPAL FINDINGS: To identify these regions, coatings of proteolytically derived or recombinant pieces of fibronectin containing modules in addition to (7)F3-(10)F3 were tested for effects on fibronectin assembly by adherent fibronectin-null fibroblasts. Pieces as large as one comprising modules (2)F3-(14)F3, which include the heparin-binding and cell adhesion domains, were not effective in supporting fibronectin assembly. Addition of module (1)F3 or the C-terminal modules to modules (2)F3-(14)F3 resulted in some activity, and addition of both (1)F3 and the C-terminal modules resulted in a construct, (1)F3-C, that best mimicked the activity of a coating of intact fibronectin. Constructs (1)F3-C V0, (1)F3-C V64, and (1)F3-C Delta(V(15)F3(10)F1) were all able to support fibronectin assembly, suggesting that (1)F3 through (11)F1 and/or (12)F1 were important for activity. Coatings in which the active parts of (1)F3-C were present in different proteins were much less active than intact (1)F3-C. CONCLUSIONS: These results suggest that (1)F3 acts together with C-terminal modules to induce display of fibronectin assembly sites on adherent cells.

Transmission of single HIV-1 genomes and dynamics of early immune escape revealed by ultra-deep sequencing.

We used ultra-deep sequencing to obtain tens of thousands of HIV-1 sequences from regions targeted by CD8+ T lymphocytes from longitudinal samples from three acutely infected subjects, and modeled viral evolution during the critical first weeks of infection. Previous studies suggested that a single virus established productive infection, but these conclusions were tempered because of limited sampling; now, we have greatly increased our confidence in this observation through modeling the observed earliest sample diversity based on vastly more extensive sampling. Conventional sequencing of HIV-1 from acute/early infection has shown different patterns of escape at different epitopes; we investigated the earliest escapes in exquisite detail. Over 3-6 weeks, ultradeep sequencing revealed that the virus explored an extraordinary array of potential escape routes in the process of evading the earliest CD8 T-lymphocyte responses--using 454 sequencing, we identified over 50 variant forms of each targeted epitope during early immune escape, while only 2-7 variants were detected in the same samples via conventional sequencing. In contrast to the diversity seen within epitopes, non-epitope regions, including the Envelope V3 region, which was sequenced as a control in each subject, displayed very low levels of variation. In early infection, in the regions sequenced, the consensus forms did not have a fitness advantage large enough to trigger reversion to consensus amino acids in the absence of immune pressure. In one subject, a genetic bottleneck was observed, with extensive diversity at the second time point narrowing to two dominant escape forms by the third time point, all within two months of infection. Traces of immune escape were observed in the earliest samples, suggesting that immune pressure is present and effective earlier than previously reported; quantifying the loss rate of the founder virus suggests a direct role for CD8 T-lymphocyte responses in viral containment after peak viremia. Dramatic shifts in the frequencies of epitope variants during the first weeks of infection revealed a complex interplay between viral fitness and immune escape.

A C-terminal motif found in the beta2-adrenergic receptor, P2Y1 receptor and cystic fibrosis transmembrane conductance regulator determines binding to the Na+/H+ exchanger regulatory factor family of PDZ proteins.

The Na+/H+ exchanger regulatory factor (NHERF) binds to the tail of the beta2-adrenergic receptor and plays a role in adrenergic regulation of Na+/H+ exchange. NHERF contains two PDZ domains, the first of which is required for its interaction with the beta2 receptor. Mutagenesis studies of the beta2 receptor tail revealed that the optimal C-terminal motif for binding to the first PDZ domain of NHERF is D-S/T-x-L, a motif distinct from those recognized by other PDZ domains. The first PDZ domain of NHERF-2, a protein that is 52% identical to NHERF and also known as E3KARP, SIP-1, and TKA-1, exhibits binding preferences very similar to those of the first PDZ domain of NHERF. The delineation of the preferred binding motif for the first PDZ domain of the NHERF family of proteins allows for predictions for other proteins that may interact with NHERF or NHERF-2. For example, as would be predicted from the beta2 receptor tail mutagenesis studies, NHERF binds to the tail of the purinergic P2Y1 receptor, a seven-transmembrane receptor with an intracellular C-terminal tail ending in D-T-S-L. NHERF also binds to the tail of the cystic fibrosis transmembrane conductance regulator, which ends in D-T-R-L. Because the preferred binding motif of the first PDZ domain of the NHERF family of proteins is found at the C termini of a variety of intracellular proteins, NHERF and NHERF-2 may be multifunctional adaptor proteins involved in many previously unsuspected aspects of intracellular signaling.

Elucidation of hepatitis C virus transmission and early diversification by single genome sequencing.

A precise molecular identification of transmitted hepatitis C virus (HCV) genomes could illuminate key aspects of transmission biology, immunopathogenesis and natural history. We used single genome sequencing of 2,922 half or quarter genomes from plasma viral RNA to identify transmitted/founder (T/F) viruses in 17 subjects with acute community-acquired HCV infection. Sequences from 13 of 17 acute subjects, but none of 14 chronic controls, exhibited one or more discrete low diversity viral lineages. Sequences within each lineage generally revealed a star-like phylogeny of mutations that coalesced to unambiguous T/F viral genomes. Numbers of transmitted viruses leading to productive clinical infection were estimated to range from 1 to 37 or more (median = 4). Four acutely infected subjects showed a distinctly different pattern of virus diversity that deviated from a star-like phylogeny. In these cases, empirical analysis and mathematical modeling suggested high multiplicity virus transmission from individuals who themselves were acutely infected or had experienced a virus population bottleneck due to antiviral drug therapy. These results provide new quantitative and qualitative insights into HCV transmission, revealing for the first time virus-host interactions that successful vaccines or treatment interventions will need to overcome. Our findings further suggest a novel experimental strategy for identifying full-length T/F genomes for proteome-wide analyses of HCV biology and adaptation to antiviral drug or immune pressures.

HIV-1 envelope gp41 antibodies can originate from terminal ileum B cells that share cross-reactivity with commensal bacteria.

Monoclonal antibodies derived from blood plasma cells of acute HIV-1-infected individuals are predominantly targeted to the HIV Env gp41 and cross-reactive with commensal bacteria. To understand this phenomenon, we examined anti-HIV responses in ileum B cells using recombinant antibody technology and probed their relationship to commensal bacteria. The dominant ileum B cell response was to Env gp41. Remarkably, a majority (82%) of the ileum anti-gp41 antibodies cross-reacted with commensal bacteria, and of those, 43% showed non-HIV-1 antigen polyreactivity. Pyrosequencing revealed shared HIV-1 antibody clonal lineages between ileum and blood. Mutated immunoglobulin G antibodies cross-reactive with both Env gp41 and microbiota could also be isolated from the ileum of HIV-1 uninfected individuals. Thus, the gp41 commensal bacterial antigen cross-reactive antibodies originate in the intestine, and the gp41 Env response in HIV-1 infection can be derived from a preinfection memory B cell pool triggered by commensal bacteria that cross-react with Env.

Hybrid error correction and de novo assembly of single-molecule sequencing reads.

Single-molecule sequencing instruments can generate multikilobase sequences with the potential to greatly improve genome and transcriptome assembly. However, the error rates of single-molecule reads are high, which has limited their use thus far to resequencing bacteria. To address this limitation, we introduce a correction algorithm and assembly strategy that uses short, high-fidelity sequences to correct the error in single-molecule sequences. We demonstrate the utility of this approach on reads generated by a PacBio RS instrument from phage, prokaryotic and eukaryotic whole genomes, including the previously unsequenced genome of the parrot Melopsittacus undulatus, as well as for RNA-Seq reads of the corn (Zea mays) transcriptome. Our long-read correction achieves >99.9% base-call accuracy, leading to substantially better assemblies than current sequencing strategies: in the best example, the median contig size was quintupled relative to high-coverage, second-generation assemblies. Greater gains are predicted if read lengths continue to increase, including the prospect of single-contig bacterial chromosome assembly.

Sequencing three crocodilian genomes to illuminate the evolution of archosaurs and amniotes.

The International Crocodilian Genomes Working Group (ICGWG) will sequence and assemble the American alligator (Alligator mississippiensis), saltwater crocodile (Crocodylus porosus) and Indian gharial (Gavialis gangeticus) genomes. The status of these projects and our planned analyses are described.

High-coverage sequencing and annotated assemblies of the budgerigar genome.

BACKGROUND: Parrots belong to a group of behaviorally advanced vertebrates and have an advanced ability of vocal learning relative to other vocal-learning birds. They can imitate human speech, synchronize their body movements to a rhythmic beat, and understand complex concepts of referential meaning to sounds. However, little is known about the genetics of these traits. Elucidating the genetic bases would require whole genome sequencing and a robust assembly of a parrot genome. FINDINGS: We present a genomic resource for the budgerigar, an Australian Parakeet (Melopsittacus undulatus) -- the most widely studied parrot species in neuroscience and behavior. We present genomic sequence data that includes over 300× raw read coverage from multiple sequencing technologies and chromosome optical maps from a single male animal. The reads and optical maps were used to create three hybrid assemblies representing some of the largest genomic scaffolds to date for a bird; two of which were annotated based on similarities to reference sets of non-redundant human, zebra finch and chicken proteins, and budgerigar transcriptome sequence assemblies. The sequence reads for this project were in part generated and used for both the Assemblathon 2 competition and the first de novo assembly of a giga-scale vertebrate genome utilizing PacBio single-molecule sequencing. CONCLUSIONS: Across several quality metrics, these budgerigar assemblies are comparable to or better than the chicken and zebra finch genome assemblies built from traditional Sanger sequencing reads, and are sufficient to analyze regions that are difficult to sequence and assemble, including those not yet assembled in prior bird genomes, and promoter regions of genes differentially regulated in vocal learning brain regions. This work provides valuable data and material for genome technology development and for investigating the genomics of complex behavioral traits.

Whole genome sequencing identifies circulating Beijing-lineage Mycobacterium tuberculosis strains in Guatemala and an associated urban outbreak.

Limited data are available regarding the molecular epidemiology of Mycobacterium tuberculosis (Mtb) strains circulating in Guatemala. Beijing-lineage Mtb strains have gained prevalence worldwide and are associated with increased virulence and drug resistance, but there have been only a few cases reported in Central America. Here we report the first whole genome sequencing of Central American Beijing-lineage strains of Mtb. We find that multiple Beijing-lineage strains, derived from independent founding events, are currently circulating in Guatemala, but overall still represent a relatively small proportion of disease burden. Finally, we identify a specific Beijing-lineage outbreak centered on a poor neighborhood in Guatemala City.

Guidelines on artificial nutrition versus hydration in terminal cancer patients

Whether a terminally ill cancer patient should be actively fed or simply hydrated through subcutaneous or intravenous infusion of isotonic fluids is a matter of ongoing controversy among clinicians involved in the care of these patients. Under the auspices of the European Association for Palliative Care, a committee of experts developed guidelines to help clinicians make a reasonable decision on what type of nutritional support should be provided on a case-by-case basis. It was acknowledged that part of the controversy related to the definition of the terminal cancer patient, since this is a heterogeneous group of patients with different needs, expectations, and potential for a medical intervention. A major difficulty is the prediction of life expectancy and the patient's likely response to vigorous nutritional support. In an attempt to reach a decision on the type of treatment support (artificial nutrition vs. hydration) which would best meet the needs and expectations of the patient, we propose a three-step process: Step I: define the eight key elements necessary to reach a decision: Step II: make the decision; and Step III: reevaluate the patient and the proposed treatment at specified intervals. Step I involves assessing the patient concerning the following: 1) oncological/clinical condition; 2) symptoms; 3) expected length of survival; 4) hydration and nutritional status; 5) spontaneous or voluntary nutrient intake; 6) psychological profile; 7) gut function and potential route of administration; and 8) need for special services based on type of nutritional support prescribed. Step II involves the overall assessment of pros and cons, based on information determined in Step I, in order to reach an appropriate decision based on a well-defined end point (i.e. improvement of quality of life; maintaining patient survival; attaining rehydration). Step III involves the periodic reevaluation of the decision made in Step II based on the proposed goal and the attained result.

Solution structure of the N-terminal transactivation domain of ERM modified by SUMO-1.

ERM is a member of the PEA3 group of the Ets transcription factor family that plays important roles in development and tumorigenesis. The PEA3s share an N-terminal transactivation domain (TADn) whose activity is inhibited by small ubiquitin-like modifier (SUMO). However, the consequences of sumoylation and its underlying molecular mechanism remain unclear. The domain structure of ERM TADn alone or modified by SUMO-1 was analyzed using small-angle X-ray scattering (SAXS). Low resolution shapes determined ab initio from the scattering data indicated an elongated shape and an unstructured conformation of TADn in solution. Covalent attachment of SUMO-1 does not perturb the structure of TADn as indicated by the linear arrangement of the SUMO moiety with respect to TADn. Thus, ERM belongs to the growing family of proteins that contain intrinsically unstructured regions. The flexible nature of TADn may be instrumental for ERM recognition and binding to diverse molecular partners.