PanTools v3: functional annotation, classification and phylogenomics.

SUMMARY: The ever-increasing number of sequenced genomes necessitates the development of pangenomic approaches for comparative genomics. Introduced in 2016, PanTools is a platform that allows pangenome construction, homology grouping and pangenomic read mapping. The use of graph database technology makes PanTools versatile, applicable from small viral genomes like SARS-CoV-2 up to large plant or animal genomes like tomato or human. Here, we present our third major update to PanTools that enables the integration of functional annotations and provides both gene-level analyses and phylogenetics. AVAILABILITY AND IMPLEMENTATION: PanTools is implemented in Java 8 and released under the GNU GPLv3 license. Software and documentation are available at https://git.wur.nl/bioinformatics/pantools. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

How to Get Better: Taking Notes Mediates the Effect of a Video Tutorial on Number Series.

In recent studies, test-score increases have been shown in rule-based intelligence test tasks, such as number series, after watching a corresponding video tutorial. An open question remains regarding the mechanisms involved. Specifically, taking notes to describe the relations between numbers might be linked to test scores, and is hypothesized to mediate the effect of a number series video tutorial on number series test scores. Therefore, an experimental group (EG) watching a number series video tutorial (n = 58) was compared with a control group (CG) watching an irrelevant tutorial (n = 52) before working on number series items. Results revealed higher number series scores in the EG than the CG (d = .48), more items with provided notes in the EG than in the CG (d = .41), and substantial correlations between the number of items with notes and the number series sum scores in both groups (EG: r = .66; CG: r = .75). The effect of the video tutorial on the number series sum score was mediated by the number of items with notes (indirect effect = 3.41, SE = 1.74). Theoretical and practical implications as well as future research directions are discussed.

The Pectobacterium pangenome, with a focus on Pectobacterium brasiliense, shows a robust core and extensive exchange of genes from a shared gene pool.

BACKGROUND: Bacterial plant pathogens of the Pectobacterium genus are responsible for a wide spectrum of diseases in plants, including important crops such as potato, tomato, lettuce, and banana. Investigation of the genetic diversity underlying virulence and host specificity can be performed at genome level by using a comprehensive comparative approach called pangenomics. A pangenomic approach, using newly developed functionalities in PanTools, was applied to analyze the complex phylogeny of the Pectobacterium genus. We specifically used the pangenome to investigate genetic differences between virulent and avirulent strains of P. brasiliense, a potato blackleg causing species dominantly present in Western Europe. RESULTS: Here we generated a multilevel pangenome for Pectobacterium, comprising 197 strains across 19 species, including type strains, with a focus on P. brasiliense. The extensive phylogenetic analysis of the Pectobacterium genus showed robust distinct clades, with most detail provided by 452,388 parsimony-informative single-nucleotide polymorphisms identified in single-copy orthologs. The average Pectobacterium genome consists of 47% core genes, 1% unique genes, and 52% accessory genes. Using the pangenome, we zoomed in on differences between virulent and avirulent P. brasiliense strains and identified 86 genes associated to virulent strains. We found that the organization of genes is highly structured and linked with gene conservation, function, and transcriptional orientation. CONCLUSION: The pangenome analysis demonstrates that evolution in Pectobacteria is a highly dynamic process, including gene acquisitions partly in clusters, genome rearrangements, and loss of genes. Pectobacterium species are typically not characterized by a set of species-specific genes, but instead present themselves using new gene combinations from the shared gene pool. A multilevel pangenomic approach, fusing DNA, protein, biological function, taxonomic group, and phenotypes, facilitates studies in a flexible taxonomic context.

Meiotic recombination profiling of interspecific hybrid F1 tomato pollen by linked read sequencing.

Genome wide screening of pooled pollen samples from a single interspecific F1 hybrid obtained from a cross between tomato, Solanum lycopersicum and its wild relative, Solanum pimpinellifolium using linked read sequencing of the haploid nuclei, allowed profiling of the crossover (CO) and gene conversion (GC) landscape. We observed a striking overlap between cold regions of CO in the male gametes and our previously established F6 recombinant inbred lines (RILs) population. COs were overrepresented in non-coding regions in the gene promoter and 5'UTR regions of genes. Poly-A/T and AT rich motifs were found enriched in 1 kb promoter regions flanking the CO sites. Non-crossover associated allelic and ectopic GCs were detected in most chromosomes, confirming that besides CO, GC represents also a source for genetic diversity and genome plasticity in tomato. Furthermore, we identified processed break junctions pointing at the involvement of both homology directed and non-homology directed repair pathways, suggesting a recombination machinery in tomato that is more complex than currently anticipated.

PERGOLA: fast and deterministic linkage mapping of polyploids.

BACKGROUND: A large share of agriculturally and horticulturally important plant species are polyploid. Linkage maps are used to locate associations between genes and traits by breeders and geneticists. Linkage map creation for polyploid species is not supported by standard tools. We want to overcome this limitation and validate our results with simulation studies. RESULTS: We developed PERGOLA, a deterministic and heuristic method that addresses this problem. We show that it creates correct linkage groups, marker orders and distances for simulated and real datasets. We compare it to existing tools and demonstrate that it overcomes limitations in ploidy and outperforms them in computational time and mapping accuracy. We represent linkage maps as dendrograms and show that this has advantages in the comparison of different maps. CONCLUSIONS: PERGOLA can be used successfully to calculate linkage maps for diploid and polyploid species and outperforms existing tools.

Advantages of continuous genotype values over genotype classes for GWAS in higher polyploids: a comparative study in hexaploid chrysanthemum.

BACKGROUND: Association studies are an essential part of modern plant breeding, but are limited for polyploid crops. The increased number of possible genotype classes complicates the differentiation between them. Available methods are limited with respect to the ploidy level or data producing technologies. While genotype classification is an established noise reduction step in diploids, it gains complexity with increasing ploidy levels. Eventually, the errors produced by misclassifications exceed the benefits of genotype classes. Alternatively, continuous genotype values can be used for association analysis in higher polyploids. We associated continuous genotypes to three different traits and compared the results to the output of the genotype caller SuperMASSA. Linear, Bayesian and partial least squares regression were applied, to determine if the use of continuous genotypes is limited to a specific method. A disease, a flowering and a growth trait with h (2) of 0.51, 0.78 and 0.91 were associated with a hexaploid chrysanthemum genotypes. The data set consisted of 55,825 probes and 228 samples. RESULTS: We were able to detect associating probes using continuous genotypes for multiple traits, using different regression methods. The identified probe sets were overlapping, but not identical between the methods. Baysian regression was the most restrictive method, resulting in ten probes for one trait and none for the others. Linear and partial least squares regression led to numerous associating probes. Association based on genotype classes resulted in similar values, but missed several significant probes. A simulation study was used to successfully validate the number of associating markers. CONCLUSIONS: Association of various phenotypic traits with continuous genotypes is successful with both uni- and multivariate regression methods. Genotype calling does not improve the association and shows no advantages in this study. Instead, use of continuous genotypes simplifies the analysis, saves computational time and results more potential markers.

A Staphylococcus aureus TIR domain protein virulence factor blocks TLR2-mediated NF-κB signaling.

Signaling through Toll-like receptors (TLRs), crucial molecules in the induction of host defense responses, requires adaptor proteins that contain a Toll/interleukin-1 receptor (TIR) domain. The pathogen Staphylococcus aureus produces several innate immune-evasion molecules that interfere with the host's innate immune response. A database search analysis suggested the presence of a gene encoding a homologue of the human TIR domain in S. aureus MSSA476 which was named staphylococcal TIR domain protein (TirS). Ectopic expression of TirS in human embryonic kidney, macrophage and keratinocyte cell lines interfered with signaling through TLR2, including MyD88 and TIRAP, NF-κB and/or mitogen-activated protein kinase pathways. Moreover, the presence of TirS reduced the levels of cytokines MCP-1 and G-CSF secreted in response to S. aureus. The effects on NF-κB pathway were confirmed using S. aureus MSSA476 wild type, an isogenic mutant MSSA476ΔtirS, and complemented MSSA476ΔtirS +pTirS in a Transwell system where bacteria and host cells were physically separated. Finally, in a systematic mouse infection model, TirS promoted bacterial accumulation in several organs 4 days postinfection. The results of this study reveal a new S. aureus virulence factor that can interfere with PAMP-induced innate immune signaling in vitro and bacterial survival in vivo.

MYCN-regulated miRNA-92 inhibits secretion of the tumor suppressor DICKKOPF-3 (DKK3) in neuroblastoma.

The MYCN oncogene is frequently amplified in neuroblastoma. It is one of the most consistent markers of bad prognosis for this disease. Dickkopf-3 (DKK3) is a secreted protein of the DKK family of Wnt regulators. It functions as a tumor suppressor in a range of cancers, including neuroblastoma. MYCN was recently found to downregulate DKK3 mRNA. In this study, we show that MYCN knockdown in MYCN-amplified (MNA) neuroblastoma cell lines increases secretion of endogenous DKK3 to the culture media. MicroRNAs (miRNAs) are ∼20 nt long single-stranded RNA molecules that downregulate messenger RNAs by targeting the 3' untranslated region (3'UTR). Many miRNAs regulate genes involved in the pathogenesis of cancer and are extensively deregulated in different tumors. Using miRNA target prediction software, we found several MYCN-regulated miRNAs that could target the 3'UTR sequence of DKK3, including mir-92a, mir-92b and let-7e. Luciferase expression from a reporter vector containing the DKK3-3'UTR was decreased when this construct was cotransfected with mir-92a, mir-92b or let-7e in HEK293 cells. Mutation of the mir-92 seed sequence in the 3'UTR completely rescued the observed decrease in reporter expression when cotransfected with mir-92a and mir-92b. Antagomir and miRNA-mimic transfections in neuroblastoma cell lines confirmed that DKK3 secretion to the culture media is regulated by mir-92. Consistent with reports from other cancers, we found DKK3 to be expressed in the endothelium of primary neuroblastoma samples and to be absent in tumors with MYCN amplification. Our data demonstrate that MYCN-regulated miRNAs are able to modulate the expression of the tumor suppressor DKK3 in neuroblastoma.

Conditional expression of retrovirally delivered anti-MYCN shRNA as an in vitro model system to study neuronal differentiation in MYCN-amplified neuroblastoma.

BACKGROUND: Neuroblastoma is a childhood cancer derived from immature cells of the sympathetic nervous system. The disease is clinically heterogeneous, ranging from neuronal differentiated benign ganglioneuromas to aggressive metastatic tumours with poor prognosis. Amplification of the MYCN oncogene is a well established poor prognostic factor found in up to 40% of high risk neuroblastomas.Using neuroblastoma cell lines to study neuronal differentiation in vitro is now well established. Several protocols, including exposure to various agents and growth factors, will differentiate neuroblastoma cell lines into neuron-like cells. These cells are characterized by a neuronal morphology with long extensively branched neurites and expression of several neurospecific markers. RESULTS: In this study we use retrovirally delivered inducible short-hairpin RNA (shRNA) modules to knock down MYCN expression in MYCN-amplified (MNA) neuroblastoma cell lines. By addition of the inducer doxycycline, we show that the Kelly and SK-N-BE(2) neuroblastoma cell lines efficiently differentiate into neuron-like cells with an extensive network of neurites. These cells are further characterized by increased expression of the neuronal differentiation markers NFL and GAP43. In addition, we show that induced expression of retrovirally delivered anti-MYCN shRNA inhibits cell proliferation by increasing the fraction of MNA neuroblastoma cells in the G1 phase of the cell cycle and that the clonogenic growth potential of these cells was also dramatically reduced. CONCLUSION: We have developed an efficient MYCN-knockdown in vitro model system to study neuronal differentiation in MNA neuroblastomas.

Inhibition of mir-21, which is up-regulated during MYCN knockdown-mediated differentiation, does not prevent differentiation of neuroblastoma cells.

BACKGROUND: Neuroblastoma is a malignant childhood tumour arising from precursor cells of the sympathetic nervous system. Genomic amplification of the MYCN oncogene is associated with dismal prognosis. For this group of high-risk tumours, the induction of tumour cell differentiation is part of current treatment protocols. MicroRNAs (miRNAs) are small non-coding RNA molecules that effectively reduce the translation of target mRNAs. MiRNAs play an important role in cell proliferation, apoptosis, differentiation and cancer. In this study, we investigated the role of N-myc on miRNA expression in MYCN-amplified neuroblastoma. We performed a miRNA profiling study on SK-N-BE (2) cells, and determined differentially expressed miRNAs during differentiation initiated by MYCN knockdown, using anti-MYCN short-hairpin RNA (shRNA) technology. RESULTS: Microarray analyses revealed 23 miRNAs differentially expressed during the MYCN knockdown-mediated neuronal differentiation of MNA neuroblastoma cells. The expression changes were bidirectional, with 11 and 12 miRNAs being up- and down-regulated, respectively. Among the down-regulated miRNAs, we found several members of the mir-17 family of miRNAs. Mir-21, an established oncomir in a variety of cancer types, became strongly up-regulated upon MYCN knockdown and the subsequent differentiation. Neither overexpression of mir-21 in the high-MYCN neuroblastoma cells, nor repression of increased mir-21 levels during MYCN knockdown-mediated differentiation had any significant effects on cell differentiation or proliferation. CONCLUSIONS: We describe a subset of miRNAs that were altered during the N-myc deprived differentiation of MYCN-amplified neuroblastoma cells. In this context, N-myc acts as both an activator and suppressor of miRNA expression. Mir-21 was up-regulated during cell differentiation, but inhibition of mir-21 did not prevent this process. We were unable to establish a role for this miRNA during differentiation and proliferation of the two neuroblastoma cell lines used in this study.

Presynaptic and postsynaptic origin of multicomponent extracellular waveforms at the endbulb of Held-spherical bushy cell synapse.

Extracellular signals from the endbulb of Held-spherical bushy cell (SBC) synapse exhibit up to three component waves ('P', 'A' and 'B'). Signals lacking the third component (B) are frequently observed but as the origin of each of the components is uncertain, interpretation of this lack of B has been controversial: is it a failure to release transmitter or a failure to generate or propagate an action potential? Our aim was to determine the origin of each component. We combined single- and multiunit in vitro methods in Mongolian gerbils and Wistar rats and used pharmacological tools to modulate glutamate receptors or voltage-gated sodium channels. Simultaneous extra- and intracellular recordings from single SBCs demonstrated a presynaptic origin of the P-component, consistent with data obtained with multielectrode array recordings of local field potentials. The later components (A and B) correspond to the excitatory postsynaptic potential (EPSP) and action potential of the SBC, respectively. These results allow a clear interpretation of in vivo extracellular signals. We conclude that action potential failures occurring at the endbulb-SBC synaptic junction largely reflect failures of the EPSP to trigger an action potential and not failures of synaptic transmission. The data provide the basis for future investigation of convergence of excitatory and inhibitory inputs in modulating transmission at a fully functional neuronal system using physiological stimulation.

Integrating gene expression and GO classification for PCA by preclustering.

BACKGROUND: Gene expression data can be analyzed by summarizing groups of individual gene expression profiles based on GO annotation information. The mean expression profile per group can then be used to identify interesting GO categories in relation to the experimental settings. However, the expression profiles present in GO classes are often heterogeneous, i.e., there are several different expression profiles within one class. As a result, important experimental findings can be obscured because the summarizing profile does not seem to be of interest. We propose to tackle this problem by finding homogeneous subclasses within GO categories: preclustering. RESULTS: Two microarray datasets are analyzed. First, a selection of genes from a well-known Saccharomyces cerevisiae dataset is used. The GO class "cell wall organization and biogenesis" is shown as a specific example. After preclustering, this term can be associated with different phases in the cell cycle, where it could not be associated with a specific phase previously. Second, a dataset of differentiation of human Mesenchymal Stem Cells (MSC) into osteoblasts is used. For this dataset results are shown in which the GO term "skeletal development" is a specific example of a heterogeneous GO class for which better associations can be made after preclustering. The Intra Cluster Correlation (ICC), a measure of cluster tightness, is applied to identify relevant clusters. CONCLUSIONS: We show that this method leads to an improved interpretability of results in Principal Component Analysis.

Osteo-transcriptomics of human mesenchymal stem cells: accelerated gene expression and osteoblast differentiation induced by vitamin D reveals c-MYC as an enhancer of BMP2-induced osteogenesis.

Bone marrow-derived human mesenchymal stem cells (hMSCs) have the in vitro capacity to differentiate into osteoblasts, chondrocytes or adipocytes, depending on the applied stimulus. In order to identify novel regulators of osteogenesis in hMSCs, osteo-transcriptomics was performed whereby differentiation induced by dexamethasone (DEX), DEX+ bone morphogenetic protein 2 (BMP2), and DEX+ Vitamin D(3) (1,25(OH)(2)D(3)) was studied over a course of 12 days. Microarray analysis revealed that 2095 genes were significantly regulated by DEX+ 1,25(OH)(2)D(3), of which 961 showed accelerated expression kinetics compared to treatment by DEX alone. The majority of these genes were accelerated 24-48 h after onset of osteogenic treatment. Gene ontology (GO) analysis of these 1,25(OH)(2)D(3)-accelerated genes indicated their involvement in biological processes related to cellular differentiation and cell cycle regulation. When compared to cells treated with DEX or DEX+BMP2, treatment with DEX+ 1,25(OH)(2)D(3) clearly accelerated osteoprogenitor commitment and osteoblast maturation, as measured by alkaline phosphatase (ALP) activity and calcification of the matrix. Cell cycle progression, as observed after initial growth arrest, was not significantly accelerated by 1,25(OH)(2)D(3) and was not required for onset and progression of osteogenesis. However, expression of c-Myc was accelerated by 1,25(OH)(2)D(3), and binding sites for c-MYC were enriched in promoters of genes accelerated by 1,25(OH)(2)D(3). Lentiviral overexpression of c-MYC strongly promoted DEX+ BMP2-induced osteoblast differentiation and matrix maturation. In conclusion, our studies show for the first time that 1,25(OH)(2)D(3) strongly accelerates expression of genes involved in differentiation of hMSCs and, moreover, identify c-MYC as a novel regulator of osteogenesis.

Genome-wide mRNA expression analysis of hepatic adaptation to high-fat diets reveals switch from an inflammatory to steatotic transcriptional program.

BACKGROUND: Excessive exposure to dietary fats is an important factor in the initiation of obesity and metabolic syndrome associated pathologies. The cellular processes associated with the onset and progression of diet-induced metabolic syndrome are insufficiently understood. PRINCIPAL FINDINGS: To identify the mechanisms underlying the pathological changes associated with short and long-term exposure to excess dietary fat, hepatic gene expression of ApoE3Leiden mice fed chow and two types of high-fat (HF) diets was monitored using microarrays during a 16-week period. A functional characterization of 1663 HF-responsive genes reveals perturbations in lipid, cholesterol and oxidative metabolism, immune and inflammatory responses and stress-related pathways. The major changes in gene expression take place during the early (day 3) and late (week 12) phases of HF feeding. This is also associated with characteristic opposite regulation of many HF-affected pathways between these two phases. The most prominent switch occurs in the expression of inflammatory/immune pathways (early activation, late repression) and lipogenic/adipogenic pathways (early repression, late activation). Transcriptional network analysis identifies NF-kappaB, NEMO, Akt, PPARgamma and SREBP1 as the key controllers of these processes and suggests that direct regulatory interactions between these factors may govern the transition from early (stressed, inflammatory) to late (pathological, steatotic) hepatic adaptation to HF feeding. This transition observed by hepatic gene expression analysis is confirmed by expression of inflammatory proteins in plasma and the late increase in hepatic triglyceride content. In addition, the genes most predictive of fat accumulation in liver during 16-week high-fat feeding period are uncovered by regression analysis of hepatic gene expression and triglyceride levels. CONCLUSIONS: The transition from an inflammatory to a steatotic transcriptional program, possibly driven by the reciprocal activation of NF-kappaB and PPARgamma regulators, emerges as the principal signature of the hepatic adaptation to excess dietary fat. These findings may be of essential interest for devising new strategies aiming to prevent the progression of high-fat diet induced pathologies.

The hangover gene negatively regulates bouton addition at the Drosophila neuromuscular junction.

The synaptic growth of neurons during the development and adult life of an animal is a very dynamic and highly regulated process. During larval development in Drosophila new boutons and branches are added at the glutamatergic neuromuscular junction (NMJ) until a balance between neuronal activity and morphological structures is reached. Analysis of several Drosophila mutants suggest that bouton number and size might be regulated by separate signaling processes [Budnik, V., 1996. Synapse maturation and structural plasticity at Drosophila neuromuscular junctions. Curr. Opin. Neurobiol. 6, 858-867.]. Here we show a new role for Hangover as a negative regulator of bouton number at the NMJ. The hangover gene (hang) encodes a nuclear zinc finger protein. It has a function in neuronal plasticity mediating ethanol tolerance, a behavior that develops upon previous experience with ethanol. hang(AE10) mutants have more boutons and an extended synaptic span. Moreover, Hang expression in the motoneuron is required for the regulation of bouton number and the overall length of muscle innervation. However, the increase in bouton number does not correlate with a change in synaptic transmission, suggesting a mechanism independent from neuronal activity leads to the surplus of synaptic boutons. In contrast, we find that expression levels of the cell adhesion molecule Fasciclin II (FASII) are reduced in the hang mutant. This finding suggests that the increase in bouton number in hang mutants is caused by a reduction in FASII expression, thus, linking the regulation of nuclear gene expression with the addition of boutons at the NMJ regulated by cell adhesion molecules.

Comparison of RNAi efficiency mediated by tetracycline-responsive H1 and U6 promoter variants in mammalian cell lines.

Conditional expression of short hairpin RNAs (shRNAs) to knock down target genes is a powerful tool to study gene function. The most common inducible expression systems are based on tetracycline-regulated RNA polymerase III promoters. During the last years, several tetracycline-inducible U6 and H1 promoter variants have been reported in different experimental settings showing variable efficiencies. In this study, we compare the most common variants of these promoters in several mammalian cell lines. For all cell lines tested, we find that several inducible U6 and H1 promoters containing single tetracycline operator (tetO) sequences show high-transcriptional background in the non-induced state. Promoter variants containing two tetO sequences show tight suppression of transcription in the non-induced state, and high tet responsiveness and high gene knockdown efficiency upon induction in all cell lines tested. We report a variant of the H1 promoter containing two O2-type tetO sequences flanking the TATA box that shows little transcriptional background in the non-induced state and up to 90% target knockdown when the inducer molecule (dox-doxycycline) is added. This inducible system for RNAi-based gene silencing is a good candidate for use both in basic research on gene function and for potential therapeutic applications.

Quantifying desorption of saturated hydrocarbons from silicon with quantum calculations and scanning tunneling microscopy.

Electron stimulated desorption of cyclopentene from the Si(100)-(2 x 1) surface is studied experimentally with cryogenic UHV STM and theoretically with transport, electronic structure, and dynamical calculations. Unexpectedly for a saturated hydrocarbon on silicon, desorption is observed at bias magnitudes as low as 2.5 V, albeit the desorption yields are a factor of 500 to 1000 lower than previously reported for unsaturated molecules on silicon. The low threshold voltage for desorption is attributed to hybridization of the molecule with the silicon surface, which results in low-lying ionic resonances within 2-3 eV of the Fermi level. These resonances are long-lived, spatially localized, and displaced in equilibrium with respect to the neutral state. This study highlights the importance of nuclear dynamics in silicon-based molecular electronics and suggests new guidelines for the control of such dynamics.

RNA reprogramming of alpha-mannosidase mRNA sequences in vitro by myxomycete group IC1 and IE ribozymes.

Trans-splicing group I ribozymes have been introduced in order to mediate RNA reprogramming (including RNA repair) of therapeutically relevant RNA transcripts. Efficient RNA reprogramming depends on the appropriate efficiency of the reaction, and several attempts, including optimization of target recognition and ribozyme catalysis, have been performed. In most studies, the Tetrahymena group IC1 ribozyme has been applied. Here we investigate the potential of group IC1 and group IE intron ribozymes, derived from the myxomycetes Didymium and Fuligo, in addition to the Tetrahymena ribozyme, for RNA reprogramming of a mutated alpha-mannosidase mRNA sequence. Randomized internal guide sequences were introduced for all four ribozymes and used to select accessible sites within isolated mutant alpha-mannosidase mRNA from mammalian COS-7 cells. Two accessible sites common to all the group I ribozymes were identified and further investigated in RNA reprogramming by trans-splicing analyses. All the myxomycete ribozymes performed the trans-splicing reaction with high fidelity, resulting in the conversion of mutated alpha-mannosidase RNA into wild-type sequence. RNA protection analysis revealed that the myxomycete ribozymes perform trans-splicing at approximately similar efficiencies as the Tetrahymena ribozyme. Interestingly, the relative efficiency among the ribozymes tested correlates with structural features of the P4-P6-folding domain, consistent with the fact that efficient folding is essential for group I intron trans-splicing.

Distinct gene expression patterns in a tamoxifen-sensitive human mammary carcinoma xenograft and its tamoxifen-resistant subline MaCa 3366/TAM.

The reasons why human mammary tumors become resistant to tamoxifen therapy are mainly unknown. Changes in gene expression may occur as cells acquire resistance to antiestrogens. We therefore undertook a comparative gene expression analysis of tamoxifen-sensitive and tamoxifen-resistant human breast cancer in vivo models using Affymetrix oligonucleotide arrays to analyze differential gene expression. Total RNAs from the tamoxifen-sensitive patient-derived mammary carcinoma xenograft MaCa 3366 and the tamoxifen-resistant model MaCa 3366/TAM were hybridized to Affymetrix HuGeneFL and to Hu95Av2 arrays. Pairwise comparisons and clustering algorithms were applied to identify differentially expressed genes and patterns of gene expression. As revealed by cluster analysis, the tamoxifen-sensitive and the tamoxifen-resistant breast carcinomas differed regarding their gene expression pattern. More than 100 transcripts are changed in abundance in MaCa 3366/TAM as compared with MaCa 3366. Among the genes that are differentially expressed in the tamoxifen-resistant tumors, there are several IFN-inducible and estrogen-responsive genes, and genes known to be involved in breast carcinogenesis. The genes neuronatin (NNAT) and bone marrow stem cell antigen 2 (BST2) were sharply up-regulated in MaCa 3366/TAM. The differential expression of four genes (NNAT, BST2, IGFBP5, and BCAS1) was confirmed by Taqman PCR. Our results provide the starting point for deriving markers for tamoxifen resistance by differential gene expression profiling in a human breast cancer model of acquired tamoxifen resistance. Finally, genes whose expression profiles are distinctly changed between the two xenograft lines will be further evaluated as potential targets for diagnostic or therapeutic approaches of tamoxifen-resistant breast cancer.

Comparison of proteomic and genomic analyses of the human breast cancer cell line T47D and the antiestrogen-resistant derivative T47D-r.

In search of novel mechanisms leading to the development of antiestrogen-resistance in human breast tumors, we analyzed differences in the gene and protein expression pattern of the human breast carcinoma cell line T47D and its derivative T47D-r, which is resistant toward the pure antiestrogen ZM 182780 (Faslodex trade mark, fulvestrant). Affymetrix DNA chip hybridizations on the commercially available HuGeneFL and Hu95A arrays were carried out in parallel to the proteomics analysis where the total cellular protein content of T47D or T47D-r was separated on two-dimensional gels. Thirty-eight proteins were found to be reproducibly up- or down-regulated more than 2-fold in T47D-r versus T47D in the proteomics analysis. Comparison with differential mRNA analysis revealed that 19 of these were up- or down-regulated in parallel with the corresponding mRNA molecules, among which are the protease cathepsin D, the GTPases Rab11a and MxA, and the secreted protein hAG-2. For 11 proteins, the corresponding mRNA was not found to be differentially expressed, and for eight proteins an inverse regulation was found at the mRNA level. In summary, mRNA expression data, when combined with proteomic information, provide a more detailed picture of how breast cancer cells are altered in their antiestrogen-resistant compared with the antiestrogen-sensitive state.