WUSCHEL-related homeobox family genes in rice control lateral root primordium size.

The development of a plastic root system is essential for stable crop production under variable environments. Rice plants have two types of lateral roots (LRs): S-type (short and thin) and L-type (long, thick, and capable of further branching). LR types are determined at the primordium stage, with a larger primordium size in L-types than S-types. Despite the importance of LR types for rice adaptability to variable water conditions, molecular mechanisms underlying the primordium size control of LRs are unknown. Here, we show that two WUSCHEL-related homeobox (WOX) genes have opposing roles in controlling LR primordium (LRP) size in rice. Root tip excision on seminal roots induced L-type LR formation with wider primordia formed from an early developmental stage. QHB/OsWOX5 was isolated as a causative gene of a mutant that is defective in S-type LR formation but produces more L-type LRs than wild-type (WT) plants following root tip excision. A transcriptome analysis revealed that OsWOX10 is highly up-regulated in L-type LRPs. OsWOX10 overexpression in LRPs increased the LR diameter in an expression-dependent manner. Conversely, the mutation in OsWOX10 decreased the L-type LR diameter under mild drought conditions. The qhb mutants had higher OsWOX10 expression than WT after root tip excision. A yeast one-hybrid assay revealed that the transcriptional repressive activity of QHB was lost in qhb mutants. An electrophoresis mobility shift assay revealed that OsWOX10 is a potential target of QHB. These data suggest that QHB represses LR diameter increase, repressing OsWOX10 Our findings could help improve root system plasticity under variable environments.

Nascent MSKIK peptide cancels ribosomal stalling by arrest peptides in Escherichia coli.

The insertion of the DNA sequence encoding SKIK peptide adjacent to the M start codon of a difficult-to-express protein enhances protein production in Escherichia coli. In this report, we reveal that the increased production of the SKIK-tagged protein is not due to codon usage of the SKIK sequence. Furthermore, we found that insertion of SKIK or MSKIK just before the SecM arrest peptide (FSTPVWISQAQGIRAGP), which causes ribosomal stalling on mRNA, greatly increased the production of the protein containing the SecM arrest peptide in the E. coli-reconstituted cell-free protein synthesis system (PURE system). A similar translation enhancement phenomenon by MSKIK was observed for the CmlA leader peptide, a ribosome arrest peptide, whose arrest is induced by chloramphenicol. These results strongly suggest that the nascent MSKIK peptide prevents or releases ribosomal stalling immediately following its generation during the translation process, resulting in an increase of protein production.

Robo2 contains a cryptic binding site for neural EGFL-like (NELL) protein 1/2.

The signaling pathways that are mediated by Slit ligands and their Roundabout (Robo) family of receptors play multifunctional roles in the development of the nervous system and other organs. A recent study identified neural epidermal growth factor-like (NEL)-like 2 (NELL2) as a novel ligand for Robo3. In this study, we carried out a comprehensive analysis of the interaction between NELL1 and the Robo family of receptors and demonstrated that Robo2 contains a cryptic binding site for both NELL1 and NELL2. NELL1/2 binds to the first fibronectin type III (FNIII) domain of Robo2 but not to intact Robo2. Mutation analysis revealed that several amino acids within the first FNIII domain are critical for NELL1 binding to Robo2 but not to Robo1. The Robo2 deletion mutants without the fourth immunoglobulin domain and single amino acid substitution mutants that can influence the architecture of the ectodomain facilitated binding to NELL1/2. Acidic conditions increased the binding affinity of Robo2 for NELL1. These results suggest that Robo2 functions as a receptor for NELL1/2, particularly under circumstances where Robo2 undergoes proteolytic digestion. If this is not the case, conformational changes of the ectodomain of Robo2 may unmask the binding site for NELL1/2.

Increasing cell-free gene expression yields from linear templates in Escherichia coli and Vibrio natriegens extracts by using DNA-binding proteins.

In crude extract-based cell-free protein synthesis (CFPS), DNA templates are transcribed and translated into functional proteins. Although linear expression templates (LETs) are less laborious and expensive to generate, plasmid templates are often desired over polymerase chain reaction-generated LETs due to increased stability and protection against exonucleases present in the extract of the reaction. Here we demonstrate that addition of a double stranded DNA-binding protein to the CFPS reaction, termed single-chain Cro protein (scCro), achieves terminal protection of LETs. This CroP-LET (scCro-based protection of LET) method effectively increases superfolder green fluorescent protein (sfGFP) expression levels from LETs in Escherichia coli CFPS reactions by sixfold. Our yields are comparable to other strategies that provide chemical and enzymatic DNA stabilization in E. coli CFPS. Notably, we also report that the CroP-LET method successfully enhanced yields in CFPS platforms derived from nonmodel organisms. Our results show that CroP-LET increased sfGFP yields by 18-fold in the Vibrio natriegens CFPS platform. With the fast-expanding applications of CFPS platforms, this method provides a practical and generalizable solution to protect linear expression DNA templates.

Ultra-high-throughput analysis of functional biomolecules using in vitro selection and bioinformatics.

Functional analysis of biomolecules, including nucleic acids and proteins, is important for understanding biological mechanisms in living cells such as gene expression and metabolism. To analyze diverse biomolecular functions, large-scale screening systems for biomolecules have been developed for various applications such as to improve enzyme activity and identify target binding molecules. One of these systems, the Bead Display system, utilizes emulsion technology and is a powerful tool for rapidly screening functional nucleic acids or proteins in vitro. Furthermore, an analytical pipeline that consists of genomic systematic evolution of ligands by exponential enrichment (gSELEX)-Seq, gene expression analysis, and bioinformatics was shown to be a robust platform for comprehensively identifying genes regulated by a transcription factor. This review provides an overview of the biomolecular screening methods developed to date.

Comprehensive investigation of the gene expression system regulated by an Aspergillus oryzae transcription factor XlnR using integrated mining of gSELEX-Seq and microarray data.

BACKGROUND: Transcription factors (TFs) specifically bind to DNA sequences and control the expression of target genes. AoXlnR is a key TF involved in the expression of xylanolytic and cellulolytic enzymes in the filamentous fungi, Aspergillus oryzae. Genomic SELEX-Seq (gSELEX-Seq) can reveal the in vitro binding sites of a TF in a genome. To date, the gene expression network controlled by AoXlnR in A. oryzae is not fully explored. In this study, the data from gSELEX-Seq analysis and data mining were applied toward a comprehensive investigation of the AoXlnR-regulated transcriptional network in A. oryzae. RESULTS: Around 2000 promoters were selected as AoXlnR-binding DNAs using gSELEX-Seq, consequently identifying the genes downstream of them. On the other hand, 72 differentially expressed genes (DEGs) related to AoXlnR had been determined by microarray analysis. The intersecting set of genes, that were found using the gSELEX-Seq and the microarray analysis, had 51 genes. Further, the canonical AoXlnR-binding motifs, 5'-GGCT(A/G) A-3', were successfully identified in gSELEX-Seq. The motif numbers in each promoter of the DEGs and differential expression levels were correlated by in silico analysis. The analysis showed that the presence of both 5'-GGCTAA-3' and 5'-GGCTGA-3' motif has significantly high correlation with the differential expression levels of the genes. CONCLUSIONS: Genes regulated directly by AoXlnR were identified by integrated mining of data obtained from gSELEX-Seq and microarray. The data mining of the promoters of differentially expressed genes revealed the close relation between the presence of the AoXlnR-binding motifs and the expression levels of the downstream genes. The knowledge obtained in this study can contribute greatly to the elucidation of AoXlnR-mediated cellulose and xylan metabolic network in A. oryzae. The pipeline, which is based on integrated mining of data consisting of both in vitro characterization of the DNA-binding sites and TF phenotype, can be a robust platform for comprehensive analysis of the gene expression network via the TFs.

Leucine depletion extends the lifespans of leucine-auxotrophic fission yeast by inducing Ecl1 family genes via the transcription factor Fil1.

Many studies show that lifespans of various model organisms can be extended by limiting the quantities of nutrients that are necessary for proliferation. In Schizosaccharomyces pombe, the Ecl1 family genes have been associated with lifespan control and are necessary for cell responses to nutrient depletion, but their functions and mechanisms of action remain uncharacterized. Herein, we show that leucine depletion extends the chronological lifespan (CLS) of leucine-auxotrophic cells. Furthermore, depletion of leucine extended CLS and caused cell miniaturization and cell cycle arrest at the G1 phase, and all of these processes depended on Ecl1 family genes. Although depletion of leucine raises the expression of ecl1+ by about 100-fold in leucine-auxotrophic cells, these conditions did not affect ecl1+ expression in leucine-auxotrophic fil1 mutants that were isolated in deletion set screens using 79 mutants disrupting a transcription factor. Fil1 is a GATA-type zinc finger transcription factor that reportedly binds directly to the upstream regions of ecl1+ and ecl2+. Accordingly, we suggest that Ecl1 family genes are induced in response to environmental stresses, such as oxidative stress and heat stress, or by nutritional depletion of nitrogen or sulfur sources or the amino acid leucine. We also propose that these genes play important roles in the maintenance of cell survival until conditions that favor proliferation are restored.

A simple, real-time assay of horseradish peroxidase using biolayer interferometry.

Horseradish peroxidase (HRP) isoenzyme C1a is one of the most widely used enzymes for various analytical methods in bioscience research and medical fields. In these fields, real-time monitoring of HRP activity is highly desirable because the utility of HRP as a reporter enzyme would be expanded. In this study, we developed a simple assay system enabling real-time monitoring of HRP activity by using biolayer interferometry (BLI). The HRP activity was quantitatively detected on a BLI sensor chip by tracing a binding response of tyramide, a substrate of HRP, onto an immobilized protein. This system could be applied to analyses related to oxidase activity, as well as to the functional analysis of recombinant HRP.

Spatial arrangement of proteins using scCro-tag: application for an in situ enzymatic microbead assay.

In natural systems, various metabolic reactions are often spatially organized to increase enzyme activity and specificity. Thus, by spatially arranging enzyme molecules in synthetic systems to imitate these natural systems, it is possible to promote a high rate of enzymatic turnover. In this present study, a normal and mutant form of the scCro DNA-binding protein were shown to bind orthogonally to specific recognition sequences under appropriate conditions. Furthermore, these DNA-binding tags were used to establish an enzyme assay system based on the spatial arrangement of transglutaminase and its substrate at the molecular level. Together, the results of the present study suggest that the scCro-tag may be a powerful tool to facilitate the synthetic spatial arrangement of proteins on a DNA ligand.

Rapid and cost-effective epitope mapping using PURE ribosome display coupled with next-generation sequencing and bioinformatics.

A novel, efficient and cost-effective approach for epitope identification of an antibody has been developed using a ribosome display platform. This platform, known as PURE ribosome display, utilizes an Escherichia coli-based reconstituted cell-free protein synthesis system (PURE system). It stabilizes the mRNA-ribosome-peptide complex via a ribosome-arrest peptide sequence. This system was complemented by next-generation sequencing (NGS) and an algorithm for analyzing binding epitopes. To showcase the effectiveness of this method, selection conditions were refined using the anti-PA tag monoclonal antibody with the PA tag peptide as a model. Subsequently, a random peptide library was constructed using 10 NNK triplet oligonucleotides via the PURE ribosome display. The resulting random peptide library-ribosome-mRNA complex was selected using a commercially available anti-HA (YPYDVPDYA) tag monoclonal antibody, followed by NGS and bioinformatic analysis. Our approach successfully identified the DVPDY sequence as an epitope within the hemagglutinin amino acid sequence, which was then experimentally validated. This platform provided a valuable tool for investigating continuous epitopes in antibodies.

Construction of transcript regulation mechanism prediction models based on binding motif environment of transcription factor AoXlnR in Aspergillus oryzae.

DNA-binding transcription factors (TFs) play a central role in transcriptional regulation mechanisms, mainly through their specific binding to target sites on the genome and regulation of the expression of downstream genes. Therefore, a comprehensive analysis of the function of these TFs will lead to the understanding of various biological mechanisms. However, the functions of TFs in vivo are diverse and complicated, and the identified binding sites on the genome are not necessarily involved in the regulation of downstream gene expression. In this study, we investigated whether DNA structural information around the binding site of TFs can be used to predict the involvement of the binding site in the regulation of the expression of genes located downstream of the binding site. Specifically, we calculated the structural parameters based on the DNA shape around the DNA binding motif located upstream of the gene whose expression is directly regulated by one TF AoXlnR from Aspergillus oryzae, and showed that the presence or absence of expression regulation can be predicted from the sequence information with high accuracy ([Formula: see text]-1.0) by machine learning incorporating these parameters.

The Identification of a Target Gene of the Transcription Factor KojR and Elucidation of Its Role in Carbon Metabolism for Kojic Acid Biosynthesis in Aspergillus oryzae.

DNA-binding transcription factors are broadly characterized as proteins that bind to specific sequences within genomic DNA and modulate the expression of downstream genes. This study focused on KojR, a transcription factor involved in the metabolism of kojic acid, which is an organic acid synthesized in Aspergillus oryzae and is known for its tyrosinase-inhibitory properties. However, the regulatory mechanism underlying KojR-mediated kojic acid synthesis remains unclear. Hence, we aimed to obtain a comprehensive identification of KojR-associated genes using genomic systematic evolution of ligands by exponential enrichment with high-throughput DNA sequencing (gSELEX-Seq) and RNA-Seq. During the genome-wide exploration of KojR-binding sites via gSELEX-Seq and identification of KojR-dependent differentially expressed genes (DEGs) using RNA-Seq, we confirmed that KojR preferentially binds to 5'-CGGCTAATGCGG-3', and KojR directly regulates kojT, as was previously reported. We also observed that kojA expression, which may be controlled by KojR, was significantly reduced in a ΔkojR strain. Notably, no binding of KojR to the kojA promoter region was detected. Furthermore, certain KojR-dependent DEGs identified in the present study were associated with enzymes implicated in the carbon metabolic pathway of A. oryzae. This strongly indicates that KojR plays a central role in carbon metabolism in A. oryzae.

FcRY is a key molecule controlling maternal blood IgY transfer to yolks during egg development in avian species.

Maternal immunoglobulin transfer plays a key role in conferring passive immunity to neonates. Maternal blood immunoglobulin Y (IgY) in avian species is transported to newly-hatched chicks in two steps: 1) IgY is transported from the maternal circulation to the yolk of maturing oocytes, 2) the IgY deposited in yolk is transported to the circulation of the embryo via the yolk sac membrane. An IgY-Fc receptor, FcRY, is involved in the second step, but the mechanism of the first step is still unclear. We determined whether FcRY was also the basis for maternal blood IgY transfer to the yolk in the first step during egg development. Immunohistochemistry revealed that FcRY was expressed in the capillary endothelial cells in the internal theca layer of the ovarian follicle. Substitution of the amino acid residue in Fc region of IgY substantially changed the transport efficiency of IgY into egg yolks when intravenously-injected into laying quail; the G365A mutant had a high transport efficiency, but the Y363A mutant lacked transport ability. Binding analyses of IgY mutants to FcRY indicated that the mutant with a high transport efficiency (G365A) had a strong binding activity to FcRY; the mutants with a low transport efficiency (G365D, N408A) had a weak binding activity to FcRY. One exception, the Y363A mutant had a remarkably strong binding affinity to FcRY, with a small dissociation rate. The injection of neutralizing FcRY antibodies in laying quail markedly reduced IgY uptake into egg yolks. The neutralization also showed that FcRY was engaged in prolongation of half-life of IgY in the blood; FcRY is therefore a multifunctional receptor that controls avian immunity. The pattern of the transport of the IgY mutants from the maternal blood to the egg yolk was found to be identical to that from the fertilized egg yolk to the newly-hatched chick blood circulation, via the yolk sac membrane. FcRY is therefore a critical IgY receptor that regulates the IgY uptake from the maternal blood circulation into the yolk of avian species, further indicating that the two steps of maternal-newly-hatched IgY transfer are controlled by a single receptor.

AtMYB50 regulates root cell elongation by upregulating PECTIN METHYLESTERASE INHIBITOR 8 in Arabidopsis thaliana.

Plant root development involves multiple signal transduction pathways. Notably, phytohormones like auxin and cytokinin are well characterized for their molecular mechanisms of action. Reactive oxygen species (ROS) serve as crucial signaling molecules in controlling root development. The transcription factor, UPBEAT1 (UPB1) is responsible for maintaining ROS homeostasis at the root tip, influencing the transition from cell proliferation to differentiation. While UPB1 directly regulates peroxidase expression to control ROS homeostasis, it targets genes other than peroxidases, suggesting its involvement in root growth through non-ROS signals. Our investigation focused on the transcription factor MYB50, a direct target of UPB1, in Arabidopsis thaliana. By analyzing multiple fluorescent proteins and conducting RNA-seq and ChIP-seq, we unraveled a step in the MYB50 regulatory gene network. This analysis, in conjunction with the UPB1 regulatory network, demonstrated that MYB50 directly regulates the expression of PECTIN METHYLESTERASE INHIBITOR 8 (PMEI8). Overexpressing PMEI8, similar to the MYB50, resulted in reduced mature cell length. These findings establish MYB50 as a regulator of root growth within the UPB1 gene regulatory network. Our study presents a model involving transcriptional regulation by MYB50 in the UPB1 regulated root growth system and sheds light on cell elongation via pectin modification.

The draft genome sequence of the Japanese rhinoceros beetle Trypoxylus dichotomus septentrionalis towards an understanding of horn formation.

The Japanese rhinoceros beetle Trypoxylus dichotomus is a giant beetle with distinctive exaggerated horns present on the head and prothoracic regions of the male. T. dichotomus has been used as a research model in various fields such as evolutionary developmental biology, ecology, ethology, biomimetics, and drug discovery. In this study, de novo assembly of 615 Mb, representing 80% of the genome estimated by flow cytometry, was obtained using the 10 × Chromium platform. The scaffold N50 length of the genome assembly was 8.02 Mb, with repetitive elements predicted to comprise 49.5% of the assembly. In total, 23,987 protein-coding genes were predicted in the genome. In addition, de novo assembly of the mitochondrial genome yielded a contig of 20,217 bp. We also analyzed the transcriptome by generating 16 RNA-seq libraries from a variety of tissues of both sexes and developmental stages, which allowed us to identify 13 co-expressed gene modules. We focused on the genes related to horn formation and obtained new insights into the evolution of the gene repertoire and sexual dimorphism as exemplified by the sex-specific splicing pattern of the doublesex gene. This genomic information will be an excellent resource for further functional and evolutionary analyses, including the evolutionary origin and genetic regulation of beetle horns and the molecular mechanisms underlying sexual dimorphism.

Comprehensive analysis of the DNA-binding specificity of an Aspergillus nidulans transcription factor, AmyR, using a bead display system.

The in vitro DNA binding profile of Aspergillus nidulans transcription factor AmyR was analyzed by a novel approach employing a genetic library of beads and flow cytometry analysis. An artificial library with 22 randomized nucleotides was constructed and subjected to a protein-DNA binding reaction with MalE-tagged AmyR. DNA fragments with potential AmyR-binding sites were labeled with fluorescence-conjugated antibody to be enriched by flow cytometry through 5 rounds of successive selection. Finally, a binding motif with a single CGG triplet was obtained from DNA fragments showing weak AmyR binding, while another motif with dual CGG triplets was discovered with stronger binding fragments. An informative motif, CGGNNNTTTNTCGG, was found to exist only in the promoter region of highly AmyR-dependent genes. These results suggest that this system is a powerful tool for the rapid and comprehensive analysis of the binding preferences of transcription factors.

Frequency and Clinicopathological Characteristics of Patients With KRAS/BRAF Double-Mutant Colorectal Cancer: An In Silico Study.

KRAS and BRAF mutations are currently thought to be mutually exclusive as their co-occurrence is extremely rare. Therefore, clinicopathological and molecular characteristics of colorectal carcinoma with KRAS/BRAF double mutations are unclear. We aimed to investigate the frequency and clinicopathological characteristics of double-mutant colorectal carcinoma and its differences from KRAS/BRAF single-mutant colorectal carcinoma using bioinformatics tools. We estimated the KRAS/BRAF double mutation frequency in the whole exon and coding sequences via bioinformatic analyses of three datasets from cBioPortal. We compared the clinicopathological characteristics, microsatellite instability status, BRAF classification, and tumor mutation burden of patients harboring the double mutants with those of patients harboring KRAS or BRAF single mutations. We integrated three large datasets and found that the frequency of the KRAS/BRAF double mutation in the dataset was 1.2% (29/2347). The double mutation occurred more frequently in males, with a slightly higher occurrence in the right side of the colon. Sex, histological type, histological grade, microsatellite instability, and tumor mutation burden of the patients harboring KRAS-mutant, BRAF-mutant, and double-mutant colorectal carcinoma varied significantly. The frequency of double-mutant colorectal carcinoma was 60 times higher than that previously reported. Significantly fewer double-mutant colorectal carcinoma cases were classified as BRAF class 1 and more were classified as unknown. Our findings indicate that the biological characteristics of double-mutant tumors are different from those of single-mutant tumors.

Bioinformatics in bioscience and bioengineering: Recent advances, applications, and perspectives.

Recent advances have led to the emergence of highly comprehensive and analytical approaches, such as omics analysis and high-resolution, time-resolved bioimaging analysis. These technologies have made it possible to obtain vast data from a single measurement. Subsequently, large datasets have pioneered the data-driven approach, an alternative to the traditional hypothesis-testing system, for researchers. However, processing, interpreting, and elucidating enormous datasets is no longer possible without computation. Bioinformatics is a field that has developed over long periods, intending to understand biological phenomena using methods collected from information science and statistics, thus solving this proposed research challenge. This review presents the latest methodologies and applications in sequencing, imaging, and mass spectrometry that were developed using bioinformatics. We presented the features of individual techniques and outlines in each part, avoiding the use of complex algorithms and formulas to allow beginning researchers to understand an overview. In the section on sequencing, we focused on comparative genomic, transcriptomic, and bacterial microbiome analyses, which are frequently used as applications of next-generation sequencing. Bioinformatic methods for handling sequence data and case studies were described. In the section on imaging, we introduced the analytical methods and microscopy imaging informatics techniques used in animal cell biology and plant physiology. We introduce informatics technologies for maximizing the value of measured data, including predicting the structure of unknown molecules and untargeted analysis in the section on mass spectrometry. Finally, we discuss the future outlook of this field. We anticipate that this review will assist biologists in using bioinformatics more effectively.

Auxin Distribution in Lateral Root Primordium Development Affects the Size and Lateral Root Diameter of Rice.

Lateral roots (LRs) occupy a large part of the root system and play a central role in plant water and nutrient uptake. Monocot plants, such as rice, produce two types of LRs: the S-type (short and thin) and the L-type (long, thick, and capable of further branching). Because of the ability to produce higher-order branches, the L-type LR formation contributes to efficient root system expansion. Auxin plays a major role in regulating the root system development, but its involvement in developing different types of LRs is largely unknown. Here, we show that auxin distribution is involved in regulating LR diameter. Dynamin-related protein (DRP) genes were isolated as causative genes of the mutants with increased L-type LR number and diameter than wild-type (WT). In the drp mutants, reduced endocytic activity was detected in rice protoplast and LRs with a decreased OsPIN1b-GFP endocytosis in the protoplast. Analysis of auxin distribution using auxin-responsive promoter DR5 revealed the upregulated auxin signaling in L-type LR primordia (LRP) of the WT and the mutants. The application of polar auxin transport inhibitors enhanced the effect of exogenous auxin to increase LR diameter with upregulated auxin signaling in the basal part of LRP. Inducible repression of auxin signaling in the mOsIAA3-GR system suppressed the increase in LR diameter after root tip excision, suggesting a positive role of auxin signaling in LR diameter increase. A positive regulator of LR diameter, OsWOX10, was auxin-inducible and upregulated in the drp mutants more than the WT, and revealed as a potential target of ARF transcriptional activator. Therefore, auxin signaling upregulation in LRP, especially at the basal part, induces OsWOX10 expression, increasing LR diameter.

Comprehensive analysis of transglutaminase substrate preference by cDNA display coupled with next-generation sequencing and bioinformatics.

cDNA display is an in vitro display technology based on a covalent linkage between a protein and its corresponding mRNA/cDNA, widely used for the selection of proteins and peptides from large libraries (1012) in a high throughput manner, based on their binding affinity. Here, we developed a platform using cDNA display and next-generation sequencing (NGS) for rapid and comprehensive substrate profiling of transglutaminase 2 (TG2), an enzyme crosslinking glutamine and lysine residues in proteins. After screening and selection of the control peptide library randomized at the reactive glutamine, a combinatorial library of displayed peptides randomized at positions - 1, + 1, + 2, and + 3 from the reactive glutamine was screened followed by NGS and bioinformatic analysis, which indicated a strong preference of TG2 towards peptides with glutamine at position - 1 (Gln-Gln motif), and isoleucine or valine at position + 3. The highly enriched peptides indeed contained the indicated sequence and showed a higher reactivity as TG2 substrates than the peptide previously selected by phage display, thus representing the novel candidate peptide probes for TG2 research. Furthermore, the obtained information on substrate profiling can be used to identify potential TG2 protein targets. This platform will be further used for the substrate profiling of other TG isozymes, as well as for the selection and evolution of larger biomolecules.