Search NASA⌕ Search

SEARCH · Search NASA

Results for “Sequence Alignment”

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 469 records · Page 26

BSSD 2019 Performance Metric Q4

The LANL SFA in Terrestrial Microbial Carbon Cycling aims to inform climate modeling and enable carbon management in terrestrial ecosystems by discovering widespread biological processes that control carbon storage and release in temperate biome soils. To achieve these goals, computational approaches are essential to analyze increasingly complex ‘omics’ data from microbial communities in environmental samples. The SFA continues to leverage advances in ‘omics' measurement technology and computation tools to decipher microbial community function in environmental samples. The computational approach is aligned with the recent shift in the SFA’s research approach. When the SFA began 10 years ago, computational approaches principally addressed “who is there” with extremely limited insight into function. Accessible ‘omics' data for complex soil communities were primarily amplicon DNA sequences of taxonomic or functional marker genes obtained by PCR, cloning, and Sanger sequencing. Computational tools for community analysis were nascent. Over the past decade, genome resources (DOE-Integrated Microbial Genomics database) and data acquisition increased 1000-fold, common computational pipelines for many analysis tasks have arisen, and the variety of ‘omics’ measurements has expanded to routinely include shotgun metagenomics, metatranscriptomics, and versions of metabolomics. The SFA is capitalizing on these advances while also developing new computational techniques to decipher how communities with different patterns of carbon cycling function from the species to ecosystem level.

54 ENVIRONMENTAL SCIENCES↗

Evolution of intermittent filaments in the scrape-off layer of NSTX

Filamentary structures naturally arise from background turbulence in the scrape-off layer (SOL) of plasmas, leading to significant particle and heat transport that can degrade overall plasma confinement. This enhanced transport can contribute to unacceptably high heat loads on plasma-facing components. As such, understanding the physics of SOL plasma filaments is critical for predicting and mitigating their effects in future fusion devices. On the National Spherical Torus Experiment (NSTX), plasma filaments—commonly referred to as “blobs”—were investigated using the gas-puff imaging (GPI) diagnostic in the edge and SOL regions. The analysis involved identifying, segmenting, and tracking the characteristic contours of the blobs in each frame of the GPI video sequences. Their evolution was characterized through shape descriptors, velocity, and angular velocity derived from their contour coordinates. The results indicate that as the blob area increases, their shapes become more concave and less circular, suggesting reduced structural stability in larger blobs. This result aligns with previous theoretical results where it was shown that larger blobs are more susceptible to instabilities [Krasheninnikov et al., J. Plasma Phys. 74, 679–717 (2008) and D'Ippolito et al., Phys. Plasmas 18, 060501 (2011)]. A positive correlation was observed between radial velocity and radial position, suggesting radially outward acceleration of the filaments, potentially driven by decreasing viscous drag toward the far SOL. Interestingly, blobs in background SOL turbulence exhibited minimal spinning in contrast to filaments originating from edge localized modes, which show substantial rotation during their paths [Lampert et al., Phys. Plasmas 29, 102502 (2022)]. Statistical analysis of the solidity and total curvature shape descriptors, along with their temporal evolution, revealed relatively broad, near-Gaussian distributions. This suggests that blob morphology is strongly influenced by stochastic turbulent processes in the surrounding plasma environment. Blob parameters were also compared with bulk plasma and radial profile measurements. Notable trends were found between blob rotation and poloidal velocity with collisionality and line-integrated density. These findings contribute to a deeper understanding of blob dynamics and provide valuable insights for refining SOL turbulence models.

Covariance and correlation↗

Gaia: An AI-enabled genomic context–aware platform for protein sequence annotation

Protein sequence similarity search is fundamental to biology research, but current methods are typically not able to consider crucial genomic context information indicative of protein function, especially in microbial systems. Here, we present Gaia (Genomic AI Annotator), a sequence annotation platform that enables rapid, context-aware protein sequence search across genomic datasets. Gaia leverages gLM2, a mixed-modality genomic language model trained on both amino acid sequences and their genomic neighborhoods to generate embeddings that integrate sequence-structure-context information. This approach allows for the identification of functionally and/or evolutionarily related genes that are found in conserved genomic contexts, which may be missed by traditional sequence- or structure-based search alone. Gaia enables real-time search of a curated database comprising more than 85 million protein clusters from 131,744 microbial genomes. We compare the homolog retrieval performance of Gaia search against other embedding and alignment-based approaches. We provide Gaia as a web-based, freely available tool.

Jha, Nishant↗

Gene regulatory networks controlling carbohydrate selective deconstruction pathways in fungi

The work performed at Clark University (PI, David Hibbett) was a subcomponent of the research lead by Jonathan Schilling at University of Minnesota (award number DE-SC0019427). The overall goal of this project was to address how fungi metabolize lignocellulose. Specifically, we targeted a carbohydrate-selective group of fungi known as 'brown rot' fungi that are well-aligned to be harness for bioconversion to biofuels and other products. These organisms orchestrate a two-step mechanism that initiates decomposition using oxidative radical-based chemistry and follows with enzymatic attack. We have helped characterize that this mechanism occurs in this sequence, but currently lack the fundamental regulatory information needed to harness these fungi for application. Thus, the major goal of this project was to use systems biology to create an integrated model of brown rot wood decay that couples fungal gene regulation with the metabolite feedbacks used to coordinate a two-step oxidative-hydrolytic mechanism. The research at Clark University focused on phylogenomic analyses of data produced by other project partners. Hence, the Clark component of the research was staged to initiate somewhat after that in other partner laboratories. The major Clark-based accomplishment was to advise on selection of fungal species for analyses of brown rot decay chemistry, and to interpret results in a phylogenetic context. The Clark University group also performed complementary analyses of gene families encoding plant cell wall-degrading enzymes in the the basidiomycete genus Lentinula, as well as analyses of substrate-specific gene expression and RNA editing in brown rot fungi.

59 BASIC BIOLOGICAL SCIENCES↗

A YOUNG STELLAR CLUSTER WITHIN THE RCW41 H II REGION: DEEP NIR PHOTOMETRY AND OPTICAL/NIR POLARIMETRY

The RCW41 star-forming region is embedded within the Vela Molecular Ridge, hosting a massive stellar cluster surrounded by a conspicuous H II region. Understanding the role of interstellar magnetic fields and studying the newborn stellar population is crucial to building a consistent picture of the physical processes acting on this kind of environment. We carried out a detailed study of the interstellar polarization toward RCW41 with data from an optical and near-infrared polarimetric survey. Additionally, deep near-infrared images from the 3.5 meter New Technology Telescope were used to study the photometric properties of the embedded young stellar cluster, revealing several YSO candidates. By using a set of pre-main-sequence isochrones, a mean cluster age in the range 2.5-5.0 million years was determined, and evidence of sequential star formation was revealed. An abrupt decrease in R-band polarization degree was noticed toward the central ionized area, probably due to low grain alignment efficiency caused by the turbulent environment and/or the weak intensity of magnetic fields. The distortion of magnetic field lines exhibits dual behavior, with the mean orientation outside the area approximately following the borders of the star-forming region and directed radially toward the cluster inside the ionized area, in agreement with simulations of expanding H II regions. The spectral dependence of polarization allowed a meaningful determination of the total-to-selective extinction ratio by fittings of the Serkowski relation. Furthermore, a large rotation of polarization angle as a function of wavelength was detected toward several embedded stars.

79 ASTRONOMY AND ASTROPHYSICS↗

Unpinning the skyrmion lattice in MnSi: Effect of substitutional disorder

By employing magnetization and small angle neutron scattering measurements, we have investigated the behavior of the skyrmion lattice (SKL) and the helical order in MnSi 0.992 Ga 0.008 . Our results indicate that the order of the SKL is sensitive to the orientation of an applied magnetic field with respect to the crystal lattice and to variations in the sequence of small temperature and applied magnetic field changes. The disorder caused by the substitution of the heavier element Ga for Si is sufficient to reduce the pinning of the SKL to the underlying crystalline lattice, reducing the propensity for the SKL to be aligned with the crystal lattice. This tendency is most evident when the applied field is not well oriented with respect to the high symmetry axes of the crystal resulting in disorder in the long range SKL while maintaining sharp short range (radial) order. We have also investigated the effect of substituting heavier elements into MnSi on the reorientation process of the helical domains with field cycling in MnSi 0.992 Ga 0.008 and Mn 0.985 Ir 0.015 Si A comparison of the reorientation process in these materials with field reduction indicates that the substitution of heavier elements on either Mn or Si sites creates a higher energy barrier for the reorientation of the helical order and for the formation of domains.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Genomic Encyclopedia of Bacterial and Archaeal Type Strains, Phase III: the genomes of soil and plant-associated and newly described type strains

The Genomic Encyclopedia of Bacteria and Archaea (GEBA) project was launched by the JGI in 2007 as a pilot project to sequence about 250 bacterial and archaeal genomes of elevated phylogenetic diversity. Here in this paper, we propose to extend this approach to type strains of prokaryotes associated with soil or plants and their close relatives as well as type strains from newly described species. Understanding the microbiology of soil and plants is critical to many DOE mission areas, such as biofuel production from biomass, biogeochemistry, and carbon cycling. We are also targeting type strains of novel species while they are being described. Since 2006, about 630 new species have been described per year, many of which are closely aligned to DOE areas of interest in soil, agriculture, degradation of pollutants, biofuel production, biogeochemical transformation, and biodiversity

59 BASIC BIOLOGICAL SCIENCES↗

Effect of PbI{sub 2} deposition rate on two-step PVD/CVD all-vacuum prepared perovskite

In this work we fabricate all-vacuum processed methyl ammonium lead halide perovskite by a sequence of physical vapour deposition of PbI{sub 2} and chemical vapour deposition (CVD) of CH{sub 3}NH{sub 3}I under a static atmosphere. We demonstrate that for higher deposition rate the (001) planes of PbI{sub 2} film show a higher degree of alignment parallel to the sample's surface. From X-ray diffraction data of the resulted perovskite film we derive that the intercalation rate of CH{sub 3}NH{sub 3}I is fostered for PbI{sub 2} films with higher degree of (001) planes alignment. The stoichiometry of the produced perovskite film is also studied by Hard X-ray photoelectron spectroscopy measurements. Complete all-vacuum perovskite solar cells were fabricated on glass/ITO substrates coated by an ultra-thin (5 nm) Zn-phthalocyanine film as hole selective layer. A dependence of residual PbI{sub 2} on the solar cells performance is displayed, while photovoltaic devices with efficiency up to η=11.6% were achieved. - Graphical abstract: A two-step PVD/CVD processed perovskite film with the CVD intercalation rate of CH{sub 3}NCH{sub 3} molecules been fostered by increasing the PVD rate of PbI{sub 2} and prolonging the CVD time. - Highlights: • A simple PVD/CVD process for perovskite film production. • Increased PVD rate yields better alignment of the PbI{sub 2} (001) crystallite planes. • CH{sub 3}NH{sub 3}I intercalation process fostered by increased PbI{sub 2} PVD rate. • Stoichiometric CH{sub 3}NH{sub 3}PbI{sub 3} suitable as absorber in photovoltaic applications • Reduced PbI{sub 2} residue at the bottom of CH{sub 3}NH{sub 3}PbI{sub 3} improves device performance.

36 MATERIALS SCIENCE↗

Electrostatic carrier doping of GdTiO{sub 3}/SrTiO{sub 3} interfaces

Heterostructures and superlattices consisting of a prototype Mott insulator, GdTiO{sub 3}, and the band insulator SrTiO{sub 3} are grown by molecular beam epitaxy and show intrinsic electronic reconstruction, approximately 1/2 electron per surface unit cell at each GdTiO{sub 3}/SrTiO{sub 3} interface. The sheet carrier densities in all structures containing more than one unit cell of SrTiO{sub 3} are independent of layer thicknesses and growth sequences, indicating that the mobile carriers are in a high concentration, two-dimensional electron gas bound to the interface. These carrier densities closely meet the electrostatic requirements for compensating the fixed charge at these polar interfaces. Based on the experimental results, insights into interfacial band alignments, charge distribution, and the influence of different electrostatic boundary conditions are obtained.

36 MATERIALS SCIENCE↗

BiG-SCAPE 2.0 and BiG-SLiCE 2.0: scalable, accurate and interactive sequence clustering of metabolic gene clusters

Microbial metabolic gene clusters encode the biosynthesis or catabolism of metabolites that facilitate ecological specialization, mediate microbiome interactions and constitute a major source of medicines and crop protection agents. Here, we present BiG-SCAPE and BiG-SLiCE 2.0, next-generation methods that facilitate scalable, accurate and interactive gene cluster analyses. BiG-SCAPE 2.0 updates its classification, alignment methods, and visualizations, enabling more accurate analysis, up to 8x faster runtimes and halved memory requirements. BiG-SLiCE 2.0 updates its distance metric, pHMM database, and classification logic, resulting in increased sensitivity nearing that of BiG-SCAPE. Analysis of 260,630 biosynthetic gene clusters from publicly available genomes reveals that both tools generate concurring estimates of gene cluster diversity, thus providing significantly extended methodological support for recent evidence indicating that the vast majority of natural product diversity remains unexplored. Together, these updates will facilitate global genome mining efforts for natural product discovery and microbiome analyses scalable with current data sizes.

Draisma, Arjan [Wageningen University & Research (↗

Differentiable hybrid neural network approach for enhancing reactor dynamics simulations

Reactor dynamics simulations provide essential insights into the time-dependent behavior of nuclear reactors under various operating conditions. However, high-fidelity simulations can be computationally intensive, requiring significant computational resources. Here, to address this challenge, this study employs a differentiable hybrid model that utilizes neural networks as a corrector to enhance the performance of a low-fidelity simulation, aligning its predictions with those of a high-fidelity simulation. Low-fidelity and high-fidelity simulations were obtained by adjusting the mesh size in the System Dynamics Analysis Tool. The differentiable hybrid model was trained in two approaches: time-step-wise and sequence-wise. It was then applied to simulate various transients in a molten salt reactor. Its performance was evaluated by comparing its responses to transients against those of the high-fidelity simulation. An additional approach was performed using a data-driven model to correct the low-fidelity simulation. In comparison, the differentiable hybrid model showed significant improvements in transient prediction, effectively addressing the limitations of the low-fidelity simulations. The results highlighted the robustness of the differentiable hybrid model in both training approaches. It delivered simulations that were at least 3.8 times faster than high-fidelity models. In the time-step-wise approach, it achieved at least a 39% improvement in accuracy. In the sequence-wise approach, it showed at least an 81% accuracy improvement over the full transient. This approach offers a promising path for improving computational efficiency without compromising accuracy in nuclear reactor simulations, making it suitable for real-time digital twin applications.

42 - ENGINEERING↗

Vaccine potential of ESAT-6 protein fused with consensus CD4+ T-cell epitopes of PE/PPE proteins against highly pathogenic Mycobacterium tuberculosis strain HN878

Pro-Glu/Pro-Pro-Glu (PE/PPE) family proteins in Mycobacterium tuberculosis (Mtb) are contributors to pathogenesis and immune evasion. These proteins have a unique structure in which the sequence is conserved. We investigated the vaccine potential of ESAT-6 fused with consensus CD4{sup +} T-cell epitopes of PE/PPE proteins against highly pathogenic Mtb strain HN878 in a murine model. We selected consensus CD4{sup +} T-cell epitopes of PE/PPE proteins by multiple alignments, investigated their IFN-γ response during Mtb infection, and produced their fused ESAT-6 vaccine antigens. Our results showed an increased immune response in PE/PPE peptide -ESAT-6 fusion protein immunization group compared to ESAT-6 only immunization group. After challenge with Mtb strain HN878, we observed that induced CD4{sup +} T-cells secreted double-positive cytokine IL-2{sup +}/IFN-γ{sup +}, which is considered to be associated with protective T-cell immunity. Additionally, lower numbers of colony-forming units were observed in the spleen of fusion protein immunization groups than in those of single ESAT-6 group. Therefore, conjugation of consensus CD4{sup +} T-cell epitopes in N terminus of PE/PPE to vaccine antigens could potentially increase the protective efficacy of subunit vaccine.

60 APPLIED LIFE SCIENCES↗

P finder: genomic and metagenomic annotation of RNase P RNA gene (rnpB)

Abstract Background The rnpB gene encodes for an essential catalytic RNA (RNase P). Like other essential RNAs, RNase P’s sequence is highly variable. However, unlike other essential RNAs (i.e. tRNA, 16 S, 6 S,...) its structure is also variable with at least 5 distinct structure types observed in prokaryotes. This structural variability makes it labor intensive and challenging to create and maintain covariance models for the detection of RNase P RNA in genomic and metagenomic sequences. The lack of a facile and rapid annotation algorithm has led to the rnpB gene being the most grossly under annotated essential gene in completed prokaryotic genomes with only a 24% annotation rate. Here we describe the coupling of the largest RNase P RNA database with the local alignment scoring algorithm to create the most sensitive and rapid prokaryote rnpB gene identification and annotation algorithm to date. Results Of the 2772 completed microbial genomes downloaded from GenBank only 665 genomes had an annotated rnpB gene. We applied P Finder to these genomes and were able to identify 2733 or nearly 99% of the 2772 microbial genomes examined. From these results four new rnpB genes that encode the minimal T-type P RNase P RNAs were identified computationally for the first time. In addition, only the second C-type RNase P RNA was identified in Sphaerobacter thermophilus . Of special note, no RNase P RNAs were detected in several obligate endosymbionts of sap sucking insects suggesting a novel evolutionary adaptation. Conclusions The coupling of the largest RNase P RNA database and associated structure class identification with the P Finder algorithm is both sensitive and rapid, yielding high quality results to aid researchers annotating either genomic or metagenomic data. It is the only algorithm to date that can identify challenging RNAse P classes such as C-type and the minimal T-type RNase P RNAs. P Finder is written in C# and has a user-friendly GUI that can run on multiple 64-bit windows platforms (Windows Vista/7/8/10). P Finder is free available for download at https://github.com/JChristopherEllis/P-Finder as well as a small sample RNase P RNA file for testing.

59 BASIC BIOLOGICAL SCIENCES↗

Using Calibrated Sodium Data for Preliminary Validation of the SRT Code for Advanced Reactors

Various types of non-light water reactors are currently engaged in the U.S. licensing process. Because of inherent differences compared with well-established large light water reactors, appropriate assessment tools are needed. Specifically, source term analysis, which determines environmental dose impacts from potential accident scenarios, is a crucial part of design and licensing. The U.S. Nuclear Regulatory Commission has emphasized the importance of mechanistic source term analysis for advanced reactor deployments. To align with these needs, Argonne National Laboratory has developed the Simplified Radionuclide Transport (SRT) source term analysis code for metal fuel Sodium-cooled Fast Reactors (SFRs) and microreactors. SRT conducts time-dependent radionuclide transport and retention in SFRs for core and ex-core radionuclide source accident sequences. The main objective of SRT is to provide rapid sensitivity and uncertainty analyses, incorporating parametric uncertainties and summarizing probabilistic results. As part of the code validation process, a study focused on the bubble scrubbing module was performed using an experiment recently carried out by the University of Wisconsin-Madison. Based on the analysis, the modeling approach in SRT provides accurate results for small and large aerosols, while slight underprediction of radionuclide aerosol removal are observed for medium sized aerosols. However, the deviation is minor, considering the highly uncertain phenomenon and range of results, and is in the conservative direction. In addition, uncertainty information derived from the experiments is further implemented, reflecting the actual span of parameters, which leads to enhanced agreement with code predictions. The results demonstrate that SRT provides reasonable predictions for the bubble scrubbing process in sodium pool.

Kam, Dong Hoon↗

Crystal Structure of the HEAT Domain from the Pre-mRNA Processing Factor Symplekin

The majority of eukaryotic pre-mRNAs are processed by 3'-end cleavage and polyadenylation, although in metazoa the replication-dependent histone mRNAs are processed by 3'-end cleavage but not polyadenylation. The macromolecular complex responsible for processing both canonical and histone pre-mRNAs contains the {approx} 1160-residue protein Symplekin. Secondary-structural prediction algorithms identified putative HEAT domains in the 300 N-terminal residues of all Symplekins of known sequence. The structure and dynamics of this domain were investigated to begin elucidating the role Symplekin plays in mRNA maturation. The crystal structure of the Drosophila melanogaster Symplekin HEAT domain was determined to 2.4 {angstrom} resolution with single-wavelength anomalous dispersion phasing methods. The structure exhibits five canonical HEAT repeats along with an extended 31-amino-acid loop (loop 8) between the fourth and fifth repeat that is conserved within closely related Symplekin sequences. Molecular dynamics simulations of this domain show that the presence of loop 8 dampens correlated and anticorrelated motion in the HEAT domain, therefore providing a neutral surface for potential protein-protein interactions. HEAT domains are often employed for such macromolecular contacts. The Symplekin HEAT region not only structurally aligns with several established scaffolding proteins, but also has been reported to contact proteins essential for regulating 3'-end processing. Together, these data support the conclusion that the Symplekin HEAT domain serves as a scaffold for protein-protein interactions essential to the mRNA maturation process.

36 MATERIALS SCIENCE↗

Complex Fault Geometry of the 2020 M ww 6.5 Monte Cristo Range, Nevada, Earthquake Sequence

On 15 May 2020 an M ww 6.5 earthquake occurred beneath the Monte Cristo Range in the Mina Deflection region of western Nevada. Rapid deployment of eight temporary seismic stations enabled detailed analysis of its productive and slowly decaying aftershock sequence (p=0.8), which included ~18,000 autodetected events in 3.5 months. Double-difference, waveform-based relative relocation of 16,714 earthquakes reveals a complex network of faults, many of which cross the inferred 35-km-long east–northeast-striking, left-lateral mainshock rupture. Seismicity aligns with left-lateral, right-lateral, and normal mechanism moment tensors of 128 of the largest earthquakes. The mainshock occurred near the middle of the aftershock zone at the intersection of two distinct zones of seismicity. In the western section, numerous subparallel, shallow, north-northeast-striking faults form a broad flower-structure-like fault mesh that coalesces at depth into a near-vertical, left-lateral fault. We infer the near-vertical fault to be a region of significant slip in the mainshock and an eastward extension of the left-lateral Candelaria fault. Near the mainshock hypocenter, seismicity occurs on a northeast-striking, west-dipping structure that extends north from the eastern Columbus Salt Marsh normal fault. Together, these two intersecting structures bound the Columbus Salt Marsh tectonic basin. East of this intersection and the mainshock hypocenter, seismicity occurs in a narrow, near-vertical, east-northeast-striking fault zone through to its eastern terminus. At the eastern end, the aftershock zone broadens and extends northwest toward the southern extension of the northwest-striking, right-lateral Petrified Springs fault system. The eastern section hosts significantly fewer aftershocks than the western section, but has more moment release. We infer that shallow aftershocks throughout the system highlight fault-fracture meshes that connect mapped fault systems at depth. Here, comparing earthquake data with surface ruptures and a simple geodetic fault model sheds light on the complexity of this recent M 6.5 Walker Lane earthquake.

58 GEOSCIENCES↗

Structure and Function of the N-Terminal Domain of the Vesicular Stomatitis Virus RNA Polymerase

Viruses have various mechanisms to duplicate their genomes and produce virus-specific mRNAs. Negative-strand RNA viruses encode their own polymerases to perform each of these processes. For the nonsegmented negative-strand RNA viruses, the polymerase is comprised of the large polymerase subunit (L) and the phosphoprotein (P). L proteins from members of the Rhabdoviridae, Paramyxoviridae, and Filoviridae share sequence and predicted secondary structure homology. Here, we present the structure of the N-terminal domain (conserved region I) of the L protein from a rhabdovirus, vesicular stomatitis virus, at 1.8-Å resolution. The strictly and strongly conserved residues in this domain cluster in a single area of the protein. Serial mutation of these residues shows that many of the amino acids are essential for viral transcription but not for mRNA capping. Three-dimensional alignments show that this domain shares structural homology with polymerases from other viral families, including segmented negative-strand RNA and double-stranded RNA (dsRNA) viruses. Negative-strand RNA viruses include a diverse set of viral families that infect animals and plants, causing serious illness and economic impact. Furthermore, the members of this group of viruses share a set of functionally conserved proteins that are essential to their replication cycle. Among this set of proteins is the viral polymerase, which performs a unique set of reactions to produce genome- and subgenome-length RNA transcripts. In this article, we study the polymerase of vesicular stomatitis virus, a member of the rhabdoviruses, which has served in the past as a model to study negative-strand RNA virus replication. We have identified a site in the N-terminal domain of the polymerase that is essential to viral transcription and that shares sequence homology with members of the paramyxoviruses and the filoviruses. Newly identified sites such as that described here could prove to be useful targets in the design of new therapeutics against negative-strand RNA viruses.

59 BASIC BIOLOGICAL SCIENCES↗

Transcriptome sequencing and differential gene expression analysis in Viola yedoensis Makino (Fam. Violaceae) responsive to cadmium (Cd) pollution

Viola yedoensis Makino is an important Chinese traditional medicine plant adapted to cadmium (Cd) pollution regions. Illumina sequencing technology was used to sequence the transcriptome of V. yedoensis Makino. We sequenced Cd-treated (VIYCd) and untreated (VIYCK) samples of V. yedoensis, and obtained 100,410,834 and 83,587,676 high quality reads, respectively. After de novo assembly and quantitative assessment, 109,800 unigenes were finally generated with an average length of 661 bp. We then obtained functional annotations by aligning unigenes with public protein databases including NR, NT, SwissProt, KEGG and COG. In addition, 892 differentially expressed genes (DEGs) were investigated between the two libraries of untreated (VIYCK) and Cd-treated (VIYCd) plants. Moreover, 15 randomly selected DEGs were further validated with qRT-PCR and the results were highly accordant with the Solexa analysis. This study firstly generated a successful global analysis of the V. yedoensis transcriptome and it will provide for further studies on gene expression, genomics, and functional genomics in Violaceae. - Highlights: • A de novo assembly generated 109,800 unigenes and 5,4479 of them were annotated. • 31,285 could be classified into 26 COG categories. • 263 biosynthesis pathways were predicted and classified into five categories. • 892 DEGs were detected and 15 of them were validated by qRT-PCR.

60 APPLIED LIFE SCIENCES↗