Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sequence analysis”

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.

At least 307 records · Page 17

Structure of Type IIL Restriction-Modification Enzyme MmeI in Complex with DNA Has Implications for Engineering New Specificities

The creation of restriction enzymes with programmable DNA-binding and -cleavage specificities has long been a goal of modern biology. The recently discovered Type IIL MmeI family of restriction-and-modification (RM) enzymes that possess a shared target recognition domain provides a framework for engineering such new specificities. However, a lack of structural information on Type IIL enzymes has limited the repertoire that can be rationally engineered. We report here a crystal structure of MmeI in complex with its DNA substrate and an S-adenosylmethionine analog (Sinefungin). The structure uncovers for the first time the interactions that underlie MmeI-DNA recognition and methylation (5’-TCCRAC-3’; R = purine) and provides a molecular basis for changing specificity at four of the six base pairs of the recognition sequence (5’-TCCRAC-3’). Surprisingly, the enzyme is resilient to specificity changes at the first position of the recognition sequence (5’-TCCRAC-3’). Collectively, the structure provides a basis for engineering further derivatives of MmeI and delineates which base pairs of the recognition sequence are more amenable to alterations than others.

59 BASIC BIOLOGICAL SCIENCES↗

A1 Signature pattern and phylogenetic analysis of full-length env genes in 20 hemophiliacs infected with Korean subclade of HIV-1 subtype B

Twenty patients with hemophilia (HP) were diagnosed with HIV-1, one to two years after exposure to domestic clotting factor IX (DCF), in 1989–90. Previous analysis of pol and vif genes confirmed that HIV-1 transmission to 20 HPs occurred through infusion of DCF. In this study, we determined full-length env gene sequences (about 2,550 bp) in 21 HPs, three plasma donors whose plasma were used in DCF production, and 45 local controls. The env gene from frozen stored sera obtained 1–3 years after diagnosis as well as from samples collected years after infection, was amplified by RT-PCR, and subjected to direct sequencing. Phylogenetic analysis revealed that all sequences were subtype B, with 109 sequences from 64 patients (20 HPs, 3 serum donors, 41 local controls, LCs) belonging to the Korean subclade (KSB) and 13 sequences from 5 patients (1HP infected outside Korea, 4 LCs) not in the KSB. Sequences of the env gene from donors O and P plus the 20 HPs comprised two subclusters within KSB. In addition, signature pattern analysis revealed signature nucleotides at 45 positions between the HPs and LCs (P < 0.05). Surprisingly, specific signature nucleotides positions were conserved 100% in clusters O and P only (at 4 and 1 positions, respectively) with none in LCs. Within both clusters, sequence identity was high. These results are consistent with our previous conclusion that these 20 HPs were infected with viruses in the DCF used for treatment. In addition, we found that there are 25 signature pattern residues in amino acids in env gene of KSB distinct from subtype B.

59 BASIC BIOLOGICAL SCIENCES↗

Unusual modifications of protein biomarkers expressed by plasmid, prophage, and bacterial host of pathogenic Escherichia coli identified using top‐down proteomic analysis

Rationale Pathogenic bacteria often carry prophage (bacterial viruses) and plasmids (small circular pieces of DNA) that may harbor toxin, antibacterial, and antibiotic resistance genes. Proteomic characterization of pathogenic bacteria should include the identification of host proteins and proteins produced by prophage and plasmid genomes. Methods Protein biomarkers of two strains of Shiga toxin–producingEscherichia coli(STEC) were identified using antibiotic induction, matrix‐assisted laser desorption/ionization tandem time‐of‐flight (MALDI‐TOF‐TOF) tandem mass spectrometry (MS/MS) with post‐source decay (PSD), top‐down proteomic (TDP) analysis, and plasmid sequencing. Alphafold2 was also used to compare predicted in silico structures of the identified proteins to prominent fragment ions generated using MS/MS‐PSD. Strain samples were also analyzed with and without chemical reduction treatment to detect the attachment of pendant groups bound by thioester or disulfide bonds. Results Shiga toxin was detected and/or identified in both STEC strains. For the first time, we also identified the osmotically inducible protein (OsmY) whose sequence unexpectedly had two forms: a full and a truncated sequence. The truncated OsmY terminates in the middle of an α‐helix as determined by Alphafold2. A plasmid‐encoded colicin immunity protein was also identified with and without attachment of an unidentified cysteine‐bound pendant group (~307 Da). Plasmid sequencing confirmed top‐down analysis and the identification of a promoter upstream of the immunity gene that is activated by antibiotic induction, that is, SOS box. Conclusions TDP analysis, coupled with other techniques (e.g., antibiotic induction, chemical reduction, plasmid sequencing, and in silico protein modeling), is a powerful tool to identify proteins (and their modifications), including prophage‐ and plasmid‐encoded proteins, produced by pathogenic microorganisms.

Biochemistry & Molecular Biology↗

Bioinformatic analysis of neurotropic HIV envelope sequences identifies polymorphisms in the gp120 bridging sheet that increase macrophage-tropism through enhanced interactions with CCR5

Macrophages express low levels of the CD4 receptor compared to T-cells. Macrophage-tropic HIV strains replicating in brain of untreated patients with HIV-associated dementia (HAD) express Envs that are adapted to overcome this restriction through mechanisms that are poorly understood. Here, bioinformatic analysis of env sequence datasets together with functional studies identified polymorphisms in the β3 strand of the HIV gp120 bridging sheet that increase M-tropism. D197, which results in loss of an N-glycan located near the HIV Env trimer apex, was detected in brain in some HAD patients, while position 200 was estimated to be under positive selection. D197 and T/V200 increased fusion and infection of cells expressing low CD4 by enhancing gp120 binding to CCR5. These results identify polymorphisms in the HIV gp120 bridging sheet that overcome the restriction to macrophage infection imposed by low CD4 through enhanced gp120–CCR5 interactions, thereby promoting infection of brain and other macrophage-rich tissues. - Highlights: • We analyze HIV Env sequences and identify amino acids in beta 3 of the gp120 bridging sheet that enhance macrophage tropism. • These amino acids at positions 197 and 200 are present in brain of some patients with HIV-associated dementia. • D197 results in loss of a glycan near the HIV Env trimer apex, which may increase exposure of V3. • These variants may promote infection of macrophages in the brain by enhancing gp120–CCR5 interactions.

60 APPLIED LIFE SCIENCES↗

Genetic Analysis of the Full-Length gag Gene from the Earliest Korean Subclade B of HIV-1: An Outbreak among Korean Hemophiliacs

We determined the earliest full-length HIV-1 gag gene sequences in 110 patients with HIV-1, including 20 hemophiliacs (HPs) and 90 local controls (LCs). The gag gene from stored sera was amplified using RT-PCR, and was subjected to direct sequencing. Phylogenetic analysis indicated that 94 and 16 sequences belonged to the Korean subclade of HIV-1 subtype B (KSB) and subtype B, respectively. A total of 12 signature pattern amino acids were found within the KSB, distinct from the worldwide consensus of subtype B. Within the KSB, the gag gene sequences from donors O and P and those from the 20 HPs comprised two subclusters. In particular, sequences from donor O strongly clustered with those of eight HPs. Moreover, signature pattern analysis indicated that 14 signature nucleotides were shared between the HPs and LCs within KSB (p < 0.01). Among the 14 nucleotides, positions 9 and 5 belonged to clusters O and P, respectively. In conclusion, signature pattern analysis for the gag gene revealed 12 signature pattern residues within the KSB and also confirmed the previous conclusion that the 20 HPs were infected with viruses due to incompletely inactivated clotting factor IX. This study is the first genetic analysis of the HIV-1 gag gene in Korea.

59 BASIC BIOLOGICAL SCIENCES↗

MerCat: a versatile k-mer counter and diversity estimator for database-independent property analysis obtained from metagenomic and/or metatranscriptomic sequencing data

MerCat is a parallel, highly scalable and modular property software package for robust analysis of features in next-generation sequencing data. MerCat inputs include assembled contigs and raw sequence reads from any platform resulting in feature abundance counts tables. MerCat allows for direct analysis of data properties without reference sequence database dependency commonly used by search tools such as BLAST and/or DIAMOND for compositional analysis of whole community shotgun sequencing (e.g. metagenomes and metatranscriptomes).

metabolomics↗

Breathing dynamics based parameter sensitivity analysis of hetero-polymeric DNA

We study the parameter sensitivity of hetero-polymeric DNA within the purview of DNA breathing dynamics. The degree of correlation between the mean bubble size and the model parameters is estimated for this purpose for three different DNA sequences. The analysis leads us to a better understanding of the sequence dependent nature of the breathing dynamics of hetero-polymeric DNA. Out of the 14 model parameters for DNA stability in the statistical Poland-Scheraga approach, the hydrogen bond interaction ε{sub hb}(AT) for an AT base pair and the ring factor ξ turn out to be the most sensitive parameters. In addition, the stacking interaction ε{sub st}(TA-TA) for an TA-TA nearest neighbor pair of base-pairs is found to be the most sensitive one among all stacking interactions. Moreover, we also establish that the nature of stacking interaction has a deciding effect on the DNA breathing dynamics, not the number of times a particular stacking interaction appears in a sequence. We show that the sensitivity analysis can be used as an effective measure to guide a stochastic optimization technique to find the kinetic rate constants related to the dynamics as opposed to the case where the rate constants are measured using the conventional unbiased way of optimization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

National Prevalence of Caprine Prion Protein Genetic Variability at Codons 146, 211, and 222 in Goat Herds in the United States

Scrapie is a neurodegenerative disease that impacts sheep and goats, characterized by gradual and progressive changes in neurological function. Recent research shows that the scrapie incubation period is significantly influenced by specific variations in amino acids within the prion protein gene (PRNP). The objective of this study was to estimate the national prevalence of caprine PRNP genetic variability at codons 146, 211, and 222 in goat populations across the United States. A total of 3052 blood, ear tissue, and brain tissue samples were collected from goats from 50 states. The participating states were categorized into four Veterinary Service (VS) district regions. The samples underwent DNA extraction, and the PRNP variants corresponding to codons 146, 211, and 222 were amplified and sequenced. The analysis of PRNP variants, when compared to the PRNP reference sequence, revealed seven alleles in twelve genotypes. The homozygous 146NN, 211RR, and 222QQ alleles, which have been linked to an increased risk of scrapie, were found to be the most prevalent among all the goats. The heterozygous 222QK, 211RQ, 146SD, 146ND, and 146NS alleles and the homozygous 222KK, 146SS, and 146DD alleles, known to be associated with reduced scrapie susceptibility and a prolonged incubation period after experimental challenge, were found in 1.098% (222QK), 2.33% (211RQ), 0.58% (146SD), 3.13% (146ND), 20.68% (146NS), 0.005% (222KK), 3.31% (146SS), and 0.67% (146DD) of goats, respectively. The 222QK allele was found most frequently in goats tested from the east (VS District 1, 1.59%) and southwest (VS District 4, 1.08%) regions, whereas the 211RQ allele was found most often in goats tested from the Midwest (VS District 2, 8.03%) and east (VS District 1, 6.53%) regions. The 146NS allele was found most frequently in goats tested from the northwest (VS District 3, 29.02%) and southwest (VS District 4, 20.69%) regions. Our results showed that the prevalence of less susceptible genotypes at PRNP codon 146 may be sufficient to use genetic susceptibility testing in some herds. This may reduce the number of goats removed as part of a herd clean-up plan and may promote the selective breeding goats for less susceptible alleles in high-risk herds at the national level.

59 BASIC BIOLOGICAL SCIENCES↗

DNA sequence determinants controlling affinity, stability and shape of DNA complexes bound by the nucleoid protein Fis

The abundant Fis nucleoid protein selectively binds poorly related DNA sequences with high affinities to regulate diverse DNA reactions. Fis binds DNA primarily through DNA backbone contacts and selects target sites by reading conformational properties of DNA sequences, most prominently intrinsic minor groove widths. High-affinity binding requires Fis-stabilized DNA conformational changes that vary depending on DNA sequence. In order to better understand the molecular basis for high affinity site recognition, we analyzed the effects of DNA sequence within and flanking the core Fis binding site on binding affinity and DNA structure. X-ray crystal structures of Fis-DNA complexes containing variable sequences in the noncontacted center of the binding site or variations within the major groove interfaces show that the DNA can adapt to the Fis dimer surface asymmetrically. We show that the presence and position of pyrimidine-purine base steps within the major groove interfaces affect both local DNA bending and minor groove compression to modulate affinities and lifetimes of Fis-DNA complexes. Sequences flanking the core binding site also modulate complex affinities, lifetimes, and the degree of local and global Fis-induced DNA bending. In particular, a G immediately upstream of the 15 bp core sequence inhibits binding and bending, and A-tracts within the flanking base pairs increase both complex lifetimes and global DNA curvatures. Taken together, our observations support a revised DNA motif specifying high-affinity Fis binding and highlight the range of conformations that Fis-bound DNA can adopt. Lastly, the affinities and DNA conformations of individual Fis-DNA complexes are likely to be tailored to their context-specific biological functions.

59 BASIC BIOLOGICAL SCIENCES↗

Comprehensive analysis of differentially expressed mRNA and circRNA in Ankylosing spondylitis patients’ platelets

Highlights: • 4996 mRNAs and 2942 circRNAs were identified differentially expressed in platelets of AS patients. • Several platelet-derived immune mediators were upregulated to regulate the interaction between platelets and immune cells. • Two downregulated circRNA were identified, and the corresponding circRNA-miRNA-mRNA regulatory network was constructed. • mRNAs and circRNAs found herein could be used as diagnostic biomarkers or new therapeutic targets. Ankylosing spondylitis (AS) is a chronic inflammatory disease significantly decreasing the quality of life. Platelets play an important and active role in the development of AS. Accumulating evidence demonstrated platelets contain diverse RNA repository inherited from megakaryocytes or microvesicles. Platelet RNAs are dynamically affected by pathological conditions and could be used as diagnostic or prognostic biomarkers. However, the role of the platelet RNAs in AS is elusive. In this study, we compared mRNA and circRNA profiles in platelets between AS patients and healthy controls using RNA sequencing and bioinformatic analysis, and found 4996 mRNAs and 2942 circRNAs were differently expressed. The significantly over-expressed mRNAs in AS patients are involved in platelet activity, gap junction, focal adhesion, rap1 and toll and Imd signaling pathway. The previous identified platelet-derived immune mediators such as P2Y1, P2Y12, PF4, GPIbα, CD40L, ICAM2, CCL5 (RANTES), TGF-β (TGF-β1 and TGF-β2) and PDGF (PDGFB and PDGFA) are also included in these over expressed mRNAs, implying these factors may trigger inflammatory cascades and promote the development of AS. Additionally, we found two down-regulated circRNA (circPTPN22 and circFCHSD2) from the intersection analyses of platelets and spinal ligament tissues of AS patients. The circRNA-miRNA-mRNA regulatory network of these two circRNAs was constructed, and the target mRNAs were enriched in Th17 cell differentiation, inflammatory bowel disease, cell adhesion molecules, cytokine-cytokine receptor interaction, Jak-STAT and Wnt signaling pathway, all these pathways participate in the bone remodeling and pro/anti-inflammatory immune regulation in AS. Then, qRT-PCR was performed to validate the expression of selected key mRNAs and circRNAs and the results demonstrated that the expression levels of P2Y12, GPIbα, circPTPN22 and circFCHSD2 were consistent with the sequencing analysis. In addition, the high expression of five predicted miRNAs interacting with circPTPN22 and circFCHSD2 were also detected in AS by qRT-PCR. Taken together, our study presents a comprehensive overview of mRNAs and circRNAs in platelets in AS patients and offers new insight into the mechanisms of platelet involving in the pathogenesis of AS. The mRNAs and circRNAs identified in this study may serve as candidates for diagnosis and targeted treatment of AS.

60 APPLIED LIFE SCIENCES↗

Molecular evolution of human coronavirus-NL63, -229E, -HKU1 and -OC43 in hospitalized children in China

Human coronaviruses (HCoVs) HCoV-NL63, HCoV-229E, HCoV-HKU1 and HCoV-OC43 have been circulated in the human population worldwide, and they are associated with a broad range of respiratory diseases with varying severity. However, there are neither effective therapeutic drugs nor licensed vaccines available for the treatment and prevention of infections by the four HCoVs. In this study, we collected nasopharyngeal aspirates of children hospitalized for respiratory tract infection in China during 2014–2018 and conducted next-generation sequencing. Sequences of four HCoVs were then selected for an in-depth analysis. Genome sequences of 2 HCoV-NL63, 8 HCoV-229E, 2 HCoV-HKU1, and 6 HCoV-OC43 were obtained. Based on the full-length S gene, a strong temporal signal was found in HCoV-229E and the molecular evolutionary rate was 6 × 10 −4 substitutions/site/year. Based on the maximum-likelihood (ML) phylogenetic tree of complete S gene, we designated H78 as a new sub-genotype C2 of HCoV-HKU1, and the obtained P43 sequence was grouped into the reported novel genotype K of HCoV-OC43 circulating in Guangzhou, China. Based on the complete genome, potential recombination events were found to occur as two phenomena, namely intraspecies and interspecies. Moreover, we observed two amino acid substitutions in the S1 subunit of obtained HCoV-NL63 (G534V) and HCoV-HKU1 (H512R), while residues 534 and 512 are important for the binding of angiotensin-converting enzyme 2 and neutralizing antibodies, respectively. Our findings might provide a clue for the molecular evolution of the four HCoVs and help in the early diagnosis, treatment and prevention of broad-spectrum HCoV infection.

Shao, Nan↗

Evolution of substrate specificity in bacterial AA10 lytic polysaccharide monooxygenases

Background: Understanding the diversity of lignocellulose-degrading enzymes in nature will provide insights for the improvement of cellulolytic enzyme cocktails used in the biofuels industry. Two families of enzymes, fungal AA9 and bacterial AA10, have recently been characterized as crystalline cellulose or chitin-cleaving lytic polysaccharide monooxygenases (LPMOs). Here we analyze the sequences, structures, and evolution of LPMOs to understand the factors that may influence substrate specificity both within and between these enzyme families. Results: Comparative analysis of sequences, solved structures, and homology models from AA9 and AA10 LPMO families demonstrated that, although these two LPMO families are highly conserved, structurally they have minimal sequence similarity outside the active site residues. Phylogenetic analysis of the AA10 family identified clades with putative chitinolytic and cellulolytic activities. Estimation of the rate of synonymous versus non-synonymous substitutions (dN/dS) within two major AA10 subclades showed distinct selective pressures between putative cellulolytic genes (subclade A) and CBP21-like chitinolytic genes (subclade D). Estimation of site-specific selection demonstrated that changes in the active sites were strongly negatively selected in all subclades. Furthermore, all codons in the subclade D had dN/dS values of less than 0.7, whereas codons in the cellulolytic subclade had dN/dS values of greater than 1.5. Positively selected codons were enriched at sites localized on the surface of the protein adjacent to the active site. Conclusions: The structural similarity but absence of significant sequence similarity between AA9 and AA10 families suggests that these enzyme families share an ancient ancestral protein. Combined analysis of amino acid sites under Darwinian selection and structural homology modeling identified a subclade of AA10 with diversifying selection at different surfaces, potentially used for cellulose-binding and protein-protein interactions. Together, these data indicate that AA10 LPMOs are under selection to change their function, which may optimize cellulolytic activity. This work provides a phylogenetic basis for identifying and classifying additional cellulolytic or chitinolytic LPMOs.

09 BIOMASS FUELS↗

Data for High Yield Production of 3-Hydroxypropionic Acid Using Issatchenkia orientalis

Biomanufacturing provides a more sustainable alternative to fossil-based chemical manufacturing. 3-Hydroxypropionic acid (3HP) is a top Department of Energy value-added chemical and precursor to bioplastics, yet cost-effective microbial production remains elusive. Here, we establish the acid-tolerant yeast Issatchenkia orientalis as a robust host for low-pH 3HP biosynthesis. Genome-scale modeling identifies the β-alanine pathway as optimal, offering the highest theoretical yield and lowest oxygen requirement. Thermodynamic analysis confirms its favorability under acidic conditions. Using sequence similarity network analysis, we discover highly active aspartate 1-decarboxylase (PAND), β-alanine-pyruvate aminotransferase (BAPAT), and 3HP dehydrogenase (YDFG), which significantly improve the pathway efficiency. Next, to further elevate the production, pathway optimization through multi-copy PAND integration, byproduct elimination (knockouts of pyruvate decarboxylase and glycerol-3-phosphate dehydrogenase), and reinforcement of aspartate flux by overexpression of pyruvate carboxylase and aspartate amino transferase improves the titer to 29 g/L in shake flasks. Fed-batch fermentation at pH 4 with low-cost corn steep liquor medium further increases the production to 92 g/L with 0.7 g/g yield and 0.55 g/L/h productivity. Techno-economic analysis indicates that such performance could potentially enable a financially viable process for sustainable acrylic acid production. This work establishes I. orientalis as a next-generation platform for cost-effective 3HP production and paves the way toward industrial commercialization.

Bioproducts↗

Sequence, structure prediction, and epitope analysis of the polymorphic membrane protein family in Chlamydia trachomatis

The polymorphic membrane proteins (Pmps) are a family of autotransporters that play an important role in infection, adhesion and immunity in Chlamydia trachomatis. Here we show that the characteristic GGA(I,L,V) and FxxN tetrapeptide repeats fit into a larger repeat sequence, which correspond to the coils of a large beta-helical domain in high quality structure predictions. Analysis of the protein using structure prediction algorithms provided novel insight to the chlamydial Pmp family of proteins. While the tetrapeptide motifs themselves are predicted to play a structural role in folding and close stacking of the beta-helical backbone of the passenger domain, we found many of the interesting features of Pmps are localized to the side loops jutting out from the beta helix including protease cleavage, host cell adhesion, and B-cell epitopes; while T-cell epitopes are predominantly found in the beta-helix itself. This analysis more accurately defines the Pmp family of Chlamydia and may better inform rational vaccine design and functional studies.

59 BASIC BIOLOGICAL SCIENCES↗

Physical, resource supply, and biological controls on nutrient processing along the river continuum

Nutrient impairment has led to damages to US surface and groundwater systems in excess of 100 billion dollars per year. Therefore, there is a strong need to develop methods to predict the transport, uptake, and export of nutrients along fluvial networks. We present results that are based on a data-driven mechanistic understanding of three factors that largely control nutrient uptake and export: 1) interactions between transport-related processes (mass transfer to metabolically active zones), 2) resource supply dynamics (nutrient concentration, stoichiometric constraints, etc.), and 3) biological controls (microbial community structure and function). Our results were generated from column experiments conducted along the Jemez River-Rio Grande continuum, which spans four orders of magnitude in mean annual discharge, more than 2000 m in altitude, and more than 500 km of stream longitude. Two resource supply injections were performed on each of the columns, i.e., a nitrate only addition, followed by a stoichiometrically ‘balanced’ 106Carbon:16Nitrogen:1Phosphorus addition. We quantified NO3-N uptake kinetics while constraining three variables: stream order, sediment type and type of injection (N vs stoichiometrically ‘balanced’ C:N:P). Following the laboratory nutrient uptake experiments, the columns were destructively sampled and the contents were homogenized to collect subsamples for DNA sequencing. Amplicon analysis was carried out as described by the Earth Microbiome Protocol for 16s and ITS sequencing.

54 ENVIRONMENTAL SCIENCES↗

Sequence Design of Random Heteropolymers as Protein Mimics

Random heteropolymers (RHPs) have been computationally designed and experimentally shown to recapitulate protein-like phase behavior and function. However, unlike proteins, RHP sequences are only statistically defined and cannot be sequenced. Recent developments in reversible-deactivation radical polymerization allowed simulated polymer sequences based on the well-established Mayo–Lewis equation to more accurately reflect ground-truth sequences that are experimentally synthesized. This led to opportunities to perform bioinformatics-inspired analysis on simulated sequences to guide the design, synthesis, and interpretation of RHPs. We compared batches on the order of 10000 simulated RHP sequences that vary by synthetically controllable and measurable RHP characteristics such as chemical heterogeneity and average degree of polymerization. Our analysis spans across 3 levels: segments along a single chain, sequences within a batch, and batch-averaged statistics. We discuss simulator fidelity and highlight the importance of robust segment definition. Examples are presented that demonstrate the use of simulated sequence analysis for in-silico iterative design to mimic protein hydrophobic/hydrophilic segment distributions in RHPs and compare RHP and protein sequence segments to explain experimental results of RHPs that mimic protein function. To facilitate the community use of this workflow, the simulator and analysis modules have been made available through an open source toolkit, the RHPapp.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A high-resolution single-molecule sequencing-based Arabidopsis transcriptome using novel methods of Iso-seq analysis

Accurate and comprehensive annotation of transcript sequences is essential for transcript quantification and differential gene and transcript expression analysis. Single-molecule long-read sequencing technologies provide improved integrity of transcript structures including alternative splicing, and transcription start and polyadenylation sites. However, accuracy is significantly affected by sequencing errors, mRNA degradation, or incomplete cDNA synthesis. We present a new and comprehensive Arabidopsis thaliana Reference Transcript Dataset 3 (AtRTD3). AtRTD3 contains over 169,000 transcripts—twice that of the best current Arabidopsis transcriptome and including over 1500 novel genes. Seventy-eight percent of transcripts are from Iso-seq with accurately defined splice junctions and transcription start and end sites. We develop novel methods to determine splice junctions and transcription start and end sites accurately. Mismatch profiles around splice junctions provide a powerful feature to distinguish correct splice junctions and remove false splice junctions. Stratified approaches identify high-confidence transcription start and end sites and remove fragmentary transcripts due to degradation. AtRTD3 is a major improvement over existing transcriptomes as demonstrated by analysis of an Arabidopsis cold response RNA-seq time-series. AtRTD3 provides higher resolution of transcript expression profiling and identifies cold-induced differential transcription start and polyadenylation site usage. AtRTD3 is the most comprehensive Arabidopsis transcriptome currently. It improves the precision of differential gene and transcript expression, differential alternative splicing, and transcription start/end site usage analysis from RNA-seq data. The novel methods for identifying accurate splice junctions and transcription start/end sites are widely applicable and will improve single-molecule sequencing analysis from any species.

transcription start and end sites↗