The DDB1–DCAF1–Vpr–UNG2 crystal structure reveals how HIV-1 Vpr steers human UNG2 toward destruction
Explore the source record for details and available documents.
Wu, Ying↗
SEARCH · Engineering Papers
Results for “DNA REPLICATION”
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Structural basis of G-quadruplex unfolding by the DEAH/RHA helicase DHX36
Guanine-rich nucleic acid sequences challenge the replication, transcription, and translation machinery by spontaneously folding into G-quadruplexes, the unfolding of which requires forces greater than most polymerases can exert. Eukaryotic cells contain numerous helicases that can unfold G-quadruplexes. The molecular basis of the recognition and unfolding of G-quadruplexes by helicases remains poorly understood. DHX36 (also known as RHAU and G4R1), a member of the DEAH/RHA family of helicases, binds both DNA and RNA G-quadruplexes with extremely high affinity, is consistently found bound to G-quadruplexes in cells, and is a major source of G-quadruplex unfolding activity in HeLa cell lysates. DHX36 is a multi-functional helicase that has been implicated in G-quadruplex-mediated transcriptional and post-transcriptional regulation, and is essential for heart development, haematopoiesis, and embryogenesis in mice. Here we report the co-crystal structure of bovine DHX36 bound to a DNA with a G-quadruplex and a 3' single-stranded DNA segment. Furthermore, we show that the N-terminal DHX36-specific motif folds into a DNA-binding-induced α-helix that, together with the OB-fold-like subdomain, selectively binds parallel G-quadruplexes. Comparison with unliganded and ATP-analogue-bound DHX36 structures, together with single-molecule fluorescence resonance energy transfer (FRET) analysis, suggests that G-quadruplex binding alone induces rearrangements of the helicase core; by pulling on the single-stranded DNA tail, these rearrangements drive G-quadruplex unfolding one residue at a time.
59 BASIC BIOLOGICAL SCIENCES↗
Toward three-dimensional DNA industrial nanorobots
Nanoscale industrial robots have potential as manufacturing platforms and are capable of automatically performing repetitive tasks to handle and produce nanomaterials with consistent precision and accuracy. We demonstrate a DNA industrial nanorobot that fabricates a three-dimensional (3D), optically active chiral structure from optically inactive parts. By making use of externally controlled temperature and ultraviolet (UV) light, our programmable robot, ~100 nanometers in size, grabs different parts, positions and aligns them so that they can be welded, releases the construct, and returns to its original configuration ready for its next operation. Our robot can also self-replicate its 3D structure and functions, surpassing single-step templating (restricted to two dimensions) by using folding to access the third dimension and more degrees of freedom. Our introduction of multiple-axis precise folding and positioning as a tool/technology for nanomanufacturing will open the door to more complex and useful nano- and microdevices.
Robotics↗
Identification of the cAMP response element that controls transcriptional activation of the insulin-like growth factor-I gene by prostaglandin E2 in osteoblasts
Insulin-like growth factor-I (IGF-I), a multifunctional growth factor, plays a key role in skeletal growth and can enhance bone cell replication and differentiation. We previously showed that prostaglandin E2 (PGE2) and other agents that increase cAMP activated IGF-I gene transcription in primary rat osteoblast cultures through promoter 1 (P1), the major IGF-I promoter, and found that transcriptional induction was mediated by protein kinase A. We now have identified a short segment of P1 that is essential for full hormonal regulation and have characterized inducible DNA-protein interactions involving this site. Transient transfections of IGF-I P1 reporter genes into primary rat osteoblasts showed that the 328-base pair untranslated region of exon 1 was required for a full 5.3-fold response to PGE2; mutation in a previously footprinted site, HS3D (base pairs +193 to +215), reduced induction by 65%. PGE2 stimulated nuclear protein binding to HS3D. Binding, as determined by gel mobility shift assay, was not seen in nuclear extracts from untreated osteoblast cultures, was detected within 2 h of PGE2 treatment, and was maximal by 4 h. This DNA-protein interaction was not observed in cytoplasmic extracts from PGE2-treated cultures, indicating nuclear localization of the protein kinase A-activated factor(s). Activation of this factor was not blocked by cycloheximide (Chx), and Chx did not impair stimulation of IGF-I gene expression by PGE2. In contrast, binding to a consensus cAMP response element (CRE; 5'-TGACGTCA-3') from the rat somatostatin gene was not modulated by PGE2 or Chx. Competition gel mobility shift analysis using mutated DNA probes identified 5'-CGCAATCG-3' as the minimal sequence needed for inducible binding. All modified IGF-I P1 promoterreporter genes with mutations within this CRE sequence also showed a diminished functional response to PGE2. These results identify the CRE within the 5'-untranslated region of IGF-I exon 1 that is required for hormonal activation of IGF-I gene transcription by cAMP in osteoblasts.
NASA Discipline Musculoskeletal↗
Pooled PPIseq: Screening the SARS-CoV-2 and human interface with a scalable multiplexed protein-protein interaction assay platform
Protein-Protein Interactions (PPIs) are a key interface between virus and host, and these interactions are important to both viral reprogramming of the host and to host restriction of viral infection. In particular, viral-host PPI networks can be used to further our understanding of the molecular mechanisms of tissue specificity, host range, and virulence. At higher scales, viral-host PPI screening could also be used to screen for small-molecule antivirals that interfere with essential viral-host interactions, or to explore how the PPI networks between interacting viral and host genomes co-evolve. Current high-throughput PPI assays have screened entire viral-host PPI networks. However, these studies are time consuming, often require specialized equipment, and are difficult to further scale. Here, we develop methods that make larger-scale viral-host PPI screening more accessible. This approach combines the mDHFR split-tag reporter with the iSeq2 interaction-barcoding system to permit massively-multiplexed PPI quantification by simple pooled engineering of barcoded constructs, integration of these constructs into budding yeast, and fitness measurements by pooled cell competitions and barcode-sequencing. We applied this method to screen for PPIs between SARS-CoV-2 proteins and human proteins, screening in triplicate >180,000 ORF-ORF combinations represented by >1,000,000 barcoded lineages. Our results complement previous screens by identifying 74 putative PPIs, including interactions between ORF7A with the taste receptors TAS2R41 and TAS2R7, and between NSP4 with the transmembrane KDELR2 and KDELR3. We show that this PPI screening method is highly scalable, enabling larger studies aimed at generating a broad understanding of how viral effector proteins converge on cellular targets to effect replication.
60 APPLIED LIFE SCIENCES↗
Defensive symbiosis against giant viruses in amoebae
Protists are important regulators of microbial communities and key components in food webs with impact on nutrient cycling and ecosystem functioning. In turn, their activity is shaped by diverse intracellular parasites, including bacterial symbionts and viruses. Yet, bacteria–virus interactions within protists are poorly understood. Here, we studied the role of bacterial symbionts of free-living amoebae in the establishment of infections with nucleocytoplasmic large DNA viruses (Nucleocytoviricota). To investigate these interactions in a system that would also be relevant in nature, we first isolated and characterized a giant virus (Viennavirus, family Marseilleviridae) and a sympatric potential Acanthamoeba host infected with bacterial symbionts. Subsequently, coinfection experiments were carried out, using the fresh environmental isolates as well as additional amoeba laboratory strains. Employing fluorescence in situ hybridization and qPCR, we show that the bacterial symbiont, identified as Parachlamydia acanthamoebae, represses the replication of the sympatric Viennavirus in both recent environmental isolates as well as Acanthamoeba laboratory strains. In the presence of the symbiont, virions are still taken up, but viral factory maturation is inhibited, leading to survival of the amoeba host. The symbiont also suppressed the replication of the more complex Acanthamoeba polyphaga mimivirus and Tupanvirus deep ocean (Mimiviridae). Our work provides an example of an intracellular bacterial symbiont protecting a protist host against virus infections. The impact of virus–symbiont interactions on microbial population dynamics and eventually ecosystem processes requires further attention.
59 BASIC BIOLOGICAL SCIENCES↗
Measurement Error and Resolution in Quantitative Stable Isotope Probing: Implications for Experimental Design
Quantitative stable isotope probing (qSIP) estimates isotope tracer incorporation into DNA of individual microbes and can link microbial biodiversity and biogeochemistry in complex communities. As with any quantitative estimation technique, qSIP involves measurement error, and a fuller understanding of error, precision, and statistical power benefits qSIP experimental design and data interpretation. We used several qSIP data sets—from soil and seawater microbiomes—to evaluate how variance in isotope incorporation estimates depends on organism abundance and resolution of the density fractionation scheme. We assessed statistical power for replicated qSIP studies, plus sensitivity and specificity for unreplicated designs. As a taxon’s abundance increases, the variance of its weighted mean density declines. Nine fractions appear to be a reasonable trade-off between cost and precision for most qSIP applications. Increasing the number of density fractions beyond that reduces variance, although the magnitude of this benefit declines with additional fractions. Our analysis suggests that, if a taxon has an isotope enrichment of 10 atom% excess, there is a 60% chance that this will be detected as significantly different from zero (with alpha 0.1). With five replicates, isotope enrichment of 5 atom% could be detected with power (0.6) and alpha (0.1). Finally, we illustrate the importance of internal standards, which can help to calibrate per sample conversions of %GC to mean weighted density. These results should benefit researchers designing future SIP experiments and provide a useful reference for metagenomic SIP applications where both financial and computational limitations constrain experimental scope.
59 BASIC BIOLOGICAL SCIENCES↗
CCAAT/enhancer-binding protein delta activates insulin-like growth factor-I gene transcription in osteoblasts. Identification of a novel cyclic AMP signaling pathway in bone
Insulin-like growth factor-I (IGF-I) plays a key role in skeletal growth by stimulating bone cell replication and differentiation. We previously showed that prostaglandin E2 (PGE2) and other cAMP-activating agents enhanced IGF-I gene transcription in cultured primary rat osteoblasts through promoter 1, the major IGF-I promoter, and identified a short segment of the promoter, termed HS3D, that was essential for hormonal regulation of IGF-I gene expression. We now demonstrate that CCAAT/enhancer-binding protein (C/EBP) delta is a major component of a PGE2-stimulated DNA-protein complex involving HS3D and find that C/EBPdelta transactivates IGF-I promoter 1 through this site. Competition gel shift studies first indicated that a core C/EBP half-site (GCAAT) was required for binding of a labeled HS3D oligomer to osteoblast nuclear proteins. Southwestern blotting and UV-cross-linking studies showed that the HS3D probe recognized a approximately 35-kDa nuclear protein, and antibody supershift assays indicated that C/EBPdelta comprised most of the PGE2-activated gel-shifted complex. C/EBPdelta was detected by Western immunoblotting in osteoblast nuclear extracts after treatment of cells with PGE2. An HS3D oligonucleotide competed effectively with a high affinity C/EBP site from the rat albumin gene for binding to osteoblast nuclear proteins. Co-transfection of osteoblast cell cultures with a C/EBPdelta expression plasmid enhanced basal and PGE2-activated IGF-I promoter 1-luciferase activity but did not stimulate a reporter gene lacking an HS3D site. By contrast, an expression plasmid for the related protein, C/EBPbeta, did not alter basal IGF-I gene activity but did increase the response to PGE2. In osteoblasts and in COS-7 cells, C/EBPdelta, but not C/EBPbeta, transactivated a reporter gene containing four tandem copies of HS3D fused to a minimal promoter; neither transcription factor stimulated a gene with four copies of an HS3D mutant that was unable to bind osteoblast nuclear proteins. These results identify C/EBPdelta as a hormonally activated inducer of IGF-I gene transcription in osteoblasts and show that the HS3D element within IGF-I promoter 1 is a high affinity binding site for this protein.
NASA Discipline Musculoskeletal↗
Expression and purification of affinity-tagged variants of the photochemical reaction center from Heliobacterium modesticaldum
The heliobacterial photochemical reaction center (HbRC) from the chlorophototrophic Firmicutes bacterium Heliobacterium modesticaldum is the only homodimeric type I RC whose structure is known. Using genetic techniques recently established in our lab, we have developed a rapid heterologous expression system for the HbRC core polypeptide PshA. Our system relies on rescue of the non-chlorophototrophic ΔpshA::cbp2p - aph3 strain of Hbt. modesticaldum by expression of a heterologous pshA gene from a replicating shuttle vector. Furthermore, we constructed two tagged variants of PshA, one with an N-terminal octahistidine tag and one with an internal hexahistidine tag, which facilitate rapid purification of pure, active HbRC cores in milligram quantities. We constructed a suite of shuttle vectors bearing untagged or tagged versions of pshA driven by various promoters. Surprisingly, we found that the eno and gapDH_2 promoters from Clostridium thermocellum drive better expression of pshA than fragments of DNA derived from the region upstream of the pshA locus on the Hbt. modesticaldum genome. This “ pshA rescue” strategy additionally provided a useful window into how Hbt. modesticaldum regulates pigment synthesis and growth rate when chlorophototrophic output decreases.
59 BASIC BIOLOGICAL SCIENCES↗
Mechanistic insights into the role of Hop2-Mnd1 in meiotic homologous DNA pairing
The Hop2–Mnd1 complex functions with the DMC1 recombinase in meiotic recombination. Hop2–Mnd1 stabilizes the DMC1-single-stranded DNA (ssDNA) filament and promotes the capture of the double stranded DNA partner by the recombinase filament to assemble the synaptic complex. Herein, we define the action mechanism of Hop2–Mnd1 in DMC1-mediated recombination. Small angle X-ray scattering analysis and electron microscopy reveal that the heterodimeric Hop2–Mnd1 is a V-shaped molecule. We show that the protein complex harbors three distinct DNA binding sites, and determine their functional relevance. Specifically, the N-terminal double-stranded DNA binding functions of Hop2 and Mnd1 co-operate to mediate synaptic complex assembly, whereas ssDNA binding by the Hop2 C-terminus helps stabilize the DMC1-ssDNA filament. A model of the Hop2-Mnd1-DMC1-ssDNA ensemble is proposed to explain how it mediates homologous DNA pairing in meiotic recombination.
59 BASIC BIOLOGICAL SCIENCES↗
Subclinical Shed of Infectious Varicella zoster Virus in Astronauts
Aerosol borne varicella zoster virus (VZV) enters the nasopharynx and replicates in tonsillar T-cells, resulting in viremia and varicella (chickenpox). Virus then becomes latent in cranial nerve, dorsal root and autonomic nervous system ganglia along the entire neuraxis (1). Decades later, as cell-mediated immunity to VZV declines (4), latent VZV can reactivate to produce zoster (shingles). Infectious VZV is present in patients with varicella or zoster, but shed of infectious virus in the absence of disease has not been shown. We previously detected VZV DNA in saliva of astronauts during and shortly after spaceflight, suggesting stress induced subclinical virus reactivation (3). We show here that VZV DNA as well as infectious virus in present in astronaut saliva. VZV DNA was detected in saliva during and after a 13-day spaceflight in 2 of 3 astronauts (Fig. panel A). Ten days before liftoff, there was a rise in serum anti-VZV antibody in subjects 1 and 2, consistent with virus reactivation. In subject 3, VZV DNA was not detected in saliva, and there was no rise in anti-VZV antibody titer. Subject 3 may have been protected from virus reactivation by having zoster <10 years ago, which provides a boost in cell-medicated immunity to VZV (2). No VZV DNA was detected in astronaut saliva months before spaceflight, or in saliva of 10 age/sex-matched healthy control subjects sampled on alternate days for 3 weeks (88 saliva samples). Saliva taken 2-6 days after landing from all 3 subjects was cultured on human fetal lung cells (Fig. panel B). Infectious VZV was recovered from saliva of subjects 1 and 2 on the second day after landing. Virus specificity was confirmed by antibody staining and DNA analysis which showed it to be VZV of European descent, common in the US (5). Further, both antibody staining and DNA PCR demonstrated that no HSV-1 was detected in any infected culture. This is the first report of infectious VZV shedding in the absence of clinical disease. Spaceflight presents a uniquely stressful environment which includes physical isolation and confinement, anxiety, sleep deprivation, as well as exposure to increased radiation and microgravity. It is interesting that in our study, VZV and not HSV-1 reactivation was detected, since stress-induced HSV-1 reactivation has been reported (6). Future studies are needed to determine the specific inducer of VZV reactivation.
Cohrs, Randall J.↗
Defining the genetic components of callus formation: A GWAS approach
A characteristic feature of plant cells is the ability to form callus from parenchyma cells in response to biotic and abiotic stimuli. Tissue culture propagation of recalcitrant plant species and genetic engineering for desired phenotypes typically depends on efficient in vitro callus generation. Callus formation is under genetic regulation, and consequently, a molecular understanding of this process underlies successful generation for propagation materials and/or introduction of genetic elements in experimental or industrial applications. Herein, we identified 11 genetic loci significantly associated with callus formation in Populus trichocarpa using a genome-wide association study (GWAS) approach. Eight of the 11 significant gene associations were consistent across biological replications, exceeding a chromosome-wide–log10 (p) = 4.46 [p = 3.47E-05] Bonferroni-adjusted significance threshold. These eight genes were used as hub genes in a high-resolution co-expression network analysis to gain insight into the genome-wide basis of callus formation. A network of positively and negatively co-expressed genes, including several transcription factors, was identified. As proof-of-principle, a transient protoplast assay confirmed the negative regulation of a Chloroplast Nucleoid DNA-binding-related gene (Potri.018G014800) by the LEC2 transcription factor. Many of the candidate genes and co-expressed genes were 1) linked to cell division and cell cycling in plants and 2) showed homology to tumor and cancer-related genes in humans. The GWAS approach based on a high-resolution marker set, and the ability to manipulate targets genes in vitro, provided a catalog of high-confidence genes linked to callus formation that can serve as an important resource for successful manipulation of model and non-model plant species, and likewise, suggests a robust method of discovering common homologous functions across organisms.
59 BASIC BIOLOGICAL SCIENCES↗
Successful Validation of RNA Purification and Quantitative Real-Time PCR Analysis of Gene Expression on the International Space Station
The NASA Ames WetLab-2 system was developed to offer new on-orbit gene expression analysis capabilities to ISS researchers and can be used to conduct on-orbit RNA isolation and quantitative real time PCR (RT-qPCR) analysis of gene expression from a wide range of biological samples ranging from microbes to mammalian tissues. On orbit validation included three quantitative PCR (qPCR) runs using an E. coli genomic DNA template pre-loaded at three different concentrations. The flight Ct values for the DNA standards showed no statistically significant differences relative to ground controls although there was increased noise in Ct curves, likely due to microgravity-related bubble retention in the optical windows. RNA was successfully purified from both E. coli and mouse liver samples and successfully generated singleplex, duplex and triplex data although with higher standard deviations than ground controls, also likely due to bubbles. Using volunteer science activities, a potential bubble reduction strategy was tested and resulted in smooth amplification curves and tighter Cts between replicates. The WetLab-2 validation experiment demonstrates a novel molecular biology workbench on ISS which allows scientists to purify and stabilize RNA, and to conduct RT-qPCR analyses on-orbit with rapid results. This novel ability is an important step towards utilizing ISS as a National Laboratory facility with the capability to conduct and adjust science experiments in real time without sample return, and opens new possibilities for rapid medical diagnostics and biological environmental monitoring on ISS.
microgravity↗
Antiretroviral therapy reveals triphasic decay of intact SIV genomes and persistence of ancestral variants
The decay kinetics of HIV-1-infected cells are critical to understand virus persistence. We evaluated the frequency of simian immunodeficiency virus (SIV)-infected cells for 4 years of antiretroviral therapy (ART). The intact proviral DNA assay (IPDA) and an assay for hypermutated proviruses revealed short- and long-term infected cell dynamics in macaques starting ART ~1 year after infection. Intact SIV genomes in circulating CD4+T cells showed triphasic decay with an initial phase slower than the decay of the plasma virus, a second phase faster than the second phase decay of intact HIV-1, and a stable third phase reached after 1.6–2.9 years. Hypermutated proviruses showed bi- or mono-phasic decay, reflecting different selective pressures. Viruses replicating at ART initiation had mutations conferring antibody escape. With time on ART, viruses with fewer mutations became more prominent, reflecting decay of variants replicating at ART initiation. Collectively, these findings confirm ART efficacy and indicate that cells enter the reservoir throughout untreated infection.
60 APPLIED LIFE SCIENCES↗
A genetic toolbox for metabolic engineering of Issatchenkia orientalis
The nonconventional yeast Issatchenkia orientalis can grow under highly acidic conditions and has been explored for production of various organic acids. Yet, its broader application is hampered by the lack of efficient genetic tools to enable sophisticated metabolic manipulations. We recently constructed an episomal plasmid based on the autonomously replicating sequence (ARS) from Saccharomyces cerevisiae (ScARS) in I. orientalis and developed a CRISPR/Cas9 system for multiplexed gene deletions. Here we report three additional genetic tools including: (1) identification of a 0.8 kb centromere-like (CEN-L) sequence from the I. orientalis genome by using bioinformatics and functional screening; (2) discovery and characterization of a set of constitutive promoters and terminators under different culture conditions by using RNA-Seq analysis and a fluorescent reporter; and (3) development of a rapid and efficient in vivo DNA assembly method in I. orientalis, which exhibited ~100% fidelity when assembling a 7 kb-plasmid from seven DNA fragments ranging from 0.7 kb to 1.7 kb. As evidence of concept, we used these genetic tools to rapidly construct a functional xylose utilization pathway in I. orientalis.
59 BASIC BIOLOGICAL SCIENCES↗
Data from A Genetic Toolbox for Metabolic Engineering of Issatchenkia orientalis
The nonconventional yeast Issatchenkia orientalis can grow under highly acidic conditions and has been explored for production of various organic acids. However, its broader application is hampered by the lack of efficient genetic tools to enable sophisticated metabolic manipulations. We recently constructed an episomal plasmid based on the autonomously replicating sequence (ARS) from Saccharomyces cerevisiae (ScARS) in I. orientalis and developed a CRISPR/Cas9 system for multiplexed gene deletions. Here we report three additional genetic tools including: (1) identification of a 0.8 kb centromere-like (CEN-L) sequence from the I. orientalis genome by using bioinformatics and functional screening; (2) discovery and characterization of a set of constitutive promoters and terminators under different culture conditions by using RNA-Seq analysis and a fluorescent reporter; and (3) development of a rapid and efficient in vivo DNA assembly method in I. orientalis , which exhibited ~100% fidelity when assembling a 7 kb-plasmid from seven DNA fragments ranging from 0.7 kb to 1.7 kb. As proof of concept, we used these genetic tools to rapidly construct a functional xylose utilization pathway in I. orientalis .
Conversion↗
Human APOBEC3 Induced Mutation of Human Immunodeficiency Virus Type-1 Contributes to Adaptation and Evolution in Natural Infection
Human APOBEC3 proteins are cytidine deaminases that contribute broadly to innate immunity through the control of exogenous retrovirus replication and endogenous retroelement retrotransposition. As an intrinsic antiretroviral defense mechanism, APOBEC3 proteins induce extensive guanosine-to-adenosine (G-to-A) mutagenesis and inhibit synthesis of nascent human immunodeficiency virus-type 1 (HIV-1) cDNA. Human APOBEC3 proteins have additionally been proposed to induce infrequent, potentially non-lethal G-to-A mutations that make subtle contributions to sequence diversification of the viral genome and adaptation though acquisition of beneficial mutations. Using single-cycle HIV-1 infections in culture and highly parallel DNA sequencing, we defined trinucleotide contexts of the edited sites for APOBEC3D, APOBEC3F, APOBEC3G, and APOBEC3H. We then compared these APOBEC3 editing contexts with the patterns of G-to-A mutations in HIV-1 DNA in cells obtained sequentially from ten patients with primary HIV-1 infection. Viral substitutions were highest in the preferred trinucleotide contexts of the edited sites for the APOBEC3 deaminases. Consistent with the effects of immune selection, amino acid changes accumulated at the APOBEC3 editing contexts located within human leukocyte antigen (HLA)- appropriate epitopes that are known or predicted to enable peptide binding. Thus, APOBEC3 activity may induce mutations that influence the genetic diversity and adaptation of the HIV-1 population in natural infection.
59 BASIC BIOLOGICAL SCIENCES↗
Scalable Parallel Methods for Analyzing Metagenomics Data at Extreme Scale
The field of bioinformatics and computational biology is currently experiencing a data revolution. The exciting prospect of making fundamental biological discoveries is fueling the rapid development and deployment of numerous cost-effective, high-throughput next-generation sequencing technologies. The result is that the DNA and protein sequence repositories are being bombarded with new sequence information. Databases are continuing to report a Moore’s law-like growth trajectory in their database sizes, roughly doubling every 18 months. In what seems to be a paradigm-shift, individual projects are now capable of generating billions of raw sequence data that need to be analyzed in the presence of already annotated sequence information. While it is clear that data-driven methods, such as sequencing homology detection, are becoming the mainstay in the field of computational life sciences, the algorithmic advancements essential for implementing complex data analytics at scale have mostly lagged behind. Sequence homology detection is central to a number of bioinformatics applications including genome sequencing and protein family characterization. Given millions of sequences, the goal is to identify all pairs of sequences that are highly similar (or “homologous”) on the basis of alignment criteria. While there are optimal alignment algorithms to compute pairwise homology, their deployment for large-scale is currently not feasible; instead, heuristic methods are used at the expense of quality. In this dissertation, we present the design and evaluation of a parallel implementation for conducting optimal homology detection on distributed memory supercomputers. Our approach uses a combination of techniques from asynchronous load balancing (viz. work stealing, dynamic task counters), data replication, and exact-matching filters to achieve homology detection at scale. Results for a collection of 2.56M sequences show parallel efficiencies of ~75-100% on up to 8K cores, representing a time-to-solution of 33 seconds. We extend this work with a detailed analysis of single-node sequence alignment performance using the latest CPU vector instruction set extensions. Preliminary results reveal that current sequence alignment algorithms are unable to fully utilize widening vector registers.
59 BASIC BIOLOGICAL SCIENCES↗