Directed evolution of proteins by exon shuffling

A bird's-eye view: exploration of the flavin-containing monooxygenase superfamily in common wheat

The Flavin Monooxygenase (FMO) gene superfamily in plants is involved in various processes most widely documented for its involvement in auxin biosynthesis, specialized metabolite biosynthesis, and plant microbial defense signaling. The roles of FMOs in defense signaling and disease resistance have recently come into focus as they may present opportunities to increase immune responses in plants including leading to systemic acquired resistance, but are not well characterized. We present a comprehensive catalogue of FMOs found in genomes across vascular plants and explore, in depth, 170 wheat TaFMO genes for sequence architecture, cis-acting regulatory elements, and changes due to Transposable Element insertions. A molecular phylogeny separates TaFMOs into three clades (A, B, and C) for which we further report gene duplication patterns, and differential rates of homoeologue expansion and retention among TaFMO subclades. We discuss Clade B TaFMOs where gene expansion is similarly seen in other cereal genomes. Transcriptome data from various studies point towards involvement of subclade B2 TaFMOs in disease responses against both biotrophic and necrotrophic pathogens, substantiated by promoter element analysis. We hypothesize that certain TaFMOs are responsive to both abiotic and biotic stresses, providing potential targets for enhancing disease resistance, plant yield and other important agronomic traits. Altogether, FMOs in wheat and other crop plants present an untapped resource to be exploited for improving the quality of crops.

A systematic screen for co-option of transposable elements across the fungal kingdom

How novel protein functions are acquired is a central question in molecular biology. Key paths to novelty include gene duplications, recombination or horizontal acquisition. Transposable elements (TEs) are increasingly recognized as a major source of novel domain-encoding sequences. However, the impact of TE coding sequences on the evolution of the proteome remains understudied. Here, we analyzed 1237 genomes spanning the phylogenetic breadth of the fungal kingdom. We scanned proteomes for evidence of co-occurrence of TE-derived domains along with other conventional protein functional domains. We detected more than 13,000 predicted proteins containing potentially TE-derived domain, of which 825 were identified in more than five genomes, indicating that many host-TE fusions may have persisted over long evolutionary time scales. We used the phylogenetic context to identify the origin and retention of individual TE-derived domains. The most common TE-derived domains are helicases derived from Academ, Kolobok or Helitron. We found putative TE co-options at a higher rate in genomes of the Saccharomycotina, providing an unexpected source of protein novelty in these generally TE depleted genomes. We investigated in detail a candidate host-TE fusion with a heterochromatic transcriptional silencing function that may play a role in TE and gene regulation in ascomycetes. The affected gene underwent multiple full or partial losses within the phylum. Overall, our work establishes a kingdom-wide view of putative host-TE fusions and facilitates systematic investigations of candidate fusion proteins.

Mosquitoes Possess Specialized Cuticular Proteins That Are Evolutionarily Related to the Elastic Protein Resilin

Resilin is an elastic protein that is vital to insects' vigorous movement. Canonical resilin proteins possess the R&R Consensus, a chitin-binding domain conserved in a family of cuticular proteins, and highly repetitive sequences conferring elastic properties. In the malaria vector mosquito, Anopheles gambiae, however, a cuticular protein has been found that has an R&R Consensus resembling that of resilin but lacks the repetitive sequences (here, we call it resilin-related or resilin-r). The relationship between resilin-r and resilin was unclear. It was also unknown whether resilin-r is conserved in mosquitoes. In this paper, phylogenetic and structural analyses were performed to reveal the relationship of resilin homologous proteins from holometabolous insects. Their chitin-binding abilities were also assessed. A resilin-r was found in each mosquito species, and these proteins constitute a clade with resilin from other insects based on the R&R Consensus sequences, indicating an evolutionary relationship between resilin-r and resilin. The resilin-r showed chitin-binding activity as same as resilin, but had distinct structural features from resilin, suggesting that it plays specialized roles in the mosquito cuticle. Another resilin-like protein was found to exist in each holometabolous insect that possesses resilin-like repetitive sequences but lacks the R&R Consensus. These results suggest that similar evolutionary events occurred to create resilin-r and resilin-like proteins.

Mutexa: A Computational Ecosystem for Intelligent Protein Engineering

Protein engineering holds immense promise in shaping the future of biomedicine and biotechnology. This Review focuses on our ongoing development of Mutexa, a computational ecosystem designed to enable "intelligent protein engineering". In this vision, researchers will seamlessly acquire sequences of protein variants with desired functions as biocatalysts, therapeutic peptides, and diagnostic proteins through a finely-tuned computational machine, akin to Amazon Alexa's role as a versatile virtual assistant. The technical foundation of Mutexa has been established through the development of a database that combines and relates enzyme structures and their respective functions (e.g., IntEnzyDB), workflow software packages that enable high-throughput protein modeling (e.g., EnzyHTP and LassoHTP), and scoring functions that map the sequence-structure-function relationship of proteins (e.g., EnzyKR and DeepLasso). We will showcase the applications of these tools in benchmarking the convergence conditions of enzyme functional descriptors across mutants, investigating protein electrostatics and cavity distributions in SAM-dependent methyltransferases, and understanding the role of nonelectrostatic dynamic effects in enzyme catalysis. Finally, we will conclude by addressing the future steps and fundamental challenges in our endeavor to develop new Mutexa applications that assist the identification of beneficial mutants in protein engineering.

Origin and Evolution of RAS Membrane Targeting

KRAS, HRAS and NRAS proto-oncogenes belong to a family of 40 highly homologous genes, which in turn are a subset of a superfamily of >160 genes encoding small GTPases. RAS proteins consist of a globular G-domain (aa1-166) and a 22-23 aa unstructured hypervariable region (HVR) that mediates membrane targeting. The evolutionary origins of the RAS isoforms, their HVRs and alternative splicing of the KRAS locus has not been explored. We found that KRAS is basal to the RAS proto-oncogene family and its duplication generated HRAS in the common ancestor of vertebrates. In a second round of duplication HRAS generated NRAS and KRAS generated an additional RAS gene we have designated KRASBL, absent in mammals and birds. KRAS4A arose through a duplication and insertion of the 4th exon of NRAS into the 3rd intron of KRAS. We found evolutionary conservation of a short polybasic region (PBR1) in HRAS, NRAS and KRAS4A, a second polybasic region (PBR2) in KRAS4A, two neutralized basic residues (NB) and a serine in KRAS4B and KRASBL, and a modification of the CaaX motif in vertebrates with farnesyl rather than geranylgeranyl polyisoprene lipids, suggesting that a less hydrophobic membrane anchor is critical to RAS protein function. The persistence of four RAS isoforms through >400 million years of evolution argues strongly for differential function.

Identification of a novel amphioxus leucine-rich repeat receptor involved in phagocytosis reveals a role for Slit2-N-type LRR in bacterial elimination

The basal chordate amphioxus is a model for tracing the origin and evolution of vertebrate immunity. To explore the evolution of immunoreceptor signaling pathways, we searched the associated receptors of the amphioxus B. belcheri (Bb) homolog of immunoreceptor signaling adaptor protein Grb2. Mass-spectrum analysis of BbGrb2 immunoprecipitates from B. belcheri intestine lysates revealed a folate receptor (FR) domain- and leucine-rich repeat (LRR)-containing protein (FrLRR). Sequence and structural analysis showed that FrLRR is a membrane protein with a predicted curved solenoid structure. The N-terminal Fr domain contains very few folate-binding sites; the following LRR region is a Slit2-type LRR, and a GPI-anchored site was predicted at the C-terminus. RT-PCR analysis showed FrLRR is a transcription-mediated fusion gene of BbFR-like and BbSlit2-N-like genes. Genomic DNA structure analysis implied the B. belcheri FrLRR gene locus and the corresponding locus in B. floridae might be generated by exon shuffling of a Slit2-N-like gene into an FR gene. RT-qPCR, immunostaining and immunoblot results showed that FrLRR was primarily distributed in B. belcheri intestinal tissue. We further demonstrated that FrLRR localized to the cell membrane and lysosomes. Functionally, FrLRR mediated and promoted bacteria-binding and phagocytosis, and FrLRR antibody blocking or Grb2 knockdown inhibited FrLRR-mediated phagocytosis. Interestingly, we found that human Slit2-N (hSlit2-N) also mediated direct bacteria-binding and phagocytosis which was inhibited by Slit2-N antibody blocking or Grb2 knockdown. Together, these results indicate FrLRR and hSlit2-N may function as phagocytotic-receptors to promote phagocytosis through Grb2, implying the Slit2-N-type-LRR-containing proteins play a role in bacterial binding and elimination.

In-silico identification and characterization of O-methyltransferase gene family in peanut (Arachis hypogaea L.) reveals their putative roles in development and stress tolerance

Cultivated peanut (Arachis hypogaea) is a leading protein and oil-providing crop and food source in many countries. At the same time, it is affected by a number of biotic and abiotic stresses. O-methyltransferases (OMTs) play important roles in secondary metabolism, biotic and abiotic stress tolerance. However, the OMT genes have not been comprehensively analyzed in peanut. In this study, we performed a genome-wide investigation of A. hypogaea OMT genes (AhOMTs). Gene structure, motifs distribution, phylogenetic history, genome collinearity and duplication of AhOMTs were studied in detail. Promoter cis-elements, protein-protein interactions, and micro-RNAs targeting AhOMTs were also predicted. We also comprehensively studied their expression in different tissues and under different stresses. We identified 116 OMT genes in the genome of cultivated peanut. Phylogenetically, AhOMTs were divided into three groups. Tandem and segmental duplication events played a role in the evolution of AhOMTs, and purifying selection pressure drove the duplication process. AhOMT promoters were enriched in several key cis-elements involved in growth and development, hormones, light, and defense-related activities. Micro-RNAs from 12 different families targeted 35 AhOMTs. GO enrichment analysis indicated that AhOMTs are highly enriched in transferase and catalytic activities, cellular metabolic and biosynthesis processes. Transcriptome datasets revealed that AhOMTs possessed varying expression levels in different tissues and under hormones, water, and temperature stress. Expression profiling based on qRT-PCR results also supported the transcriptome results. This study provides the theoretical basis for further work on the biological roles of AhOMT genes for developmental and stress responses.

Unveiling the Impact of Gene Presence/Absence Variation in Driving Inter-Individual Sequence Diversity within the CRP-I Gene Family in Mytilus spp

Mussels (Mytilus spp.) tolerate infections much better than other species living in the same marine coastal environment thanks to a highly efficient innate immune system, which exploits a remarkable diversification of effector molecules involved in mucosal and humoral responses. Among these, antimicrobial peptides (AMPs) are subjected to massive gene presence/absence variation (PAV), endowing each individual with a potentially unique repertoire of defense molecules. The unavailability of a chromosome-scale assembly has so far prevented a comprehensive evaluation of the genomic arrangement of AMP-encoding loci, preventing an accurate ascertainment of the orthology/paralogy relationships among sequence variants. Here, we characterized the CRP-I gene cluster in the blue mussel Mytilus edulis, which includes about 50 paralogous genes and pseudogenes, mostly packed in a small genomic region within chromosome 5. We further reported the occurrence of widespread PAV within this family in the Mytilus species complex and provided evidence that CRP-I peptides likely adopt a knottin fold. We functionally characterized the synthetic peptide sCRP-I H1, assessing the presence of biological activities consistent with other knottins, revealing that mussel CRP-I peptides are unlikely to act as antimicrobial agents or protease inhibitors, even though they may be used as defense molecules against infections from eukaryotic parasites.

Convergence in determining enzyme functional descriptors across Kemp eliminase variants

Molecular simulations have been extensively employed to accelerate biocatalytic discoveries. Enzyme functional descriptors derived from molecular simulations have been leveraged to guide the search for beneficial enzyme mutants. However, the ideal active-site region size for computing the descriptors over multiple enzyme variants remains untested. Here, we conducted convergence tests for dynamics-derived and electrostatic descriptors on 18 Kemp eliminase variants across six active-site regions with various boundary distances to the substrate. The tested descriptors include the root-mean-square deviation of the active-site region, the solvent accessible surface area ratio between the substrate and active site, and the projection of the electric field (EF) on the breaking C-H bond. All descriptors were evaluated using molecular mechanics methods. To understand the effects of electronic structure, the EF was also evaluated using quantum mechanics/molecular mechanics methods. The descriptor values were computed for 18 Kemp eliminase variants. Spearman correlation matrices were used to determine the region size condition under which further expansion of the region boundary does not substantially change the ranking of descriptor values. We observed that protein dynamics-derived descriptors, including RMSD_{active_site} and SASA_ratio, converge at a distance cutoff of 5 Å from the substrate. The electrostatic descriptor, EF_C-H, converges at 6 Å using molecular mechanics methods with truncated enzyme models and 4 Å using quantum mechanics/molecular mechanics methods with whole enzyme model. This study serves as a future reference to determine descriptors for predictive modeling of enzyme engineering.

A patent-based consideration of latest platforms in the art of directed evolution: a decade long untold story

Directed (or in vitro) evolution of proteins and metabolic pathways requires tools for creating genetic diversity and identifying protein variants with new or improved functional properties. Besides simplicity, reliability, speed, versatility, universal applicability and economy of the technique, the new science of synthetic biology requires improved means for construction of smart and high-quality mutant libraries to better navigate the sequence diversity. In vitro CRISPR/Cas9-mediated mutagenic (ICM) system and machine-learning (ML)-assisted approaches to directed evolution are now in the field to achieve the goal. This review describes the gene diversification strategies, screening and selection methods, in silico (computer-aided), Cas9-mediated and ML-based approaches to mutagenesis, developed especially in the last decade, and their patent position. The objective behind is to emphasize researchers the need for noting which mutagenesis, screening or selection method is patented and then selecting a suitable restriction-free approach to sequence diversity. Techniques and evolved products subject to patent rights need commercial license if their use is for purposes other than private or experimental research.

Overexpression of enzymes in glycolysis and energy metabolic pathways to enhance coenzyme Q10 production in Rhodobacter sphaeroides VK-2-3

We subjected the components of the glycolysis and energy metabolism pathways of Rhodobacter sphaeroides (R. sphaeroides) to metabolic engineering to improve the titer and yield of coenzyme Q10 (CoQ10). Phosphofructokinase (PFK), cyclic adenylate-dependent protein kinase (PKAC), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), and adenosine triphosphate hydrolase (KdpC) were overexpressed in R. sphaeroides VK-2-3 (VK-2-3). The strains were labeled R. sphaeroides PFK (RS.PFK), RS.PKAC, RS.PFK–PKAC, RS.KdpC, RS.GAPDH, and RS.KdpC–GAPDH. Results showed that the CoQ10 titers of RS.PFK, RS.PKAC, and RS.PFK–PKAC were 300.96 ± 0.87, 405.94 ± 4.77, and 379.94 ± 0.42 mg/l, respectively. The CoQ10 titers of RS.PFK and VK-2-3 were not significantly different; however, those for RS.PKAC and RS.PFK–PKAC were 13 and 6% higher than that of VK-2-3, respectively. Further, the titers of RS.KdpC, RS.GAPDH, and RS.KdpC–GAPDH were 360.17 ± 0.39, 409.79 ± 0.76, and 359.87 ± 1.14 mg/l, respectively. The titers of RS.KdpC and RS.KdpC–GAPDH were not significantly different from that for VK-2-3, whereas that for RS.GAPDH was 14% higher than that of VK-2-3. Finally, when the cultures of RS.GAPDH and VK-2-3 were scaled up in 5-L fermenters, the CoQ10 titers and RS.GAPDH yields increased by 44.3 and 37.8%, respectively, compared with VK-2-3.To the best of our knowledge, the glycolysis pathway of R. sphaeroides was studied for the first time in this study. We genetically modified the components of the energy metabolism pathway to obtain the strain with high yield of CoQ10 mutant RS.GAPDH. The findings of this study can serve as a basis for future studies involving metabolic engineering of CoQ10-producing strains.

Evolutionary Insight into Immunothrombosis as a Healing Mechanism

Both invertebrates and vertebrates possess a cluster of immediate and local wound-sealing, pathogen-killing, and tissue healing responses known as immunoclotting and immunothrombosis, respectively, to cope with two life-threatening emergencies, namely, bleeding and microbial invasion. Despite their convergence in function, immunoclotting and immunothrombosis are deployed by different blood cells and intravascular multidomain proteins. In vertebrates, these proteins share some domains with intrinsic chemical affinities useful in generating cooperative networks such as pathogen and damage pattern recognition molecules. Moreover, many of the proteins involved in coagulation and fibrinolysis in humans are multifunctional molecules playing roles in other processes from inflammation to healing and beyond. In our modern society, however, the interaction of activated intravascular allosteric proteins with one another and with blood cells entails vulnerabilities posing a biological paradox: intravascular proteins that locally operate as tissue repair enhancers can nevertheless generate pathogenic processes by acting systemically. In this manuscript, we contextualize and frame the coagulation system and hemostasis through an evolutionary time scale, illustrating their role as dual players in the defense against exsanguination and pathogens while significantly influencing wound healing.

Single molecule, long-read Apoer2 sequencing identifies conserved and species-specific splicing patterns

Apolipoprotein E receptor 2 (Apoer2) is a synaptic receptor in the brain that binds disease-relevant ligand Apolipoprotein E (Apoe) and is highly alternatively spliced. We examined alternative splicing (AS) of conserved Apoer2 exons across vertebrate species and identified gain of exons in mammals encoding functional domains such as the cytoplasmic and furin inserts, and loss of an exon in primates encoding the eighth LDLa repeat, likely altering receptor surface levels and ligand-binding specificity. We utilized single molecule, long-read RNA sequencing to profile full-length Apoer2 isoforms and identified 68 and 48 unique full-length Apoer2 transcripts in the mouse and human cerebral cortex, respectively. Furthermore, we identified two exons encoding protein functional domains, the third EGF-precursor like repeat and glycosylation domain, that are tandemly skipped specifically in mouse. Our study provides new insight into Apoer2 isoform complexity in the vertebrate brain and highlights species-specific differences in splicing decisions that support functional diversity.

Tandem Exon Duplications Expanding the Alternative Splicing Repertoire

Tandem exon duplications play an important role in the evolution of eukaryotic genes, providing a generic mechanism for adaptive regulation of protein function. In recent studies, tandem exon duplications have been linked to mutually exclusive exon choice, a pattern of alternative splicing in which one and only one exon from a group of tandemly arranged exons is included in the mature transcript. Here, we revisit the problem of identifying tandem exon duplications in eukaryotic genomes using bioinformatic methods and show that tandemly duplicated exons are abundant not only in the coding parts, but also in the untranslated regions. We present a number of remarkable examples of tandem exon duplications, identify unannotated duplicated exons, and provide statistical support for their expression using large panels of RNA-seq experiments.

Forty Years of Directed Evolution and its Continuously Evolving Technology Toolbox - A Review of the Patent Landscape

Generating functional protein variants with novel or improved characteristics has been a goal of the biotechnology industry and life sciences, for decades. Rational design and directed evolution are two major pathways to achieve the desired ends. Whilst rational protein design approach has made substantial progress, the idea of using a method based on cycles of mutagenesis and natural selection to develop novel binding proteins, enzymes and structures has attracted great attention. Laboratory evolution of proteins/enzymes requires new tools and analytical approaches to create genetic diversity and identifying variants with desired traits. In this pursuit, construction of sufficiently large libraries of target molecules to search for improved variants and the need for new protocols to alter the properties of target molecules has been a continuing challenge in the directed evolution experiments. This review will discuss the in vivo and in vitro gene diversification tools, library screening or selection approaches, and artificial intelligence/machine-learning-based strategies to mutagenesis developed in the last forty years to accelerate the natural process of evolution in creating new functional protein variants, optimization of microbial strains and transformation of enzymes into industrial machines. Analyzing patent position over these techniques and mechanisms also constitutes an integral and distinctive part of this review. The aim is to provide an up-to-date resource/technology toolbox for research-based and pharmaceutical companies to discover the boundaries of competitor's intellectual property (IP) portfolio, their freedom-to-operate in the relevant IP landscape, and the need for patent due diligence analysis to rule out whether use of a particular patented mutagenesis method, library screening/selection technique falls outside the safe harbor of experimental use exemption. While so doing, we have referred to some recent cases that emphasize the significance of selecting a suitable gene diversification strategy in directed evolution experiments. This article is protected by copyright. All rights reserved.

Improved RAD51 binders through motif shuffling based on the modularity of BRC repeats

Exchanges of protein sequence modules support leaps in function unavailable through point mutations during evolution. Here we study the role of the two RAD51-interacting modules within the eight binding BRC repeats of BRCA2. We created 64 chimeric repeats by shuffling these modules and measured their binding to RAD51. We found that certain shuffled module combinations were stronger binders than any of the module combinations in the natural repeats. Surprisingly, the contribution from the two modules was poorly correlated with affinities of natural repeats, with a weak BRC8 repeat containing the most effective N-terminal module. The binding of the strongest chimera, BRC8-2, to RAD51 was improved by -2.4 kCal/mol compared to the strongest natural repeat, BRC4. A crystal structure of RAD51:BRC8-2 complex shows an improved interface fit and an extended β-hairpin in this repeat. BRC8-2 was shown to function in human cells, preventing the formation of nuclear RAD51 foci after ionizing radiation.

TopDomain: Exhaustive Protein Domain Boundary Metaprediction Combining Multisource Information and Deep Learning

Protein domains are independent, functional, and stable structural units of proteins. Accurate protein domain boundary prediction plays an important role in understanding protein structure and evolution, as well as for protein structure prediction. Current domain boundary prediction methods differ in terms of boundary definition, methodology, and training databases resulting in disparate performance for different proteins. We developed TopDomain, an exhaustive metapredictor, that uses deep neural networks to combine multisource information from sequence- and homology-based features of over 50 primary predictors. For this purpose, we developed a new domain boundary data set termed the TopDomain data set, in which the true annotations are informed by SCOPe annotations, structural domain parsers, human inspection, and deep learning. We benchmark TopDomain against 2484 targets with 3354 boundaries from the TopDomain test set and achieve F1 scores of 78.4% and 73.8% for multidomain boundary prediction within ±20 residues and ±10 residues of the true boundary, respectively. When examined on targets from CASP11-13 competitions, TopDomain achieves F1 scores of 47.5% and 42.8% for multidomain proteins. TopDomain significantly outperforms 15 widely used, state-of-the-art ab initio and homology-based domain boundary predictors. Finally, we implemented TopDomain_TMC, which accurately predicts whether domain parsing is necessary for the target protein.

Cathepsin F of Teladorsagia circumcincta is a recently evolved cysteine protease

Parasitic cysteine proteases are involved in parasite stage transition, invasion of host tissues, nutrient uptake, and immune evasion. The cysteine protease cathepsin F is the most abundant protein produced by fourth-stage larvae (L4) of the nematode Teladorsagia circumcincta, while its transcript is only detectable in L4 and adults. T. circumcincta cathepsin F is a recently evolved cysteine protease that does not fall clearly into either of the cathepsin L or F subfamilies. This protein exhibits characteristics of both cathepsins F and L, and its phylogenetic relationship to its closest homologs is distant, including proteins of closely related nematodes of the same subfamily.

Introns mediate post-transcriptional enhancement of nuclear gene expression in the green microalga Chlamydomonas reinhardtii

Efficient nuclear transgene expression in the green microalga Chlamydomonas reinhardtii is generally hindered by low transcription rates. Introns can increase transcript abundance by a process called Intron-Mediated Enhancement (IME) in this alga and has been broadly observed in other eukaryotes. However, the mechanisms of IME in microalgae are poorly understood. Here, we identified 33 native introns from highly expressed genes in C. reinhardtii selected from transcriptome studies as well as 13 non-native introns. We investigated their IME capacities and probed the mechanism of action by modification of splice sites, internal sequence motifs, and position within transgenes. Several introns were found to elicit strong IME and found to be broadly applicable in different expression constructs. We determined that IME in C. reinhardtii exclusively occurs from introns within transcribed ORFs regardless of the promoter and is not induced by traditional enhancers of transcription. Our results elucidate some mechanistic details of IME in C. reinhardtii, which are similar to those observed in higher plants yet underly distinctly different induction processes. Our findings narrow the focus of targets responsible for algal IME and provides evidence that introns are underestimated regulators of C. reinhardtii nuclear gene expression.

Although many genetic tools and basic transformation strategies exist for the model microalga Chlamydomonas reinhardtii, high-level genetic engineering with this organism is hindered by its inherent recalcitrance to foreign gene expression and limited knowledge of responsible expression regulators. In this work, we characterized the dynamics of 33 endogenous and 13 non-native introns and their effect on gene expression as artificial insertions into codon optimized transgenes. We found that introns from different origins have the capacity to increase gene expression rates. Intron-mediated enhancement was observed exclusively when these elements were placed in transcripts but not outside of transcribed mRNA regions. Insertion of different endogenous introns into coding sequences was found to positively affect expression rates through a synergy of additive transcription enhancement and exon length reduction, similar to those natively found in the C. reinhardtii genome. Our results indicate that intensive mRNA processing plays an underestimated role in the regulation of native gene expression in C. reinhardtii. In addition to internal sequence motifs, the location of artificially introduced introns greatly affected transgene expression levels. This work is highly valuable to the greater microalgal and synthetic biology research communities and contributes to broadening our understanding of eukaryotic intron-mediated enhancement.

Genome-wide analysis of the C3H zinc finger family reveals its functions in salt stress responses of Pyrus betulaefolia

Transcription factors regulate gene expression in response to various external and internal cues by activating or suppressing downstream genes. Significant progress has been made in identifying and characterizing the Cysteine3Histidine (C3H) gene family in several dicots and monocots. They are characterized by their signature motif of three cysteine and one histidine residues, and reportedly play important roles in regulation of plant growth, developmental processes and environmental responses. In this study, we performed genome-wide and deep analysis of putative C3H genes, and a total of 117 PbeC3H members, were identified in P. betulaefolia and classified into 12 groups. Results were supported by the gene structural characteristics and phylogenetic analysis. These genes were unevenly distributed on 17 chromosomes. The gene structures of the C3H genes were relatively complex but conserved in each group. The C3H genes experienced a WGD event that occurred in the ancestor genome of P. betulaefolia and apple before their divergence based on the synonymous substitutions (Ks) values. There were 35 and 37 pairs of paralogous genes in the P. betulaefolia and apple genome, respectively, and 87 pairs of orthologous genes between P. betulaefolia and apple were identified. Except for one orthologous pairs PbeC3H66 and MD05G1311700 which had undergone positive selection, the other C3H genes had undergone purifying selection. Expression profiles showed that high salinity stress could influence the expression level of C3H genes in P. betulaefolia. Four members were responsive to salt stress in roots, nine were responsive to salt stress in leaves and eight showed inhibited expression in leaves. Results suggested important roles of PbeC3H genes in response to salt stress and will be useful for better understanding the complex functions of the C3H genes, and will provide excellent candidates for salt-tolerance improvement.

Brain-related genes are specifically enriched with long phase 1 introns

Intronic gene regions are mostly considered in the scope of gene expression regulation, such as alternative splicing. However, relations between basic statistical properties of introns are much rarely studied in detail, despite vast available data. Particularly, little is known regarding the relationship between the intron length and the intron phase. Intron phase distribution is significantly different at different intron length thresholds. In this study, we performed GO enrichment analysis of gene sets with a particular intron phase at varying intron length thresholds using a list of 13823 orthologous human-mouse gene pairs. We found a specific group of 153 genes with phase 1 introns longer than 50 kilobases that were specifically expressed in brain, functionally related to synaptic signaling, and strongly associated with schizophrenia and other mental disorders. We propose that the prevalence of long phase 1 introns arises from the presence of the signal peptide sequence and is connected with 1–1 exon shuffling.

Targeted insertional mutagenesis libraries for deep domain insertion profiling

Domain recombination is a key principle in protein evolution and protein engineering, but inserting a donor domain into every position of a target protein is not easily experimentally accessible. Most contemporary domain insertion profiling approaches rely on DNA transposons, which are constrained by sequence bias. Here, we establish Saturated Programmable Insertion Engineering (SPINE), an unbiased, comprehensive, and targeted domain insertion library generation technique using oligo library synthesis and multi-step Golden Gate cloning. Through benchmarking to MuA transposon-mediated library generation on four ion channel genes, we demonstrate that SPINE-generated libraries are enriched for in-frame insertions, have drastically reduced sequence bias as well as near-complete and highly-redundant coverage. Unlike transposon-mediated domain insertion that was severely biased and sparse for some genes, SPINE generated high-quality libraries for all genes tested. Using the Inward Rectifier K⁺ channel Kir2.1, we validate the practical utility of SPINE by constructing and comparing domain insertion permissibility maps. SPINE is the first technology to enable saturated domain insertion profiling. SPINE could help explore the relationship between domain insertions and protein function, and how this relationship is shaped by evolutionary forces and can be engineered for biomedical applications.

Thermotogales origin scenario of eukaryogenesis

How eukaryotes were generated is an enigma of evolutionary biology. Widely accepted archaeal-origin eukaryogenesis scenarios, based on similarities of genes and related characteristics between archaea and eukaryotes, cannot explain several eukaryote-specific features of the last eukaryotic common ancestor, such as glycerol-3-phosphate-type membrane lipids, large cells and genomes, and endomembrane formation. Thermotogales spheroids, having multicopy-integrated large nucleoids and producing progeny in periplasm, may explain all of these features as well as endoplasmic reticulum-type signal cleavage sites, although they cannot divide. We hypothesize that the progeny chromosome is formed by random joining small DNAs in immature progeny, followed by reorganization by mechanisms including homologous recombination enabled with multicopy-integrated large genome (MILG). We propose that Thermotogales ancestor spheroids came to divide owing to the archaeal cell division genes horizontally transferred via virus-related particles, forming the first eukaryotic common ancestor (FECA). Referring to the hypothesis, the archaeal information-processing system would have been established in FECA by random joining DNAs excised from the MILG, which contained horizontally transferred archaeal and bacterial DNAs, followed by reorganization by the MILG-enabled homologous recombination. Thus, the large genome may have been a prerequisite, but not a consequence, of eukaryogenesis. The random joining of DNAs likely provided the basic mechanisms for eukaryotic evolution: producing the diversity by the formations of supergroups, novel genes, and introns that are involved in exon shuffling.

Genome-Wide Identification and Characterization of the NAC Transcription Factor Family in Musa Acuminata and Expression Analysis during Fruit Ripening

Banana (Musa acuminata, AAA group) is a representative climacteric fruit with essential nutrients and pleasant flavors. Control of its ripening determines both the fruit quality and the shelf life. NAC (NAM, ATAF, CUC2) proteins, as one of the largest superfamilies of transcription factors, play crucial roles in various functions, especially developmental processes. Thus, it is important to conduct a comprehensive identification and characterization of the NAC transcription factor family at the genomic level in M. acuminata. In this article, a total of 181 banana NAC genes were identified. Phylogenetic analysis indicated that NAC genes in M. acuminata, Arabidopsis, and rice were clustered into 18 groups (S1–S18), and MCScanX analysis disclosed that the evolution of MaNAC genes was promoted by segmental duplication events. Expression patterns of NAC genes during banana fruit ripening induced by ethylene were investigated using RNA-Seq data, and 10 MaNAC genes were identified as related to fruit ripening. A subcellular localization assay of selected MaNACs revealed that they were all localized to the nucleus. These results lay a good foundation for the investigation of NAC genes in banana toward the biological functions and evolution.

'Why genes in pieces?'-revisited

The alignment between the boundaries of protein domains and the boundaries of exons could provide evidence for the evolution of proteins via domain shuffling, but literature in the field has so far struggled to conclusively show this. Here, on larger data sets than previously possible, we do finally show that this phenomenon is indisputably found widely across the eukaryotic tree. In contrast, the alignment between exons and the boundaries of intrinsically disordered regions of proteins is not a general property of eukaryotes. Most interesting of all is the discovery that domain-exon alignment is much more common in recently evolved protein sequences than older ones.

Alternative splicing regulates stochastic NLRP3 activity

Leucine-rich repeat (LRR) domains are evolutionarily conserved in proteins that function in development and immunity. Here we report strict exonic modularity of LRR domains of several human gene families, which is a precondition for alternative splicing (AS). We provide evidence for AS of LRR domain within several Nod-like receptors, most prominently the inflammasome sensor NLRP3. Human NLRP3, but not mouse NLRP3, is expressed as two major isoforms, the full-length variant and a variant lacking exon 5. Moreover, NLRP3 AS is stochastically regulated, with NLRP3 ∆ exon 5 lacking the interaction surface for NEK7 and hence loss of activity. Our data thus reveals unexpected regulatory roles of AS through differential utilization of LRRs modules in vertebrate innate immunity.

Leucine-rich repeat (LRR) domains are commonly present in immune regulatory proteins. Here the authors show that LRR exonic modularity and alternative splicing of an LRR-containing protein, NLRP3, modulate the ratio of functional/afunctional NLRP3 isoforms to instill a stochastic regulation of NLRP3-mediated inflammation and innate immunity.

Differential response to Cadmium exposure by expression of a two and a three-domain metallothionein isoform in the land winkle Pomatias elegans: Valuating the marine heritage of a land snail

Through evolution, marine snails have adapted several times independently to terrestrial life. A prime example for such transitions is the adaptation to terrestrial conditions in members of the gastropod clade of Littorinoidea (Caenogastropoda). Some species of this lineage like the periwinkle (Littorina littorea), live in intertidal habitats, where they are intermittently exposed to semi-terrestrial conditions. Pomatias elegans is a close relative of Littorina littorea that has successfully colonized terrestrial habitats. Evolutionary transitions from marine to terrestrial conditions have often been fostered in marine ancestors by acquisition of physiological pre-adaptations to terrestrial life. Such pre-adaptations are based, among others, on the optimization of a wide repertoire of stress resistance mechanisms, such as the expression of metal inactivating metallothioneins (MTs). The objective of our study was to explore the Cd handling strategy in the terrestrial snail Pomatias elegans in comparison to that observed previously in Littorina littorea. After Cd exposure, the metal is accumulated mainly in the midgut gland of Pomatias elegans, in a similar way as in its marine relative. Upon Cd exposure, Pomatias elegans expresses Cd-specific MTs, as also described from Littorina littorea. In contrast to the latter species, however, the detoxification of Cd in Pomatias elegans is mediated by two different MT isoforms, one two-domain and one three-domain MT. Although the MT proteins of both species are homologous and clearly originate from one common ancestor, the three-domain MT isoform of Pomatias elegans has evolved independently from the three-domain MT of its marine counterpart, probably by addition of a third domain to the pre-existing two-domain MT. Obviously, the occurrence of homologous MT structures in both species is a hereditary character, whereas the differentiation into two distinct MT isoforms with different upregulation capacities in Pomatias elegans is an adaptive feature that probably emerged upon transition to life on land.

DNA double-strand breaks as drivers of neural genomic change, function, and disease

Early work from about two decades ago implicated DNA double-strand break (DSB) formation and repair in neuronal development. Findings emerging from recent studies of DSBs in proliferating neural progenitors and in mature, non-dividing neurons suggest important roles of DSBs in brain physiology, aging, cancer, psychiatric and neurodegenerative disorders. We provide an overview of some findings and speculate on what may lie ahead.

Challenges related to the immunogenicity of parenteral recombinant proteins: Underlying mechanisms and new approaches to overcome it

Immune response elicited by therapeutic proteins is an important safety and efficacy issue for regulatory agencies, drug manufacturers, clinicians, and patients. Administration of therapeutic proteins can potentially induce the production of anti-drug antibodies or cell-mediated immune responses. At first, it was speculated that the immunogenicity is related to the non-human origin of these proteins. Later on, it was confirmed that the human proteins may also show immunogenicity. In this review article, we will focus on a number of factors, which play crucial roles in the human protein immunogenicity. These factors are related to the patient's status (or intrinsic properties) and molecular characteristics of the therapeutic protein's (or extrinsic properties). Furthermore, we will discuss available in silico, in vitro, and in vivo methods for the prediction of sequences, which may generate an immune response following parenteral administration of these proteins. In summary, nowadays, it is possible for drug manufacturers to evaluate the risk of immunogenicity of therapeutic proteins and implement a management plan to overcome the problems prior to proceeding to human clinical trials.

Protein Engineering: Past, Present, and Future

The last decade has seen a dramatic increase in the utilization of enzymes as green and sustainable (bio)catalysts in pharmaceutical and industrial applications. This trend has to a significant degree been fueled by advances in scientists' and engineers' ability to customize native enzymes by protein engineering. A review of the literature quickly reveals the tremendous success of this approach; protein engineering has generated enzyme variants with improved catalytic activity, broadened or altered substrate specificity, as well as raised or reversed stereoselectivity. Enzymes have been tailored to retain activity at elevated temperatures and to function in the presence of organic solvents, salts and pH values far from physiological conditions. However, readers unfamiliar with the field will soon encounter the confusingly large number of experimental techniques that have been employed to accomplish these engineering feats. Herein, we use history to guide a brief overview of the major strategies for protein engineering-past, present, and future.

Genome-wide analysis of DUF221 domain-containing gene family in Oryza species and identification of its salinity stress-responsive members in rice

DUF221 domain-containing genes (DDP genes) play important roles in developmental biology, hormone signalling transduction, and responses to abiotic stress. Therefore to understand their structural and evolutionary relationship, we did a genome-wide analysis of this important gene family in rice. Further, through comparative genomics, DDP genes from Oryza sativa subsp. (indica), nine different wild species of rice and Arabidopsis were also identified. We also found an expansion of the DDP gene families in rice and Arabidopsis which is due to the segmental duplication events in some of the gene family members. In general, a highly purifying selection was found acting on all the deduced paralogous and orthologous DDP gene pairs. The data from microarray and subsequent qRT-PCR analysis revealed that although several OsDDPs were differentially regulated under salinity stress, yet OsDDP6 was upregulated at all the developmental stages in salt tolerant rice genotype, FL478. Interestingly, OsDDP6 was found to be involved in proline metabolism pathway as indicated by protein network analysis. The diverse gene structures, varied transmembrane topologies and the differential expression patterns implied the functional diversity in DDP genes. Therefore, the comprehensive evolutionary analysis of DDP genes from different Oryza species and Arabidopsis performed in this study will provide the basis for further functional validation studies vis-à-vis DDP genes of rice and other plant species.

Public antibodies to malaria antigens generated by two LAIR1 insertion modalities

We previously described two donors in whom the extracellular domain of LAIR1, a collagen-binding inhibitory receptor encoded on chromosome 191, was inserted between the V and the DJ segments of an antibody. This insertion generated, through somatic mutations, broadly reactive antibodies against RIFINs, a type of variant antigen expressed on the surface of Plasmodium falciparum-infected erythrocytes (IEs)2. To investigate how frequently such antibodies are produced in response to malaria infection, we screened plasma from two large cohorts of individuals living in malaria-endemic regions. We report that 5-10% of malaria-exposed individuals, but none of the European blood donors tested, have high levels of LAIR1-containing antibodies that dominate the response to IEs without conferring enhanced protection against febrile malaria. By analyzing the antibody-producing B cell clones at the protein, cDNA and gDNA level, we characterized additional LAIR1 insertions between the V and DJ segments and discovered a second insertion modality whereby the LAIR1 exon encoding the extracellular domain and flanking intronic sequences are inserted into the switch region. By exon shuffling, this mechanism leads to the production of bispecific antibodies in which the LAIR1 domain is precisely positioned at the elbow between the VH and CH1 domains. Additionally, in one donor the gDNA encoding the VH and CH1 domains was deleted, leading to the production of a camel-like LAIR1-containing antibody. Sequencing of the switch regions of memory B cells from European blood donors revealed frequent templated inserts originating from transcribed genes that, in rare cases, comprised exons with orientation and frame compatible with expression. Collectively, these results reveal different modalities of LAIR1 insertion that lead to public and dominant antibodies against IEs and suggest that insertion of templated DNA represents an additional mechanism of antibody diversification that can be selected in the immune response against pathogens and exploited for B cell engineering.

Genome-wide identification, functional prediction, and evolutionary analysis of the R2R3-MYB superfamily in Brassica napus

R2R3-MYB transcription factors (TFs) have been shown to play important roles in plants, including in development and in various stress conditions. Phylogenetic analysis showed the presence of 249 R2R3-MYB TFs in Brassica napus, called BnaR2R3-MYB TFs, clustered into 38 clades. BnaR2R3-MYB TFs were distributed on 19 chromosomes of B. napus. Sixteen gene clusters were identified. BnaR2R3-MYB TFs were characterized by motif prediction, gene structure analysis, and gene ontology. Evolutionary analysis revealed that BnaR2R3-MYB TFs are mainly formed as a result of whole-genome duplication. Orthologs and paralogs of BnaR2R3-MYB TFs were identified in B. napus, B. rapa, B. oleracea, and Arabidopsis thaliana using synteny-based methods. Purifying selection was pervasive within R2R3-MYB TFs. K_n/K_s values lower than 0.3 indicated that BnaR2R3-MYB TFs are being functionally converged. The role of gene conversion in the formation of BnaR2R3-MYB TFs was significant. Cis-regulatory elements in the upstream regions of BnaR2R3-MYB genes, miRNA targeting BnaR2R3MYB TFs, and post translational modifications were identified. Digital expression data revealed that BnaR2R3-MYB genes were highly expressed in the roots and under high salinity treatment after 24 h. BnaMYB21, BnaMYB141, and BnaMYB148 have been suggested for improving salt-tolerant B. napus. BnaR2R3-MYB genes were mostly up regulated on the 14th day post inoculation with Leptosphaeria biglobosa and L. maculan. BnaMYB150 is a candidate for increased tolerance to Leptospheria in B. napus.

Genome-wide analysis of the CCCH zinc finger family identifies tissue specific and stress responsive candidates in chickpea (Cicer arietinum L.)

The CCCH zinc finger is a group of proteins characterised by a typical motif consisting of three cysteine residues and one histidine residue. These proteins have been reported to play important roles in regulation of plant growth, developmental processes and environmental responses. In the present study, genome wide analysis of the CCCH zinc finger gene family was carried out in the available chickpea genome. Various bioinformatics tools were employed to predict 58 CCCH zinc finger genes in chickpea (designated CarC3H1-58), which were analysed for their physio-chemical properties. Phylogenetic analysis classified the proteins into 12 groups in which members of a particular group had similar structural organization. Further, the numbers as well as the types of CCCH motifs present in the CarC3H proteins were compared with those from Arabidopsis and Medicago truncatula. Synteny analysis revealed valuable information regarding the evolution of this gene family. Tandem and segmental duplication events were identified and their Ka/Ks values revealed that the CarC3H gene family in chickpea had undergone purifying selection. Digital, as well as real time qRT-PCR expression analysis was performed which helped in identification of several CarC3H members that expressed preferentially in specific chickpea tissues as well as during abiotic stresses (desiccation, cold, salinity). Moreover, molecular characterization of an important member CarC3H45 was carried out. This study provides comprehensive genomic information about the important CCCH zinc finger gene family in chickpea. The identified tissue specific and abiotic stress specific CCCH genes could be potential candidates for further characterization to delineate their functional roles in development and stress.

Combinatorial Library Based on Restriction Enzyme-mediated Modular Assembly

Combinatorial approaches in directed evolution were proven to be more efficient for exploring sequence space and innovating function of protein. Here, we presented the modular assembly of secondary structures (MASS) for constructing a combinatorial library. In this approach, secondary structure elements were extracted from natural existing protein. The common linkers were flanking secondary structure elements, and then secondary structure elements were digested by Hinf I restriction endonuclease that was used in the construction of combinatorial library for the first time. The digested DNA fragments were randomly ligated in the sense orientation, then in sequence to be amplified by PCR and transformation. This approach showed that different DNA fragments without homologous sequences could be randomly assembled to create significant sequence space. With the structure analysis of recombinants, it would be beneficial to the rational design, even to the design of protein de novo, and to evolve any genetic part or circuit.

Biomolecular engineering for nanobio/bionanotechnology

Biomolecular engineering can be used to purposefully manipulate biomolecules, such as peptides, proteins, nucleic acids and lipids, within the framework of the relations among their structures, functions and properties, as well as their applicability to such areas as developing novel biomaterials, biosensing, bioimaging, and clinical diagnostics and therapeutics. Nanotechnology can also be used to design and tune the sizes, shapes, properties and functionality of nanomaterials. As such, there are considerable overlaps between nanotechnology and biomolecular engineering, in that both are concerned with the structure and behavior of materials on the nanometer scale or smaller. Therefore, in combination with nanotechnology, biomolecular engineering is expected to open up new fields of nanobio/bionanotechnology and to contribute to the development of novel nanobiomaterials, nanobiodevices and nanobiosystems. This review highlights recent studies using engineered biological molecules (e.g., oligonucleotides, peptides, proteins, enzymes, polysaccharides, lipids, biological cofactors and ligands) combined with functional nanomaterials in nanobio/bionanotechnology applications, including therapeutics, diagnostics, biosensing, bioanalysis and biocatalysts. Furthermore, this review focuses on five areas of recent advances in biomolecular engineering: (a) nucleic acid engineering, (b) gene engineering, (c) protein engineering, (d) chemical and enzymatic conjugation technologies, and (e) linker engineering. Precisely engineered nanobiomaterials, nanobiodevices and nanobiosystems are anticipated to emerge as next-generation platforms for bioelectronics, biosensors, biocatalysts, molecular imaging modalities, biological actuators, and biomedical applications.

SITEX 2.0: Projections of protein functional sites on eukaryotic genes. Extension with orthologous genes

Functional sites define the diversity of protein functions and are the central object of research of the structural and functional organization of proteins. The mechanisms underlying protein functional sites emergence and their variability during evolution are distinguished by duplication, shuffling, insertion and deletion of the exons in genes. The study of the correlation between a site structure and exon structure serves as the basis for the in-depth understanding of sites organization. In this regard, the development of programming resources that allow the realization of the mutual projection of exon structure of genes and primary and tertiary structures of encoded proteins is still the actual problem. Previously, we developed the SitEx system that provides information about protein and gene sequences with mapped exon borders and protein functional sites amino acid positions. The database included information on proteins with known 3D structure. However, data with respect to orthologs was not available. Therefore, we added the projection of sites positions to the exon structures of orthologs in SitEx 2.0. We implemented a search through database using site conservation variability and site discontinuity through exon structure. Inclusion of the information on orthologs allowed to expand the possibilities of SitEx usage for solving problems regarding the analysis of the structural and functional organization of proteins. Database URL: http://www-bionet.sscc.ru/sitex/ .

miR-1207-3p regulates the androgen receptor in prostate cancer via FNDC1/fibronectin

Prostate cancer (PCa) is frequently diagnosed in men, and dysregulation of microRNAs is characteristic of many cancers. MicroRNA-1207-3p is encoded at the non-protein coding gene locus PVT1 on the 8q24 human chromosomal region, an established PCa susceptibility locus. However, the role of microRNA-1207-3p in PCa is unclear. We discovered that microRNA-1207-3p is significantly underexpressed in PCa cell lines in comparison to normal prostate epithelial cells. Increased expression of microRNA-1207-3p in PCa cells significantly inhibits proliferation, migration, and induces apoptosis via direct molecular targeting of FNDC1, a protein which contains a conserved protein domain of fibronectin (FN1). FNDC1, FN1, and the androgen receptor (AR) are significantly overexpressed in PCa cell lines and human PCa, and positively correlate with aggressive PCa. Prostate tumor FN1 expression in patients that experienced PCa-specific death is significantly higher than in patients that remained alive. Furthermore, FNDC1, FN1 and AR are concomitantly overexpressed in metastatic PCa. Consequently, these studies have revealed a novel microRNA-1207-3p/FNDC1/FN1/AR regulatory pathway in PCa.

Predictive model of thrombospondin-1 and vascular endothelial growth factor in breast tumor tissue

Angiogenesis, the formation of new blood capillaries from pre-existing vessels, is a hallmark of cancer. Thus far, strategies for reducing tumor angiogenesis have focused on inhibiting pro-angiogenic factors, and less is known about the therapeutic effects of mimicking the actions of angiogenesis inhibitors. Thrombospondin-1 (TSP1) is an important endogenous inhibitor of angiogenesis that has been investigated as an anti-angiogenic agent. TSP1 impedes the growth of new blood vessels in many ways, including crosstalk with pro-angiogenic factors. Owing to the complexity of TSP1 signaling, a predictive systems biology model would provide quantitative understanding of the angiogenic balance in tumor tissue. Therefore, we have developed a molecular-detailed, mechanistic model of TSP1 and vascular endothelial growth factor (VEGF), a promoter of angiogenesis, in breast tumor tissue. The model predicts the distribution of the angiogenic factors in tumor tissue, revealing that TSP1 is primarily in an inactive, cleaved form owing to the action of proteases, rather than bound to its cellular receptors or to VEGF. The model also predicts the effects of enhancing TSP1’s interactions with its receptors and with VEGF. To provide additional predictions that can guide the development of new anti-angiogenic drugs, we simulate administration of exogenous TSP1 mimetics that bind specific targets. The model predicts that the CD47-binding TSP1 mimetic markedly decreases the ratio of receptor-bound VEGF to receptor-bound TSP1, in favor of anti-angiogenesis. Thus, we have established a model that provides a quantitative framework to study the response to TSP1 mimetics.

Facile Construction of Random Gene Mutagenesis Library for Directed Evolution Without the Use of Restriction Enzyme in Escherichia coli

A foolproof protocol was developed for the construction of mutant DNA library for directed protein evolution. First, a library of linear mutant gene was generated by error-prone PCR or molecular shuffling, and a linear vector backbone was prepared by high-fidelity PCR. Second, the amplified insert and vector fragments were assembled by overlap-extension PCR with a pair of 5'-phosphorylated primers. Third, full-length linear plasmids with phosphorylated 5'-ends were self-ligated with T4 ligase, yielding circular plasmids encoding mutant variants suitable for high-efficiency transformation. Self-made competent Escherichia coli BL21(DE3) showed a transformation efficiency of 2.4 × 10(5) cfu/µg of the self-ligated circular plasmid. Using this method, three mutants of mCherry fluorescent protein were found to alter their colors and fluorescent intensities under visible and UV lights, respectively. Also, one mutant of 6-phosphorogluconate dehydrogenase from a thermophilic bacterium Moorella thermoacetica was found to show the 3.5-fold improved catalytic efficiency (kcat /Km ) on NAD(+) as compared to the wild-type. This protocol is DNA-sequence independent, and does not require restriction enzymes, special E. coli host, or labor-intensive optimization. In addition, this protocol can be used for subcloning the relatively long DNA sequences into any position of plasmids.

Functional Role of γ-Crystallin N in the Auditory Hindbrain

γ-crystallins are major components of the vertebrate lens but show expression in other tissues as well. Their extralenticular functions remain so far unclear. Here, we explored such roles in the rodent superior olivary complex in which previous analysis demonstrated developmentally regulated expression of Crygd, Cryge and Crygn. Immunohistochemistry with novel antibodies against Crygd/e and Crygn indicate that expression of Crygd/e was moderate and varied between the perinatal superior olivary complex of mice, rats, and gerbils. Crygn-immunoreactivity was more robust and consistently highest in the medial nucleus of the trapezoid body, but also present in other nuclei of the superior olivary complex. To analyze the function of Crygn in the auditory hindbrain, we used a Crygn allele with a floxed exon 2. Upon pairing with Egr2::Cre mice, exon 2, encoding the first two greek key motifs of Crygn, was deleted in the developing auditory hindbrain. Anatomical analysis of these mice revealed a 20% volume reduction in the medial nucleus of the trapezoid body and a 7% reduction in the lateral superior olive at postnatal day 25. This was due to cell loss between postnatal days 4 and 25, whereas cell size was unaffected. Auditory brainstem responses showed normal threshold but a significant increase in the amplitude of wave IV. Crygn is hence required for postmigratory survival and proper function of auditory hindbrain neurons. These results ascertain for the first time an essential extralenticular role for γ-crystallins in vivo.

Deletion of amelotin exons 3-6 is associated with amelogenesis imperfecta

Amelogenesis imperfecta (AI) is a heterogeneous group of genetic conditions that result in defective dental enamel formation. Amelotin (AMTN) is a secreted protein thought to act as a promoter of matrix mineralization in the final stage of enamel development, and is strongly expressed, almost exclusively, in maturation stage ameloblasts. Amtn overexpression and Amtn knockout mouse models have defective enamel with no other associated phenotypes, highlighting AMTN as an excellent candidate gene for human AI. However, no AMTN mutations have yet been associated with human AI. Using whole exome sequencing, we identified an 8,678 bp heterozygous genomic deletion encompassing exons 3-6 of AMTN in a Costa Rican family segregating dominant hypomineralised AI. The deletion corresponds to an in-frame deletion of 92 amino acids, shortening the protein from 209 to 117 residues. Exfoliated primary teeth from an affected family member had enamel that was of a lower mineral density compared to control enamel and exhibited structural defects at least some of which appeared to be associated with organic material as evidenced using elemental analysis. This study demonstrates for the first time that AMTN mutations cause non-syndromic human AI and explores the human phenotype, comparing it with that of mice with disrupted Amtn function.

Nothing in Evolution Makes Sense Except in the Light of Genomics: Read-Write Genome Evolution as an Active Biological Process

The 21st century genomics-based analysis of evolutionary variation reveals a number of novel features impossible to predict when Dobzhansky and other evolutionary biologists formulated the neo-Darwinian Modern Synthesis in the middle of the last century. These include three distinct realms of cell evolution; symbiogenetic fusions forming eukaryotic cells with multiple genome compartments; horizontal organelle, virus and DNA transfers; functional organization of proteins as systems of interacting domains subject to rapid evolution by exon shuffling and exonization; distributed genome networks integrated by mobile repetitive regulatory signals; and regulation of multicellular development by non-coding lncRNAs containing repetitive sequence components. Rather than single gene traits, all phenotypes involve coordinated activity by multiple interacting cell molecules. Genomes contain abundant and functional repetitive components in addition to the unique coding sequences envisaged in the early days of molecular biology. Combinatorial coding, plus the biochemical abilities cells possess to rearrange DNA molecules, constitute a powerful toolbox for adaptive genome rewriting. That is, cells possess "Read-Write Genomes" they alter by numerous biochemical processes capable of rapidly restructuring cellular DNA molecules. Rather than viewing genome evolution as a series of accidental modifications, we can now study it as a complex biological process of active self-modification.

Recombination of chl-fus gene (Plastid Origin) downstream of hop: a locus of chromosomal instability

Background

The co-chaperone Hop [heat shock protein (HSP) organizing protein] has been shown to act as an adaptor for protein folding and maturation, in concert with Hsp70 and Hsp90. The hop gene is of eukaryotic origin. Likewise, the chloroplast elongation factor G (cEF-G) catalyzes the translocation step in chloroplast protein synthesis. The chl-fus gene, which encodes the cEF-G protein, is of plastid origin. Both proteins, Hop and cEF-G, derived from domain duplications. It was demonstrated that the nuclear chl-fus gene locates in opposite orientation to a hop gene in Glycine max. We explored 53 available plant genomes from Chlorophyta to higher plants, to determine whether the chl-fus gene was transferred directly downstream of the primordial hop in the proto-eukaryote host cell. Since both genes came from exon/module duplication events, we wanted to explore the involvement of introns in the early origin and the ensuing evolutionary changes in gene structure.

Results

We reconstructed the evolutionary history of the two convergent plant genes, on the basis of their gene structure, microsynteny and microcolinearity, from 53 plant nuclear genomes. Despite a high degree (72 %) of microcolinearity among vascular plants, our results demonstrate that their adjacency was a product of chromosomal rearrangements. Based on predicted exon − intron structures, we inferred the molecular events giving rise to the current form of genes. Therefore, we propose a simple model of exon/module shuffling by intronic recombinations in which phase-0 introns were essential for domain duplication, and a phase-1 intron for transit peptide recruiting. Finally, we demonstrate a natural susceptibility of the intergenic region to recombine or delete, seriously threatening the integrity of the chl-fus gene for the future.

Conclusions

Our results are consistent with the interpretation that the chl-fus gene was transferred from the chloroplast to a chromosome different from that of hop, in the primitive photosynthetic eukaryote, and much later before the appearance of angiosperms, it was recombined downstream of hop. Exon/module shuffling mediated by symmetric intron phases (i.e., phase-0 introns) was essential for gene evolution. The intergenic region is prone to recombine, risking the integrity of both genes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1780-1) contains supplementary material, which is available to authorized users.

Methods for the directed evolution of proteins

Directed evolution has proved to be an effective strategy for improving or altering the activity of biomolecules for industrial, research and therapeutic applications. The evolution of proteins in the laboratory requires methods for generating genetic diversity and for identifying protein variants with desired properties. This Review describes some of the tools used to diversify genes, as well as informative examples of screening and selection methods that identify or isolate evolved proteins. We highlight recent cases in which directed evolution generated enzymatic activities and substrate specificities not known to exist in nature.