Strong Codon Usage Bias Research Articles

Comparative mitogenomic analysis provides evolutionary insights into Formica (Hymenoptera: Formicidae).

Formica is a large genus in the family Formicidae with high diversity in its distribution, morphology, and physiology. To better understand evolutionary characteristics of Formica, the complete mitochondrial genomes (mitogenomes) of two Formica species were determined and a comparative mitogenomic analysis for this genus was performed. The two newly sequenced Formica mitogenomes each included 37 typical mitochondrial genes and a large non-coding region (putative control region), as observed in other Formica mitogenomes. Base composition, gene order, codon usage, and tRNA secondary structure were well conserved among Formica species, whereas diversity in sequence size and structural characteristics was observed in control regions. We also observed several conserved motifs in the intergenic spacer regions. These conserved genomic features may be related to mitochondrial function and their highly conserved physiological constraints, while the diversity of the control regions may be associated with adaptive evolution among heterogenous habitats. A negative AT-skew value on the majority chain was presented in each of Formica mitogenomes, indicating a reversal of strand asymmetry in base composition. Strong codon usage bias was observed in Formica mitogenomes, which was predominantly determined by nucleotide composition. All 13 mitochondrial protein-coding genes of Formica species exhibited molecular signatures of purifying selection, as indicated by the ratio of non-synonymous substitutions to synonymous substitutions being less than 1 for each protein-coding gene. Phylogenetic analyses based on mitogenomic data obtained fairly consistent phylogenetic relationships, except for two Formica species that had unstable phylogenetic positions, indicating mitogenomic data are useful for constructing phylogenies of ants. Beyond characterizing two additional Formica mitogenomes, this study also provided some key evolutionary insights into Formica.

Comparative Study on Codon Usage Patterns across Chloroplast Genomes of Eighteen Taraxacum Species

This study investigates codon usage bias within the chloroplast genomes of 18 Taraxacum species, focusing on the base composition and various metrics including GC content, Relative Synonymous Codon Usage (RSCU), Effective Number of Codons (ENc), and GC3s. Our analysis revealed a pronounced preference for A/T-ending codons across Taraxacum species, with GC content across the first, second, and third positions of the codons (GC1, GC2, GC3) and the average GC content consistently below 50%. A detailed examination using the RSCU metric identified 29 commonly preferred A/T-ending codons, indicating a strong codon usage bias towards these endings. Specifically, the codon for leucine (UUA) emerged as highly preferred, while the codon for serine (AGC) was least favored. Through the ENc–GC3s plot analysis, we explored the forces shaping this bias, finding evidence that both mutation pressure and natural selection significantly influence codon preference, with most coding sequences showing weak bias. The PR2 plot analysis further confirmed the role of these factors by demonstrating a higher frequency of T over A and C over G at the third codon position, pointing towards a mutation bias complemented by natural selection. Collectively, our findings highlight a consistent pattern of codon usage bias in the chloroplast genomes of Taraxacum species, influenced by a combination of mutation pressure, natural selection, and possibly other unidentified factors.

Frequency of change determines effectiveness of microbial response strategies

Nature challenges microbes with change at different frequencies and demands an effective response for survival. Here, we used controlled laboratory experiments to investigate the effectiveness of different response strategies, such as post-translational modification, transcriptional regulation, and specialized versus adaptable metabolisms. For this, we inoculated replicated chemostats with an enrichment culture obtained from sulfidic stream microbiomes 16 weeks prior. The chemostats were submitted to alternatingly oxic and anoxic conditions at three frequencies, with periods of 1, 4 and 16 days. The microbial response was recorded with 16S rRNA gene amplicon sequencing, shotgun metagenomics, transcriptomics and proteomics. Metagenomics resolved provisional genomes of all abundant bacterial populations, mainly affiliated with Proteobacteria and Bacteroidetes. Almost all these populations maintained a steady growth rate under both redox conditions at all three frequencies of change. Our results supported three conclusions: (1) Oscillating oxic/anoxic conditions selected for generalistic species, rather than species specializing in only a single condition. (2) A high frequency of change selected for strong codon usage bias. (3) Alignment of transcriptomes and proteomes required multiple generations and was dependent on a low frequency of change.

Analysis of Codon Usage Bias in Xyloglucan Endotransglycosylase (XET) Genes

Xyloglucan endotransglycosylase (XET) genes are widely distributed in most plants, but the codon usage bias of XET genes has remained uncharacterized. Thus, we analyzed the codon usage bias using 4500 codons of 20 XET genes to elucidate the genetic and evolutionary patterns. Phylogenetic and hierarchical cluster analyses revealed that the 20 XET genes belonged to two groups. The closer the genetic distance, the more similar the codon usage preference. The codon usage bias of most XET genes was weak, but there was also some codon usage bias. AGA, AGG, AUC, and GUG were the top four codons (RSCU > 1.5) in the 20 XET genes. CitXET had a stronger codon usage bias, and there were eight optimal codons of CitXET (i.e., AGA, AUU, UCU, CUU, CCA, GCU, GUU, and AAA). The RSCU values underwent a correspondence analysis. The two main factors affecting codon usage bias (i.e., Axes 1 and 2) accounted for 54.8% and 17.6% of the total variation, respectively. Multiple correspondence analysis revealed that XET genes were widely distributed, with Group 1 genes being closer to Axis 1 than Group 2 genes, which were closer to Axis 2. Codons with A/U at the third codon position were distributed closer to Axis 1 than codons with G/C at the third codon position. PgXET, ZmXET, VlXET, VrXET, and PcXET were biased toward codons ending with G/C. In contrast, CitXET, DpXET, and BrpXET were strongly biased toward codons ending with A/U, indicating that these XET genes have a strong codon usage bias. Translational selection and base composition (especially A and U at the third codon position), followed by mutation pressure and natural selection, may be the most important factors affecting codon usage of 20 XET genes. These results may be useful in clarifying the codon usage bias of XET genes and the relevant evolutionary characteristics.

Genome dynamics, codon usage patterns and influencing factors in Aeromonas hydrophila phages

In the present study, genome characteristics and codon usage patterns of 44 Aeromonas hydrophila phages were studied. Phage genomes varied from 30.8 to 262.0 kb with mean±SD and median values of 111.3 ± 81.4 kb and 79.4 kb, respectively. Though the great variation in phage GC contents (35.1–62.2%) was observed, GC contents of all phages (except two phages) were significantly less than the GC content (62.4 ± 5.6%) of the host. The effective number of codons (ENC) values of phage genes ranged from 27.7 to 61 with a mean±SD value of 47.4 ± 6.8. Out of a total 5773 phage genes, 207 (3.6%), 3,528 (61.1%) and 2,012 (34.9%) genes had strong (ENC < 35), moderate (35 < ENC < 50) and low (ENC ≥ 50) codon usage bias, respectively. During relative synonymous codon usage (RSCU) analysis, shared usage of preferred codons was also observed between the phages and host. During codon adaptation index (CAI) analysis, 1028 (17.8%) phage genes showed significant adaptation towards the host. Among these genes, 797 (78.0%) genes encoded hypothetical proteins or proteins of unknown function; whereas 118 (12%) genes encoded the phage structural and packaging proteins. Segregation of ENC, RSCU and CAI analysis results based on genome size also indicated that codon usage bias was more prominent in phages with small genomes. Correlation, neutrality and GC3 versus ENC analyzes indicated a more dominant role of natural selection in shaping the codon usage patterns of A. hydrophila phages. The findings of the current study could be useful from evolutionary and host-pathogen interaction perspectives leading to efficient utilization of phages for therapeutic and other applications.

Codon usage of host-specific P genotypes (VP4) in group A rotavirus

BackgroundGroup A rotavirus (RVA) is a common causative agent of acute gastroenteritis in infants and young children worldwide. RVA P genotypes, determined by VP4 sequences, have been confirmed to infect humans and animals. However, their codon usage patterns that are essential to obtain insights into the viral evolution, host adaptability, and genetic characterization remained unclear, especially across animal hosts.ResultsWe performed a comprehensive codon usage analysis of eight host-specific RVA P genotypes, including human RVA (P[4] and P[8]), porcine RVA (P[13] and P[23]), and zoonotic RVA (P[1], P[6], P[7] and P[19]), based on 233 VP4 complete coding sequences. Nucleotide composition, relative synonymous codon usage (RSCU), and effective number of codons (ENC) were calculated. Principal component analysis (PCA) based on RSCU values was used to explore the codon usage patterns of different RVA P genotypes. In addition, mutation pressure and natural selection were identified by using ENC-plot, parity rule 2 plot, and neutrality plot analyses. All VP4 sequences preferred using A/U nucleotides (A: 0.354-0.377, U: 0.267-0.314) than G/C nucleotides across genotypes. Similarly, majority of commonly used synonymous codons were likely to end with A/U nucleotides (A: 9/18-12/18, U: 6/18-9/18). In PCA, human, porcine, and zoonotic genotypes clustered separately in terms of RSCU values, indicating the host-specific codon usage patterns; however, porcine and zoonotic genotypes were partly overlapped. Human genotypes, P[4] and P[8], had stronger codon usage bias, as indicated by more over-represented codons and lower ENC, compared to porcine and zoonotic genotypes. Moreover, natural selection was determined to be a predominant driver in shaping the codon usage bias across the eight P genotypes. In addition, mutation pressure contributed to the codon usage bias of human genotypes.ConclusionsOur study identified a strong codon usage bias of human RVA P genotypes attributable to both natural selection and mutation pressure, whereas similar codon usage bias between porcine and zoonotic genotypes predominantly attributable to natural selection. It further suggests possible cross-species transmission. Therefore, it warrants further surveillance of RVA P genotypes for early identification of zoonotic infection.

Whole genome sequencing and comparative genomic analyses of Pseudomonas aeruginosa strain isolated from arable soil reveal novel insights into heavy metal resistance and codon biology

Elevated concentration of non-essential persistent heavy metals and metalloids in the soil is detrimental to essential soil microbes and plants, resulting in diminished diversity and biomass. Thus, isolation, screening, and whole genomic analysis of potent strains of bacteria from arable lands with inherent capabilities of heavy metal resistance and plant growth promotion hold the key for bio remedial applications. This study is an attempt to do the same. In this study, a potent strain of Pseudomonas aeruginosa was isolated from paddy fields, followed by metabolic profiling using FTIR, metal uptake analysis employing ICP-MS, whole genome sequencing and comparative codon usage analysis. ICP-MS study provided insights into a high degree of Cd uptake during the exponential phase of growth under cumulative metal stress to Cd, Zn and Co, which was further corroborated by the detection ofcadAgene along with czcCBA operon in the genome upon performing whole-genome sequencing. This potent strain ofPseudomonas aeruginosaalso harboured genes, such ascopA, chrA, znuA, mgtE, corA, and others conferring resistance against different heavy metals, such as Cd, Zn, Co, Cu, Cr, etc. A comparative codon usage bias analysis at the genomic and genic level, whereby several heavy metal resistant genes were considered in the backdrop of two housekeeping genes among 40Pseudomonasspp. indicated the presence of a relatively strong codon usage bias in the studied strain. With this work, an effort was made to explore heavy metal-resistant bacteria (isolated from arable soil) and whole genome sequence analysis to get insight into metal resistance for future bio remedial applications.

Mitochondrial genome characterization of Gryllodes sigillatus (Orthoptera: Gryllidae) and its phylogenetic implications

Gryllodes sigillatus is a cricket widely distributed throughout the world. In this study, we reported the first complete mitogenome sequence of Genus Gryllodes and inferred its phylogeny. The mitogenome of G. sigillatus was 16,369 bp and consisted of a control region and a typical set of 37 genes. It was AT-rich with strong codon usage bias and possessed a gene arrangement of trnE-trnS1-trnN. Phylogenetic analysis indicated G. sigillatus was sister species to Velarifictorus hemelytrus, together belonging to the Family Gryllidae. Our findings would contribute to understanding mitogenomic evolution and phylogeny of Ensifera.

Mitogenome evolution in ladybirds: Potential association with dietary adaptation.

Dietary shifts can alter the relative availability of different nutrients and are therefore associated with metabolic adaptation in animals. The Coccinellidae (ladybirds) exhibits three major types of feeding habits and provides a useful model to study the effects of dietary changes on the evolution of mitogenomes, which encode proteins directly involved in energy metabolism. Here, mitogenomes of three coccinellid species were newly sequenced. These data were combined with other ten previously sequenced coccinellid mitogenomes to explore the relationship between mitogenome evolution and diets. Our results indicate that mitogenomic data can be effectively used to resolve phylogenetic relationships of Coccinellidae. Strong codon usage bias in coccinellid mitogenomes was predominantly determined by nucleotide composition. The 13 mitochondrial protein‐coding genes (PCGs) globally evolved under negative constraints, with some PCGs showing a stronger purifying selection. Six PCGs (nad3, nad4L, and nad5 from Complex I; cox1 and cox3 from Complex IV; and atp6 from Complex V) displayed signs of positive selection. Of these, adaptive changes in cox3 were potentially associated with metabolic differences resulting from dietary shifts in Coccinellidae. Our results provide insights into the adaptive evolution of coccinellid mitogenomes in response to both dietary shifts and other life history traits.

Codon Usage in the Iflaviridae Family Is Not Diverse Though the Family Members Are Isolated from Diverse Host Taxa.

All iflavirus members belong to the unique genus, Iflavirus, of the family, Iflaviridae. The host taxa and sequence identities of these viruses are diverse. A codon usage bias, maintained by a balance between selection, mutation, and genetic drift, exists in a wide variety of organisms. We characterized the codon usage patterns of 44 iflavirus genomes that were isolated from the classes, Insecta, Arachnida, Mammalia, and Malacostraca. Iflaviruses lack a strong codon usage bias when they are evaluated using an effective number of codons. The odds ratios of the majority of dinucleotides are within the normal range. However, the dinucleotides at the 1st–2nd codon positions are more biased than those at the 2nd–3rd codon positions. Plots of effective numbers of codons, relative neutrality analysis, and PR2 bias analysis all indicate that selection pressure dominates mutations in shaping codon usage patterns in the family, Iflaviridae. When these viruses were grouped into their host taxa, we found that the indices, including the nucleotide composition, effective number of codons, relative synonymous codon usage, and the influencing factors behind the codon usage patterns, all show that there are non-significant differences between the six host-taxa-groups. Our results disagree with our assumption that diverse viruses should possess diverse codon usage patterns, suggesting that the nucleotide composition and codon usage in the family, Iflaviridae, are not host taxa-specific signatures.

Nonoptimal codon usage influences protein structure in intrinsically disordered regions.

Synonymous codons are not used with equal frequencies in most genomes. Codon usage has been proposed to play a role in regulating translation kinetics and co-translational protein folding. The relationship between codon usage and protein structures and the in vivo role of codon usage in eukaryotic protein folding is not clear. Here, we show that there is a strong codon usage bias in the filamentous fungus Neurospora. Importantly, we found genome-wide correlations between codon choices and predicted protein secondary structures: Nonoptimal codons are preferentially used in intrinsically disordered regions, and more optimal codons are used in structured domains. The functional importance of such correlations in vivo was confirmed by structure-based codon manipulation of codons in the Neurospora circadian clock gene frequency (frq). The codon optimization of the predicted disordered, but not well-structured regions of FRQ impairs clock function and altered FRQ structures. Furthermore, the correlations between codon usage and protein disorder tendency are conserved in other eukaryotes. Together, these results suggest that codon choices and protein structures co-evolve to ensure proper protein folding in eukaryotic organisms.

Capturing the variant surface glycoprotein repertoire (the VSGnome) of Trypanosoma brucei Lister 427

Trypanosoma brucei evades the adaptive immune response through the expression of antigenically distinct Variant Surface Glycoprotein (VSG) coats. To understand the progression and mechanisms of VSG switching, and to identify the VSGs expressed in populations of trypanosomes, it is desirable to predetermine the available repertoire of VSG genes (the ‘VSGnome’). To date, the catalog of VSG genes present in any strain is far from complete and the majority of current information regarding VSGs is derived from the TREU927 strain that is not commonly used as an experimental model. We have assembled, annotated and analyzed 2563 distinct and previously unsequenced genes encoding complete and partial VSGs of the widely used Lister 427 strain of T. brucei. Around 80% of the VSGnome consists of incomplete genes or pseudogenes. Read-depth analysis demonstrated that most VSGs exist as single copies, but 360 exist as two or more indistinguishable copies. The assembled regions include five functional metacyclic VSG expression sites. One third of minichromosome sub-telomeres contain a VSG (64–67 VSGs on ∼96 minichromosomes), of which 85% appear to be functionally competent. The minichromosomal repertoire is very dynamic, differing among clones of the same strain. Few VSGs are unique along their entire length: frequent recombination events are likely to have shaped (and to continue to shape) the repertoire. In spite of their low sequence conservation and short window of expression, VSGs show evidence of purifying selection, with ∼40% of non-synonymous mutations being removed from the population. VSGs show a strong codon-usage bias that is distinct from that of any other group of trypanosome genes. VSG sequences are generally very divergent between Lister 427 and TREU927 strains of T. brucei, but those that are highly similar are not found in ‘protected’ genomic environments, but may reflect genetic exchange among populations.

Complete mitochondrial genome of the Yellownose skate: Zearaja chilensis (Rajiformes, Rajidae)

The complete sequence of mitochondrial DNA of a Yellownose skate, Zearaja chilensis was determined for the first time. It is 16,909 bp in length covering 2 rRNA, 22 tRNA and 13 protein coding genes with the identical gene order and structure as those of other Rajidae species. The nucleotide of L-strand is composed of low G (14.3%), and slightly high A + T (58.9%) nucleotides. The strong codon usage bias against the use of G (6.0%) is found at the third codon positions. Twelve of the 13 protein coding genes use ATG as the start codon while COX1 starts with GTG. As for the stop codon, only ND4 shows an incomplete stop codon TA. This is the first report of the mitogenome for a species in the genus Zearaja, providing a valuable source of genetic information on the evolution of the family Rajidae and the genus Zearaja as well as for establishment of a sustainble fishery management plan of the species.

Complete mitochondrial genome of the brownbanded bamboo shark Chiloscyllium punctatum

The complete mitochondrial genome of the brownbanded bamboo shark Chiloscyllium punctatum was obtained for the first time in this study. The mitogenome is 16,703 bp in length, containing 13 protein-coding genes, 2 rRNA genes, 22 tRNA genes and one control region. The gene order and direction of transcription of C. punctatum is identical to most vertebrates. The strong codon usage bias against the G was found in the third codon positions of H-strand-encoded protein genes. The tRNA-Ser2 lost dihydrouridine arm and could not be fold into the typical cloverleaf secondary structure. In the control region, the termination-associated sequence domain, the central conserved domain and the conserved sequence blocks domain were identified.

English

Ubiquitin which has the function of selective protein degradation may play an important role in the regulation of insect growth and development. The coding sequence of an ubiquitin gene from the larvae of the rice stem borer, Chilo suppressalis Walker (Lepidoptera: Pyralidae) named CsUB (GenBank Accession No. GU238420) was cloned by RT-PCR and sequenced in this study, with primers according to the sequences of ubiquitin genes from Homo sapiens, Drosophila melanogaster and Lepidopteran insects. Sequence analysis showed that the length of the coding sequence is 228 bp, encoding 76 amino acids with calculated molecular weight of 8.50 kDa and the theoretical isoeletric point of 5.26. Signal sequence and transmembrane domain had not been found. Multiple sequence alignment indicated that CsUB gene sequence with other known gene sequences of invertebrates and vertebrates had a high degree of homology (more than 72% similarity) and a shorter genetic distance (lower than 0.360). During the genetic diversity analysis, the total of 104 polymorphic sites was detected from 18 ubiquitin gene sequences and 18 haplotypes were sorted. Abundant genetic diversity and strong codon usage bias were found by the haplotype diversity (1.000), average number of nucleotide differences (47.475), nucleotide diversity (0.20866), effective number of codons (44.526), codon bias index (0.559) and scaled Chi-square (0.779). The predicated secondary structure composition of CsUB protein had about 32.89% extended strands, 36.84% random colis, 15.79% alpha helixes and 14.47% beta turns. Subcellular localization analysis showed that CsUB protein of cytoplasm, cell nucleus, mitochondrion, cell skeleton and plasma membrane occupied about 47.80, 26.10, 17.40, 4.30 and 4.30%, respectively. Sequence, homology and structural analysis confirmed that CsUB gene was highly conserved during evolution and belonged to ubiquitin gene family. The results might provide some fundamental data for further studies on expressed characteristics and physiological functions of CsUB gene. Key words : Chilo suppressalis Walker, ubiquitin, gene cloning, bioinformatics.

Codon usage roles in human papillomavirus

Human papillomavirus (HPV) genomes, similar to other virus genomes, frequently have a G + C content significantly different from their host species. The HPV genomes show a strong codon usage bias to 18 codons, with 14 showing T at the third position amongst degenerately encoded amino acids. The codon usage pattern in HPV genome plays an important role, which regulates low or non-translational expression of the viral capsid genes and results in very weak protein expression of oncogenes in a wide range of mammalian cells. Codon modification has been proved to be a powerful technology to overcome the translational blockage and weak expression of both HPV capsid genes and oncogenes in different expression systems. Furthermore, keratinocytes are the host cells of HPV infection; the codon usage in HPV capsid genes matches available aminoacyl-tRNAs in differentiated keratinocytes to modulate their protein expression. HPV DNA vaccines with codon optimization have been shown to have higher immunogenicity and induce both strong cellular and humoral responses in animal models, which may be a promising form of therapeutic HPV vaccines.

Expression of papillomavirus L1 proteins regulated by authentic gene codon usage is favoured in G2/M-like cells in differentiating keratinocytes

We investigated whether differentiation-dependent expression of papillomavirus (PV) L1 genes is influenced by the cell cycle state in keratinocytes (KCs) grown in vitro or in vivo. In primary keratinocytes, flow cytometry revealed a clear shift from predominantly G0/G1 to G2/M cells from day 1 to day 7, with a three-fold increase in G2/M-like cells in day 7 keratinocytes that showed approximately 50% of the cells expressed a terminal differentiation marker involucrin. The correlation between the levels of the L1 proteins expressed from authentic (Nat) L1 genes of HPV6b and BPV1 and the frequencies of the G2/M-like KCs was significantly positive, while in contrast, a significantly negative correlation in the levels of L1 proteins expressed from codon-modified (Mod) L1 genes of HPV6b and BPV1 with the frequencies of the G2/M-like KCs was observed. Experiments using cell cycle arrest reagents (all-trans retinoic acid (RA) and colchicine) confirmed that L1 proteins expressed from PV Nat L1 genes were facilitated in G2/M-like KCs upon differentiation. Using immunofluorescence microscopy, it appears that L1 proteins from PV Nat L1 genes were co-expressed with cyclin B1, while the L1 proteins expressed from PV Mod L1 genes were preferentially associated with cyclin D2 in KCs in vitro and in mouse skin. Our results demonstrate that (1) expression of the L1 proteins from Nat L1 genes of HPV6b and BPV1 that have strong codon usage bias with A or T at codon third position dependent on KC differentiation is favoured by the G2/M-like environment and (2) codon modifications can alter the cell differentiation-dependent and cell cycle-associated patterns of expression of the PV L1 proteins in KCs.

Unusual composition of a yeast chromosome arm is associated with its delayed replication

The 11.3-Mb genome of the yeast Lachancea (Saccharomyces) kluyveri displays an intriguing compositional heterogeneity: a region of approximately 1 Mb, covering almost the whole left arm of chromosome C (C-left), has an average GC content of 52.9%, which is significantly higher than the 40.4% global GC content of the rest of the genome. This region contains the MAT locus, which remains normal in composition. The excess of GC base pairs affects both coding and noncoding sequences, and thus is not due to selective pressure acting on protein sequences. It leads to a strong codon usage bias and alters the amino acid composition of the 457 proteins encoded on C-left that do not show obvious bias for functional categories, or the presence of paralogs or orthologs of essential genes of Saccharomyces cerevisiae. They share significant synteny conservation with other species of the Saccharomycetaceae, and phylogenetic analysis indicates that C-left originates from a Lachancea species. In contrast, there is a complete absence of transposable elements in C-left, whereas 18 elements per megabase are distributed across the rest of the genome. Comparative hybridization of synchronized cells using high-density genome arrays reveals that C-left is replicated later during S phase than the rest of the genome. Two possible primary causes of this major compositional heterogeneity are discussed: an ancient hybridization of two related species with very distinct GC composition, or an intrinsic mechanism, possibly associated with the loss of the silent cassettes from C-left that progressively increased the GC content and generated the delayed replication of this chromosomal arm.

Consistent and contrasting properties of lineage-specific genes in the apicomplexan parasites Plasmodium and Theileria

BackgroundLineage-specific genes, the genes that are restricted to a limited subset of related organisms, may be important in adaptation. In parasitic organisms, lineage-specific gene products are possible targets for vaccine development or therapeutics when these genes are absent from the host genome.ResultsIn this study, we utilized comparative approaches based on a phylogenetic framework to characterize lineage-specific genes in the parasitic protozoan phylum Apicomplexa. Genes from species in two major apicomplexan genera, Plasmodium and Theileria, were categorized into six levels of lineage specificity based on a nine-species phylogeny. In both genera, lineage-specific genes tend to have a higher level of sequence divergence among sister species. In addition, species-specific genes possess a strong codon usage bias compared to other genes in the genome. We found that a large number of genus- or species-specific genes are putative surface antigens that may be involved in host-parasite interactions. Interestingly, the two parasite lineages exhibit several notable differences. In Plasmodium, the (G + C) content at the third codon position increases with lineage specificity while Theileria shows the opposite trend. Surface antigens in Plasmodium are species-specific and mainly located in sub-telomeric regions. In contrast, surface antigens in Theileria are conserved at the genus level and distributed across the entire lengths of chromosomes.ConclusionOur results provide further support for the model that gene duplication followed by rapid divergence is a major mechanism for generating lineage-specific genes. The result that many lineage-specific genes are putative surface antigens supports the hypothesis that lineage-specific genes could be important in parasite adaptation. The contrasting properties between the lineage-specific genes in two major apicomplexan genera indicate that the mechanisms of generating lineage-specific genes and the subsequent evolutionary fates can differ between related parasite lineages. Future studies that focus on improving functional annotation of parasite genomes and collection of genetic variation data at within- and between-species levels will be important in facilitating our understanding of parasite adaptation and natural selection.

Detecting Key Structural Features within Highly Recombined Genes

Many microorganisms exhibit high levels of intragenic recombination following horizontal gene transfer events. Furthermore, many microbial genes are subject to strong diversifying selection as part of the pathogenic process. A multiple sequence alignment is an essential starting point for many of the tools that provide fundamental insights on gene structure and evolution, such as phylogenetics; however, an accurate alignment is not always possible to attain. In this study, a new analytic approach was developed in order to better quantify the genetic organization of highly diversified genes whose alleles do not align. This BLAST-based method, denoted BLAST Miner, employs an iterative process that places short segments of highly similar sequence into discrete datasets that are designated “modules.” The relative positions of modules along the length of the genes, and their frequency of occurrence, are used to identify sequence duplications, insertions, and rearrangements. Partial alleles of sof from Streptococcus pyogenes, encoding a surface protein under host immune selection, were analyzed for module content. High-frequency Modules 6 and 13 were identified and examined in depth. Nucleotide sequences corresponding to both modules contain numerous duplications and inverted repeats, whereby many codons form palindromic pairs. Combined with evidence for a strong codon usage bias, data suggest that Module 6 and 13 sequences are under selection to preserve their nucleic acid secondary structure. The concentration of overlapping tandem and inverted repeats within a small region of DNA is highly suggestive of a mechanistic role for Module 6 and 13 sequences in promoting aberrant recombination. Analysis of pbp2X alleles from Streptococcus pneumoniae, encoding cell wall enzymes that confer antibiotic resistance, supports the broad applicability of this tool in deciphering the genetic organization of highly recombined genes. BLAST Miner shares with phylogenetics the important predictive quality that leads to the generation of testable hypotheses based on sequence data.