Characterization of a catalase-peroxidase variant (L333V-KatG) identified in an INH-resistant Mycobacterium tuberculosis clinical isolate

Mycobacterium tuberculosis catalase-peroxidase (Mt-KatG) is a bifunctional heme-dependent enzyme that has been shown to activate isoniazid (INH), the widely used antibiotic against tuberculosis (TB). The L333V-KatG variant has been associated with INH resistance in clinical M. tuberculosis isolates from Mexico. To understand better the mechanisms of INH activation, its catalytic properties (catalase, peroxidase, and IN-NAD formation) and crystal structure were compared with those of the wild-type enzyme (WT-KatG). The rate of IN-NAD formation mediated by WT-KatG was 23% greater than L333V-KatG when INH concentration is varied. In contrast to WT-KatG, the crystal structure of the L333V-KatG variant has a perhydroxy modification of the indole nitrogen of W107 from MYW adduct. L333V-KatG shows most of the active site residues in a similar position to WT-KatG; only R418 is in the R-conformation instead of the double R and Y conformation present in WT-KatG. L333V-KatG shows a small displacement respect to WT-KatG in the helix from R385 to L404 towards the mutation site, an increase in length of the coordination bond between H270 and heme Fe, and a longer H-bond between proximal D381 and W321, compared to WT-KatG; these small displacements could explain the altered redox potential of the heme, and result in a less active and stable enzyme.

A bioinformatics pipeline for Mycobacterium tuberculosis sequencing that cleans contaminant reads from sputum samples

Next-Generation Sequencing (NGS) is widely used to investigate genomic variation. In several studies, the genetic variation of Mycobacterium tuberculosis has been analyzed in sputum samples without previous culture, using target enrichment methodologies for NGS. Alignments obtained by different programs generally map the sequences under default parameters, and from these results, it is assumed that only Mycobacterium reads will be obtained. However, variants of interest microorganism in clinical samples can be confused with a vast collection of reads from other bacteria, viruses, and human DNA. Currently, there are no standardized pipelines, and the cleaning success is never verified since there is a lack of rigorous controls to identify and remove reads from other sputum-microorganisms genetically similar to M. tuberculosis. Therefore, we designed a bioinformatic pipeline to process NGS data from sputum samples, including several filters and quality control points to identify and eliminate non-M. tuberculosis reads to obtain a reliable genetic variant report. Our proposal uses the SURPI software as a taxonomic classifier to filter input sequences and perform a mapping that provides the highest percentage of Mycobacterium reads, minimizing the reads from other microorganisms. We then use the filtered sequences to perform variant calling with the GATK software, ensuring the mapping quality, realignment, recalibration, hard-filtering, and post-filter to increase the reliability of the reported variants. Using default mapping parameters, we identified reads of contaminant bacteria, such as Streptococcus, Rhotia, Actinomyces, and Veillonella. Our final mapping strategy allowed a sequence identity of 97.8% between the input reads and the whole M. tuberculosis reference genome H37Rv using a genomic edit distance of three, thus removing 98.8% of the off-target sequences with a Mycobacterium reads loss of 1.7%. Finally, more than 200 unreliable genetic variants were removed during the variant calling, increasing the report’s reliability.

Altered Microbial Composition of Drug-Sensitive and Drug-Resistant TB Patients Compared with Healthy Volunteers

Mycobacterium tuberculosis infection has three discernible outcomes: active tuberculosis, latent tuberculosis, or clearance of the bacterium. The outcome of the infection depends on the interaction of the bacterium, the immune system, and the microbiome of the host. The current study uses 16S rRNA sequencing to determine the diversity and composition of the respiratory microbiome of drug-resistant and drug-sensitive tuberculosis patients as well as healthy volunteers. Tuberculosis patients exhibited increased microbial diversity and differentially abundant bacteria than healthy volunteers. Compositional differences were also observed when comparing drug-sensitive or -resistant tuberculosis patients. Finally, we defined and assessed the differences in the core sputum microbiota between tuberculosis patients and healthy volunteers. Our observations collectively suggest that in sputum, Mycobacterium tuberculosis infection is related to altered bacterial diversity and compositional differences of core members of the microbiome, with potential implications for the bacterial pulmonary ecosystem's stability and function.

Sex Dimorphic Changes in Trh Gene Methylation and Thyroid-Axis Response to Energy Demands in Maternally Separated Rats

The hypothalamus-pituitary-thyroid (HPT) axis regulates energy balance through the pleiotropic action of thyroid hormones. HPT basal activity and stimulation by cold or voluntary exercise are repressed by previous chronic stress in adults. Maternal separation (MS) modifies HPT basal activity; we thus studied the response of the axis to energy demands and analyzed possible epigenetic changes on Trh promoter. Nonhandled (NH) or MS male Wistar rats were cold exposed 1 h at adulthood; Trh expression in the hypothalamic paraventricular nucleus (PVN) and serum thyrotropin (TSH) concentration were increased only in NH rats. Two weeks of voluntary exercise decreased fat mass and increased Trh expression, and thyroid hormones concentration changed proportionally to running distance in NH male rats and MS male rats. Although NH females ran more than MS and much more than males, exercise decreased body weight and fat mass only in NH rats with no change on any parameter of the HPT axis but increased Pomc expression in arcuate-nucleus of NH and Npy in MS females. Overall, the methylation pattern of PVN Trh gene promoter was similar in NH males and females; MS modified methylation of specific CpG sites, a thyroid hormone receptor (THR)-binding site present after the initiation site was hypomethylated in MS males; in MS females, the THR binding site of the proximal promoter (site 4) and 2 sites in the first intron were hypermethylated. Our studies showed that, in a sex-dimorphic manner, MS blunted the responses of HPT axis to energy demands in adult animals and caused methylation changes on Trh promoter that could alter T3 feedback.

Transporters, TBC1D4, and ARID5B Variants to Explain Glycated Hemoglobin Variability in Patients with Type 2 Diabetes

INTRODUCTION

Genetic variants could aid in predicting antidiabetic drug response by associating them with markers of glucose control, such as glycated hemoglobin (HbA1c). However, pharmacogenetic implementation for antidiabetics is still under development, as the list of actionable markers is being populated and validated. This study explores potential associations between genetic variants and plasma levels of HbA1c in 100 patients under treatment with metformin.

METHODS

HbA1c was measured in a clinical chemistry analyzer (Roche), genotyping was performed in an Illumina-GSA array and data were analyzed using PLINK. Association and prediction models were developed using R and a 10-fold cross-validation approach.

RESULTS

We identified genetic variants on SLC47A1, SLC28A1, ABCG2, TBC1D4, and ARID5B that can explain up to 55% of the interindividual variability of HbA1c plasma levels in diabetic patients under treatment. Variants on SLC47A1, SLC28A1, and ABCG2 likely impact the pharmacokinetics (PK) of metformin, while the role of the two latter can be related to insulin resistance and regulation of adipogenesis.

CONCLUSIONS

Our results confirm previous genetic associations and point to previously unassociated gene variants for metformin PK and glucose control.

Association of FAAH p.Pro129Thr and COMT p.Ala72Ser with schizophrenia and comorbid substance use through next-generation sequencing: an exploratory analysis

OBJECTIVE

Individuals with schizophrenia and substance use disorders have a poor prognosis and increased psychiatric symptoms. The present study aimed to explore the association of 106 genes in individuals with schizophrenia and comorbid substance use through a next-generation sequencing (NGS) analysis and different in silico algorithms.

METHODS

We included 105 individuals diagnosed with schizophrenia and a family history of schizophrenia, of whom 49 (46.67%) presented comorbid substance use. Using NGS, we sequenced 106 genes previously associated with schizophrenia. Logistic regression models were used to assess differences in allele frequencies, and a generalized gene-set analysis was performed at the gene level. Functional annotations were performed using different algorithms and databases.

RESULTS

We identified a total of 3,109 variants, of which 25 were associated with schizophrenia and comorbid substance use and were located in regulatory and coding regions. We found low-frequency variants in COMT p.Ala72Ser, independently of p.Val158Met, that were associated with substance use. The endocannabinoid functional variant FAAH p.Pro129Thr was also associated with substance use.

CONCLUSIONS

Genetic variants of genes related to dopaminergic and cannabinoid neurotransmitter systems were associated with comorbid substance use in schizophrenia. Nevertheless, more studies with larger sample sizes are needed to confirm our findings.

Whole genome variation in 27 Mexican indigenous populations, demographic and biomedical insights

There has been limited study of Native American whole genome diversity to date, which impairs effective implementation of personalized medicine and a detailed description of its demographic history. Here we report high coverage whole genome sequencing of 76 unrelated individuals, from 27 indigenous groups across Mexico, with more than 97% average Native American ancestry. On average, each individual has 3.26 million Single Nucleotide Variants and short indels, that together comprise a catalog of 9,737,152 variants, 44,118 of which are novel. We report 497 common Single Nucleotide Variants (with allele frequency > 5%) mapped to drug responses and 316,577 in enhancer or promoter elements; interestingly we found some of these enhancer variants in PPARG, a nuclear receptor involved in highly prevalent health problems in Mexican population, such as obesity, diabetes, and insulin resistance. By detecting signals of positive selection we report 24 enriched key pathways under selection, most of them related to immune mechanisms. No missense variants in ACE2, the receptor responsible for the entry of the SARS CoV-2 virus, were found in any individual. Population genomics and phylogenetic analyses demonstrated stratification in a Northern-Central-Southern axis, with major substructure in the Central region. The Seri, a northern group with the most genetic divergence in our study, showed a distinctive genomic context with the most novel variants, and the most population specific genotypes. Genome-wide analysis showed that the average haplotype blocks are longer in Native Mexicans than in other world populations. With this dataset we describe previously undetected population level variation in Native Mexicans, helping to reduce the gap in genomic data representation of such groups.

Pharmacogenomics: Current Actionable Variants

Pharmacogenomics (PGx), one of the several tools of precision medicine, has been slowly implemented in the clinic during the past decades. This process generally starts with direct and indirect genotype-phenotype associations of gene variants and drug efficacy, or adverse drug reactions, followed by replication and validation studies. Institutional efforts led by the PGx Research Network, The PGx Knowledge Base, and The Clinical Pharmacogenetics Implementation Consortium, mine all available data for further validation or research in additional populations. This data mining gives rise to a detailed classification of over 200 druggene pairs which, with enough documentation, may become part of a publishable guideline to aid clinicians in drug selection and dosing using genetics. The US Food and Drug Administration utilizes these guidelines to issue warnings and recommendations for specific drugs and their cautioning serves clinicians and pharmacists worldwide. Here, we aim to discuss the steps of this process and list existing actionable drug-gene pairs. Moreover, we describe the current status of PGx knowledge in populations from Mexico for actionable variants on the 19 genes listed by present PGx guidelines affecting 47 drugs. Our review collects current allele frequency information for these actionable variants, lists gaps of PGx information for relevant markers, and highlights the importance of continuing PGx research in Native and Mestizo populations.

Pharmacogenomics: Current Actionable Variants

Pharmacogenomics (PGx), one of the several tools of precision medicine, has been slowly implemented in the clinic during the past decades. This process generally starts with direct and indirect genotype-phenotype associations of gene variants and drug efficacy, or adverse drug reactions, followed by replication and validation studies. Institutional efforts led by the PGx Research Network, The PGx Knowledge Base, and The Clinical Pharmacogenetics Implementation Consortium, mine all available data for further validation or research in additional populations. This data mining gives rise to a detailed classification of over 200 druggene pairs which, with enough documentation, may become part of a publishable guideline to aid clinicians in drug selection and dosing using genetics. The US Food and Drug Administration utilizes these guidelines to issue warnings and recommendations for specific drugs and their cautioning serves clinicians and pharmacists worldwide. Here, we aim to discuss the steps of this process and list existing actionable drug-gene pairs. Moreover, we describe the current status of PGx knowledge in populations from Mexico for actionable variants on the 19 genes listed by present PGx guidelines affecting 47 drugs. Our review collects current allele frequency information for these actionable variants, lists gaps of PGx information for relevant markers, and highlights the importance of continuing PGx research in Native and Mestizo populations.

Genomic diversity of prevalent Staphylococcus epidermidis multidrug-resistant strains isolated from a Children's Hospital in México City in an eight-years survey

Staphylococcus epidermidis is a human commensal and pathogen worldwide distributed. In this work, we surveyed for multi-resistant S. epidermidis strains in eight years at a children's health-care unit in México City. Multidrug-resistant S. epidermidis were present in all years of the study, including resistance to methicillin, beta-lactams, fluoroquinolones, and macrolides. To understand the genetic basis of antibiotic resistance and its association with virulence and gene exchange, we sequenced the genomes of 17 S. epidermidis isolates. Whole-genome nucleotide identities between all the pairs of S. epidermidis strains were about 97% to 99%. We inferred a clonal structure and eight Multilocus Sequence Types (MLSTs) in the S. epidermidis sequenced collection. The profile of virulence includes genes involved in biofilm formation and phenol-soluble modulins (PSMs). Half of the S. epidermidis analyzed lacked the ica operon for biofilm formation. Likely, they are commensal S. epidermidis strains but multi-antibiotic resistant. Uneven distribution of insertion sequences, phages, and CRISPR-Cas immunity phage systems suggest frequent horizontal gene transfer. Rates of recombination between S. epidermidis strains were more prevalent than the mutation rate and affected the whole genome. Therefore, the multidrug resistance, independently of the pathogenic traits, might explain the persistence of specific highly adapted S. epidermidis clonal lineages in nosocomial settings.

Variation in Actionable Pharmacogenetic Markers in Natives and Mestizos From Mexico

The identification and characterization of pharmacogenetic variants in Latin American populations is still an ongoing endeavor. Here, we investigated SNVs on genes listed by the Pharmacogenomics Knowledge Base in 1284 Mestizos and 94 Natives from Mexico. Five institutional cohorts with NGS data were retrieved from different research projects at INMEGEN, sequencing files were filtered for 55 pharmacogenes present in all cohorts to identify novel and known variation. Bioinformatic tools VEP, PROVEAN, and FATHMM were used to assess, in silico, the functional impact of this variation. Next, we focused on 17 genes with actionable variants that have been clinically implemented. Allele frequencies were compared with major continental groups and differences discussed in the scope of a pharmacogenomic impact. We observed a wide genetic variability for known and novel SNVs, the largest variation was on UGT1A > ACE > COMT > ABCB1 and the lowest on APOE and NAT2. Although with allele frequencies around 1%, novel variation was observed in 16 of 17 PGKB genes. In Natives we identified 59 variants and 58 in Mestizos. Several genes did not show novel variation, on CYP2B6, CYP2D6, and CYP3A4 in Natives; and APOE, UGT1A, and VKORC1 in Mestizos. Similarities in allele frequency, comparing major continental groups for VIP pharmacogenes, hint towards a comparable PGx for drugs metabolized by UGT1A1, DPYD, ABCB1, CBR3, COMT, and TPMT; in contrast to variants on CYP3A5 and CYP2B6 for which significant MAF differences were identified. Our observations offer some discernment into the extent of pharmacogenetic variation registered up-to-date in Mexicans and contribute to quantitatively dissect actionable pharmacogenetic variants in Natives and Mestizos.

Whole genomic sequencing as a tool for diagnosis of drug and multidrug-resistance tuberculosis in an endemic region in Mexico

Background

Whole genome sequencing (WGS) has been proposed as a tool for diagnosing drug resistance in tuberculosis. However, reports of its effectiveness in endemic countries with important numbers of drug resistance are scarce. The goal of this study was to evaluate the effectiveness of this procedure in isolates from a tuberculosis endemic region in Mexico.

Methods

WGS analysis was performed in 81 tuberculosis positive clinical isolates with a known phenotypic profile of resistance against first-line drugs (isoniazid, rifampin, ethambutol, pyrazinamide and streptomycin). Mutations related to drug resistance were identified for each isolate; drug resistant genotypes were predicted and compared with the phenotypic profile. Genotypes and transmission clusters based on genetic distances were also characterized.

Findings

Prediction by WGS analysis of resistance against isoniazid, rifampicin, ethambutol, pyrazinamide and streptomycin showed sensitivity values of 84%, 96%, 71%, 75% and 29%, while specificity values were 100%, 94%, 90%, 90% and 98%, respectively. Prediction of multidrug resistance showed a sensitivity of 89% and specificity of 97%. Moreover, WGS analysis revealed polymorphisms related to second-line drug resistance, enabling classification of eight and two clinical isolates as pre- and extreme drug-resistant cases, respectively. Lastly, four lineages were identified in the population (L1, L2, L3 and L4). The most frequent of these was L4, which included 90% (77) of the isolates. Six transmission clusters were identified; the most frequent was TC6, which included 13 isolates with a L4.1.1 and a predominantly multidrug-resistant condition.

Conclusions

The results illustrate the utility of WGS for establishing the potential for prediction of resistance against first and second line drugs in isolates of tuberculosis from the region. They also demonstrate the feasibility of this procedure for use as a tool to support the epidemiological surveillance of drug- and multidrug-resistant tuberculosis.

Exome sequencing of 20,791 cases of type 2 diabetes and 24,440 controls

Protein-coding genetic variants that strongly affect disease risk can yield relevant clues to disease pathogenesis. Here we report exome-sequencing analyses of 20,791 individuals with type 2 diabetes (T2D) and 24,440 non-diabetic control participants from 5 ancestries. We identify gene-level associations of rare variants (with minor allele frequencies of less than 0.5%) in 4 genes at exome-wide significance, including a series of more than 30 SLC30A8 alleles that conveys protection against T2D, and in 12 gene sets, including those corresponding to T2D drug targets (P = 6.1 × 10-3) and candidate genes from knockout mice (P = 5.2 × 10-3). Within our study, the strongest T2D gene-level signals for rare variants explain at most 25% of the heritability of the strongest common single-variant signals, and the gene-level effect sizes of the rare variants that we observed in established T2D drug targets will require 75,000-185,000 sequenced cases to achieve exome-wide significance. We propose a method to interpret these modest rare-variant associations and to incorporate these associations into future target or gene prioritization efforts.

Microbiome-MX 2018: microbiota and microbiome opportunities in Mexico, a megadiverse country

A weekly conference series paired with lectures entitled "Microbiome-MX: exploring the Microbiota and Microbiome Research in Mexico" was organized to provide a multidisciplinary overview of the most recent research done in Mexico using high-throughput sequencing. Scientists and postgraduate students from several disciplines such as microbiology, bioinformatics, virology, immunology, nutrition, and medical genomics gathered to discuss state of the art in each of their respective subjects of expertise, as well as advances, applications and new opportunities on microbiota/microbiome research. In particular, high-throughput sequencing is a crucial tool to understand the challenges of a megadiverse developing country as Mexico, and moreover to know the scientific capital and capabilities available for collaboration. The conference series addressed three main topics important for Mexico: i) the complex role of microbiota in health and prevalent diseases such as obesity, diabetes, inflammatory bowel disease, tuberculosis, HIV, autoimmune diseases and gastric cancer; ii) the use of local, traditional and prehispanic products as pre/probiotics to modulate the microbiota and improve human health; and iii) the impact of the microbiota in shaping the biodiversity of economically important terrestrial and marine ecosystems. Herein, we summarize the contributions that Mexican microbiota/microbiome research is making to the global trends, describing the highlights of the conferences and lectures, rather than a review of the state-of-the-art of this research. This meeting report also presents the efforts of a multidisciplinary group of scientist to encourage collaborations and bringing this research field closer for younger generations.

EXPLORATORY ANALYSIS OF RARE AND NOVEL VARIANTS IN MEXICAN PATIENTS DIAGNOSED WITH SCHIZOPHRENIA AND DEMENTIA

BACKGROUND

Schizophrenia (SCZ) and dementia, often related, are two of the most common neuropsychiatric diseases; epidemiological studies have shown that SCZ patients present a 2-fold increased risk for dementia compared to non-schizophrenic individuals. We explored the presence of rare and novel damaging gene variants in patients diagnosed with late-onset dementia of Alzheimer's type (DAT) or SCZ.

METHODS

We included 7 DAT and 12 SCZ patients and performed high-depth targeted sequencing of 184 genes.

RESULTS

We found novel and rare damaging variants in 18 genes in these Mexican patients. Carriers of these variants showed extreme phenotypes, including, treatment-resistant SCZ or cognitive decline. Furthermore, we found a variation on ABCC1 as a possible link between psychosis and cognitive impairment.

DISCUSSION

As an exploratory analysis, we report some interesting variations that should be corroborated in larger sample size studies.

Mexican Carriers of the HNF1A p.E508K Variant Do Not Experience an Enhanced Response to Sulfonylureas

OBJECTIVE

To assess whether an ethnic-specific variant (p.E508K) in the maturity-onset diabetes of the young (MODY) gene hepatocyte nuclear factor-1α (HNF1A) found in Mexicans is associated with higher sensitivity to sulfonylureas, as documented in patients with MODY3.

RESEARCH DESIGN AND METHODS

We recruited 96 participants (46 variant carriers and 50 age- and sex-matched noncarriers). Response to glipizide (one 2.5-5.0-mg dose), metformin (four 500-mg doses), and an oral glucose challenge was evaluated using a previously validated protocol. Glucose and insulin levels and their areas under the curve (AUCs) were compared between groups.

RESULTS

Carriers of the p.E508K variant had a lower maximum insulin peak during the glipizide challenge as compared with noncarriers with diabetes (P < 0.05). Also, carriers had a lower insulin response after the oral glucose challenge. Following an oral glucose tolerance test in the presence of metformin, carriers of the p.E508K variant with diabetes had a lower maximum insulin peak and total and incremental insulin AUC value as compared with noncarriers with diabetes (P < 0.05). A similar but nonsignificant trend was seen in participants without type 2 diabetes.

CONCLUSIONS

Carriers of variant p.E508K in HNF1A have a reduced insulin response rather than the increased sensitivity to sulfonylureas seen in patients with MODY3.

Identification and characterization of Taenia solium enolase as a plasminogen-binding protein

The larval stage of Taenia solium (cysticerci) is the causal agent of human and swine cysticercosis. When ingested by the host, T. solium eggs are activated and hatch in the intestine, releasing oncospheres that migrate to various tissues and evolve into cysticerci. Plasminogen (Plg) receptor proteins have been reported to play a role in migration processes for several pathogens. This work is aimed to identify Plg-binding proteins in T. solium cysticerci and determine whether T. solium recombinant enolase (rTsEnoA) is capable of specifically binding and activating human Plg. To identify Plg-binding proteins, a 2D-SDS-PAGE ligand blotting was performed, and recognized spots were identified by MS/MS. Seven proteins from T. solium cysticerci were found capable of binding Plg: fascicilin-1, fasciclin-2, enolase, MAPK, annexin, actin, and cytosolic malate dehydrogenase. To determine whether rTsEnoA binds human Plg, a ligand blotting was performed and the results were confirmed by ELISA both in the presence and absence of εACA, a competitive Plg inhibitor. Finally, rTsEnoA-bound Plg was activated to plasmin in the presence of tPA. To better understand the evolution of enolase isoforms in T. solium, a phylogenetic inference analysis including 75 enolase amino acid sequences was conducted. The origin of flatworm enolase isoforms, except for Eno4, is independent of their vertebrate counterparts. Therefore, herein we propose to designate tapeworm protein isoforms as A, B, C, and 4. In conclusion, recombinant enolase showed a strong plasminogen binding and activating activity in vitro. T. solium enolase could play a role in parasite invasion along with other plasminogen-binding proteins.

A novel approach for human whole transcriptome analysis based on absolute gene expression of microarray data

Background

In spite of the emergence of RNA sequencing (RNA-seq), microarrays remain in widespread use for gene expression analysis in the clinic. There are over 767,000 RNA microarrays from human samples in public repositories, which are an invaluable resource for biomedical research and personalized medicine. The absolute gene expression analysis allows the transcriptome profiling of all expressed genes under a specific biological condition without the need of a reference sample. However, the background fluorescence represents a challenge to determine the absolute gene expression in microarrays. Given that the Y chromosome is absent in female subjects, we used it as a new approach for absolute gene expression analysis in which the fluorescence of the Y chromosome genes of female subjects was used as the background fluorescence for all the probes in the microarray. This fluorescence was used to establish an absolute gene expression threshold, allowing the differentiation between expressed and non-expressed genes in microarrays.

Methods

We extracted the RNA from 16 children leukocyte samples (nine males and seven females, ages 6-10 years). An Affymetrix Gene Chip Human Gene 1.0 ST Array was carried out for each sample and the fluorescence of 124 genes of the Y chromosome was used to calculate the absolute gene expression threshold. After that, several expressed and non-expressed genes according to our absolute gene expression threshold were compared against the expression obtained using real-time quantitative polymerase chain reaction (RT-qPCR).

Results

From the 124 genes of the Y chromosome, three genes (DDX3Y, TXLNG2P and EIF1AY) that displayed significant differences between sexes were used to calculate the absolute gene expression threshold. Using this threshold, we selected 13 expressed and non-expressed genes and confirmed their expression level by RT-qPCR. Then, we selected the top 5% most expressed genes and found that several KEGG pathways were significantly enriched. Interestingly, these pathways were related to the typical functions of leukocytes cells, such as antigen processing and presentation and natural killer cell mediated cytotoxicity. We also applied this method to obtain the absolute gene expression threshold in already published microarray data of liver cells, where the top 5% expressed genes showed an enrichment of typical KEGG pathways for liver cells. Our results suggest that the three selected genes of the Y chromosome can be used to calculate an absolute gene expression threshold, allowing a transcriptome profiling of microarray data without the need of an additional reference experiment.

Discussion

Our approach based on the establishment of a threshold for absolute gene expression analysis will allow a new way to analyze thousands of microarrays from public databases. This allows the study of different human diseases without the need of having additional samples for relative expression experiments.

Experimental and Theoretical Approaches To Investigate the Immunogenicity of Taenia solium-Derived KE7 Antigen

Taenia solium cysticercosis, a parasitic disease that affects human health in various regions of the world, is preventable by vaccination. Both the 97-amino-acid-long KETc7 peptide and its carboxyl-terminal, 18-amino-acid-long sequence (GK-1) are found in Taenia crassiceps Both peptides have proven protective capacity against cysticercosis and are part of the highly conserved, cestode-native, 264-amino-acid long protein KE7. KE7 belongs to a ubiquitously distributed family of proteins associated with membrane processes and may participate in several vital cell pathways. The aim of this study was to identify the T. solium KE7 (TsKE7) full-length protein and to determine its immunogenic properties. Recombinant TsKE7 (rTsKE7) was expressed in Escherichia coli Rosetta2 cells and used to obtain mouse polyclonal antibodies. Anti-rTsKE7 antibodies detected the expected native protein among the 350 spots developed from T. solium cyst vesicular fluid in a mass spectrometry-coupled immune proteomic analysis. These antibodies were then used to screen a phage-displayed 7-random-peptide library to map B-cell epitopes. The recognized phages displayed 9 peptides, with the consensus motif Y(F/Y)PS sequence, which includes YYYPS (named GK-1M, for being a GK-1 mimotope), exactly matching a part of GK-1. GK-1M was recognized by 58% of serum samples from cysticercotic pigs with 100% specificity but induced weak protection against murine cysticercosis. In silico analysis revealed a universal T-cell epitope(s) in native TsKE7 potentially capable of stimulating cytotoxic T lymphocytes and helper T lymphocytes under different major histocompatibility complex class I and class II mouse haplotypes. Altogether, these results provide a rationale for the efficacy of the KETc7, rTsKE7, and GK-1 peptides as vaccines.

Quantitative multiplexed proteomics of Taenia solium cysts obtained from the skeletal muscle and central nervous system of pigs

In human and porcine cysticercosis caused by the tapeworm Taenia solium, the larval stage (cysts) can infest several tissues including the central nervous system (CNS) and the skeletal muscles (SM). The cyst’s proteomics changes associated with the tissue localization in the host tissues have been poorly studied. Quantitative multiplexed proteomics has the power to evaluate global proteome changes in response to different conditions. Here, using a TMT-multiplexed strategy we identified and quantified over 4,200 proteins in cysts obtained from the SM and CNS of pigs, of which 891 were host proteins. To our knowledge, this is the most extensive intermixing of host and parasite proteins reported for tapeworm infections.Several antigens in cysticercosis, i.e., GP50, paramyosin and a calcium-binding protein were enriched in skeletal muscle cysts. Our results suggested the occurrence of tissue-enriched antigen that could be useful in the improvement of the immunodiagnosis for cysticercosis. Using several algorithms for epitope detection, we selected 42 highly antigenic proteins enriched for each tissue localization of the cysts. Taking into account the fold changes and the antigen/epitope contents, we selected 10 proteins and produced synthetic peptides from the best epitopes. Nine peptides were recognized by serum antibodies of cysticercotic pigs, suggesting that those peptides are antigens. Mixtures of peptides derived from SM and CNS cysts yielded better results than mixtures of peptides derived from a single tissue location, however the identification of the ‘optimal’ tissue-enriched antigens remains to be discovered. Through machine learning technologies, we determined that a reliable immunodiagnostic test for porcine cysticercosis required at least five different antigenic determinants.

Human and porcine cysticercosis caused by Taenia solium is a parasite disease still endemic in developing countries. The cysts can be located in different host tissues, including different organs of the central nervous system and the skeletal muscles. The molecular mechanisms associated with the tissue localization of the cysts are not well understood. Here, we described the proteome changes of the cysts obtained from different host tissues from infected pigs using quantitative multiplex proteomics. We explored the diversity of host proteins identified in the cyst’s protein extracts and we also explored the immune-localization of several host-related proteins within the cysts, and propose their possible function. We identified several proteins and antigens enriched for a given tissue localization. Several synthetic peptides designed from these tissue-enriched antigens were tested trough ELISA. Using a combination of peptide mixtures and machine learning technologies we were able to distinguish non cysticercotic and cysticercotic pig’s sera. The tissue-enriched proteins/antigens could be useful for the development of improved immuno-diagnostic tests capable of discriminate the tissue-localization of the cysts.

Prediction of atorvastatin plasmatic concentrations in healthy volunteers using integrated pharmacogenetics sequencing

AIM

To use variants found by next-generation sequencing to predict atorvastatin plasmatic concentration profiles (AUC) in healthy volunteers.

SUBJECTS & METHODS

A total of 60 healthy Mexican volunteers were enrolled in this study. We used variants with a predicted functional effect across 20 genes involved in atorvastatin metabolism to construct a regression model using a support vector approach with a radial basis function kernel to predict AUC refining it afterwards in order to explain a greater extent of the variance.

RESULTS

The final support vector regression model using 60 variants (including six novel variants) explained 94.52% of the variance in atorvastatin AUC.

CONCLUSION

An integrated analysis of several genes known to intervene in the different steps of metabolism is required to predict atorvastatin's AUC.

Combining metagenomics, metatranscriptomics and viromics to explore novel microbial interactions: towards a systems-level understanding of human microbiome

The advances in experimental methods and the development of high performance bioinformatic tools have substantially improved our understanding of microbial communities associated with human niches. Many studies have documented that changes in microbial abundance and composition of the human microbiome is associated with human health and diseased state. The majority of research on human microbiome is typically focused in the analysis of one level of biological information, i.e., metagenomics or metatranscriptomics. In this review, we describe some of the different experimental and bioinformatic strategies applied to analyze the 16S rRNA gene profiling and shotgun sequencing data of the human microbiome. We also discuss how some of the recent insights in the combination of metagenomics, metatranscriptomics and viromics can provide more detailed description on the interactions between microorganisms and viruses in oral and gut microbiomes. Recent studies on viromics have begun to gain importance due to the potential involvement of viruses in microbial dysbiosis. In addition, metatranscriptomic combined with metagenomic analysis have shown that a substantial fraction of microbial transcripts can be differentially regulated relative to their microbial genomic abundances. Thus, understanding the molecular interactions in the microbiome using the combination of metagenomics, metatranscriptomics and viromics is one of the main challenges towards a system level understanding of human microbiome.

Genome analysis of Excretory/Secretory proteins in Taenia solium reveals their Abundance of Antigenic Regions (AAR)

Excretory/Secretory (ES) proteins play an important role in the host-parasite interactions. Experimental identification of ES proteins is time-consuming and expensive. Alternative bioinformatics approaches are cost-effective and can be used to prioritize the experimental analysis of therapeutic targets for parasitic diseases. Here we predicted and functionally annotated the ES proteins in T. solium genome using an integration of bioinformatics tools. Additionally, we developed a novel measurement to evaluate the potential antigenicity of T. solium secretome using sequence length and number of antigenic regions of ES proteins. This measurement was formalized as the Abundance of Antigenic Regions (AAR) value. AAR value for secretome showed a similar value to that obtained for a set of experimentally determined antigenic proteins and was different to the calculated value for the non-ES proteins of T. solium genome. Furthermore, we calculated the AAR values for known helminth secretomes and they were similar to that obtained for T. solium. The results reveal the utility of AAR value as a novel genomic measurement to evaluate the potential antigenicity of secretomes. This comprehensive analysis of T. solium secretome provides functional information for future experimental studies, including the identification of novel ES proteins of therapeutic, diagnosis and immunological interest.

Human genetics. The genetics of Mexico recapitulates Native American substructure and affects biomedical traits

Mexico harbors great cultural and ethnic diversity, yet fine-scale patterns of human genome-wide variation from this region remain largely uncharacterized. We studied genomic variation within Mexico from over 1000 individuals representing 20 indigenous and 11 mestizo populations. We found striking genetic stratification among indigenous populations within Mexico at varying degrees of geographic isolation. Some groups were as differentiated as Europeans are from East Asians. Pre-Columbian genetic substructure is recapitulated in the indigenous ancestry of admixed mestizo individuals across the country. Furthermore, two independently phenotyped cohorts of Mexicans and Mexican Americans showed a significant association between subcontinental ancestry and lung function. Thus, accounting for fine-scale ancestry patterns is critical for medical and population genetic studies within Mexico, in Mexican-descent populations, and likely in many other populations worldwide.

Association of a low-frequency variant in HNF1A with type 2 diabetes in a Latino population

IMPORTANCE

Latino populations have one of the highest prevalences of type 2 diabetes worldwide.

OBJECTIVES

To investigate the association between rare protein-coding genetic variants and prevalence of type 2 diabetes in a large Latino population and to explore potential molecular and physiological mechanisms for the observed relationships.

DESIGN, SETTING, AND PARTICIPANTS

Whole-exome sequencing was performed on DNA samples from 3756 Mexican and US Latino individuals (1794 with type 2 diabetes and 1962 without diabetes) recruited from 1993 to 2013. One variant was further tested for allele frequency and association with type 2 diabetes in large multiethnic data sets of 14,276 participants and characterized in experimental assays.

MAIN OUTCOME AND MEASURES

Prevalence of type 2 diabetes. Secondary outcomes included age of onset, body mass index, and effect on protein function.

RESULTS

A single rare missense variant (c.1522G>A [p.E508K]) was associated with type 2 diabetes prevalence (odds ratio [OR], 5.48; 95% CI, 2.83-10.61; P = 4.4 × 10(-7)) in hepatocyte nuclear factor 1-α (HNF1A), the gene responsible for maturity onset diabetes of the young type 3 (MODY3). This variant was observed in 0.36% of participants without type 2 diabetes and 2.1% of participants with it. In multiethnic replication data sets, the p.E508K variant was seen only in Latino patients (n = 1443 with type 2 diabetes and 1673 without it) and was associated with type 2 diabetes (OR, 4.16; 95% CI, 1.75-9.92; P = .0013). In experimental assays, HNF-1A protein encoding the p.E508K mutant demonstrated reduced transactivation activity of its target promoter compared with a wild-type protein. In our data, carriers and noncarriers of the p.E508K mutation with type 2 diabetes had no significant differences in compared clinical characteristics, including age at onset. The mean (SD) age for carriers was 45.3 years (11.2) vs 47.5 years (11.5) for noncarriers (P = .49) and the mean (SD) BMI for carriers was 28.2 (5.5) vs 29.3 (5.3) for noncarriers (P = .19).

CONCLUSIONS AND RELEVANCE

Using whole-exome sequencing, we identified a single low-frequency variant in the MODY3-causing gene HNF1A that is associated with type 2 diabetes in Latino populations and may affect protein function. This finding may have implications for screening and therapeutic modification in this population, but additional studies are required.

Alternative splice variants in TIM barrel proteins from human genome correlate with the structural and evolutionary modularity of this versatile protein fold

After the surprisingly low number of genes identified in the human genome, alternative splicing emerged as a major mechanism to generate protein diversity in higher eukaryotes. However, it is still not known if its prevalence along the genome evolution has contributed to the overall functional protein diversity or if it simply reflects splicing noise. The (βα)₈ barrel or TIM barrel is one of the most frequent, versatile, and ancient fold encountered among enzymes. Here, we analyze the structural modifications present in TIM barrel proteins from the human genome product of alternative splicing events. We found that 87% of all splicing events involved deletions; most of these events resulted in protein fragments that corresponded to the (βα)₂, (βα)₄, (βα)₅, (βα)₆, and (βα)₇ subdomains of TIM barrels. Because approximately 7% of all the splicing events involved internal β-strand substitutions, we decided, based on the genomic data, to design β-strand and α-helix substitutions in a well-studied TIM barrel enzyme. The biochemical characterization of one of the chimeric variants suggests that some of the splice variants in the human genome with β-strand substitutions may be evolving novel functions via either the oligomeric state or substrate specificity. We provide results of how the splice variants represent subdomains that correlate with the independently folding and evolving structural units previously reported. This work is the first to observe a link between the structural features of the barrel and a recurrent genetic mechanism. Our results suggest that it is reasonable to expect that a sizeable fraction of splice variants found in the human genome represent structurally viable functional proteins. Our data provide additional support for the hypothesis of the origin of the TIM barrel fold through the assembly of smaller subdomains. We suggest a model of how nature explores new proteins through alternative splicing as a mechanism to diversify the proteins encoded in the human genome.

Mutant Escherichia coli penicillin acylase with enhanced stability at alkaline pH

Increased stability at alkaline pH should be a valuable attribute for the utilization of penicillin acylase in bioreactors employed to convert penicillins into 6-aminopenicillanic acid, a precursor of semisynthetic penicillins. In these systems, base is added for pH control, which results in local alkaline conditions that promote enzyme inactivation. Hydrolysis and synthesis reactions are also pH dependent. Here, we report work in which the gene coding for Escherichia coli penicillin acylase was subjected to oligonucleotide-directed random mutagenesis at regions coding for amino acids predicted to be at the surface of the enzyme. The resulting mutant library, cloned in E. coli, was screened by a filter paper assay of the colonies for the presence of penicillin acylase activity with enhanced stability at alkaline pH. Characterization of one of the selected clones revealed the presence of a mutation, Trp431-Arg, which would presumably alter the surface charge of the protein. In vitro experiments demonstrated a near twofold increase in the half-life of the mutant enzyme when stored at pH 8.5 as compared with the wild-type enzyme, with a comparable specific activity at several pH values. In general, the mutant displayed increased stability toward the basic side in the pH-stability profile. (c) 1995 John Wiley & Sons, Inc.

Exploring the Structure-Function Loop Adaptability of a (β/α)(8)-Barrel Enzyme through Loop Swapping and Hinge Variability

Evolution of proteins involves sequence changes that are frequently localized at loop regions, revealing their important role in natural evolution. However, the development of strategies to understand and imitate such events constitutes a challenge to design novel enzymes in the laboratory. In this study, we show how to adapt loop swapping as semiautonomous units of functional groups in an enzyme with the (β/α)(8)-barrel and how this functional adaptation can be measured in vivo. To mimic the natural mechanism providing loop variability in antibodies, we developed an overlap PCR strategy. This includes introduction of sequence diversity at two hinge residues, which connect the new loops with the rest of the protein scaffold, and we demonstrate that this is necessary for a successful exploration of functional sequence space. This design allowed us to explore the sequence requirements to functional adaptation of each loop replacement that may not be sampled otherwise. Libraries generated following this strategy were evaluated in terms of their folding competence and their functional proficiency, an observation that was formalized as a Structure-Function Loop Adaptability value. Molecular details about the function and structure of some variants were obtained by enzyme kinetics and circular dichroism. This strategy yields functional variants that retain the original activity at higher frequencies, suggesting a new strategy for protein engineering that incorporates a more divergent sequence exploration beyond that limited to point mutations. We discuss how this approach may provide insights into the mechanism of enzyme evolution and function.

Identification and analysis of residues contained on beta --> alpha loops of the dual-substrate (beta alpha)8 phosphoribosyl isomerase A specific for its phosphoribosyl anthranilate isomerase activity

A good model to experimentally explore evolutionary hypothesis related to enzyme function is the ancient-like dual-substrate (beta alpha)(8) phosphoribosyl isomerase A (PriA), which takes part in both histidine and tryptophan biosynthesis in Streptomyces coelicolor and related organisms. In this study, we determined the Michaelis-Menten enzyme kinetics for both isomerase activities in wild-type PriA from S. coelicolor and in selected single-residue monofunctional mutants, identified after Escherichia coli in vivo complementation experiments. Structural and functional analyses of a hitherto unnoticed residue contained on the functionally important beta --> alpha loop 5, namely, Arg(139), which was postulated on structural grounds to be important for the dual-substrate specificity of PriA, is presented for the first time. Indeed, enzyme kinetics analyses done on the mutant variants PriA_Ser(81)Thr and PriA_Arg(139)Asn showed that these residues, which are contained on beta --> alpha loops and in close proximity to the N-terminal phosphate-binding site, are essential solely for the phosphoribosyl anthranilate isomerase activity of PriA. Moreover, analysis of the X-ray crystallographic structure of PriA_Arg(139)Asn elucidated at 1.95 A herein strongly implicates the occurrence of conformational changes in this beta --> alpha loop as a major structural feature related to the evolution of the dual-substrate specificity of PriA. It is suggested that PriA has evolved by tuning a fine energetic balance that allows the sufficient degree of structural flexibility needed for accommodating two topologically dissimilar substrates--within a bifunctional and thus highly constrained active site--without compromising its structural stability.

Protein design through systematic catalytic loop exchange in the (beta/alpha)8 fold

Protein engineering by directed evolution has proven effective in achieving various functional modifications, but the well-established protocols for the introduction of variability, typically limited to random point mutations, seriously restrict the scope of the approach. In an attempt to overcome this limitation, we sought to explore variant libraries with richer diversity at regions recognized as functionally important through an exchange of natural components, thus combining design with combinatorial diversity. With this approach, we expected to maintain interactions important for protein stability while directing the introduction of variability to areas important for catalysis. Our strategy consisted in loop exchange over a (beta/alpha)(8) fold. Phosphoribosylanthranilate isomerase was chosen as scaffold, and we investigated its tolerance to loop exchange by fusing variant libraries to the chloramphenicol acetyl transferase coding gene as an in vivo folding reporter. We replaced loops 2, 4, and 6 of phosphoribosylanthranilate isomerase with loops of varied types and sizes from enzymes sharing the same fold. To allow for a better structural fit, saturation mutagenesis was adopted at two amino acid positions preceding the exchanged loop. Our results showed that 30% to 90% of the generated mutants in the different libraries were folded. Some variants were selected for further characterization after removal of chloramphenicol acetyl transferase gene, and their stability was studied by circular dichroism and fluorescence spectroscopy. The sequences of 545 clones show that the introduction of variability at "hinges" connecting the loops with the scaffold exhibited a noticeable effect on the appearance of folded proteins. Also, we observed that each position accepted foreign loops of different sizes and sequences. We believe our work provides the basis of a general method of exchanging variably sized loops within the (beta/alpha)(8) fold, affording a novel starting point for the screening of novel activities as well as modest diversions from an original activity.

Growth rate of a non-fermentative Escherichia coli strain is influenced by NAD+ regeneration

By complementing a non-fermentative Escherichia coli (ldhA (-) pflB (-)) strain with the recombinant Zymomonas mobilis ethanol pathway (pdc, adhB), we evaluated the effect of different levels of enzymatic activity on growth rate demonstrating that there is a direct relationship between anaerobic growth rate and the total specific activity of pyruvate decarboxylase, which is the limiting enzyme of this specific fermentative NAD(+) regenerating pathway. This relationship was proved to be useful to establish a selection strategy based on growth rate for the analysis of lctE libraries, which encode lactate dehydrogenase from Bacillus subtilis.

The effect of amino acid deletions and substitutions in the longest loop of GFP

Background

The effect of single and multiple amino acid substitutions in the green fluorescent protein (GFP) from Aequorea victoria has been extensively explored, yielding several proteins of diverse spectral properties. However, the role of amino acid deletions in this protein -as with most proteins- is still unknown, due to the technical difficulties involved in generating combinatorial in-phase amino acid deletions on a target region.

Results

In this study, the region I129-L142 of superglo GFP (sgGFP), corresponding to the longest loop of the protein and located far away from the central chromophore, was subjected to a random amino acid deletion approach, employing an in-house recently developed mutagenesis method termed Codon-Based Random Deletion (COBARDE). Only two mutants out of 16384 possible variant proteins retained fluorescence: sgGFP-Δ I129 and sgGFP-Δ D130. Interestingly, both mutants were thermosensitive and at 30°C sgGFP-Δ D130 was more fluorescent than the parent protein. In contrast with deletions, substitutions of single amino acids from residues F131 to L142 were well tolerated. The substitution analysis revealed a particular importance of residues F131, G135, I137, L138, H140 and L142 for the stability of the protein.

Conclusion

The behavior of GFP variants with both amino acid deletions and substitutions demonstrate that this loop is playing an important structural role in GFP folding. Some of the amino acids which tolerated any substitution but no deletion are simply acting as "spacers" to localize important residues in the protein structure.

Engineering Transport Protein Function: Theoretical and Technical Considerations Using the Sugar-Transporting Phosphotransferase System of Escherichia coli as a Model System

Potential experimental approaches for developing and applying protein-engineering protocols to transmembrane transport systems are described. We specifically consider procedures designed to alter protein function. These procedures are designed for the specific purposes of (1) changing protein interaction specificities and (2) changing a protein's catalytic function. We use sugar-transporting bacterial phosphotransferase systems as model systems to illustrate the proposed approaches. These and other similar procedures are likely to prove to be of utility for biotechnological manipulation of proteins as well as for elucidating potential evolutionary pathways taken for the appearance of novel functions within a protein family.

An endangered oasis of aquatic microbial biodiversity in the Chihuahuan desert

The Cuatro Cienegas basin in the Chihuahuan desert is a system of springs, streams, and pools. These ecosystems support >70 endemic species and abundant living stromatolites and other microbial communities, representing a desert oasis of high biodiversity. Here, we combine data from molecular microbiology and geology to document the microbial biodiversity of this unique environment. Ten water samples from locations within the Cuatro Cienegas basin and two neighboring valleys as well as three samples of wet sediments were analyzed. The phylogeny of prokaryotic populations in the samples was determined by characterizing cultured organisms and by PCR amplification and sequencing of 16S rRNA genes from total community DNA. The composition of microbial communities was also assessed by determining profiles of terminal restriction site polymorphisms of 16S rRNA genes in total community DNA. There were 250 different phylotypes among the 350 cultivated strains. Ninety-eight partial 16S rRNA gene sequences were obtained and classified. The clones represented 38 unique phylotypes from ten major lineages of Bacteria and one of Archaea. Unexpectedly, 50% of the phylotypes were most closely related to marine taxa, even though these environments have not been in contact with the ocean for tens of millions of years. Furthermore, terminal restriction site polymorphism profiles and geological data suggest that the aquatic ecosystems of Cuatro Cienegas are hydrologically interconnected with adjacent valleys recently targeted for agricultural intensification. The findings underscore the conservation value of desert aquatic ecosystems and the urgent need for study and preservation of freshwater microbial communities.

The genome project of Taenia solium

We have constituted a consortium of key laboratories at the National Autonomous University of Mexico to carry out a genomic project for Taenia solium. This project will provide powerful resources for the study of taeniasis/cysticercosis, and, in conjunction with the Echinococcus granulosus and Echinococcus multilocularis genome project of expressed sequence tags (ESTs), will mark the advent of genomics for cestode parasites. Our project is planned in two consecutive stages. The first stage is being carried out to determine some basic parameters of the T. solium genome. Afterwards, we will evaluate the best strategy for the second stage, a full blown genome project. We have estimated the T. solium genome size by two different approaches: cytofluorometry on isolated cyton nuclei, as well as a probabilistic calculation based on approximately 2000 sequenced genomic clones, approximately 3000 ESTs, resulting in size estimates of 270 and 251 Mb, respectively. In terms of sequencing, our goal for the first stage is to characterize several thousand EST's (from adult worm and cysticerci cDNA libraries) and genomic clones. Results obtained so far from about 16,000 sequenced ESTs from the adult stage, show that only about 40% of the T. solium coding sequences have a previously sequenced homologue. Many of the best hits are found with mammalian genes, especially with humans. However, 1.5% of the hits lack homologues in humans, making these genes immediate candidates for investigation on pharmaco-therapy, diagnostics and vaccination. Most T. solium ESTs are related to gene regulation, and signal transduction. Other important functions are housekeeping, metabolism, cell division, cytoskeleton, proteases, vacuolar transport, hormone response, and extracellular matrix activities. Preliminary results also suggest that the genome of T. solium is not highly repetitive.

Generation of variability by in vivo recombination of halves of a (beta/alpha)8 barrel protein

Similar to what has been achieved with nucleic acids, directed evolution of proteins would be greatly facilitated by the availability of large libraries and efficient selection methods. So far, host cell transformation efficiency has been a bottleneck, practically limiting libraries to sizes less than 10(9). One way to circumvent this problem has been implemented with antibody systems, where contribution to the binding site is provided by two different polypeptides (light and heavy chains). The central concept is the construction of binary systems in which the gene from the two chains are separated by a cre-lox recombinase recognition site, packaged in a phage, and subsequently introduced, by multiple infection, into a recombinase expressing cell [Sblattero D, Bradbury A. Nat Biotechnol 2000;18(1):75-80]. Here, we describe the development of a system which applies the same concept to a single-domain enzyme, the cytoplasmic (beta/alpha)8 barrel protein phosphoribosyl anthranilate isomerase (PRAI) from E. coli. For that purpose, we identified the site at which a loop containing the recognition sequence for cre-lox recombinase could be inserted yielding a functional enzyme. We evaluated the effect of this insertion on the capability of the engineered gene to complement a trp F-E. coli strain and the efficiency of the system to recover the original sequence from an abundance of non-functional mutant genes.

Improving random mutagenesis by purification of the oligonucleotide variants

An oligonucleotide-based mutagenesis method is presented where, contrary to most classical mutagenic approaches, preselection of the variants is performed at the oligonucleotide level to avoid cloning of non-desired sequences. The method relies on the generation of differentially phosphate-protected oligonucleotides. Protection of the phosphates is accomplished by substoichiometric incorporation of an Fmoc-protected and n-propyl-protected trinucleotide phosphoramidite during ordinary oligonucleotide assembly. Instead of the alkali-labile beta-cyanoethyl group introduced in ordinary assembly, the trinucleotide introduces the alkali-stable n-propyl group. As a result, single mutants carry three ionic phosphates less than the wild-type sequence, double mutants carry six ionic phosphates less and so on. This difference in ionic ratio enables separation of the variants by conventional polyacrilamide gel electrophoresis. In the exemplified library described herein, two sub-populations containing mainly triple and quadruple mutants were selected out of five possible sub-populations.

Metabolic engineering and protein directed evolution increase the yield of L-phenylalanine synthesized from glucose in Escherichia coli

L-phenylalanine (L-Phe) is an aromatic amino acid with diverse commercial applications. Technologies for industrial microbial synthesis of L-Phe using glucose as a starting raw material currently achieve a relatively low conversion yield (Y(Phe/Glc)). The purpose of this work was to study the effect of PTS (phosphotransferase transport system) inactivation and overexpression of different versions of feedback inhibition resistant chorismate mutase-prephenate dehydratase (CM-PDT) on the yield (Y(Phe/Glc)) and productivity of L-Phe synthesized from glucose. The E. coli JM101 strain and its mutant derivative PB12 (PTS(-)Glc(+) phenotype) were used as hosts. PB12 has an inactive PTS, but is capable of transporting and phosphorylating glucose by using an alternative system constituted by galactose permease (GalP) and glucokinase activities (Glk). JM101 and PB12 were transformed with three plasmids, harboring genes that encode for a feedback inhibition resistant DAHP synthase (aroG(fbr)), transketolase (tktA) and either a truncated CM-PDT (pheA(fbr)) or its derived evolved genes (pheA(ev1) or pheA(ev2)). Resting-cells experiments with these engineered strains showed that JM101 and PB12 strains expressing either pheA(ev1) or pheA(ev2) genes produced l-Phe from glucose with Y(Phe/Glc) of 0.21 and 0.33 g/g, corresponding to 38 and 60% of the maximum theoretical yield (0.55 g/g), respectively. In addition, in both engineered strains the reached q(Phe) high levels of 40 mg/g-dcw.h. The metabolic engineering strategy followed in this work, including a strain with an inactive PTS, resulted in a positive impact over the Y(Phe/Glc), enhancing it nearly 57% compared with its PTS(+) counterpart. This is the first report wherein PTS inactivation was a successful strategy to improve the Y(Phe/Glc).

Combinatorial codon-based amino acid substitutions

Twenty Fmoc-protected trinucleotide phosphoramidites representing a complete set of codons for the natural amino acids were chemically synthesized for the first time. A pool of these reagents was incorporated into oligonucleotides at substoichiometric levels to generate two libraries of variants that randomly carry either few or many codon replacements on a region encoding nine amino acids of the bacterial enzyme TEM-1 beta-lactamase. Assembly of the libraries was performed in a completely automated mode through a simple modification of ordinary protocols. This technology eliminates codon redundancy, stop codons and enables complete exploration of sequence space for single, double and triple mutations throughout a protein region spanning several residues. Sequence analysis of many non-selected clones revealed a good incorporation of the trinucleotides, producing combinations of mutations quite different from those obtained using conventional degenerate oligonucleotides. Ceftazidime-selection experiments yielded several never before reported variants containing novel amino acid combinations in the beta-lactamase omega loop region.

Protein evolution by codon-based random deletions

A method to delete in-phase codons throughout a defined target region of a gene has been developed. This approach, named the codon-based random deletion (COBARDE) method, is able to delete complete codons in a random and combinatorial mode. Robustness, automation and fine-tuning of the mutagenesis rate are essential characteristics of the method, which is based on the assembly of oligonucleotides and on the use of two transient orthogonal protecting groups during the chemical synthesis. The performance of the method for protein function evolution was demonstrated by changing the substrate specificity of TEM-1 beta-lactamase. Functional ceftazidime-resistant beta-lactamase variants containing several deleted residues inside the catalytically important omega-loop region were found. The results show that the COBARDE method is a useful new molecular tool to access previously unexplorable sequence space.

Production of a fully functional, permuted single-chain penicillin G acylase

Penicillin G acylase (PGA) is a heterodimeric enzyme synthesized as a single-polypeptide precursor that undergoes an autocatalytic processing to remove an internal spacer peptide to produce the active enzyme. We constructed a single-chain PGA not dependent on autoproteolytic processing. The mature sequence of the beta-domain was expressed as the N terminus of a new polypeptide, connected by a random tetra-peptide to the alpha-domain, to afford a permuted protein. We found several active enzymes among variants differing in their linker peptides. Protein expression analysis showed that the functional single-chain variants were produced when using a Sec-dependent leader peptide, or when expressed inside the bacterial cytoplasm. Active-site titration experiments showed that the single-chain proteins displayed similar k(cat) values to the ones obtained with the wild-type enzyme. Interestingly, the single-chain proteins also displayed close to 100% of functional active sites compared to 40% to 70% functional yield usually obtained with the heterodimeric protein.

In vivo fragment complementation of a (β/α)8barrel protein: generation of variability by recombination

The high representation of the TIM barrel as a scaffold for enzymatic proteins makes it an interesting model for protein engineering. Based on previous reports of folding mechanisms of TIM barrels that suggest an independent folding unit formed by six (beta/alpha) subunits, we interrupted the gene of phosphoribosylanthranilate isomerase (PRAI) from Escherichia coli at three different positions to yield fragments with different combinations of (beta/alpha) subunits. When these constructions were expressed as polycistrons in a TrpF-E. coli strain, complementation of the function only occurred with fragments beta1-alpha4 and beta5-alpha8, demonstrating that (beta/alpha)(4) subunits are stable enough to survive in vivo conditions and to assemble to yield a functional enzyme. The expression of these fragments in a separated plasmid/phagemid system to complement the function gave a slower complementation in the TrpF-E. coli strain; this was overcome by introducing extra secondary elements to the structure that reinforce their interaction.

Alpha-amylase from Bacillus licheniformis mutants near to the catalytic site: effects on hydrolytic and transglycosylation activity

The alpha-amylase from Bacillus licheniformis is the most widely used enzyme in the starch industry owing to its hyperthermostability, converting starch to medium-sized oligosaccharides. Based on sequence alignment of homologous amylases, we found a semi-conserved sequence pattern near the active site between transglycosidic and hydrolytic amylases, which suggested that hydrophobicity may play a role in modifying the transglycosylation/hydrolysis ratio. Based on this analysis, we replaced residue Val286 by Phe and Tyr in Bacillus licheniformis alpha-amylase. Surprisingly, the two resultant mutant enzymes, Val286Phe and Val286Tyr, showed two different behaviors. Val286Tyr mutant was 5-fold more active for hydrolysis of starch than the wild-type enzyme. In contrast, the Val286Phe mutant, differing only by one hydroxyl group, was 3-fold less hydrolytic than the wild-type enzyme and apparently had a higher transglycosylation/hydrolysis ratio. These results are discussed in terms of affinity of subsites, hydrophobicity and electrostatic environment in the active site. The engineered enzyme reported here may represent an attractive alternative for the starch transformation industries as it affords direct and substantial material savings and requires no process modifications.