Biased Amino Acid Research Articles

Insect Mouthpart Transcriptome Unveils Extension of Cuticular Protein Repertoire and Complex Organization.

SummaryInsects have developed intriguing cuticles with very specific structures and functions, including microstructures governing their interactions with transmitted microbes, such as in aphid mouthparts harboring virus receptors within such microstructures. Here, we provide the first transcriptome analysis of an insect mouthpart cuticle (“retort organs” [ROs], the stylets' precursors). This analysis defined stylets as a complex composite material. The retort transcriptome also allowed us to propose an algorithmic definition of a new cuticular protein (CP) family with low complexity and biased amino acid composition. Finally, we identified a differentially expressed gene encoding a pyrokinin (PK) neuropeptide precursor and characterizing the mandibular glands. Injection of three predicted synthetic peptides PK1/2/3 into aphids prior to ecdysis caused a molt-specific phenotype with altered head formation. Our study provides the most complete description to date of the potential protein composition of aphid stylets, which should improve the understanding of the transmission of stylet-borne viruses.

Ragp: Pipeline for mining of plant hydroxyproline-rich glycoproteins with implementation in R.

Hydroxyproline-rich glycoproteins (HRGPs) are one of the most complex families of macromolecules found in plants, due to the diversity of glycans decorating the protein backbone, as well as the heterogeneity of the protein backbones. While this diversity is responsible for a wide array of physiological functions associated with HRGPs, it hinders attempts for homology based identification. Current approaches, based on identifying sequences with characteristic motifs and biased amino acid composition, are limited to prototypical sequences. Ragp is an R package for mining and analysis of HRGPs, with emphasis on arabinogalactan proteins. The ragp filtering pipeline exploits one of the HRGPs key features, the presence of hydroxyprolines which represent glycosylation sites. Main package features include prediction of proline hydroxylation sites, amino acid motif and bias analyses, efficient communication with web servers for prediction of N-terminal signal peptides, glycosylphosphatidylinositol modification sites and disordered regions and the ability to annotate sequences through hmmscan and subsequent GO enrichment, based on predicted Pfam domains. As such, ragp extends R's rich ecosystem for high-throughput sequence data analyses. The ragp R package is available under the MIT Open Source license and is freely available to download from GitHub at: https://github.com/missuse/ragp.

Bacillus lacisalsi sp. nov., a moderately haloalkaliphilic bacterium isolated from a saline-alkaline lake.

An alkaliphilic and moderately halophilic strain, designated YSP-3T, characterised by optimal growth at pH 9.0 and at 8.0% (w/v) NaCl, was isolated from Yangshapao Lake, Jilin Province, China. Cells of this strain is Gram-positive, straight rods and form a central or sub-terminal ellipsoidal endospore. Phylogenetic analysis based on 16S rRNA gene sequences indicated that it was grouped in the genus Bacillus with Bacillus aurantiacus K1-5T and Bacillus populi FJAT-45347T as the close relative (97.5 and 97.2% 16S rRNA gene sequence similarity, respectively). Genomic relatedness between strain YSP-3T and its close relative was evaluated using average nucleotide identity, digital DNA-DNA hybridization and average amino acid identity with the values of 70.3-85.1%, 19.7-20.1% and 71.5-71.6%, respectively. Comparative genomics analysis showed that strain YSP-3T has distinct amino acid bias and significantly differences from foreign invasion events during evolution relative to the reference strains. Cell-wall peptidoglycan contains meso-diaminopimelic acid. The predominant polar lipids are phosphatidylethanolamine, diphosphatidylglycerol, phosphatidylglycerol. The predominant quinone is menaquinone-7. The major fatty acids of strain YSP-3T are anteiso-C15:0, iso-C15:0, iso-C16:0, anteiso-C17:0 and Iso-C14:0. DNA G + C content of strain YSP-3T is 48.3mol%. Based on genomics analysis, physiological, biochemical and chemotaxonomic data, strain YSP-3T represent a novel species, for which the name Bacillus lacisalsi sp. nov. is proposed. The type strain is YSP-3T ( = ACCC 60365T = KCTC 33934T).

Universal principles of membrane protein assembly, composition and evolution.

Structural diversity in α-helical membrane proteins (MP) arises from variations in helix-helix crossings and contacts that may bias amino acid usage. Here, we reveal systematic changes in transmembrane amino acid frequencies (f) as a function of the number of helices (n). For eukarya, breaks in f(n) trends of packing (Ala, Gly and Pro), polar, and hydrophobic residues identify different MP assembly principles for 2≤n≤7, 8≤n≤12 and n≥13. In bacteria, the first f break already occurs after n = 6 in correlation to an earlier n peak in MP size distribution and dominance of packing over polar interactions. In contrast to the later n brackets, the integration levels of helix bundles continuously increased in the first, most populous brackets indicating the formation of single structural units (domains). The larger first bracket of eukarya relates to a balance of polar and packing interactions that enlarges helix-helix combinatorial possibilities (MP diversity). Between the evolutionary old, packing and new, polar residues f anti-correlations extend over all biological taxa, broadly ordering them according to evolutionary history and allowing f estimates for the earliest forms of life. Next to evolutionary history, the amino acid composition of MP is determined by size (n), proteome diversity, and effective amino acid cost.

Sequence Versus Composition: What Prescribes IDP Biophysical Properties?

Intrinsically disordered proteins (IDPs) represent a distinct class of proteins and are distinguished from globular proteins by conformational plasticity, high evolvability and a broad functional repertoire. Some of their properties are reminiscent of early proteins, but their abundance in eukaryotes, functional properties and compositional bias suggest that IDPs appeared at later evolutionary stages. The spectrum of IDP properties and their determinants are still not well defined. This study compares rudimentary physicochemical properties of IDPs and globular proteins using bioinformatic analysis on the level of their native sequences and random sequence permutations, addressing the contributions of composition versus sequence as determinants of the properties. IDPs have, on average, lower predicted secondary structure contents and aggregation propensities and biased amino acid compositions. However, our study shows that IDPs exhibit a broad range of these properties. Induced fold IDPs exhibit very similar compositions and secondary structure/aggregation propensities to globular proteins, and can be distinguished from unfoldable IDPs based on analysis of these sequence properties. While amino acid composition seems to be a major determinant of aggregation and secondary structure propensities, sequence randomization does not result in dramatic changes to these properties, but for both IDPs and globular proteins seems to fine-tune the tradeoff between folding and aggregation.

Low cationicity is important for systemic in vivo efficacy of database-derived peptides against drug-resistant Gram-positive pathogens

As bacterial resistance to traditional antibiotics continues to emerge, new alternatives are urgently needed. Antimicrobial peptides (AMPs) are important candidates. However, how AMPs are designed with in vivo efficacy is poorly understood. Our study was designed to understand structural moieties of cationic peptides that would lead to their successful use as antibacterial agents. In contrast to the common perception, serum binding and peptide stability were not the major reasons for in vivo failure in our studies. Rather, our systematic study of a series of peptides with varying lysines revealed the significance of low cationicity for systemic in vivo efficacy against Gram-positive pathogens. We propose that peptides with biased amino acid compositions are not favored to associate with multiple host factors and are more likely to show in vivo efficacy. Thus, our results uncover a useful design strategy for developing potent peptides against multidrug-resistant pathogens.

Codon selection reduces GC content bias in nucleic acids encoding for intrinsically disordered proteins.

Protein-coding nucleic acids exhibit composition and codon biases between sequences coding for intrinsically disordered regions (IDRs) and those coding for structured regions. IDRs are regions of proteins that are folding self-insufficient and which function without the prerequisite of folded structure. Several authors have investigated composition bias or codon selection in regions encoding for IDRs, primarily in Eukaryota, and concluded that elevated GC content is the result of the biased amino acid composition of IDRs. We substantively extend previous work by examining GC content in regions encoding IDRs, from 44 species in Eukaryota, Archaea, and Bacteria, spanning a wide range of GC content. We confirm that regions coding for IDRs show a significantly elevated GC content, even across all domains of life. Although this is largely attributable to the amino acid composition bias of IDRs, we show that this bias is independent of the overall GC content and, most importantly, we are the first to observe that GC content bias in IDRs is significantly different than expected from IDR amino acid composition alone. We empirically find compensatory codon selection that reduces the observed GC content bias in IDRs. This selection is dependent on the overall GC content of the organism. The codon selection bias manifests as use of infrequent, AT-rich codons in encoding IDRs. Further, we find these relationships to be independent of the intrinsic disorder prediction method used, and independent of estimated translation efficiency. These observations are consistent with the previous work, and we speculate on whether the observed biases are causal or symptomatic of other driving forces.

Analysis of the Preserved Amino Acid Bias in Peptide Profiles of Iron Age Teeth from a Tropical Environment Enable Sexing of Individuals Using Amelogenin MRM

The first dental proteomic profile of Iron Age individuals (ca. 2000-1000 years B.P.), collected from the site of Long Long Rak rock shelter in northwest Thailand is described. A bias toward the preservation of the positively charged aromatic, and polar amino acids is observed. It is evident that the 212 proteins identified (2 peptide, FDR <1%) comprise a palimpsest of alterations that occurred both ante-mortem and post-mortem. Conservation of amino acids within the taphonomically resistant crystalline matrix enabled the identification of both X and Y chromosome linked amelogenin peptides. A novel multiple reaction monitoring method using the sex specific amelogenin protein isoforms is described and indicate the teeth are of male origin. Functional analysis shows an enrichment of pathways associated with metabolic disorders and shows a capacity for harboring these conditions prior to death. Stable isotope analysis using carbon isotopes highlights the strongly C3 based (≈80%) diet of the Long Long Rak cemetery people, which probably comprised rice combined with protein from freshwater fish among other food items. The combination of proteomics and stable isotope analysis provides a complementary strategy for assessing the demography, diet, lifestyle, and possible diseases experienced by ancient populations.

Solution Ensemble of the C-Terminal Domain from the Transcription Factor Pdx1 Resembles an Excluded Volume Polymer.

The pancreatic and duodenal homeobox 1 (Pdx1) is an essential pancreatic transcription factor. The C-terminal intrinsically disordered domain of Pdx1 (Pdx1-C) has a heavily biased amino acid composition; most notably, 18 of 83 residues are proline, including a hexaproline cluster near the middle of the chain. For these reasons, Pdx1-C is an attractive target for structure characterization, given the availability of suitable methods. To determine the solution ensembles of disordered proteins, we have developed a suite of 13C direct-detect NMR experiments that provide high spectral quality, even in the presence of strong proline enrichment. Here, we have extended our suite of NMR experiments to include four new pulse programs designed to record backbone residual dipolar couplings in a 13C,15N-CON detection format. Using our NMR strategy, in combination with small-angle X-ray scattering measurements and Monte Carlo simulations, we have determined that Pdx1-C is extended in solution, with a radius of gyration and internal scaling similar to that of an excluded volume polymer, and a subtle tendency toward a collapsed structure to the N-terminal side of the hexaproline sequence. This structure leaves Pdx1-C exposed for interactions with trans-regulatory co-factors that contribute with Pdx1 to transcription control in the cell.

Comparative functional analysis of proteins containing low-complexity predicted amyloid regions.

BackgroundIn both prokaryotic and eukaryotic proteins, repeated occurrence of a single or a group of few amino acids are found. These regions are termed as low complexity regions (LCRs). It has been observed that amino acid bias in LCR is directly linked to their uncontrolled expansion and amyloid formation. But a comparative analysis of the behavior of LCR based on their constituent amino acids and their association with amyloidogenic propensity is not available.MethodsFirstly we grouped all LCRs on the basis of their composition: homo-polymers, positively charged amino acids, negatively charged amino acids, polar amino acids and hydrophobic amino acids. We analyzed the compositional pattern of LCRs in each group and their propensity to form amyloids. The functional characteristics of proteins containing different groups of LCRs were explored using DAVID. In addition, we also analyzed the classes, pathways and functions of human proteins that form amyloids in LCRs.ResultsAmong homopolymeric LCRs, the most common was Gln repeats. LCRs composed of repeats of Met and aromatic amino acids were amongst the least occurring. The results revealed that LCRs composed of negatively charged and polar amino acids were more common in comparison to LCRs formed by positively charged and hydrophobic amino acids. We also noted that generally proteins with LCRs were involved in transcription but those with Gly repeats were associated to translational activities. Our analysis suggests that proteins in which LCR is composed of hydrophobic residues are more prone toward amyloid formation. We also found that the human proteins with amyloid forming LCRs were generally involved in binding and catalytic activity.DiscussionThe presented analysis summarizes the most common and least occurring LCRs in proteins. Our results show that though repeats of Gln are the most abundant but Asn repeats make longest stretch of low complexity. The results showed that potential of LCRs to form amyloids varies with their amino acid composition.

Three new Diplozoidae mitogenomes expose unusual compositional biases within the Monogenea class: implications for phylogenetic studies

BackgroundAs the topologies produced by previous molecular and morphological studies were contradictory and unstable (polytomy), evolutionary relationships within the Diplozoidae family and the Monogenea class (controversial relationships among the Discocotylinea, Microcotylinea and Gastrocotylinea suborders) remain unresolved. Complete mitogenomes carry a relatively large amount of information, sufficient to provide a much higher phylogenetic resolution than traditionally used morphological traits and/or single molecular markers. However, their implementation is hampered by the scarcity of available monogenean mitogenomes. Therefore, we sequenced and characterized mitogenomes belonging to three Diplozoidae family species, and conducted comparative genomic and phylogenomic analyses for the entire Monogenea class.ResultsTaxonomic identification was inconclusive, so two of the species were identified merely to the genus level. The complete mitogenomes of Sindiplozoon sp. and Eudiplozoon sp. are 14,334 bp and 15,239 bp in size, respectively. Paradiplozoon opsariichthydis (15,385 bp) is incomplete: an approximately 2000 bp-long gap within a non-coding region could not be sequenced. Each genome contains the standard 36 genes (atp8 is missing). G + T content and the degree of GC- and AT-skews of these three mitogenome (and their individual elements) were higher than in other monogeneans. nad2, atp6 and nad6 were the most variable PCGs, whereas cox1, nad1 and cytb were the most conserved. Mitochondrial phylogenomics analysis, conducted using concatenated amino acid sequences of all PCGs, indicates that evolutionary relationships of the three genera are: (Eudiplozoon, (Paradiplozoon, Sindiplozoon)); and of the three suborders: (Discocotylinea, (Microcotylinea, Gastrocotylinea)). These intergeneric relationships were also supported by the skewness and principal component analyses.ConclusionsOur results show that nad2, atp6 and nad6 (fast-evolving) would be better candidates than cox1 (slow-evolving) for species identification and population genetics studies in Diplozoidae. Nucleotide bias and codon and amino acid usage patterns of the three diplozoid mitogenomes are more similar to cestodes and trematodes than to other monogenean flatworms. This unusual mutational bias was reflected in disproportionately long branches in the phylogram. Our study offsets the scarcity of molecular data for the subclass Polyopisthocotylea to some extent, and might provide important new insights into the evolutionary history of the three genera and three suborders.

Rapid and Error-Free Site-Directed Mutagenesis by a PCR-Free In Vitro CRISPR/Cas9-Mediated Mutagenic System.

The quality and efficiency of any PCR-based mutagenesis technique may not be optimal due to, among other things, amino acid bias, which means that the development of efficient PCR-free methods is desirable. Here, we present a highly efficient in vitro CRISPR/Cas9-mediated mutagenic (ICM) system that allows rapid construction of designed mutants in a PCR-free manner. First, it involves plasmid digestion by utilizing a complex of Cas9 with specific single guide RNA (sgRNA) followed by degradation with T5 exonuclease to generate a 15 nt homologous region. Second, primers containing the desired mutations are annealed to form the double-stranded DNA fragments, which are then ligated into the linearized plasmid. In theory, neither the size of the target plasmid nor the unavailable restriction enzyme site poses any problems that may arise in traditional techniques. In this study, single and multiple site-directed mutagenesis were successfully performed even for a large size plasmid (up to 9.0 kb). Moreover, a PCR-free site-saturation mutagenesis library on single site and two adjacent sites of a green fluorescent protein was also generated with promising results. This demonstrates the great potential of the ICM system for creating high-quality mutant libraries in directed evolution as an alternative to PCR-based saturation mutagenesis, thus facilitating research on synthetic biology.

Plasmid-Based One-Pot Saturation Mutagenesis and Robot-Based Automated Screening for Protein Engineering.

We evaluated a method for protein engineering using plasmid-based one-pot saturation mutagenesis and robot-based automated screening. When the biases in nucleotides and amino acids were assessed for a loss-of-function point mutation in green fluorescent protein, the ratios of gain-of-function mutants were not significantly different from the expected values for the primers among the three different suppliers. However, deep sequencing analysis revealed that the ratios of nucleotides in the primers were highly biased among the suppliers. Biases for NNB were less severe than for NNN. We applied this method to screen a fusion protein of two chitinases, ChiA and ChiB (ChiAB). Three NNB codons as well as tyrosine and serine (X1YSX2X3) were inserted to modify the surface structure of ChiAB. We observed significant amino acid bias at the X3 position in water-soluble, active ChiAB-X1YSX2X3 mutants. Examination of the crystal structure of one active mutant, ChiAB-FYSFV, revealed that the X3 residue plays an important role in structure stabilization.

GC bias lead to increased small amino acids and random coils of proteins in cold-water fishes

BackgroundTemperature adaptation of biological molecules is fundamental in evolutionary studies but remains unsolved. Fishes living in cold water are adapted to low temperatures through adaptive modification of their biological molecules, which enables their functioning in extreme cold. To study nucleotide and amino acid preference in cold-water fishes, we investigated the substitution asymmetry of codons and amino acids in protein-coding DNA sequences between cold-water fishes and tropical fishes., The former includes two Antarctic fishes, Dissostichus mawsoni (Antarctic toothfish), Gymnodraco acuticeps (Antarctic dragonfish), and two temperate fishes, Gadus morhua (Atlantic cod) and Gasterosteus aculeatus (stickleback), and the latter includes three tropical fishes, including Danio rerio (zebrafish), Oreochromis niloticus (Nile tilapia) and Xiphophorus maculatus (Platyfish).ResultsCold-water fishes showed preference for Guanines and cytosines (GCs) in both synonymous and nonsynonymous codon substitution when compared with tropical fishes. Amino acids coded by GC-rich codons are favored in the temperate fishes, while those coded by AT-rich codons are disfavored. Similar trends were discovered in Antarctic fishes but were statistically weaker. The preference of GC rich codons in nonsynonymous substitution tends to increase ratio of small amino acid in proteins, which was demonstrated by biased small amino acid substitutions in the cold-water species when compared with the tropical species, especially in the temperate species. Prediction and comparison of secondary structure of the proteomes showed that frequency of random coils are significantly larger in the cold-water fish proteomes than those of the tropical fishes.ConclusionsOur results suggested that natural selection in cold temperature might favor biased GC content in the coding DNA sequences, which lead to increased frequency of small amino acids and consequently increased random coils in the proteomes of cold-water fishes.

The unique evolutionary pattern of the Hydroxyproline-rich glycoproteins superfamily in Chinese white pear (Pyrus bretschneideri)

BackgroundThe hydroxyproline-rich glycoprotein (HRGP) superfamily, comprising three families (arabinogalactan-proteins, AGPs; extensins, EXTs; proline-rich proteins, PRPs), is a class of proline-rich proteins that exhibit high diversity and are involved in many aspects of plant biology.ResultsIn this study, 838 HRGPs were identified from Chinese white pear (Pyrus bretschneideri) by searching for biased amino acid composition and conserved motifs. 405 HRGPs were derived from whole genome duplication (WGD) events which is suggested to be the major force of driving HRGPs expansion and the recent WGD event shared by apple and pear generated most duplicated HRGPs in pear. This duplication event drived the structural variation of the HRGPs encoding hydroxyproline (Hyp)-rich motifs. The rate of HRGPs evolution mainly impacted the Hyp-rich motifs even in chimeric HRGPs. During the evolution of 53 PRPs that are also typified by 7-deoxyloganetin glucosyltransferase-like genes, the duplication from PRP to non-PRP was indirectly modified by positive selection. These results suggested that the rate of HRGP evolution mainly influenced the Hyp-rich motifs even in chimeric HRGPs. The expression divergence of HRGPs was higher than that of other commonly duplicated genes. In pear pistil, 601 HRGPs exhibited expression, while in pear pollen, 285 HRGPs were expressed. The qPCR results revealed that Pbr036330.1 and Pbr010506.1 showed different expression profile in self-incompatibility of pear pistil.ConclusionsThe researches indicated that WGD events was the main duplication type during the evolution of HRGPs, and the highly variable Hyp-motifs might be accountable for the expansion, evolution and expression divergence of HRGPs and that this divergence may be responsible for the gain of new functions in plants.

Beating Bias in the Directed Evolution of Proteins: Combining High-Fidelity on-Chip Solid-Phase Gene Synthesis with Efficient Gene Assembly for Combinatorial Library Construction.

Saturation mutagenesis (SM) constitutes a widely used technique in the directed evolution of selective enzymes as catalysts in organic chemistry and in the manipulation of metabolic paths and genomes, but the quality of the libraries is far from optimal due to the inherent amino acid bias. Herein, it is shown how this fundamental problem can be solved by applying high-fidelity solid-phase chemical gene synthesis on silicon chips followed by efficient gene assembly. Limonene epoxide hydrolase was chosen as the catalyst in the model desymmetrization of cyclohexene oxide with the stereoselective formation of (R,R)- and (S,S)-cyclohexane-1,2-diol. A traditional combinatorial PCR-based SM library, produced by simultaneous randomization at several residues by using a reduced amino acid alphabet, and the respective synthetic library were constructed and compared. Statistical analysis at the DNA level with massive sequencing demonstrates that, in the synthetic approach, 97 % of the theoretically possible DNA mutants are formed, whereas the traditional SM library contained only about 50 %. Screening at the protein level also showed the superiority of the synthetic library; many highly (R,R)- and (S,S)-selective variants being discovered are not found in the traditional SM library. With the prices of synthetic genes decreasing, this approach may point the way to future directed evolution.

The hnRNP-like Nab3 termination factor can employ heterologous prion-like domains in place of its own essential low complexity domain.

Many RNA-binding proteins possess domains with a biased amino acid content. A common property of these low complexity domains (LCDs) is that they assemble into an ordered amyloid form, juxtaposing RNA recognition motifs in a subcellular compartment in which RNA metabolism is focused. Yeast Nab3 is one such protein that contains RNA-binding domains and a low complexity, glutamine/proline-rich, prion-like domain that can self-assemble. Nab3 also contains a region of structural homology to human hnRNP-C that resembles a leucine zipper which can oligomerize. Here we show that the LCD and the human hnRNP-C homology domains of Nab3 were experimentally separable, as cells were viable with either segment, but not when both were missing. In exploiting the lethality of deleting these regions of Nab3, we were able to test if heterologous prion-like domains known to assemble into amyloid, could substitute for the native sequence. Those from the hnRNP-like protein Hrp1, the canonical prion Sup35, or the epsin-related protein Ent2, could rescue viability and enable the new Nab3 chimeric protein to support transcription termination. Other low complexity domains from RNA-binding, termination-related proteins or a yeast prion, could not. As well, an unbiased genetic selection revealed a new protein sequence that could rescue the loss of Nab3’s essential domain via multimerization. This new sequence and Sup35’s prion domain could also rescue the lethal loss of Hrp1’s prion-like domain when substituted for it. This suggests there are different cross-functional classes of amyloid-forming LCDs and that appending merely any assembly-competent LCD to Nab3 does not restore function or rescue viability. The analysis has revealed the functional complexity of LCDs and provides a means by which the differing classes of LCD can be dissected and understood.

An unfolded protein-induced conformational switch activates mammalian IRE1.

The unfolded protein response (UPR) adjusts the cell's protein folding capacity in the endoplasmic reticulum (ER) according to need. IRE1 is the most conserved UPR sensor in eukaryotic cells. It has remained controversial, however, whether mammalian and yeast IRE1 use a common mechanism for ER stress sensing. Here, we show that similar to yeast, human IRE1α's ER-lumenal domain (hIRE1α LD) binds peptides with a characteristic amino acid bias. Peptides and unfolded proteins bind to hIRE1α LD's MHC-like groove and induce allosteric changes that lead to its oligomerization. Mutation of a hydrophobic patch at the oligomerization interface decoupled peptide binding to hIRE1α LD from its oligomerization, yet retained peptide-induced allosteric coupling within the domain. Importantly, impairing oligomerization of hIRE1α LD abolished IRE1's activity in living cells. Our results provide evidence for a unifying mechanism of IRE1 activation that relies on unfolded protein binding-induced oligomerization.

Codon usage bias in 5′ terminal coding sequences reveals distinct enrichment of gene functions

Codon bias at the 5′ terminal of coding sequence (CDS) is known to be distinct from the rest of the CDS. A number of events occur in this short region to regulate early translation elongation and co-translational translocation. In the genes encoding secretory proteins, there is a special signal sequence which has a higher occurrence of rare codons. In this study, we analyzed codon bias of secretory genes in several eukaryotes. The results showed that secretory genes in the species except mammals had a higher occurrence of rare codons in the 5′ terminal of CDS, and the bias was greater than the same region of non-secretory genes. GO analysis revealed that secretory genes containing rare codon clusters in different regions were responsible for various roles in gene functions. Moreover, codon bias in the region encoding the hydrophobic region of protein is similar in secretory and non-secretory genes, indicating that codon bias in secretory genes was partly influenced by amino acid bias. Rare codon clusters are found more frequently in specific regions, and continuous rare codons are not favoured probably because they will increase the probability of ribosome collision and drop-off. Based on ribosome profiling data, there is no significant difference in the average translation efficiencies between rare and optimal codons. Higher ribosomal density in the 5′ terminal may result from ribosome pausing which could be involved in different translation events. These findings collectively provided rich information on codon bias in secretory genes, which may shed light on the co-effect of codon bias, mRNA structure and tRNA abundance in translational regulations.

Insights into the Evolution of Hydroxyproline-Rich Glycoproteins from 1000 Plant Transcriptomes.

The carbohydrate-rich cell walls of land plants and algae have been the focus of much interest given the value of cell wall-based products to our current and future economies. Hydroxyproline-rich glycoproteins (HRGPs), a major group of wall glycoproteins, play important roles in plant growth and development, yet little is known about how they have evolved in parallel with the polysaccharide components of walls. We investigate the origins and evolution of the HRGP superfamily, which is commonly divided into three major multigene families: the arabinogalactan proteins (AGPs), extensins (EXTs), and proline-rich proteins. Using motif and amino acid bias, a newly developed bioinformatics pipeline, we identified HRGPs in sequences from the 1000 Plants transcriptome project (www.onekp.com). Our analyses provide new insights into the evolution of HRGPs across major evolutionary milestones, including the transition to land and the early radiation of angiosperms. Significantly, data mining reveals the origin of glycosylphosphatidylinositol (GPI)-anchored AGPs in green algae and a 3- to 4-fold increase in GPI-AGPs in liverworts and mosses. The first detection of cross-linking (CL)-EXTs is observed in bryophytes, which suggests that CL-EXTs arose though the juxtaposition of preexisting SPn EXT glycomotifs with refined Y-based motifs. We also detected the loss of CL-EXT in a few lineages, including the grass family (Poaceae), that have a cell wall composition distinct from other monocots and eudicots. A key challenge in HRGP research is tracking individual HRGPs throughout evolution. Using the 1000 Plants output, we were able to find putative orthologs of Arabidopsis pollen-specific GPI-AGPs in basal eudicots.