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.

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.

De novo design and molecular assembly of a transmembrane diporphyrin-binding protein complex.

The de novo design of membrane proteins remains difficult despite recent advances in understanding the factors that drive membrane protein folding and association. We have designed a membrane protein PRIME (PoRphyrins In MEmbrane) that positions two non-natural iron diphenylporphyrins (Fe(III)DPP's) sufficiently close to provide a multicentered pathway for transmembrane electron transfer. Computational methods previously used for the design of multiporphyrin water-soluble helical proteins were extended to this membrane target. Four helices were arranged in a D(2)-symmetrical bundle to bind two Fe(II/III) diphenylporphyrins in a bis-His geometry further stabilized by second-shell hydrogen bonds. UV-vis absorbance, CD spectroscopy, analytical ultracentrifugation, redox potentiometry, and EPR demonstrate that PRIME binds the cofactor with high affinity and specificity in the expected geometry.

Morphing low-affinity ligands into high-avidity nanoparticles by thermally triggered self-assembly of a genetically encoded polymer.

Multivalency is the increase in avidity resulting from the simultaneous interaction of multiple ligands with multiple receptors. This phenomenon, seen in antibody-antigen and virus-cell membrane interactions, is useful in designing bioinspired materials for targeted delivery of drugs or imaging agents. While increased avidity offered by multivalent targeting is attractive, it can also promote nonspecific receptor interaction in nontarget tissues, reducing the effectiveness of multivalent targeting. Here, we present a thermal targeting strategy--dynamic affinity modulation (DAM)--using elastin-like polypeptide diblock copolymers (ELP(BC)s) that self-assemble from a low-affinity to high-avidity state by a tunable thermal "switch", thereby restricting activity to the desired site of action. We used an in vitro cell binding assay to investigate the effect of the thermally triggered self-assembly of these ELP(BC)s on their receptor-mediated binding and cellular uptake. The data presented herein show that (1) ligand presentation does not disrupt ELP(BC) self-assembly; (2) both multivalent ligand presentation and upregulated receptor expression are needed for receptor-mediated interaction; (3) increased size of the hydrophobic segment of the block copolymer promotes multivalent interaction with membrane receptors, potentially due to changes in the nanoscale architecture of the micelle; and (4) nanoscale presentation of the ligand is important, as presentation of the ligand by micrometer-sized aggregates of an ELP showed a low level of binding/uptake by receptor-positive cells compared to its presentation on the corona of a micelle. These data validate the concept of thermally triggered DAM and provide rational design parameters for future applications of this technology for targeted drug delivery.

Telomere disruption results in non-random formation of de novo dicentric chromosomes involving acrocentric human chromosomes.

Genome rearrangement often produces chromosomes with two centromeres (dicentrics) that are inherently unstable because of bridge formation and breakage during cell division. However, mammalian dicentrics, and particularly those in humans, can be quite stable, usually because one centromere is functionally silenced. Molecular mechanisms of centromere inactivation are poorly understood since there are few systems to experimentally create dicentric human chromosomes. Here, we describe a human cell culture model that enriches for de novo dicentrics. We demonstrate that transient disruption of human telomere structure non-randomly produces dicentric fusions involving acrocentric chromosomes. The induced dicentrics vary in structure near fusion breakpoints and like naturally-occurring dicentrics, exhibit various inter-centromeric distances. Many functional dicentrics persist for months after formation. Even those with distantly spaced centromeres remain functionally dicentric for 20 cell generations. Other dicentrics within the population reflect centromere inactivation. In some cases, centromere inactivation occurs by an apparently epigenetic mechanism. In other dicentrics, the size of the alpha-satellite DNA array associated with CENP-A is reduced compared to the same array before dicentric formation. Extra-chromosomal fragments that contained CENP-A often appear in the same cells as dicentrics. Some of these fragments are derived from the same alpha-satellite DNA array as inactivated centromeres. Our results indicate that dicentric human chromosomes undergo alternative fates after formation. Many retain two active centromeres and are stable through multiple cell divisions. Others undergo centromere inactivation. This event occurs within a broad temporal window and can involve deletion of chromatin that marks the locus as a site for CENP-A maintenance/replenishment.

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.

The Effects of Geography, Environment and Phylogeny on Community Assembly and Gene Flow Dynamics

It is increasingly evident that evolutionary processes play a role in how ecological communities are assembled. However the extend to which evolution influences how plants respond to spatial and environmental gradients and interact with each other is less clear. In this dissertation I leverage evolutionary tools and thinking to understand how space and environment affect community composition and patterns of gene flow in a unique system of Atlantic rainforest and restinga (sandy coastal plains) habitats in Southeastern Brazil.

In chapter one I investigate how space and environment affect the population genetic structure and gene flow of Aechmea nudicaulis, a bromeliad species that co-occurs in forest and restinga habitats. I genotyped seven microsatellite loci and sequenced one chloroplast DNA region for individuals collected in 7 pairs of forest / restinga sites. Bayesian genetic clustering analyses show that populations of A. nudicaulis are geographically structured in northern and southern populations, a pattern consistent with broader scale phylogeographic dynamics of the Atlantic rainforest. On the other hand, explicit migration models based on the coalescent estimate that inter-habitat gene flow is less common than gene flow between populations in the same habitat type, despite their geographic discontinuity. I conclude that there is evidence for repeated colonization of the restingas from forest populations even though the steep environmental gradient between habitats is a stronger barrier to gene flow than geographic distance.

In chapter two I use data on 2800 individual plants finely mapped in a restinga plot and on first-year survival of 500 seedlings to understand the roles of phylogeny, functional traits and abiotic conditions in the spatial structuring of that community. I demonstrate that phylogeny is a poor predictor of functional traits in and that convergence in these traits is pervasive. In general, the community is not phylogenetically structured, with at best 14% of the plots deviating significantly from the null model. The functional traits SLA, leaf dry matter content (LDMC), and maximum height also showed no clear pattern of spatial structuring. On the other hand, leaf area is strongly overdispersed across all spatial scales. Although leaf area overdispersion would be generally taken as evidence of competition, I argue that interpretation is probably misleading. Finally, I show that seedling survival is dramatically increased when they grow shaded by an adult individual, suggesting that seedlings are being facilitated. Phylogenetic distance to their adult neighbor has no influence on rates of survival though. Taken together, these results indicate that phylogeny has very limited influence on the fine scale assembly of restinga communities.

Patterns of de novo allo B cells and antibody formation in chronic cardiac allograft rejection after alemtuzumab treatment.

Even though the etiology of chronic rejection (CR) is multifactorial, donor specific antibody (DSA) is considered to have a causal effect on CR development. Currently the antibody-mediated mechanisms during CR are poorly understood due to lack of proper animal models and tools. In a clinical setting, we previously demonstrated that induction therapy by lymphocyte depletion, using alemtuzumab (anti-human CD52), is associated with an increased incidence of serum alloantibody, C4d deposition and antibody-mediated rejection in human patients. In this study, the effects of T cell depletion in the development of antibody-mediated rejection were examined using human CD52 transgenic (CD52Tg) mice treated with alemtuzumab. Fully mismatched cardiac allografts were transplanted into alemtuzumab treated CD52Tg mice and showed no acute rejection while untreated recipients acutely rejected their grafts. However, approximately half of long-term recipients showed increased degree of vasculopathy, fibrosis and perivascular C3d depositions at posttransplant day 100. The development of CR correlated with DSA and C3d deposition in the graft. Using novel tracking tools to monitor donor-specific B cells, alloreactive B cells were shown to increase in accordance with DSA detection. The current animal model could provide a means of testing strategies to understand mechanisms and developing therapeutic approaches to prevent chronic rejection.

Enhanced de novo alloantibody and antibody-mediated injury in rhesus macaques.

Chronic allograft rejection is a major impediment to long-term transplant success. Humoral immune responses to alloantigens are a growing clinical problem in transplantation, with mounting evidence associating alloantibodies with the development of chronic rejection. Nearly a third of transplant recipients develop de novo antibodies, for which no established therapies are effective at preventing or eliminating, highlighting the need for a nonhuman primate model of antibody-mediated rejection. In this report, we demonstrate that depletion using anti-CD3 immunotoxin (IT) combined with maintenance immunosuppression that included tacrolimus with or without alefacept reliably prolonged renal allograft survival in rhesus monkeys. In these animals, a preferential skewing toward CD4 repopulation and proliferation was observed, particularly with the addition of alefacept. Furthermore, alefacept-treated animals demonstrated increased alloantibody production (100%) and morphologic features of antibody-mediated injury. In vitro, alefacept was found to enhance CD4 effector memory T cell proliferation. In conclusion, alefacept administration after depletion and with tacrolimus promotes a CD4+memory T cell and alloantibody response, with morphologic changes reflecting antibody-mediated allograft injury. Early and consistent de novo alloantibody production with associated histological changes makes this nonhuman primate model an attractive candidate for evaluating targeted therapeutics.