De Novo Sequencing Research Articles

Leaf Transcriptome Sequencing for Identifying Genic-SSR Markers and SNP Heterozygosity in Crossbred Mango Variety 'Amrapali' (Mangifera indica L.).

Mango (Mangifera indica L.) is called “king of fruits” due to its sweetness, richness of taste, diversity, large production volume and a variety of end usage. Despite its huge economic importance genomic resources in mango are scarce and genetics of useful horticultural traits are poorly understood. Here we generated deep coverage leaf RNA sequence data for mango parental varieties ‘Neelam’, ‘Dashehari’ and their hybrid ‘Amrapali’ using next generation sequencing technologies. De-novo sequence assembly generated 27,528, 20,771 and 35,182 transcripts for the three genotypes, respectively. The transcripts were further assembled into a non-redundant set of 70,057 unigenes that were used for SSR and SNP identification and annotation. Total 5,465 SSR loci were identified in 4,912 unigenes with 288 type I SSR (n ≥ 20 bp). One hundred type I SSR markers were randomly selected of which 43 yielded PCR amplicons of expected size in the first round of validation and were designated as validated genic-SSR markers. Further, 22,306 SNPs were identified by aligning high quality sequence reads of the three mango varieties to the reference unigene set, revealing significantly enhanced SNP heterozygosity in the hybrid Amrapali. The present study on leaf RNA sequencing of mango varieties and their hybrid provides useful genomic resource for genetic improvement of mango.

Recombinant production, purification and characterization of vessel dilator in E. coli

Vessel dilator is a 3.9-KDa potent anticancer peptide and a valuable candidate in the treatment of conditions such as congestive heart failure and acute renal failure amongst others. Here we report the recombinant production of vessel dilator in Escherichia coli. Three different synthetic ORF's dubbed VDI, VDII and VDIII, each encoding a trimmer of the vessel dilator peptide attached to a His tag sequence at their C- terminal, were synthesized and placed in pET21c expression vectors. The highest yield, following expression in E. coli BL21 (DE3), was recorded with VDII that carried the shortest fusion partner. Subsequent to the initial capture of the fusion protein by a Ni affinity column, the vessel dilator monomers were cleaved by trypsin treatment, and further purified to at least 90% homogeneity by anion exchange chromatography. De-novo sequencing and in vivo anticancer activity tests were used to verify the peptide sequence and its biological activity, respectively. The final yield was estimated to be approximately 15 mg of the purified vessel dilator per gram wet weight of the bacterial cells.

Genetic structure and diversity of the selfing model grass Brachypodium stacei (Poaceae) in Western Mediterranean: out of the Iberian Peninsula and into the islands.

Annual Mediterranean species of the genus Brachypodium are promising model plants for energy crops since their selfing nature and short-life cycles are an advantage in breeding programs. The false brome, B. distachyon, has already been sequenced and new genomic initiatives have triggered the de-novo genome sequencing of its close relatives such as B. stacei, a species that was until recently mistaken for B. distachyon. However, the success of these initiatives hinges on detailed knowledge about the distribution of genetic variation within and among populations for the effective use of germplasm in a breeding program. Understanding population genetic diversity and genetic structure is also an important prerequisite for designing effective experimental populations for genomic wide studies. However, population genetic data are still limited in B. stacei. We therefore selected and amplified 10 nuclear microsatellite markers to depict patterns of population structure and genetic variation among 181 individuals from 19 populations of B. stacei occurring in its predominant range, the western Mediterranean area: mainland Iberian Peninsula, continental Balearic Islands and oceanic Canary Islands. Our genetic results support the occurrence of a predominant selfing system with extremely high levels of homozygosity across the analyzed populations. Despite the low level of genetic variation found, two different genetic clusters were retrieved, one clustering all SE Iberian mainland populations and the island of Minorca and another one grouping all S Iberian mainland populations, the Canary Islands and all Majorcan populations except one that clustered with the former group. These results, together with a high sharing of alleles (89%) suggest different colonization routes from the mainland Iberian Peninsula into the islands. A recent colonization scenario could explain the relatively low levels of genetic diversity and low number of alleles found in the Canary Islands populations while older colonization events are hypothesized to explain the high genetic diversity values found in the Majorcan populations. Our study provides widely applicable information about geographical patterns of genetic variation in B. stacei. Among others, the genetic pattern and the existence of local alleles will need to be adequately reflected in the germplasm collection of B. stacei for efficient genome wide association studies.

Time-space analysis of highly pathogenic avian influenza H5N2 outbreak in the US.

BackgroundIn early 2015, highly pathogenic avian influenza H5N2 caused outbreaks in commercial poultry farms in Minnesota and neighboring states where more than 48 million birds were affected. To date, the origin and transmission pathways of HPAI H5N2 have not been conclusively established.MethodsIn this study, we analyzed forty-six samples from turkeys and their environment that were collected at different time-points of the outbreak to identify origins and within outbreak evolutionary changes. We performed de-novo whole genome sequencing from primary samples and the most recent common ancestors of the PB2, PA, HA5, M and NS segments were traced back to Japanese HPAI H5N8 isolates. These segments appeared to have diverged from the ancestor around June and November 2014.ResultsThe time to most recent common ancestor analysis for PB1, NP and NA2 segments suggest two likely possibilities of reassortant HPAI H5N2 origin - either a reassortment in Alaska area or multiple reassortments with North American low pathogenic avian influenza strains, before the HPAI H5N2 outbreak strain emerged. Within the outbreak, viruses clustered into two and three subgroups suggesting high substitution rates of 0.702x10-2 - 1.665x10-2 (subs/site/year), over the 5-month outbreak period.ConclusionsData are suggestive of a fast evolving HPAI strain within an outbreak that should be taken into consideration in developing appropriate control strategies in the future.Electronic supplementary materialThe online version of this article (doi:10.1186/s12985-016-0605-4) contains supplementary material, which is available to authorized users.

Sequence protein identification by randomized sequence database and transcriptome mass spectrometry (SPIDER-TMS): from manual to automatic application of a 'denovo sequencing' approach.

Sequence protein identification by a randomized sequence database and transcriptome mass spectrometry software package has been developed at the University of Basilicata in Potenza (Italy) and designed to facilitate the determination of the amino acid sequence of a peptide as well as an unequivocal identification of proteins in a high-throughput manner with enormous advantages of time, economical resource and expertise. The software package is a valid tool for the automation of a denovo sequencing approach, overcoming the main limits and a versatile platform useful in the proteomic field for an unequivocal identification of proteins, starting from tandem mass spectrometry data. The strength of this software is that it is a user-friendly and non-statistical approach, so protein identification can be considered unambiguous.

De Novo Sequencing and Analysis of Lemongrass Transcriptome Provide First Insights into the Essential Oil Biosynthesis of Aromatic Grasses.

Aromatic grasses of the genus Cymbopogon (Poaceae family) represent unique group of plants that produce diverse composition of monoterpene rich essential oils, which have great value in flavor, fragrance, cosmetic, and aromatherapy industries. Despite the commercial importance of these natural aromatic oils, their biosynthesis at the molecular level remains unexplored. As the first step toward understanding the essential oil biosynthesis, we performed de novo transcriptome assembly and analysis of C. flexuosus (lemongrass) by employing Illumina sequencing. Mining of transcriptome data and subsequent phylogenetic analysis led to identification of terpene synthases, pyrophosphatases, alcohol dehydrogenases, aldo-keto reductases, carotenoid cleavage dioxygenases, alcohol acetyltransferases, and aldehyde dehydrogenases, which are potentially involved in essential oil biosynthesis. Comparative essential oil profiling and mRNA expression analysis in three Cymbopogon species (C. flexuosus, aldehyde type; C. martinii, alcohol type; and C. winterianus, intermediate type) with varying essential oil composition indicated the involvement of identified candidate genes in the formation of alcohols, aldehydes, and acetates. Molecular modeling and docking further supported the role of identified protein sequences in aroma formation in Cymbopogon. Also, simple sequence repeats were found in the transcriptome with many linked to terpene pathway genes including the genes potentially involved in aroma biosynthesis. This work provides the first insights into the essential oil biosynthesis of aromatic grasses, and the identified candidate genes and markers can be a great resource for biotechnological and molecular breeding approaches to modulate the essential oil composition.

De Novo Sequencing and Resurrection of a Human Astrovirus-Neutralizing Antibody.

Monoclonal antibody (mAb) therapeutics targeting cancer, autoimmune diseases, inflammatory diseases, and infectious diseases are growing exponentially. Although numerous panels of mAbs targeting infectious disease agents have been developed, their progression into clinically useful mAbs is often hindered by the lack of sequence information and/or loss of hybridoma cells that produce them. Here we combine the power of crystallography and mass spectrometry to determine the amino acid sequence and glycosylation modification of the Fab fragment of a potent human astrovirus-neutralizing mAb. We used this information to engineer a recombinant antibody single-chain variable fragment that has the same specificity as the parent monoclonal antibody to bind to the astrovirus capsid protein. This antibody can now potentially be developed as a therapeutic and diagnostic agent.

A transcriptomic insight into the infective juvenile stage of the insect parasitic nematode, Heterorhabditis indica.

BackgroundNematodes are the most numerous animals in the soil. Insect parasitic nematodes of the genus Heterorhabditis are capable of selectively seeking, infecting and killing their insect-hosts in the soil. The infective juvenile (IJ) stage of the Heterorhabditis nematodes is analogous to Caenorhabditis elegans dauer juvenile stage, which remains in ‘arrested development’ till it finds and infects a new insect-host in the soil. H. indica is the most prevalent species of Heterorhabditis in India. To understand the genes and molecular processes that govern the biology of the IJ stage, and to create a resource to facilitate functional genomics and genetic exploration, we sequenced the transcriptome of H. indica IJs.ResultsThe de-novo sequence assembly using Velvet-Oases pipeline resulted in 13,593 unique transcripts at N50 of 1,371 bp, of which 53 % were annotated by blastx. H. indica transcripts showed higher orthology with parasitic nematodes as compared to free living nematodes. In-silico expression analysis showed 30 % of transcripts expressing with ≥100 FPKM value. All the four canonical dauer formation pathways like cGMP-PKG, insulin, dafachronic acid and TGF-β were active in the IJ stage. Several other signaling pathways were highly represented in the transcriptome. Twenty-four orthologs of C. elegans RNAi pathway effector genes were discovered in H. indica, including nrde-3 that is reported for the first time in any of the parasitic nematodes. An ortholog of C. elegans tol-1 was also identified. Further, 272 kinases belonging to 137 groups, and several previously unidentified members of important gene classes were identified.ConclusionsWe generated high-quality transcriptome sequence data from H. indica IJs for the first time. The transcripts showed high similarity with the parasitic nematodes, M. hapla, and A. suum as opposed to C. elegans, a species to which H. indica is more closely related. The high representation of transcripts from several signaling pathways in the IJs indicates that despite being a developmentally arrested stage; IJs are a hotbed of signaling and are actively interacting with their environment.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-016-2510-z) contains supplementary material, which is available to authorized users.

De Novo Sequencing Assisted Approach for Characterizing Mixture MS/MS Spectra.

Extensive research has been conducted for the computational analysis of mass spectrometry based proteomics data. However, there are still remaining challenges, among which, one particular challenge is the low identification rate of the collected spectral data. A specific contributing factor is the existence of mixture spectra in the collected MS/MS spectra which are generated by the concurrent fragmentation of multiple precursors in one sequencing attempt. The quite frequently observed mixture spectra necessitates the development of effective computational approaches to characterize those non-conventional spectral data. In this research, we proposed an approach for matching the query mixture spectra with a pair of peptide sequences acquired from the protein database by incorporating a special de novo assisted filtration strategy. The experiment results on two different datasets of MS/MS spectra containing mixed ion fragments from multiple peptides demonstrated the efficiency of the integrated filtration strategy in reducing examination space and verified the effectiveness of the proposed matching scheme as well.

A Proteomic Workflow Using High-Throughput De Novo Sequencing Towards Complementation of Genome Information for Improved Comparative Crop Science.

The proteomic study of non-model organisms, such as many crop plants, is challenging due to the lack of comprehensive genome information. Changing environmental conditions require the study and selection of adapted cultivars. Mutations, inherent to cultivars, hamper protein identification and thus considerably complicate the qualitative and quantitative comparison in large-scale systems biology approaches. With this workflow, cultivar-specific mutations are detected from high-throughput comparative MS analyses, by extracting sequence polymorphisms with de novo sequencing. Stringent criteria are suggested to filter for confidential mutations. Subsequently, these polymorphisms complement the initially used database, which is ready to use with any preferred database search algorithm. In our example, we thereby identified 26 specific mutations in two cultivars of Pisum sativum and achieved an increased number (17 %) of peptide spectrum matches.

De Novo Sequencing and Analysis of the Safflower Transcriptome to Discover Putative Genes Associated with Safflor Yellow in Carthamus tinctorius L.

Safflower (Carthamus tinctorius L.), an important traditional Chinese medicine, is cultured widely for its pharmacological effects, but little is known regarding the genes related to the metabolic regulation of the safflower’s yellow pigment. To investigate genes related to safflor yellow biosynthesis, 454 pyrosequencing of flower RNA at different developmental stages was performed, generating large databases.In this study, we analyzed 454 sequencing data from different flowering stages in safflower. In total, 1,151,324 raw reads and 1,140,594 clean reads were produced, which were assembled into 51,591 unigenes with an average length of 679 bp and a maximum length of 5109 bp. Among the unigenes, 40,139 were in the early group, 39,768 were obtained from the full group and 28,316 were detected in both samples. With the threshold of “log2 ratio ≥ 1”, there were 34,464 differentially expressed genes, of which 18,043 were up-regulated and 16,421 were down-regulated in the early flower library. Based on the annotations of the unigenes, 281 pathways were predicted. We selected 12 putative genes and analyzed their expression levels using quantitative real time-PCR. The results were consistent with the 454 sequencing results. In addition, the expression of chalcone synthase, chalcone isomerase and anthocyanidin synthase, which are involved in safflor yellow biosynthesis and safflower yellow pigment (SYP) content, were analyzed in different flowering periods, indicating that their expression levels were related to SYP synthesis. Moreover, to further confirm the results of the 454 pyrosequencing, full-length cDNA of chalcone isomerase (CHI) and anthocyanidin synthase (ANS) were cloned from safflower petal by RACE (Rapid-amplification of cDNA ends) method according to fragment of the transcriptome.

De Novo Sequencing of Peptides from Top-Down Tandem Mass Spectra.

De novo sequencing of proteins and peptides is one of the most important problems in mass spectrometry-driven proteomics. A variety of methods have been developed to accomplish this task from a set of bottom-up tandem (MS/MS) mass spectra. However, a more recently emerged top-down technology, now gaining more and more popularity, opens new perspectives for protein analysis and characterization, implying a need for efficient algorithms to process this kind of MS/MS data. Here, we describe a method that allows for the retrieval, from a set of top-down MS/MS spectra, of long and accurate sequence fragments of the proteins contained in the sample. To this end, we outline a strategy for generating high-quality sequence tags from top-down spectra, and introduce the concept of a T-Bruijn graph by adapting to the case of tags the notion of an A-Bruijn graph widely used in genomics. The output of the proposed approach represents the set of amino acid strings spelled out by optimal paths in the connected components of a T-Bruijn graph. We illustrate its performance on top-down data sets acquired from carbonic anhydrase 2 (CAH2) and the Fab region of alemtuzumab.

Recent Developments in Computational Methods for De Novo Peptide Sequencing from Tandem Mass Spectrometry (MS/MS).

Tandem mass spectrometry (MS/MS) has emerged as a major technology for peptide sequencing. Typically, there are three kinds of methods for the peptide sequencing: database searching, peptide tagging, and de novo sequencing. De novo sequencing has drawn increasing attention because of its independence from existing protein databases and potential for identifying new proteins, proteins resulting from mutations, proteins with unexpected modifications and so on. Recently, with the improvements in the accuracy of MS/MS and development of alternative fragmentation modes of MS/MS, many new de novo sequencing methods have been formulated. This paper reviews these recently developed sequencing methods including those for alternative MS/MS spectra. The paper first introduces background knowledge on peptide sequencing and mass spectrometry, and then reviews de novo peptide sequencing methods for traditional CID spectra. After that, it focuses on the recent development of de novo methods for alternative MS/MS spectra. In addition, methods using multiple spectra from the same peptide are surveyed. Finally, conclusions and some directions of future work are discussed.

CycloBranch: De Novo Sequencing of Nonribosomal Peptides from Accurate Product Ion Mass Spectra.

Nonribosomal peptides have a wide range of biological and medical applications. Their identification by tandem mass spectrometry remains a challenging task. A new open-source de novo peptide identification engine CycloBranch was developed and successfully applied in identification or detailed characterization of 11 linear, cyclic, branched, and branch-cyclic peptides. CycloBranch is based on annotated building block databases the size of which is defined by the user according to ribosomal or nonribosomal peptide origin. The current number of involved nonisobaric and isobaric building blocks is 287 and 521, respectively. Contrary to all other peptide sequencing tools utilizing either peptide libraries or peptide fragment libraries, CycloBranch represents a true de novo sequencing engine developed for accurate mass spectrometric data. It is a stand-alone and cross-platform application with a graphical and user-friendly interface; it supports mzML, mzXML, mgf, txt, and baf file formats and can be run in parallel on multiple threads. It can be downloaded for free from http://ms.biomed.cas.cz/cyclobranch/ , where the User's manual and video tutorials can be found.

ITOC2 – 021. The melanoma immune-peptidome for T-cell-based anti-tumour immunotherapies

Metastatic malignant melanoma (MMM) is a highly aggressive disease. Immune checkpoint inhibitors have proved immunotherapeutic efficacy; however, not all patients profit and severe side effects may occur. Currently, little is known about the tumour-specific targets and the tumour-reactive effector cells that are involved. Aiming to characterise the melanoma immune-peptidome, we purified HLA-I and HLA-II complexes from MMM patient tumours and accurately identified thousands of unique HLA binding peptides per tissue sample, using state-of-the-art mass spectrometry. We applied algorithms for database search in combination with de-novo sequencing and identified many novel cancer-associated peptides from known melanoma-specific antigens (PRAME, Tyrosinase, MAGEA3, PMEL, CTAG1A, PAX3 etc.). We further extended our search to post-translationally modified peptides that potentially serve as neo-antigens. After in vitro T-cell stimulation assays, we detected elevated levels of peptide-specific T-cells that displayed functional responses against peptide-pulsed target cells as well as tumour reactivity. We established this large-scale resource of the melanoma in vivo presented immune-peptidome that contains many novel targets validated by immune reactivity. These data provide the fundament for further optimisation of treatment modalities for this life-threatening disease.

“De-novo” amino acid sequence elucidation of protein G′e by combined “Top-Down” and “Bottom-Up” mass spectrometry

Mass spectrometric de-novo sequencing was applied to review the amino acid sequence of a commercially available recombinant protein G´ with great scientific and economic importance. Substantial deviations to the published amino acid sequence (Uniprot Q54181) were found by the presence of 46 additional amino acids at the N-terminus, including a so-called "His-tag" as well as an N-terminal partial α-N-gluconoylation and α-N-phosphogluconoylation, respectively. The unexpected amino acid sequence of the commercial protein G' comprised 241 amino acids and resulted in a molecular mass of 25,998.9 ± 0.2 Da for the unmodified protein. Due to the higher mass that is caused by its extended amino acid sequence compared with the original protein G' (185 amino acids), we named this protein "protein G'e." By means of mass spectrometric peptide mapping, the suggested amino acid sequence, as well as the N-terminal partial α-N-gluconoylations, was confirmed with 100% sequence coverage. After the protein G'e sequence was determined, we were able to determine the expression vector pET-28b from Novagen with the Xho I restriction enzyme cleavage site as the best option that was used for cloning and expressing the recombinant protein G'e in E. coli. A dissociation constant (K(d)) value of 9.4 nM for protein G'e was determined thermophoretically, showing that the N-terminal flanking sequence extension did not cause significant changes in the binding affinity to immunoglobulins.

Isotopic Differentiation Protocol, a Selective Extraction of C- and N-Terminal Ions in ESI-Ms/Ms De Novo Peptide Sequencing

Simple and selective extraction method of C- and N-terminal ions in ESI-MS/MS analyses, named isotopic differentiation protocol (IDP), was developed. It takes advantage of the mass-difference between the isotopologues after methyl and deuterated methyl ester formations at the C-terminal carboxylic acid or acetamide and deuterated acetamide formations at the N-terminal amines. This method distinguishes the N-terminal and C-terminal ions in principle to facilitate the de novo peptide sequencing. An application with non-ribosomal depsipeptide plipastatin revealed that a pair of the peptide and the corresponding methyl ester at the C-terminus was also applicable to IDP method giving spectrum with enough quality. Facilitating software named “MSMS Search Tool” was also developed.

Choice of next-generation sequencing pipelines.

The next-generation sequencing (NGS) technologies are revolutionary tools which have made possible achieving remarkable advances in genetics since the beginning of the twenty-first century. Thanks to the possibility to produce large amount of sequence data, these tools are going to completely substitute other high-throughput technologies. Moreover, the large applications of NGS protocols are increasing the genetic decoding of biological systems through studies of genome anatomy and gene mapping, coupled to the transcriptome pictures. The application of NGS pipelines such as (1) de-novo genomic sequencing by mate-paired and whole-genome shotgun strategies; (2) specific gene sequencing on large bacterial communities; and (3) RNA-seq methods including whole transcriptome sequencing and Serial Analysis of Gene Expression (Sage-analysis) are fundamental in the genome-wide fields like metagenomics. Recently, the availability of these advanced protocols has allowed to overcome the usual sequencing technical issues related to the mapping specificity over standard shotgun library sequencing, the detection of large structural genomes variations and bridging sequencing gaps, as well as more precise gene annotation. In this chapter we will discuss how to manage a successful NGS pipeline from the planning of sequencing projects through the choice of the platforms up to the data analysis management.

NeatFreq: reference-free data reduction and coverage normalization for De Novo sequence assembly.

BackgroundDeep shotgun sequencing on next generation sequencing (NGS) platforms has contributed significant amounts of data to enrich our understanding of genomes, transcriptomes, amplified single-cell genomes, and metagenomes. However, deep coverage variations in short-read data sets and high sequencing error rates of modern sequencers present new computational challenges in data interpretation, including mapping and de novo assembly. New lab techniques such as multiple displacement amplification (MDA) of single cells and sequence independent single primer amplification (SISPA) allow for sequencing of organisms that cannot be cultured, but generate highly variable coverage due to amplification biases.ResultsHere we introduce NeatFreq, a software tool that reduces a data set to more uniform coverage by clustering and selecting from reads binned by their median kmer frequency (RMKF) and uniqueness. Previous algorithms normalize read coverage based on RMKF, but do not include methods for the preferred selection of (1) extremely low coverage regions produced by extremely variable sequencing of random-primed products and (2) 2-sided paired-end sequences. The algorithm increases the incorporation of the most unique, lowest coverage, segments of a genome using an error-corrected data set. NeatFreq was applied to bacterial, viral plaque, and single-cell sequencing data. The algorithm showed an increase in the rate at which the most unique reads in a genome were included in the assembled consensus while also reducing the count of duplicative and erroneous contigs (strings of high confidence overlaps) in the deliverable consensus. The results obtained from conventional Overlap-Layout-Consensus (OLC) were compared to simulated multi-de Bruijn graph assembly alternatives trained for variable coverage input using sequence before and after normalization of coverage. Coverage reduction was shown to increase processing speed and reduce memory requirements when using conventional bacterial assembly algorithms.ConclusionsThe normalization of deep coverage spikes, which would otherwise inhibit consensus resolution, enables High Throughput Sequencing (HTS) assembly projects to consistently run to completion with existing assembly software. The NeatFreq software package is free, open source and available at https://github.com/bioh4x/NeatFreq.Electronic supplementary materialThe online version of this article (doi:10.1186/s12859-014-0357-3) contains supplementary material, which is available to authorized users.

Characterization and expression of sweetfish (Pleco glossus altivelis) cathepsin D.

Cathepsin D (CTSD) is a lysosomal acidic endoproteinase that plays an important role in immune response. In this study, we obtained sweetfish (Plecoglossus altivelis) CTSD (PaCTSD) via de-novo transcriptome sequencing of sweetfish macrophages. The full length cDNA sequence of PaCTSD was 1 955 bp encoding a propeptide of 397 amino acids. The deduced protein had a calculated molecular weight of 43.17 × 10³. Multiple alignment with other known CTSD amino acid sequences revealed amino acid conservation through the teleosts. Phylogenetic tree analysis showed that PaCTSD grouped tightly with other fish CTSD, and was close to that of Atlantic salmon and rainbow trout. Subsequently, PaCTSD was prokaryotically expressed and refolded by the urea gradient method on a nickel-nitrilotriacetic acid column. Enzyme activity analysis showed that PaCTSD exhibited pH-dependent proteolytic activity. Quantitative real-time PCR showed that PaCTSD mRNA was expressed in all detected tissues in healthy sweetfish. The highest expression was observed in the spleen and white blood cells, followed by liver, head-kidney, kidney, intestine, gill, and muscle. After Listonella anguillarum infection, PaCTSD transcripts were up-regulated significantly in liver, spleen, white blood cells, and head-kidney of sweetfish. In summary, PaCTSD has proteolytic activity and is closely involved in the immune response of sweetfish.