Global analysis of proline-rich tandem repeat proteins reveals broad phylogenetic diversity in plant secretomes

Identification and analysis of proline-rich proteins and hybrid proline-rich proteins super family genes from Sorghum bicolor and their expression patterns to abiotic stress and zinc stimuli

Systematic genome-wide analysis of Sorghum bicolor revealed the identification of a total of 48 homologous genes comprising 21 proline-rich proteins (PRPs) and 27 hybrid proline-rich proteins (HyPRPs). Comprehensive scrutiny of these gene homologs was conducted for gene structure, phylogenetic investigations, chromosome mapping, and subcellular localization of proteins. Promoter analysis uncovered the regions rich with phosphorous- (BIHD), ammonium-, sulfur-responsive (SURE), and iron starvation-responsive (IRO2) along with biotic, abiotic, and development-specific cis-elements. Further, PRPs exhibit more methylation and acetylation sites in comparison with HyPRPs. miRNAs have been predicted which might play a role in cleavage and translation inhibition. Several of the SbPRP genes were stimulated in a tissue-specific manner under drought, salt, heat, and cold stresses. Additionally, exposure of plants to abscisic acid (ABA) and zinc (Zn) also triggered PRP genes in a tissue-dependent way. Among them, SbPRP17 has been found upregulated markedly in all tissues irrespective of the stress imposed. The expressions of SbHyPRPs, especially SbHyPRP2, SbHyPRP6, and SbHyPRP17 were activated under all stresses in all three tissues. On the other hand, SbHyPRP8 (root only) and SbHyPRP12 (all three tissues) were highly responsive to cold stress and ABA while SbHyPRP26 was induced by drought and Zn in the stem. Taken together, this study indicates the critical roles that SbPRPs and SbHyPRPs play during diverse abiotic stress conditions and notably the plausible roles that these genes play upon exposure to zinc, the crucial micronutrient in plants.

Bioinformatic Identification of Plant Hydroxyproline-Rich Glycoproteins

Hydroxyproline-rich glycoproteins (HRGPs) are a superfamily of plant cell wall proteins that function in diverse aspects of plant growth and development. This superfamily consists of three members: arabinogalactan-proteins (AGPs), extensins (EXTs), and proline-rich proteins (PRPs). Hybrid and chimeric HRGPs also exist. A bioinformatic software program, BIO OHIO 2.0, was developed to expedite the genome-wide identification and classification of AGPs, EXTs, and PRPs based on characteristic HRGP motifs and biased amino acid compositions. This chapter explains the principles of identifying HRGPs and provides a stepwise tutorial for using the BIO OHIO 2.0 program with genomic/proteomic data. Here, as an example, the genome/proteome of the common bean (Phaseolus vulgaris) is analyzed using the BIO OHIO 2.0 program to identify and characterize its set of HRGPs.

PRPs localized to the middle lamellae are required for cortical tissue integrity in Medicago truncatula roots

A family of repetitive proline-rich proteins interact with acidic pectins and play distinct roles in legume root cell walls affecting cortical and vascular structure. A proline-rich protein (PRP) family, composed of tandemly repeated Pro-Hyp-Val-X-Lys pentapeptide motifs, is found primarily in the Leguminosae. Four distinct size classes within this family are encoded by seven tightly linked genes: MtPRP1, MtPRP2 and MtPRP3, and four nearly identical MtPRP4 genes. Promoter fusions to β-glucuronidase showed strong expression in the stele of hairy roots for all 4 PRP genes tested, with additional expression in the cortex for PRP1, PRP2 and PRP4. All except MtPRP4 are strongly expressed in non-tumorous roots, and secreted and ionically bound to root cell walls. These PRPs are absent from root epidermal cell walls, and PRP accumulation is highly localized within the walls of root cortical and vascular tissues. Within xylem tissue, PRPs are deposited in secondary thickenings where it is spatially exclusive to lignin. In newly differentiating xylem, PRPs are deposited in the regularly spaced paired-pits and pit membranes that hydraulically connect neighboring xylem elements. Hairpin-RNA knock-down constructs reducing PRP expression in Medicago truncatula hairy root tumors disrupted cortical and vascular patterning. Immunoblots showed that the knockdown tumors had potentially compensating increases in the non-targeted PRPs, all of which cross-react with the anti-PRP antibodies. However, PRP3 knockdown differed from knockdown of PRP1 and PRP2 in that it greatly reduced viability of hairy root tumors. We hypothesize that repetitive PRPs interact with acidic pectins to form block-copolymer gels that can play distinct roles in legume root cell walls.

DNA Binding to Gold Nanoparticles through the Prism of Molecular Selection: Sequence-Affinity Relation

Native DNA strongly adsorbs to citrate-coated gold nanoparticles (AuNPs). The resulting composites (DNA/AuNPs) are valuable materials in many fields, especially in biomedicine. For this reason, the process of adsorption is a focus for intensive research. In this work, DNA adsorption to gold nanoparticles was studied using a molecular selection procedure followed by high-throughput DNA sequencing. The chemically synthesized DNA library containing a central N₂₆ randomized fragment was sieved through four cycles of adsorption to AuNPs in a tree-like selection-amplification scheme (SELEX (Selective Evolution of Ligands by EXponential enrichment)). The frequencies of occurrence of specific oligomeric DNA motifs, k-mers ( k = 1-6), in the initial and selected pools were calculated. Distribution of secondary structures in the pools was analyzed. A large set of diverse A, T, and G enriched k-mers undergo a pronounced positive selection, and these sequences demonstrate faster and strong binding to the AuNPs. For facile binding, such structural motifs should be located in the loop regions of weak intramolecular complexes-hairpins with imperfect stem, or other portion of the structure, which is unpaired under selection conditions. Our data also show that, under the conditions employed in this study, cytosine is significantly depleted during the selection process, although guanine remains unchanged. These regularities were confirmed in a series of binding experiments with a set of synthetic DNA oligonucleotides. The detailed analysis of DNA binding to AuNPs shows that the sequence specificity of this interaction is low due to its nature, although the presence and the number of specific structural motifs in DNA affect both the rate of formation and the strength of the formed noncovalent associates with AuNPs.

Sequence and structural analysis of fibronectin-binding protein reveals importance of multiple intrinsic disordered tandem repeats

The location of certain amino acid sequences like repeats along the polypeptide chain is very important in the context of forming the overall shape of the protein molecule which in fact determines its function. In gram-positive bacteria, fibronectin-binding protein (FnBP) is one such repeat containing protein, and it is a cell wall-attached protein responsible for various acute infections in human. Several studies on sequence, structure, and function of fibronectin-binding regions of FnBPs were reported; however, no detailed study was carried out on the full-length protein sequence. In the present study, we have made a thorough sequence and structure analysis on FnBP_A of Staphylococcus aureus and explored the presence of dual ligand-binding ability of fibrinogen (fg)-binding region and its molecular recognition processes. Multiple sequence alignment and protein-protein docking analysis reveal the regions which are likely involved in dual ligand binding. Further analysis of docking of FnBP_A fg-binding region and fn N-terminal modules suggests that if the latter binds to the fg-binding region of FnBP_A, it would inhibit the subsequent binding of fg because of steric hindrance. The sequence analysis further suggests that the abundance of disorder promoting residue glutamic acid and dual personality (both order/disorder promoting) residue threonine in tandem repeats of FnBP_A and B proteins possibly would help the molecule to undergo a conformational change while binding with fn by β-zipper mechanism. The segment-based power spectral analysis was carried out which helps to understand the distribution of hydrophobic residues along the sequence particularly in intrinsic disordered tandem repeats. The results presented here will help to understand the role of internal repeats and intrinsic disorder in the molecular recognition process of a pathogenic cell surface protein.

Mining of Brassica-Specific Genes (BSGs) and Their Induction in Different Developmental Stages and under Plasmodiophora brassicae Stress in Brassica rapa

Orphan genes, also called lineage-specific genes (LSGs), are important for responses to biotic and abiotic stresses, and are associated with lineage-specific structures and biological functions. To date, there have been no studies investigating gene number, gene features, or gene expression patterns of orphan genes in Brassica rapa. In this study, 1540 Brassica-specific genes (BSGs) and 1824 Cruciferae-specific genes (CSGs) were identified based on the genome of Brassica rapa. The genic features analysis indicated that BSGs and CSGs possessed a lower percentage of multi-exon genes, higher GC content, and shorter gene length than evolutionary-conserved genes (ECGs). In addition, five types of BSGs were obtained and 145 out of 529 real A subgenome-specific BSGs were verified by PCR in 51 species. In silico and semi-qPCR, gene expression analysis of BSGs suggested that BSGs are expressed in various tissue and can be induced by Plasmodiophora brassicae. Moreover, an A/C subgenome-specific BSG, BSGs1, was specifically expressed during the heading stage, indicating that the gene might be associated with leafy head formation. Our results provide valuable biological information for studying the molecular function of BSGs for Brassica-specific phenotypes and biotic stress in B. rapa.

Pipeline to Identify Hydroxyproline-Rich Glycoproteins

Intrinsically disordered proteins (IDPs) are functional proteins that lack a well-defined three-dimensional structure. The study of IDPs is a rapidly growing area as the crucial biological functions of more of these proteins are uncovered. In plants, IDPs are implicated in plant stress responses, signaling, and regulatory processes. A superfamily of cell wall proteins, the hydroxyproline-rich glycoproteins (HRGPs), have characteristic features of IDPs. Their protein backbones are rich in the disordering amino acid proline, they contain repeated sequence motifs and extensive posttranslational modifications (glycosylation), and they have been implicated in many biological functions. HRGPs are evolutionarily ancient, having been isolated from the protein-rich walls of chlorophyte algae to the cellulose-rich walls of embryophytes. Examination of HRGPs in a range of plant species should provide valuable insights into how they have evolved. Commonly divided into the arabinogalactan proteins, extensins, and proline-rich proteins, in reality, a continuum of structures exists within this diverse and heterogenous superfamily. An inability to accurately classify HRGPs leads to inconsistent gene ontologies limiting the identification of HRGP classes in existing and emerging omics data sets. We present a novel and robust motif and amino acid bias (MAAB) bioinformatics pipeline to classify HRGPs into 23 descriptive subclasses. Validation of MAAB was achieved using available genomic resources and then applied to the 1000 Plants transcriptome project (www.onekp.com) data set. Significant improvement in the detection of HRGPs using multiple-k-mer transcriptome assembly methodology was observed. The MAAB pipeline is readily adaptable and can be modified to optimize the recovery of IDPs from other organisms.

Bioinformatic Identification and Analysis of Hydroxyproline-Rich Glycoproteins in Populus trichocarpa

BACKGROUND

Hydroxyproline-rich glycoproteins (HRGPs) constitute a plant cell wall protein superfamily that functions in diverse aspects of growth and development. This superfamily contains three members: the highly glycosylated arabinogalactan-proteins (AGPs), the moderately glycosylated extensins (EXTs), and the lightly glycosylated proline-rich proteins (PRPs). Chimeric and hybrid HRGPs, however, also exist. A bioinformatics approach is employed here to identify and classify AGPs, EXTs, PRPs, chimeric HRGPs, and hybrid HRGPs from the proteins predicted by the completed genome sequence of poplar (Populus trichocarpa). This bioinformatics approach is based on searching for biased amino acid compositions and for particular protein motifs associated with known HRGPs with a newly revised and improved BIO OHIO 2.0 program. Proteins detected by the program are subsequently analyzed to identify the following: 1) repeating amino acid sequences, 2) signal peptide sequences, 3) glycosylphosphatidylinositol lipid anchor addition sequences, and 4) similar HRGPs using the Basic Local Alignment Search Tool (BLAST).

RESULTS

The program was used to identify and classify 271 HRGPs from poplar including 162 AGPs, 60 EXTs, and 49 PRPs, which are each divided into various classes. This is in contrast to a previous analysis of the Arabidopsis proteome which identified 162 HRGPs consisting of 85 AGPs, 59 EXTs, and 18 PRPs. Poplar was observed to have fewer classical EXTs, to have more fasciclin-like AGPs, plastocyanin AGPs and AG peptides, and to contain a novel class of PRPs referred to as the proline-rich peptides.

CONCLUSIONS

The newly revised and improved BIO OHIO 2.0 bioinformatics program was used to identify and classify the inventory of HRGPs in poplar in order to facilitate and guide basic and applied research on plant cell walls. The newly identified poplar HRGPs can now be examined to determine their respective structural and functional roles, including their possible applications in the areas plant biofuel and natural products for medicinal or industrial uses. Additionally, other plants whose genomes are sequenced can now be examined in a similar way using this bioinformatics program which will provide insight to the evolution of the HRGP family in the plant kingdom.

Bioinformatic Identification and Analysis of Extensins in the Plant Kingdom

Extensins (EXTs) are a family of plant cell wall hydroxyproline-rich glycoproteins (HRGPs) that are implicated to play important roles in plant growth, development, and defense. Structurally, EXTs are characterized by the repeated occurrence of serine (Ser) followed by three to five prolines (Pro) residues, which are hydroxylated as hydroxyproline (Hyp) and glycosylated. Some EXTs have Tyrosine (Tyr)-X-Tyr (where X can be any amino acid) motifs that are responsible for intramolecular or intermolecular cross-linkings. EXTs can be divided into several classes: classical EXTs, short EXTs, leucine-rich repeat extensins (LRXs), proline-rich extensin-like receptor kinases (PERKs), formin-homolog EXTs (FH EXTs), chimeric EXTs, and long chimeric EXTs. To guide future research on the EXTs and understand evolutionary history of EXTs in the plant kingdom, a bioinformatics study was conducted to identify and classify EXTs from 16 fully sequenced plant genomes, including Ostreococcus lucimarinus, Chlamydomonas reinhardtii, Volvox carteri, Klebsormidium flaccidum, Physcomitrella patens, Selaginella moellendorffii, Pinus taeda, Picea abies, Brachypodium distachyon, Zea mays, Oryza sativa, Glycine max, Medicago truncatula, Brassica rapa, Solanum lycopersicum, and Solanum tuberosum, to supplement data previously obtained from Arabidopsis thaliana and Populus trichocarpa. A total of 758 EXTs were newly identified, including 87 classical EXTs, 97 short EXTs, 61 LRXs, 75 PERKs, 54 FH EXTs, 38 long chimeric EXTs, and 346 other chimeric EXTs. Several notable findings were made: (1) classical EXTs were likely derived after the terrestrialization of plants; (2) LRXs, PERKs, and FHs were derived earlier than classical EXTs; (3) monocots have few classical EXTs; (4) Eudicots have the greatest number of classical EXTs and Tyr-X-Tyr cross-linking motifs are predominantly in classical EXTs; (5) green algae have no classical EXTs but have a number of long chimeric EXTs that are absent in embryophytes. Furthermore, phylogenetic analysis was conducted of LRXs, PERKs and FH EXTs, which shed light on the evolution of three EXT classes.

Coming of age: orphan genes in plants

Sizable minorities of protein-coding genes from every sequenced eukaryotic and prokaryotic genome are unique to the species. These so-called ‘orphan genes’ may evolve de novo from non-coding sequence or be derived from older coding material. They are often associated with environmental stress responses and species-specific traits or regulatory patterns. However, difficulties in studying genes where comparative analysis is impossible, and a bias towards broadly conserved genes, have resulted in underappreciation of their importance. We review here the identification, possible origins, evolutionary trends, and functions of orphans with an emphasis on their role in plant biology. We exemplify several evolutionary trends with an analysis of Arabidopsis thaliana and present QQS as a model orphan gene.

Efficient selection of biomineralizing DNA aptamers using deep sequencing and population clustering

DNA-based information systems drive the combinatorial optimization processes of natural evolution, including the evolution of biominerals. Advances in high-throughput DNA sequencing expand the power of DNA as a potential information platform for combinatorial engineering, but many applications remain to be developed due in part to the challenge of handling large amounts of sequence data. Here we employ high-throughput sequencing and a recently developed clustering method (AutoSOME) to identify single-stranded DNA sequence families that bind specifically to ZnO semiconductor mineral surfaces. These sequences were enriched from a diverse DNA library after a single round of screening, whereas previous screening approaches typically require 5-15 rounds of enrichment for effective sequence identification. The consensus sequence of the largest cluster was poly d(T)30. This consensus sequence exhibited clear aptamer behavior and was shown to promote the synthesis of crystalline ZnO from aqueous solution at near-neutral pH. This activity is significant, as the crystalline form of this wide-bandgap semiconductor is not typically amenable to solution synthesis in this pH range. High-resolution TEM revealed that this DNA synthesis route yields ZnO nanoparticles with an amorphous-crystalline core-shell structure, suggesting that the mechanism of mineralization involves nanoscale coacervation around the DNA template. We thus demonstrate that our new method, termed Single round Enrichment of Ligands by deep Sequencing (SEL-Seq), can facilitate biomimetic synthesis of technological nanomaterials by accelerating combinatorial selection of biomolecular-mineral interactions. Moreover, by enabling direct characterization of sequence family demographics, we anticipate that SEL-Seq will enhance aptamer discovery in applications employing additional rounds of screening.

Invited review: Caseins and the casein micelle: their biological functions, structures, and behavior in foods

A typical casein micelle contains thousands of casein molecules, most of which form thermodynamically stable complexes with nanoclusters of amorphous calcium phosphate. Like many other unfolded proteins, caseins have an actual or potential tendency to assemble into toxic amyloid fibrils, particularly at the high concentrations found in milk. Fibrils do not form in milk because an alternative aggregation pathway is followed that results in formation of the casein micelle. As a result of forming micelles, nutritious milk can be secreted and stored without causing either pathological calcification or amyloidosis of the mother's mammary tissue. The ability to sequester nanoclusters of amorphous calcium phosphate in a stable complex is not unique to caseins. It has been demonstrated using a number of noncasein secreted phosphoproteins and may be of general physiological importance in preventing calcification of other biofluids and soft tissues. Thus, competent noncasein phosphoproteins have similar patterns of phosphorylation and the same type of flexible, unfolded conformation as caseins. The ability to suppress amyloid fibril formation by forming an alternative amorphous aggregate is also not unique to caseins and underlies the action of molecular chaperones such as the small heat-shock proteins. The open structure of the protein matrix of casein micelles is fragile and easily perturbed by changes in its environment. Perturbations can cause the polypeptide chains to segregate into regions of greater and lesser density. As a result, the reliable determination of the native structure of casein micelles continues to be extremely challenging. The biological functions of caseins, such as their chaperone activity, are determined by their composition and flexible conformation and by how the casein polypeptide chains interact with each other. These same properties determine how caseins behave in the manufacture of many dairy products and how they can be used as functional ingredients in other foods.

OakContigDF159.1, a reference library for studying differential gene expression in Quercus robur during controlled biotic interactions: use for quantitative transcriptomic profiling of oak roots in ectomycorrhizal symbiosis

Oaks (Quercus spp.), which are major forest trees in the northern hemisphere, host many biotic interactions, but molecular investigation of these interactions is limited by fragmentary genome data. To date, only 75 oak expressed sequence tags (ESTs) have been characterized in ectomycorrhizal (EM) symbioses. We synthesized seven beneficial and detrimental biotic interactions between microorganisms and animals and a clone (DF159) of Quercus robur. Sixteen 454 and eight Illumina cDNA libraries from leaves and roots were prepared and merged to establish a reference for RNA-Seq transcriptomic analysis of oak EMs with Piloderma croceum. Using the Mimicking Intelligent Read Assembly (MIRA) and Trinity assembler, the OakContigDF159.1 hybrid assembly, containing 65 712 contigs with a mean length of 1003 bp, was constructed, giving broad coverage of metabolic pathways. This allowed us to identify 3018 oak contigs that were differentially expressed in EMs, with genes encoding proline-rich cell wall proteins and ethylene signalling-related transcription factors showing up-regulation while auxin and defence-related genes were down-regulated. In addition to the first report of remorin expression in EMs, the extensive coverage provided by the study permitted detection of differential regulation within large gene families (nitrogen, phosphorus and sugar transporters, aquaporins). This might indicate specific mechanisms of genome regulation in oak EMs compared with other trees.