Immunogenomic profiles associated with response to neoadjuvant chemoradiotherapy in patients with rectal cancer.

BACKGROUND: Accumulating evidence suggests that radiotherapy success has an immune-associated component. The immunogenomic profiles associated with responses to chemoradiotherapy (CRT) were assessed in patients with locally advanced rectal cancer in this study. METHODS: CD8+ tumour-infiltrating lymphocyte (TIL) and stromal lymphocyte densities were assessed by immunohistochemistry using pretreatment biopsies from patients with advanced rectal cancer who had preoperative CRT. Whole-exome sequencing and gene expression microarray analysis were conducted to investigate the genomic properties associated with the response to CRT and CD8+ TIL density. Response to CRT was determined based on Dworak tumour regression grade (TRG); tumours with complete (TRG 4) or near-complete (TRG 3) regression were grouped as good responders, and those with TRG 1 as non-responders. RESULTS: Immunohistochemical examinations (275 patients) showed that pre-CRT CD8+ TIL density was associated with better response to CRT and improved recurrence-free survival, whereas pre-CRT stromal CD8+ cell density was not associated with either response to CRT or recurrence-free survival. Whole-exome sequencing (74 patients) showed that the numbers of single-nucleotide variations (SNVs) and neoantigens predicted from SNVs were higher in good responders than in non-responders, and these correlated positively with CD8+ TIL density (rS = 0·315 and rS = 0·334 respectively). Gene expression microarray (90 patients) showed that CD8A expression correlated positively with the expression of programmed cell death 1 (PDCD1) (rS = 0·264) and lymphocyte-activation gene 3 (LAG3) (rS = 0·507). CONCLUSION: Pre-CRT neoantigen-specific CD8+ T cell priming may be a key event in CRT responses where immune checkpoint molecules could be useful targets to enhance tumour regression.

ANTECEDENTES: Las evidencias existentes sugieren que el éxito de la radioterapia tiene un componente asociado con el sistema inmunitario. En este estudio se evaluaron los perfiles inmunogenómicos asociados con la respuesta a la quimiorradioterapia (chemoradiotherapy, CRT) en pacientes con cáncer de recto localmente avanzado. MÉTODOS: Las densidades de los linfocitos infiltrantes de tumor CD8+ (tumour-infiltrating lymphocyte, TIL) y de los linfocitos del estroma se evaluaron por inmunohistoquímicas en las biopsias antes del tratamiento de pacientes con cáncer de recto localmente avanzado que recibieron CRT preoperatoria. Se realizó secuenciación de todo el exoma, así como microarrays de expresión génica, para investigar las propiedades genómicas asociadas con la respuesta a la CRT y a la densidad de los TIL CD8+. La respuesta a la CRT se determinó según el grado de regresión del tumor de Dworak (tumour regression grade, TRG), agrupándose como buenos respondedores los casos de regresión tumoral completa (TRG4) o casi completa (TRG3) y como no respondedores, los casos de grado TRG1. RESULTADOS: Los exámenes inmunohistoquímicos (n = 275) mostraron que la densidad pre-CRT de TIL CD8+ se asoció con una mejor respuesta a la CRT y una mejor supervivencia libre de recidiva, aunque la densidad de células CD8+ del estroma previa a la CRT no se asoció con la respuesta a la CRT ni con la supervivencia libre de recidiva. La secuenciación de todo el exoma (n = 74) mostró que el número de variaciones de nucleótidos únicos (single nucleotide variations, SNVs) y los neoantígenos predichos a partir de los SNVs fueron mayores en los que respondieron bien que en los que no respondieron, y éstos se correlacionaron positivamente con la densidad de los TIL CD8+ (Spearman r = 0,315 y r = 0,334 respectivamente). Los microarrays de expresión génica (n = 90) mostraron que la expresión CD8A se correlacionó positivamente con la expresión del ligando de muerte programada-1 (r = 0,264) y con el antígeno linfocitario del gen 3 (r = 0,507). CONCLUSIÓN: La activación de células T CD8+ específicas para neoantígenos previa a la CRT puede ser un evento clave en la respuesta a la misma donde las moléculas del punto de control inmunitario podrían ser dianas útiles para intensificar la regresión del tumor.

Reconstruction of transcription-translation dynamics with a model of gene networks.

A transcription-translation model of gene networks and a method to reconstruct it from gene expression data are proposed. The model is a hybrid system based on the Glass network with continuous-time dynamics and logical interactions. Transcription-translation dynamics is introduced into the Glass network. The reconstruction of gene networks is reduced to the problem of estimating logical functions from binary representations of quantities of mRNAs and proteins. The reconstruction method is applied to the gene expression data of circadian rhythms. The response characteristics of the reconstructed gene network to periodic stimuli are analyzed. The results suggest the existence of a receiver gene that responds to an external signal, consistently with biological knowledge.

Host range of poliovirus is restricted to simians because of a rapid sequence change of the poliovirus receptor gene during evolution.

The host range of most poliovirus (PV) strains is restricted to simians. This host range specificity is believed to be determined by the interaction between PV and its receptor molecule. To elucidate the molecular basis of this species-specific infection of PV, we cloned orthologs of the PV receptor (PVR) gene ( pvr) as well as those of PV receptor-related genes 1 and 2 ( prr1 and prr2) from various mammalian species. These three genes are widely present in mammalian genomes including those of non-susceptible species. Comparison of the deduced amino acid sequences of PVR orthologs revealed that the NH(2)-terminal immunoglobulin-like domain (domain 1), which is the virus binding site in the human PVR, is highly variable among species, whereas that of PRR1 is highly conserved. Domain 1 of the PVR orthologs for the ring-tailed lemur and rabbit, which are not susceptible to PV, show only 51 and 61% amino acid sequence identity to that of human PVR, respectively. Chimeric PVR proteins that have the domain 1 of the ring-tailed lemur and rabbit PVRs failed to serve as receptors for PV. These results suggest that rapid changes in the domain 1 sequence during mammalian evolution determined the host range restriction of PV.

The amino acid residues affecting the activity and azole susceptibility of rat CYP51 (sterol 14-demethylase P450).

The amino acid residues affecting the function of rat sterol 14-demethylase P450 (CYP51) were examined by means of point mutation. Forty-five mutants with respect to 27 amino acid sites were constructed and expressed in Escherichia coli. Substitution of highly conserved Y131, E369, R372, or R382 decreased the expression of CYP51 protein, indicating some structural importance of these residues. Substitution of H314, T315, or S316 caused considerable effects on the catalytic activity, and T315 was identified as the "conserved threonine" of CYP51. H314 was important for maintenance of the activity of CYP51 and was a characteristic residue of this P450, because the position corresponding to this residue is occupied by an acidic amino acid in most other P450 species. A144 was identified as a residue affecting the interaction of CYP51 with ketoconazole. Substitution of A144 with I, which occupies the corresponding position in fungal CYP51, enhanced the ketoconazole susceptibility of rat CYP51 with little change in the catalytic activity, indicating an important role of this residue in determination of the ketoconazole susceptibility of CYP51. Alteration of the catalytic activity was caused by the substitution at some other sites, whereas substitution of a few highly conserved amino acids caused little alteration of the activity of CYP51.

Gene structure and chromosomal location of a human bHLH transcriptional factor DEC1 x Stra13 x SHARP-2/BHLHB2.

DEC1/BHLHB2 is a novel cAMP-inducible basic helix-loop-helix (bHLH) transcriptional factor isolated from human chondrocyte cultures by the subtraction method [Shen et al. (1997) Biochem. Biophys. Res. Commun. 236, 294--298]. DEC1 seems to be involved in controlling the proliferation/differentiation of some cell lineages. We determined the structure of the human DEC1 gene and its chromosomal locus. Phylogenetic analysis and comparison of the gene structure showed that the DEC1 protein is a member of a new subgroup of the proline bHLH protein family that diverged earlier than other proline bHLH proteins including HES, hairy and E(spl). The human DEC1 gene spans approximately 5.7 kb and contains 5 exons. The putative promoter region contains multiple GC boxes but no TATA box. A primer extension study showed multiple transcriptional initiation sites. In the 5'-flanking region of the DEC1 gene, several transcriptional factor binding sites, including a cAMP-responsive element (CRE), were found using the transcription factor database. The DEC1 gene locates at Chromosome 3p25.3--26 by the FISH method. This is the first study to determine the genomic structure of the DEC1 gene subgroup.

Identification of novel first exons in Ad4BP/SF-1 (NR5A1) gene and their tissue- and species-specific usage.

It has been demonstrated that the mammalian Ad4BP/SF-1 (NR5A1) gene is regulated precisely in sex, tissue, and developmental stage specific manners. To clarify the complex transcriptional regulation, we investigated in the present study whether the gene transcription is regulated by multiple promoters accompanied by noncoding first exons. Novel first exons (Io and Ig) were identified downstream of the already identified exon Ia. Nucleotide sequences revealed that Ia and Ig exons were well conserved, whereas Io exon was less conserved among the mouse, rat, and human genes. Interestingly, the splice donor of the mouse and human Io and human Ig exons do not satisfy the consensus sequence. Transcripts containing Ia, Io, and Ig were detected in all rat tissues examined, while the transcript containing Io was undetectable in the corresponding tissues of mice. The lack of exon Io usage in the mouse was confirmed by transient transfection assays with cultured cells. Quantitative RT-PCR analysis revealed that the transcript containing Ig exon was the main product in the pituitary but significantly less in the spleen, suggesting that the regulation of Ad4BP/SF-1 gene transcription in the pituitary and spleen is distinct from that of other tissues. The above findings, together with the structural abnormality at the splice donor site, suggest that acquisition of the multiple first exons enables the Ad4BP/SF-1 gene to be regulated differentially in different animal species and in different tissues in the same animal.

Sterol 14-demethylase P450 (CYP51) provides a breakthrough for the discussion on the evolution of cytochrome P450 gene superfamily.

Biodiversity is the most characteristic feature of cytochrome P450. Finding of CYP51 distributing widely in biological kingdoms provided breakthroughs for the discussion on the evolution and diversification of P450. Molecular phylogenetic analysis demonstrated that CYP51 appeared in the prokaryotic era and distributed into most kingdoms concomitant with phylogenetic divergence. This is the first evolutionary evidence indicating the prokaryotic origin of P450. Modification of substrate specificity of eukaryotic CYP51s occurred independently to adapt to the different sterol precursors existing in each kingdom. Formation of CYP51 variants through the mutation of active site and the selection of the advantageous ones from them were demonstrated by the emergence of azole-resistant CYP51s in Candida albicans under the environments rich in azole antifungal agents. These findings illustrate the most probable core process of P450 diversification consisting of modification of active site and selection of the resulting variants through interaction with endogenous and exogenous chemicals.

Homology-based gene structure prediction: simplified matching algorithm using a translated codon (tron) and improved accuracy by allowing for long gaps.

MOTIVATION: Locating protein-coding exons (CDSs) on a eukaryotic genomic DNA sequence is the initial and an essential step in predicting the functions of the genes embedded in that part of the genome. Accurate prediction of CDSs may be achieved by directly matching the DNA sequence with a known protein sequence or profile of a homologous family member(s). RESULTS: A new convention for encoding a DNA sequence into a series of 23 possible letters (translated codon or tron code) was devised to improve this type of analysis. Using this convention, a dynamic programming algorithm was developed to align a DNA sequence and a protein sequence or profile so that the spliced and translated sequence optimally matches the reference the same as the standard protein sequence alignment allowing for long gaps. The objective function also takes account of frameshift errors, coding potentials, and translational initiation, termination and splicing signals. This method was tested on Caenorhabditis elegans genes of known structures. The accuracy of prediction measured in terms of a correlation coefficient (CC) was about 95% at the nucleotide level for the 288 genes tested, and 97. 0% for the 170 genes whose product and closest homologue share more than 30% identical amino acids. We also propose a strategy to improve the accuracy of prediction for a set of paralogous genes by means of iterative gene prediction and reconstruction of the reference profile derived from the predicted sequences. AVAILABILITY: The source codes for the program 'aln' written in ANSI-C and the test data will be available via anonymous FTP at ftp.genome.ad.jp/pub/genomenet/saitama-cc. CONTACT: gotoh@cancer-c.pref.saitama.jp

Association of CYP1B1 genetic polymorphism with incidence to breast and lung cancer.

Cytochrome P450 1B1 (CYP1B1) participates in the metabolic activation of a number of procarcinogens including benzo[a]pyrene and the hydroxylation of 17beta-estradiol at the C-4 position. In this study, we investigated the association between CYP1B1 genetic polymorphism and breast or lung cancer incidence. The Ala-Ser polymorphism at codon 119 in presumed substrate recognition site 1 was significantly associated with the incidence of breast or squamous cell carcinoma of the lung. On the other hand, Leu-Val polymorphism at codon 432 did not show any association to the cancers. An allele containing both Ala and Leu simultaneously, comprised 75% of alleles among 315 Japanese healthy controls, was significantly inversely associated with breast cancer incidence. When expressed in a recombinant system, this CYP1B1 cDNA showed the lowest 17beta-estradiol 4-hydroxylase activity among four different variant forms of CYP1B1. Thus, inter-individual differences in activation of procarcinogens or metabolism of oestrogen originating from genetic polymorphisms of the human CYP1B1 gene may contribute to the susceptibility of human cancers.

Membrane-bound transferrin-like protein (MTf): structure, evolution and selective expression during chondrogenic differentiation of mouse embryonic cells.

Mouse membrane-bound transferrin-like protein (MTf) cDNA was cloned to examine its expression during chondrogenic differentiation in the mouse embryonic cell line ATDC5, and to analyze the phylogenetic relationships among the MTfs of four animal species and 23 other transferrin members. Phylogenetic analysis indicated that the MTf gene diverged from the common ancestor gene earlier than the genes of the other transferrins such as serum transferrin, lactoferrin and ovotransferrin, and that the divergence occurred after the divergence of vertebrates and invertebrates. MTf, as well as the other transferrins, consists of two repeated domains. The similarity between the N-terminal and the C-terminal domains of MTf is much higher than that of the other transferrins, although the five amino acid residues required for iron binding were not conserved in the C-terminal domain of MTf in contrast to the conservation of these residues in both domains of the other transferrins. Among various adult mouse tissues, MTf mRNA was expressed at the highest level in cartilage and at a moderate level in the testis. MTf mRNA was expressed only at very low levels in the brain, spleen, thymus, muscle, lung, skin and intestine, and hardly detected in the heart, kidney, stomach and liver. In cultures of the mouse ATDC5 cell line, MTf is developmentally expressed in parallel with the expression of type II collagen and aggrecan, in the pattern commensurate with the onset of chondrogenesis to form cartilage nodules. The structural characteristics and the expression pattern suggest that during development and in adult tissues, MTf has some functions that are different from those of other transferrins.

Multiple sequence alignment: algorithms and applications.

Elucidation of interrelationships among sequence, structure, function, and evolution (FESS relationships) of a family of genes or gene products is a central theme of modern molecular biology. Multiple sequence alignment has been proven to be a powerful tool for many fields of studies such as phylogenetic reconstruction, illumination of functionally important regions, and prediction of higher order structures of proteins and RNAs. However, it is far too trivial to automatically construct a multiple alignment from a set of related sequences. A variety of methods for solving this computationally difficult problem are reviewed. Several important applications of multiple alignment for elucidation of the FESS relationships are also discussed. For a long period, progressive methods have been the only practical means to solve a multiple alignment problem of appreciable size. This situation is now changing with the development of new techniques including several classes of iterative methods. Today's progress in multiple sequence alignment methods has been made by the multidisciplinary endeavors of mathematicians, computer scientists, and biologists in various fields including biophysicists in particular. The ideas are also originated from various backgrounds, pure algorithmics, statistics, thermodynamics, and others. The outcomes are now enjoyed by researchers in many fields of biological sciences. In the near future, generalized multiple alignment may play a central role in studies of FESS relationships. The organized mixture of knowledge from multiple fields will ferment to develop fruitful results which would be hard to obtain within each area. I hope this review provides a useful information resource for future development of theory and practice in this rapidly expanding area of bioinformatics.

Structure, evolution, and liver-specific expression of sterol 12alpha-hydroxylase P450 (CYP8B).

The rat CYP8B cDNA encoding sterol 12alpha-hydroxylase was cloned and sequenced. The amino acid sequence of the heme-binding region of CYP8B was close to those of CYP7A (cholesterol 7alpha-hydroxylase) and CYP7B (oxysterol 7alpha-hydroxylase). Molecular phylogenetic analysis suggests that CYP8B and the CYP7 family derive from a common ancestor. The P450s of the CYP7 and CYP8 families, except for CYP8A (prostacyclin synthase), catalyze the oxygenation of sterols from an alpha surface in the middle of the steroid skeleton. These facts suggest that CYP8B is a P450 closely linked to those of the CYP7 family. CYP8B was expressed specifically in liver. Hepatic CYP8B mRNA level and the 12alpha-hydroxylase activity were altered by cholestyramine feeding, starvation, streptozotocin-induced diabetes mellitus, and administration of clofibrate, dexamethasone or thyroxin, indicating the pretranslational regulation of CYP8B expression. The enhanced CYP8B mRNA expression in streptozotocin-induced diabetic rats was significantly decreased by insulin within 3 h of its administration. These facts demonstrate a regulatory role of insulin in CYP8B expression as a suppressor.

Formation of azole-resistant Candida albicans by mutation of sterol 14-demethylase P450.

The sterol 14-demethylase P450 (CYP51) of a fluconazole-resistant isolate of Candida albicans, DUMC136, showed reduced susceptibility to this azole but with little change in its catalytic activity. Twelve nucleotide substitutions, resulting in four amino acid changes, were identified in the DUMC136 CYP51 gene in comparison with a reported CYP51 sequence from a wild-type, fluconazole-susceptible C. albicans strain. Seven of these substitutions, including all of those causing amino acid changes, were located within a region covering one of the putative substrate recognition sites of the enzyme (SRS-1). Polymorphisms within this region were observed in several C. albicans isolates, and some were found to be CYP51 heterozygotes. Among the amino acid changes occurring in this region, only an alteration of Y132 was common among these fluconazole-resistant isolates, which suggests the importance of this residue to the fluconazole resistance of the target enzyme. DUMC136 and another fluconazole-resistant isolate were homozygotes with respect to CYP51, although the typical wild-type, fluconazole-susceptible C. albicans was a CYP51 heterozygote. These findings suggest that part of the fluconazole-resistant phenotype of C. albicans DUMC136 was acquired through a mutation-prone area of CYP51, an area which might promote the formation of fluconazole-resistant CYP51, along with a mechanism(s) which allows the formation of a homozygote of this altered CYP51 in this diploid pathogenic yeast.

How should a subarachnoid hemorrhage grading scale be determined? A combinatorial approach based solely on the Glasgow Coma Scale.

OBJECT: The purpose of this study was to present a combinatorial approach used to develop a subarachnoid hemorrhage (SAH) grading scale based on the patient's preoperative Glasgow Coma Scale (GCS) score. METHODS: There are 4094 different combinations that can be used to compress the 13 scores of the GCS into two to 12 grades. Break points, the positions in the scale in which two adjacent scores connote a significantly different outcome, are obtained by a direct comparison of the GCS and the Glasgow Outcome Scale (GOS). Guided by the break points, the number of combinations to be considered can be limited. All possible combinations are statistically analyzed with respect to intergrade differences in outcome. Single combinations, with the maximum number of grades having maximum intergrade outcome differences for each corresponding set of adjacent grades, must be selected. The authors verified the validity of this combinatorial approach by retrospectively analyzing 1398 consecutive patients with aneurysmal SAH who underwent surgery within 7 days of the last hemorrhage episode. The patients' GCS scores were assessed just before surgery and their GOS scores were estimated 6 months post-SAH. The combinatorial approach yields only one acceptable grading scale: I (GCS Score 15); II (GCS Scores 11-14); III (GCS Scores 8-10); IV (GCS Scores 4-7); and V (GCS Score 3). CONCLUSIONS: The combinatorial approach, guided by the break points, is so simple and systematic that it can be used again in the future when revision of the grading scale becomes necessary after development of new and effective treatment modalities that improve patients' overall outcome.

CYP51-like gene of Mycobacterium tuberculosis actually encodes a P450 similar to eukaryotic CYP51.

A CYP51-like gene (Z80226) of Mycobacterium tuberculosis was expressed in Escherichia coli. The product exhibited absorption spectra characteristic of P450. The expressed P450 formed a stoichiometric complex with ketoconazole, one of the specific ligands of CYP51. These findings indicate that the CYP51-like gene of M. tuberculosis actually encodes a P450 having active-site environments similar to those of CYP51, confirming the predicted orthologous nature of this gene to eukaryotic CYP51. Although eukaryotic CYP51s are membrane-binding proteins, the expressed product was accumulated only in the soluble fraction of the host cells.

Structural and phylogenetic analyses of RGD-CAP/beta ig-h3, a fasciclin-like adhesion protein expressed in chick chondrocytes.

A cDNA for RGD-CAP/beta ig-h3 was cloned from a chick embryo chondrocyte cDNA library. The deduced amino acid sequence showed that the chick RGD-CAP/beta ig-h3 is 76-77% identical with human, mouse and pig forms of the protein, and 43% identical with human and mouse osteoblast specific factor 2 (OSF2). RGD-CAP/beta ig-h3 contained four internal repeat domains and two highly conserved sequences (H1 and H2) in each repeat. Chick RGD-CAP/beta ig-h3, as well as the mammalian RGD-CAP/beta ig-h3, contained an RGD sequence, which may serve as a recognition sequence for integrins, in the fourth repeat. Database searches revealed that the H1 and H2 sequences are conserved in some secreted or membrane proteins of several species including mammals, insects, sea urchins, plants, yeast and bacteria. Phylogenetic analysis showed that a portion of the common ancestor gene for RGD-CAP/beta ig-h3 and OSF2 was duplicated to form four repeat domains before the separation of the genes followed by the divergence of vertebrate species.

Divergent structures of Caenorhabditis elegans cytochrome P450 genes suggest the frequent loss and gain of introns during the evolution of nematodes.

The Caenorhabditis elegans genome contains more than 60 cytochrome P450 (CYP) genes. The exon-intron organizations of all of the available and potentially active C. elegans CYP genes were inferred by a newly developed program for predicting protein-coding exons based on the alignment of a genomic DNA sequence and a protein profile. From the predicted amino acid sequences, all of the C. elegans CYP genes except one were classified into three groups, which were closely related to the mammalian drug-metabolizing P450 gene families CYP2, CYP3, and CYP4. The gene structures were strikingly divergent within each group; 20, 10, and 5 unique gene organizations were identified among 40, 18, and 5 genes in the CYP2-, CYP3-, and CYP4-related groups, respectively. The degrees of divergence in gene organization were strongly correlated with those in the amino acid sequences of encoding proteins, and the minimum rate of change in an intron insertion site was estimated to be about 90 times less frequent than amino acid substitutions. Parsimonious analyses suggested that frequent loss and gain of introns has occurred during the evolution of CYP genes in each group after the divergence of nematodes, arthropods, and deuterostomia. Few, if any, incidents of intron sliding were evident, and a model that did not allow intron insertions was highly inconsistent with the observations. All of these findings are explained better by the intron-late view than by the intron-early view.

Molecular cloning and expression of fatty acid alpha-hydroxylase from Sphingomonas paucimobilis.

Fatty acid alpha-hydroxylase (FAAH) catalyzes the initial reaction in alpha-oxidation of fatty acid to produce 2-hydroxy fatty acid. FAAH activity has been detected in a wide range of organisms from prokaryotes to eukaryotes. Here, we describe cloning of the FAAH gene from Sphingomonas paucimobilis, a sphingolipid- and 2-hydroxymyristic acid-rich bacterium. The isolated gene encoded 415 amino acids. A homology search revealed that amino acid sequences highly conserved in cytochrome P450 (P450) were present in FAAH. Although the heme-binding cysteine was recognizable at position 361, the consensus in the heme-binding region was modified by an insertion. Overall, FAAH has no significant identity to the known P450s. CO difference spectrum of recombinant FAAH showed the characteristic one of P450, except this peak was at 445 nm. These results suggest bacterial FAAH is a novel member of the P450 superfamily.