Molecular evolution of transferrin: evidence for positive selection in salmonids

de B. Pereira H, M. da Rocha Filho T, A. S. Machado B, A. Moret M. New distance measure for comparing protein using cellular automata image

One of the first steps in protein sequence analysis is comparing sequences to look for similarities. We propose an information theoretical distance to compare cellular automata representing protein sequences, and determine similarities. Our approach relies in a stationary Hamming distance for the evolution of the automata according to a properly chosen rule, and to build a pairwise similarity matrix and determine common ancestors among different species in a simpler and less computationally demanding computer codes when compared to other methods.

Functional Divergence of the N-Lobe and C-Lobe of Transferrin Gene in Pungitius sinensis (Amur Stickleback)

Transferrin is an important iron-binding glycosylated protein and plays key roles in iron-binding and immune response. Here, a 2037-bp open reading frame was obtained from our previous transcriptome sequencing data of Amur stickleback, which encoded a 679 amino acid putative transferrin protein harbored obvious N-lobe and C-lobe domains. The tissue-specific expression pattern showed that the transcript was detected in a variety of tissues, with the highest signal in liver. Moreover, Streptococcus iniae pathogen stimulation can increase the expression level of this transcript, implying important immune properties for organisms. Next, N-lobes and C-lobes were obtained from 45 fish species. The phylogenetic tree showed that N-lobes and C-lobes were in two different evolutionary branches, and they had different motif composition. Functional divergence indicated a higher evolutionary rate or site-specific alteration among the N-lobe and C-lobe groups. Ka/Ks value of C-lobe group was relatively higher than that of N-lobe group, indicating a faster change rate of C-lobe sequences in evolution. Moreover, some sites experiencing positive selection were also found, which may be involved in the iron- or anion-binding, pathogen resistance and diversification of transferrin protein. Differential iron-binding activity was also detected between N-lobe and C-lobe of Amur stickleback transferrin protein with Chrome Azurol S assay. Compared with the C-lobe, the N-lobe showed stronger growth inhibitory activity of Escherichia coli, implying their potential antibacterial properties. This study will give a reference for subsequent research of transferrin proteins.

Host–Pathogen Interactions of Marine Gram-Positive Bacteria

Simple Summary

Complex interactions between marine Gram-positive pathogens and fish hosts in the marine environment can result in diseases of economically important finfish, which cause economic losses in the aquaculture industry. Understanding how these pathogens interact with the fish host and generate disease will contribute to efficient prophylactic measures and treatments. To our knowledge, there are no systematic reviews on marine Gram-positive pathogens. Therefore, here we reviewed the host–pathogen interactions of marine Gram-positive pathogens from the pathogen-centric and host-centric points of view.

Abstract

Marine Gram-positive bacterial pathogens, including Renibacterium salmoninarum, Mycobacterium marinum, Nocardia seriolae, Lactococcus garvieae, and Streptococcus spp. cause economic losses in marine fish aquaculture worldwide. Comprehensive information on these pathogens and their dynamic interactions with their respective fish–host systems are critical to developing effective prophylactic measures and treatments. While much is known about bacterial virulence and fish immune response, it is necessary to synthesize the knowledge in terms of host–pathogen interactions as a centerpiece to establish a crucial connection between the intricate details of marine Gram-positive pathogens and their fish hosts. Therefore, this review provides a holistic view and discusses the different stages of the host–pathogen interactions of marine Gram-positive pathogens. Gram-positive pathogens can invade fish tissues, evade the fish defenses, proliferate in the host system, and modulate the fish immune response. Marine Gram-positive pathogens have a unique set of virulence factors that facilitate adhesion (e.g., adhesins, hemagglutination activity, sortase, and capsules), invasion (e.g., toxins, hemolysins/cytolysins, the type VII secretion system, and immune-suppressive proteins), evasion (e.g., free radical quenching, actin-based motility, and the inhibition of phagolysosomal fusion), and proliferation and survival (e.g., heme utilization and siderophore-mediated iron acquisition systems) in the fish host. After infection, the fish host initiates specific innate and adaptive immune responses according to the extracellular or intracellular mechanism of infection. Although efforts have continued to be made in understanding the complex interplay at the host–pathogen interface, integrated omics-based investigations targeting host–pathogen–marine environment interactions hold promise for future research.

The widespread presence of a family of fish virulence plasmids in Vibrio vulnificus stresses its relevance as a zoonotic pathogen linked to fish farms

Vibrio vulnificus is a pathogen of public health concern that causes either primary septicemia after ingestion of raw shellfish or secondary septicemia after wound exposure to seawater. In consequence, shellfish and seawater are considered its main reservoirs. However, there is one aspect of its biology that is systematically overlooked: its association with fish in its natural environment. This association led in 1975 to the emergence of a zoonotic clade within phylogenetic lineage 2 following successive outbreaks of vibriosis in farmed eels. Although this clade is now worldwide distributed, no new zoonotic clades were subsequently reported. In this work, we have performed phylogenetic, genomic and functional studies to show that other zoonotic clades are in fact present in 4 of the 5 lineages of the species. Further, we associate these clades, most of them previously but incompletely described, with the acquisition of a family of fish virulence plasmids containing genes essential for resistance to the immune system of certain teleosts of interest in aquaculture. Consequently, our results provide several pieces of evidence about the importance of this species as a zoonotic agent linked to fish farms, as well as on the relevance of these artificial environments acting as drivers that accelerate the evolution of the species.

FEGS: a novel feature extraction model for protein sequences and its applications

BACKGROUND

Feature extraction of protein sequences is widely used in various research areas related to protein analysis, such as protein similarity analysis and prediction of protein functions or interactions.

RESULTS

In this study, we introduce FEGS (Feature Extraction based on Graphical and Statistical features), a novel feature extraction model of protein sequences, by developing a new technique for graphical representation of protein sequences based on the physicochemical properties of amino acids and effectively employing the statistical features of protein sequences. By fusing the graphical and statistical features, FEGS transforms a protein sequence into a 578-dimensional numerical vector. When FEGS is applied to phylogenetic analysis on five protein sequence data sets, its performance is notably better than all of the other compared methods.

CONCLUSION

The FEGS method is carefully designed, which is practically powerful for extracting features of protein sequences. The current version of FEGS is developed to be user-friendly and is expected to play a crucial role in the related studies of protein sequence analyses.

Mapping sequence to feature vector using numerical representation of codons targeted to amino acids for alignment-free sequence analysis

The phylogenetic analysis based on sequence similarity targeted to real biological taxa is one of the major challenging tasks. In this paper, we propose a novel alignment-free method, CoFASA (Codon Feature based Amino acid Sequence Analyser), for similarity analysis of nucleotide sequences. At first, we assign numerical weights to the four nucleotides. We then calculate a score of each codon based on the numerical value of the constituent nucleotides, termed as degree of codons. Accordingly, we obtain the degree of each amino acid based on the degree of codons targeted towards a specific amino acid. Utilizing the degree of twenty amino acids and their relative abundance within a given sequence, we generate 20-dimensional features for every coding DNA sequence or protein sequence. We use the features for performing phylogenetic analysis of the set of candidate sequences. We use multiple protein sequences derived from Beta-globin (BG), NADH dehydrogenase subunit 5 (ND5), Transferrins (TFs), Xylanases, low identity (<40%) and high identity (⩾40%) protein sequences (encompassing 533 and 1064 protein families) for experimental assessments. We compare our results with sixteen (16) well-known methods, including both alignment-based and alignment-free methods. Various assessment indices are used, such as the Pearson correlation coefficient, RF (Robinson-Foulds) distance and ROC score for performance analysis. While comparing the performance of CoFASA with alignment-based methods (ClustalW, ClustalΩ, MAFFT, and MUSCLE), it shows very similar results. Further, CoFASA shows better performance in comparison to well-known alignment-free methods, including LZW-Kernal, jD2Stat, FFP, spaced, and AFKS-D2s in predicting taxonomic relationship among candidate taxa. Overall, we observe that the features derived by CoFASA are very much useful in isolating the sequences according to their taxonomic labels. While our method is cost-effective, at the same time, produces consistent and satisfactory outcomes.

Transferrin Identification in Sterlet (Acipenser ruthenus) Reproductive System

Transferrins are a superfamily of iron-binding proteins and are recognized as multifunctional proteins. In the present study, transcriptomic and proteomic methods were used to identify transferrins in the reproductive organs and sperm of out-of-spawning and spermiating sterlet (Acipenser ruthenus) males. The results showed that seven transferrin transcripts were identified in the transcriptome of sterlet, and these transcripts were qualified as two different transferrin genes, serotransferrin and melanotransferrin, with several isoforms present for serotransferrin. The relative abundance of serotransferrin isoforms was higher in the kidneys and Wolffian ducts in the spermiating males compared to out-of-spawning males. In addition, transferrin was immunodetected in sterlet seminal plasma, but not in sterlet spermatozoa extract. Mass spectrometry identification of transferrin in seminal plasma but not in spermatozoa corroborates immunodetection. The identification of transferrin in the reproductive organs and seminal plasma of sterlet in this study provides the potential function of transferrin during sturgeon male reproduction.

DCGR: feature extractions from protein sequences based on CGR via remodeling multiple information

BACKGROUND

Protein feature extraction plays an important role in the areas of similarity analysis of protein sequences and prediction of protein structures, functions and interactions. The feature extraction based on graphical representation is one of the most effective and efficient ways. However, most existing methods suffer limitations from their method design.

RESULTS

We introduce DCGR, a novel method for extracting features from protein sequences based on the chaos game representation, which is developed by constructing CGR curves of protein sequences according to physicochemical properties of amino acids, followed by converting the CGR curves into multi-dimensional feature vectors by using the distributions of points in CGR images. Tested on five data sets, DCGR was significantly superior to the state-of-the-art feature extraction methods.

CONCLUSION

The DCGR is practically powerful for extracting effective features from protein sequences, and therefore important in similarity analysis of protein sequences, study of protein-protein interactions and prediction of protein functions. It is freely available at https://sourceforge.net/projects/transcriptomeassembly/files/Feature%20Extraction .

Alignment-free method for DNA sequence clustering using Fuzzy integral similarity

A larger amount of sequence data in private and public databases produced by next-generation sequencing put new challenges due to limitation associated with the alignment-based method for sequence comparison. So, there is a high need for faster sequence analysis algorithms. In this study, we developed an alignment-free algorithm for faster sequence analysis. The novelty of our approach is the inclusion of fuzzy integral with Markov chain for sequence analysis in the alignment-free model. The method estimate the parameters of a Markov chain by considering the frequencies of occurrence of all possible nucleotide pairs from each DNA sequence. These estimated Markov chain parameters were used to calculate similarity among all pairwise combinations of DNA sequences based on a fuzzy integral algorithm. This matrix is used as an input for the neighbor program in the PHYLIP package for phylogenetic tree construction. Our method was tested on eight benchmark datasets and on in-house generated datasets (18 s rDNA sequences from 11 arbuscular mycorrhizal fungi (AMF) and 16 s rDNA sequences of 40 bacterial isolates from plant interior). The results indicate that the fuzzy integral algorithm is an efficient and feasible alignment-free method for sequence analysis on the genomic scale.

Alignment-free similarity analysis for protein sequences based on fuzzy integral

Sequence comparison is an essential part of modern molecular biology research. In this study, we estimated the parameters of Markov chain by considering the frequencies of occurrence of the all possible amino acid pairs from each alignment-free protein sequence. These estimated Markov chain parameters were used to calculate similarity between two protein sequences based on a fuzzy integral algorithm. For validation, our result was compared with both alignment-based (ClustalW) and alignment-free methods on six benchmark datasets. The results indicate that our developed algorithm has a better clustering performance for protein sequence comparison.

Evolutionary conservation of transferrin genomic organization and expression characterization in seven freshwater turtles

Serum transferrin (tf), encoding an iron-binding glycoprotein, has been revealed to play important roles in iron transportation and immune response, and it also has been demonstrated to be valuable for phylogenetic analysis in vertebrates. However, the evolutionary conservation, expression profiles and positive selection of transferrin genes among freshwater turtle species remain largely unclear. Here, the genomic DNA and coding sequences of transferrin genes were cloned and characterized in seven freshwater turtles including Mauremys mutica, Mauremys sinensis, Cyclemys dentate, Mauremyssi reevesi, Heosemys grandis, Trachemys scripta and Chrysemys picta. The isolated coding sequences of turtles' tf genes were 2118 bp or 2121 bp, encoding 706 or 707 amino acids. The predicted Tf proteins of turtles share high identities with M. mutica Tf, up to 91%-98% and the M. mutica Tf has the highest identity (91%) in amino acid with the Chelomia mydas Tf, the moderate with other reptiles' Tfs (65%-59%), chicken (58%), and Human Tf (∼55%), and the lowest with zebrafish Tf (41%). Additionally, tf genes were consistently composed of 17 exons and 16 introns with the same splicing sites in introns in all the turtles examined. Moreover, 12 positive selected sites were detected in these turtles' Tf and mainly distributed on the surface of transferrin protein. Importantly, it was found that transferrin genes in all turtles examined were predominantly expressed in adult liver via real-time quantitative PCR. The molecular characterizations and expression profiles of transferrin would shed new insights into understanding the conversations and divergences of transferrin genes in turtles, even in vertebrates.

Identifying anticancer peptides by using a generalized chaos game representation

We generalize chaos game representation (CGR) to higher dimensional spaces while maintaining its bijection, keeping such method sufficiently representative and mathematically rigorous compare to previous attempts. We first state and prove the asymptotic property of CGR and our generalized chaos game representation (GCGR) method. The prediction follows that the dissimilarity of sequences which possess identical subsequences but distinct positions would be lowered exponentially by the length of the identical subsequence; this effect was taking place unbeknownst to researchers. By shining a spotlight on it now, we show the effect fundamentally supports (G)CGR as a similarity measure or feature extraction technique. We develop two feature extraction techniques: GCGR-Centroid and GCGR-Variance. We use the GCGR-Centroid to analyze the similarity between protein sequences by using the datasets 9 ND5, 24 TF and 50 beta-globin proteins. We obtain consistent results compared with previous studies which proves the significance thereof. Finally, by utilizing support vector machines, we train the anticancer peptide prediction model by using both GCGR-Centroid and GCGR-Variance, and achieve a significantly higher prediction performance by employing the 3 well-studied anticancer peptide datasets.

A novel model for protein sequence similarity analysis based on spectral radius

Advances in sequencing technologies led to rapid increase in the number and diversity of biological sequences, which facilitated development in the sequence research. In this paper, we present a new method for analyzing protein sequence similarity. We calculated the spectral radii of 20 amino acids (AAs) and put forward a novel 2-D graphical representation of protein sequences. To characterize protein sequences numerically, three groups of features were extracted and related to statistical, dynamics measurements and fluctuation complexity of the sequences. With the obtained feature vector, two models utilizing Gaussian Kernel similarity and Cosine similarity were built to measure the similarity between sequences. We applied our method to analyze the similarities/dissimilarities of four data sets. Both proposed models received consistent results with improvements when compared to that obtained by the ClustalW analysis. The novel approach we present in this study may therefore benefit protein research in medical and scientific fields.

Hepatic expression of metal-related genes and gill histology in sand flathead (Platycephalus bassensis) from a metal contaminated estuary

Hepatic gene expression and gill histology were measured in sand flathead (Platycephalus bassensis) from a metal polluted estuary. The gene expression analyses were conducted on fish from two most polluted sites and a reference site. The metal-related genes were metal-regulatory transcription factor 1 (MTF1), transferrin (TF), ferriportin1 (FPN1), ferritin and metallothionein. The transcripts of MTF1, TF, and FPN1 were significantly higher in the liver of fish caught at polluted sites, suggesting these genes are potential biomarkers for environmental exposure to metal. Strong correlations were found between the transcripts of these three genes. Four types of gill lesions such as hyperplasia and lamellar fusion, epitheliocystis, telangiectasis, and deformed filament were observed in sampled fish. There was significant difference in the prevalence of epitheliocystis and telangiectasis between the fish from the polluted areas and reference area. Gill parasites were less prevalent in the flathead from polluted sites. The gill histopathological results indicated both pollutants and infections could contribute to gill lesions.

Protein Sequence Comparison Based on Physicochemical Properties and the Position-Feature Energy Matrix

We develop a novel position-feature-based model for protein sequences by employing physicochemical properties of 20 amino acids and the measure of graph energy. The method puts the emphasis on sequence order information and describes local dynamic distributions of sequences, from which one can get a characteristic B-vector. Afterwards, we apply the relative entropy to the sequences representing B-vectors to measure their similarity/dissimilarity. The numerical results obtained in this study show that the proposed methods leads to meaningful results compared with competitors such as Clustal W.

Protein sequence analysis by incorporating modified chaos game and physicochemical properties into Chou's general pseudo amino acid composition

In this contribution we introduced a novel graphical method to compare protein sequences. By mapping a protein sequence into 3D space based on codons and physicochemical properties of 20 amino acids, we are able to get a unique P-vector from the 3D curve. This approach is consistent with wobble theory of amino acids. We compute the distance between sequences by their P-vectors to measure similarities/dissimilarities among protein sequences. Finally, we use our method to analyze four datasets and get better results compared with previous approaches.

An efficient method for measuring the similarity of protein sequences

An accurate numerical descriptor for protein sequence is introduced. It is basically a set of each three successive amino acids in the sequence (triplet), starting from left to right, in addition to the distances between each two successive amino acids in the triplet such that the summation of these distances does not exceed 8. This numerical descriptor combines two features the amino acid composition and the position of each amino acid relative to the other nearby amino acids. This numerical descriptor is used to measure the similarity between protein sequences in three sets: NADH dehydrogenase subunit 5 (ND5) proteins of different species, 24 transferrin proteins from vertebrates and 12 proteins of baculoviruses. High correlation coefficient values between our results and the results of ClustalW program are obtained. These values are higher than the values obtained in many other related works.

An efficient numerical method for protein sequences similarity analysis based on a new two-dimensional graphical representation

A new two-dimensional graphical representation of protein sequences is introduced. Twenty concentric evenly spaced circles divided by n radial lines into equal divisions are selected to represent any protein sequence of length n. Each circle represents one of the different 20 amino acids, and each radial line represents a single amino acid of the protein sequence. An efficient numerical method based on the graph is proposed to measure the similarity between two protein sequences. To prove the accuracy of our approach, the method is applied to NADH dehydrogenase subunit 5 (ND5) proteins of nine different species and 24 transferrin sequences from vertebrates. High values of correlation coefficient between our results and the results of ClustalW are obtained (approximately perfect correlations). These values are higher than the values obtained in many other related works.

cDNA cloning and expression characterization of serum transferrin gene from oriental weatherfish Misgurnus anguillicaudatus

In this study, the cDNA coding serum transferrin (stf) of Misgurnus anguillicaudatus (mastf) was cloned. mastf cDNA is composed of 2326 bp with a 2007 bp open reading frame encoding 668 amino acids. The deduced Mastf protein consists of a signal peptide, two lobes (N and C-lobes) and signature motifs of transferrin (Tf) family. The results of tissue distribution indicated that mastf mRNA was predominantly expressed in the liver. The results indicate that the mastf expression increased significantly in liver, blood, spleen and head kidney after the challenge with Aeromonas sobria, acting as a positive acute protein, suggesting that mastf is related to the immune response. The cloning and expression analysis of mastf further demonstrates the evolutionary conservation of Stf and immune function in vertebrates.

Normalized feature vectors: a novel alignment-free sequence comparison method based on the numbers of adjacent amino acids

Based on all kinds of adjacent amino acids (AAA), we map each protein primary sequence into a 400 by ((L-1)) matrix (M). In addition, we further derive a normalized 400-tuple mathematical descriptors (D), which is extracted from the primary protein sequences via singular values decomposition (SVD) of the matrix. The obtained 400-D normalized feature vectors (NFVs) further facilitate our quantitative analysis of protein sequences. Using the normalized representation of the primary protein sequences, we analyze the similarity for different sequences upon two data sets: 1) ND5 sequences from nine species and 2) transferrin sequences of 24 vertebrates. We also compared the results in this study with those from other related works. These two experiments illustrate that our proposed NFV-AAA approach does perform well in the field of similarity analysis of sequence.

Phylogenetic analyses uncover a novel clade of transferrin in nonmammalian vertebrates

Transferrin is a protein super-family involved in iron transport, a central process in cellular homeostasis. Throughout the evolution of vertebrates, transferrin members have diversified into distinct subfamilies including serotransferrin, ovotransferrin, lactoferrin, melanotransferrin, the inhibitor of carbonic anhydrase, pacifastin, and the major yolk protein in sea urchin. Previous phylogenetic analyses have established the branching order of the diverse transferrin subfamilies but were mostly focused on the transferrin repertoire present in mammals. Here, we conduct a comprehensive phylogenetic analysis of transferrin protein sequences in sequenced vertebrates, placing a special focus on the less-studied nonmammalian vertebrates. Our analyses uncover a novel transferrin clade present across fish, sauropsid, and amphibian genomes but strikingly absent from mammals. Our reconstructed scenario implies that this novel class emerged through a duplication event at the vertebrate ancestor, and that it was subsequently lost in the lineage leading to mammals. We detect footprints of accelerated evolution following the duplication event, which suggest positive selection and early functional divergence of this novel clade. Interestingly, the loss of this novel class of transferrin in mammals coincided with the divergence by duplication of lactoferrin and serotransferrin in this lineage. Altogether, our results provide novel insights on the evolution of iron-binding proteins in the various vertebrate groups.

Miiuy croaker transferrin gene and evidence for positive selection events reveal different evolutionary patterns

Transferrin (TF) is a protein that plays a central role in iron metabolism. This protein is associated with the innate immune system, which is responsible for disease defense responses after bacterial infection. The clear link between TF and the immune defense mechanism has led researchers to consider TF as a candidate gene for disease resistance. In this study, the Miichthys miiuy (miiuy croaker) TF gene (MIMI-TF) was cloned and characterized. The gene structure consisted of a coding region of 2070 nucleotides divided into 17 exons, as well as a non-coding region that included 16 introns and spans 6757 nucleotides. The deduced MIMI-TF protein consisted of 689 amino acids that comprised a signal peptide and two lobes (N- and C-lobes). MIMI-TF expression was significantly up-regulated after infection with Vibrio anguillarum. A series of model tests implemented in the CODEML program showed that TF underwent a complex evolutionary process. Branch-site models revealed that vertebrate TF was vastly different from that of invertebrates, and that the TF of the ancestors of aquatic and terrestrial organisms underwent different selection pressures. The site models detected 10 positively selected sites in extant TF genes. One site was located in the cleft between the N1 and N2 domains and was expected to affect the capability of TF to bind to or release iron indirectly. In addition, eight sites were found near the TF exterior. Two of these sites, which could have evolved from the competition for iron between pathogenic bacteria and TF, were located in potential pathogen-binding domains. Our results could be used to further investigate the function of TF and the selective mechanisms involved.

A relative Lempel-Ziv complexity: Application to comparing biological sequences

One of the main tasks in biological sequence analysis is biological sequence comparison. Numerous efficient methods have been developed for sequence comparison. Traditional sequence comparison is based on sequence alignment. In this report, we propose a novel alignment-free method based on the relative Lempel-Ziv complexity to compare biological sequences. The vertebrate transferring genomes and the spike protein sequences are prepared and tested to evaluate the validity of the method. We use this method to build phylogenetic tree of two groups of the sequences. The result demonstrates that our method is powerful and efficient.

Antimicrobial mechanisms of fish leukocytes

Early activation and coordination of innate defenses are critical for effective responses against infiltrating pathogens. Rapid engagement of immune cells provides a critical first line of defense soon after pathogen infiltration. Activation leads to a well-orchestrated set of events that sees the induction and regulation of intracellular and extracellular antimicrobial defenses. An array of regulatory mediators, highly toxic soluble molecules, degradative enzymes and antimicrobial peptides provides maximal protection against a wide range of pathogens while limiting endogenous damage to host tissues. In this review we highlight recent advances in our understanding of innate cellular antimicrobial responses of teleost fish and discuss their implications to cell survival, immunomodulation and death. The evolutionary conservation of these responses is a testament to their effectiveness against pathogen infiltration and their commitment to effective maintenance of host homeostasis. Importantly, recent developments in teleost fish systems have identified novel host defense strategies that may be unique to this lower vertebrate group or may point to previously unknown innate mechanisms that also play a significant role in higher vertebrate host immunity.

Characterization of the gilthead seabream (Sparus aurata L.) transferrin gene: genomic structure, constitutive expression and SNP variation

Transferrin (Tf) is a multi-function protein with a central role in iron metabolism, and it is this function that is associated with a role in the innate immune system response. The clear link between Tf and immune defense mechanism lead to propose Tf as a candidate gene for disease resistance. In this study, genomic and cDNA sequences of Tf gene in gilthead seabream (Sparus aurata L.) (SaTf gene), were identified and characterized. SaTf gene structure consists of a coding region of 2076 nucleotides divided into 17 exons and a no-coding region that includes 16 introns and spans 5495 nucleotides. The deduced Tf protein for gilthead seabream is composed of 691 amino acids and consists of an initial peptide and two lobes (N- and C-lobes). This gene structure is similar to that of previously described Tf genes in other fish species. RT-PCR analyses carried out in different tissues and two developmental stages showed tissue-and stage-specific Tf expression in gilthead seabream. Finally, by sequencing the transferrin genomic sequences of 20 unrelated seabreams, 31 SNPs were identified. These data allowed the estimation of the frequency of nucleotide substitution in the SaTf gene as 1SNP per 253 bp. SNPs were detected in different regions of the genomic sequence but they were mainly localized in non-coding regions, specifically, SNP frequency in non-coding regions was fifteen-fold higher than within coding regions.

Trans genomic capture and sequencing of primate exomes reveals new targets of positive selection

Comparison of protein-coding DNA sequences from diverse primates can provide insight into these species' evolutionary history and uncover the molecular basis for their phenotypic differences. Currently, the number of available primate reference genomes limits these genome-wide comparisons. Here we use targeted capture methods designed for human to sequence the protein-coding regions, or exomes, of four non-human primate species (three Old World monkeys and one New World monkey). Despite average sequence divergence of up to 4% from the human sequence probes, we are able to capture ~96% of coding sequences. Using a combination of mapping and assembly techniques, we generated high-quality full-length coding sequences for each species. Both the number of nucleotide differences and the distribution of insertion and deletion (indel) lengths indicate that the quality of the assembled sequences is very high and exceeds that of most reference genomes. Using this expanded set of primate coding sequences, we performed a genome-wide scan for genes experiencing positive selection and identified a novel class of adaptively evolving genes involved in the conversion of epithelial cells in skin, hair, and nails to keratin. Interestingly, the genes we identify under positive selection also exhibit significantly increased allele frequency differences among human populations, suggesting that they play a role in both recent and long-term adaptation. We also identify several genes that have been lost on specific primate lineages, which illustrate the broad utility of this data set for other evolutionary analyses. These results demonstrate the power of second-generation sequencing in comparative genomics and greatly expand the repertoire of available primate coding sequences.

Detecting signatures of positive selection in partial sequences generated on a large scale: pitfalls, procedures and resources

Studying the actions of selection provides insight into adaptation, population divergence and gene function. Next-generation sequencing produces large amounts of partial sequences, potentially facilitating efforts to detect signatures of selection based on comparisons between synonymous (d(S)) and nonsynonymous (d(N)) substitutions, and single nucleotide polymorphism assays placed in selected genes would improve the ability to study adaptation in population surveys. However, sequences generated by these technologies are typically short. In nonmodel organisms that are a focus of evolutionary studies, the lack of a reference genome that facilitates the assembly of short sequences has limited surveys of positive selection in large numbers of genes. Here, we describe a series of steps to facilitate these surveys. We provide PERL scripts to assist data analysis, and describe the use of commonly available programs. We demonstrate these approaches in six salmon species, which have partially duplicated genomes. We recommend using multiway blast to optimize the number of alignments between partial coding sequences. Reading frames should be manually detected after alignment with sequences in Genbank using the BLASTX program. We encourage the use of a phylogenetic approach to separate orthologs from paralogs in duplicated genomes. Simple simulations on a gene known to have undergone selection in salmon species, transferrin, showed that the ability to detect selection in short sequences (<600 bp) depended on the proportion of codons under selection (1-2%) within that sequence. This relationship was less relevant in longer sequences. In this exploratory study, we detected 11 genes showing evidence of positive selection.

Phylogenetic analysis of protein sequences based on distribution of length about common sub-string

Up to now, various approaches for phylogenetic analysis have been developed. Almost all of them put stress on analyzing nucleic acid sequences or protein primary sequences. In this paper, we propose a new sequence distance for efficient reconstruction of phylogenetic trees based on the distribution of length about common sub-sequences between two sequences. We describe some applications of this method, which not only show the validity of the method, but also suggest a number of novel phylogenetic insights.

Patterns of molecular evolution and predicted function in thaumatin-like proteins of Populus trichocarpa

Some pathogenesis-related proteins (PR proteins) are subject to positive selection, while others are under negative selection. Here, we report the patterns of molecular evolution in thaumatin-like protein (TLP, PR5 protein) genes of Populus trichocarpa. Signs of positive selection were found in 20 out of 55 Populus TLPs using the likelihood ratio test and ML-based Bayesian methods. Due to the connection between the acidic cleft and the antifungal activity, the secondary structure and three-dimensional structure analyses predicted antifungal activity beta-1,3-glucanase activities in these TLPs. Moreover, the coincidence with variable basic sites in the acidic cleft and positively selected sites suggested that fungal diseases may have been the main environmental stress that drove rapid adaptive evolution in Populus.

Positive selection drives lactoferrin evolution in mammals

Lactoferrin (LF) is a member of the transferrin family that is abundantly expressed and secreted by glandular epithelial cells. The biological functions of LF involve in iron homeostasis regulation of the body and antibacterial activity. Previous studies demonstrated that it had a high cationic N-terminal domain that could interact with glycosaminoglycans, lipopolysaccharides and the bacterial virulence protein. Two anti-microbial peptides, lactoferricin (LFcin) and lactoferrampin (LFampin), were also isolated and identified in N-terminal of LF. Although the antibacterial mechanism was carefully studied, little was known about the molecular evolution of LF. In this study, we estimated the nonsynonymous-to-synonymous substitution ratios ( ) per site using maximum likelihood method to analyze the LF evolution. The results of omega > 1 and five identified positive selection sites of amino acid suggested that the evolution of LF gene was characterized by positive selection. Further study found that the positive selection sites were either located in the LF-bacteria binding region or the peptides of LFcin and LFampin, indicating that the selection pressure was related to LF-bacteria interaction. The identification of these sites may contribute to the mechanism of bacteria-LF interaction.

Allelic discrimination, three-dimensional analysis and gene expression of multiple transferrin alleles of common carp (Cyprinus carpio L.)

We cloned and sequenced four different transferrin (Tf) alleles (C, D, F and G) of European common carp (Cyprinus carpio carpio L.) and studied allelic diversity with respect to differences in sequence, constitutive transcription and three-dimensional structure. Most of the disulfide bonds were conserved between human and carp Tf, and modeling confirmed the overall conservation of the three-dimensional structure of carp Tf. While the iron-binding sites in the C-lobe of carp Tf were completely conserved, in the N-lobe the majority of iron-coordinating residues were not conserved. This may have a serious impact on the ability of carp Tf to bind iron with both the C- and N-lobe. In contrast to human Tf, we could not detect potential N-glycosylation sites in carp Tf, which does not seem to be a glycoprotein. Comparison of the cDNA of the four Tf alleles of carp indicated 21 polymorphic sites of which 13 resulted in non-synonymous changes. Allelic diversity did not seem to influence the overall conservation of carp Tf. Neither the iron binding sites nor the receptor binding of carp Tf seemed influenced by allelic diversity. Possibly, interaction with pathogen-associated receptors for Tf could be influenced by allelic diversity. Basal gene expression of Tf alleles D and G was especially high in carp liver. Although we could detect a higher transcription level of allele D than of Tf allele G in head kidney, thymus and spleen, the differences seem minor with respect to the very high transcription level in liver. Preliminary results with Tf-typed serum suggest a difference in the ability of Tf alleles D and G to modulate LPS-induced NO production in carp macrophages.

Evolutionary dynamics of transferrin in Notropis

This study reports on the development of comparative data for the transferrin (TF) gene in cyprinid fishes, focusing on the genus Notropis. While previous studies have suggested varied patterns of adaptation influencing the diversity at this gene locus both within and among species, sequence data for (TF) in Notropis exhibit limited evidence for selection. However, there are significant effects of sequence variation associated with Notropis lutipinnis, suggesting some form of diversifying selection acting among populations of this species. Overall, the gene performs well as a second locus for phylogenetic and biogeographic inference and may help improve description of the pattern and the process of diversification in Notropis.

The induction of nitric oxide response of carp macrophages by transferrin is influenced by the allelic diversity of the molecule

The central role of transferrin (Tf) as an iron transporting protein has been extended by observations that modified versions of Tf also participate in the regulation of innate immunity. We report on the isolation of two carp Tf proteins (alleles D and G) to purity using rivanol precipitation and ion-exchange chromatography, and describe the activation of head kidney-derived carp macrophages by cleaved Tf. We demonstrate the superiority of the D-type over the G-type Tf in inducing nitric oxide (NO) and confirm previous observations that full-length Tf cannot induce NO in fish macrophages. We believe that cleaved Tf fragments should be considered to be "alarmins". We discuss the possibility that parasites such as Trypanoplasma borreli cleave Tf and use Tf fragments to their advantage by modulating the NO induction in carp macrophages.

Expression analysis of selected immune-relevant genes in channel catfish during Edwardsiella ictaluri infection

Channel catfish Ictalurus punctatus were intraperitoneally challenged with the bacterium Edwardsiella ictaluri (the causative agent of enteric septicemia of catfish), and the expression of genes presumed to function in the inducible innate defense was evaluated. End-binding protein 1 (EB1), beta1-integrin, natural-resistance-associated macrophage protein (Nramp), heat shock protein 70 (Hsp70), serum amyloid P (SAP), and transferrin gene expression profiles were determined using quantitative reverse-transcriptase-polymerase chain reaction on liver, anterior kidney, spleen, and gut. Fish were subsampled at 0, 24, 48, 72, and 96 h after bacterial or phosphate-buffered-saline injection. Posterior kidney sampling demonstrated increasing bacterial counts at 24-48 h postinjection (hpi), followed by a plateau to 96 hpi. The transferrin and SAP transcripts were liver specific. The other genes were expressed in all four tissues. In bacterially infected fish, expression of EB1 (anterior kidney, spleen, and liver), Hsp70 (anterior kidney and spleen), and Nramp (spleen and gut) significantly increased by 48 hpi. Transferrin was strongly up-regulated and SAP was downregulated by 72 hpi, indicating positive and negative acute-phase reactants, respectively. The data indicate a substantial response of innate immunity effector cells by 48 hpi, followed by suppression of bacterial growth and induction of the acute-phase response. This suggests that the 48-72-hpi time frame is critical in our model for evaluating the effectiveness of innate defenses.

Evidence for two transferrin loci in the Salmo trutta genome

To determine the organization of transferrin (TF) locus in the Salmo trutta genome, partial DNA and cDNA sequencing, fluorescent in situ hybridization (FISH) and Salmo salar BAC analysis were performed. TF expression levels and copy number prediction were assessed using real-time PCR. In addition to two previously reported DNA TF variant sequences of S. trutta and Salmo marmoratus (TF1), two novel variant sequences (TF2) were revealed in both species. Variant-specific sequence tags, characterizing two variants for each TF type (TF1 and TF2), were identified in genomic clones from each of the F1 hybrids between S. trutta and S. marmoratus. These clearly documented double heterozygote status at the TF loci. The real-time PCR data showed that each of the two TF types (TF1 and TF2) existed in one copy only and that the transcription of TF2 was considerably lower compared with TF1. Using FISH, hybridization signals were observed on two medium-sized acrocentric chromosomes of S. trutta karyotype. A TF type-specific PCR followed by a restriction analysis revealed the presence of two TF loci in the majority of analysed BAC clones. It was concluded that the TF gene is duplicated in the genome of S. trutta, and that the two TF loci are located adjacent to one another on the same chromosome. The differing transcription levels of TF1 and TF2 appear to depend on the corresponding promoter activity, which at least for TF2 seems to vary between different Salmo congeners.

Evidence of positive selection on the Atlantic salmon CD3γδ gene

Atlantic salmon are typically anadromous, spending the majority of their lifetime in oceans and returning to fresh water to breed. This diversity of environments likely results in strong selective forces shaping their genome. In this paper, we present the first genomics approach to detect positive selection operating on the Salmo salar (salmon) lineage, an important aquaculture species. We identify a panel of candidate genes that may have been subject to adaptive evolution in this species. In particular, we identify a robust signature of positive selection operating on the salmon CD3gammadelta gene, which encodes one of the protein chains essential for formation of the T-cell receptor complex and for T-cell activation. Furthermore, we identified the particular codon sites that have been subject to positive selection in fish and highlight two sites flanking an important N-glycosylation site in this molecule.