Engineering Circularized mRNAs for the Production of Spider Silk Proteins.

Biomaterials offer unique properties that make them irreplaceable for next-generation applications. Fibrous proteins, such as various caterpillar silks and especially spider silk, have strength and toughness not found in human-made materials. In early studies, proteins containing long tandem repeats, such as major ampullate spidroin 1 (MaSp1) and flagelliform silk protein (FSLP), were produced using a large DNA template composed of many tandem repeats. The hierarchical DNA assembly of the DNA template is very time-consuming and labor-intensive, which makes the fibrous proteins difficult to study and engineer. In this study, we designed a circularized mRNA (cmRNA) employing the RNA cyclase ribozyme mechanism. cmRNAs encoding spider silk protein MaSp1 and FSLP were designed based on only one unit of the template sequence but provide ribosomes with a circular and infinite translation template for production of long peptides containing tandem repeats. Using this technique, cmRNAs of MaSp1 and FSLP were successfully generated with circularization efficiencies of 8.5% and 36.7%, respectively, which supported the production of recombinant MaSp1 and FSLP larger than 110 and 88 kDa, containing tens of repeat units. Western blot analysis and mass spectrometry confirmed the authenticity of MaSp1 and FSLP, which were produced at titers of 22.1 and 81.5 mg · liter-1, respectively. IMPORTANCE Spider silk is a biomaterial with superior properties. However, its heterologous expression template is hard to construct. The cmRNA technique simplifies the construction and expression strategy by proving the ribosome a circular translation template for expression of long peptides containing tandem repeats. This revolutionary technique will allow researchers to easily build, study, and experiment with any fiber proteins with sequences either from natural genes or artificial designs. We expect a significantly accelerated development of fibrous protein-based biomaterials with the cmRNA technique.

Urine Proteome Profiling Predicts Lung Cancer from Control Cases and Other Tumors.

Development of noninvasive, reliable biomarkers for lung cancer diagnosis has many clinical benefits knowing that most of lung cancer patients are diagnosed at the late stage. For this purpose, we conducted proteomic analyses of 231 human urine samples in healthy individuals (n=33), benign pulmonary diseases (n=40), lung cancer (n=33), bladder cancer (n=17), cervical cancer (n=25), colorectal cancer (n=22), esophageal cancer (n=14), and gastric cancer (n=47) patients collected from multiple medical centers. By random forest modeling, we nominated a list of urine proteins that could separate lung cancers from other cases. With a feature selection algorithm, we selected a panel of five urinary biomarkers (FTL: Ferritin light chain; MAPK1IP1L: Mitogen-Activated Protein Kinase 1 Interacting Protein 1 Like; FGB: Fibrinogen Beta Chain; RAB33B: RAB33B, Member RAS Oncogene Family; RAB15: RAB15, Member RAS Oncogene Family) and established a combinatorial model that can correctly classify the majority of lung cancer cases both in the training set (n=46) and the test sets (n=14-47 per set) with an AUC ranging from 0.8747 to 0.9853. A combination of five urinary biomarkers not only discriminates lung cancer patients from control groups but also differentiates lung cancer from other common tumors. The biomarker panel and the predictive model, when validated by more samples in a multi-center setting, may be used as an auxiliary diagnostic tool along with imaging technology for lung cancer detection.

Proof-of-Concept Workflow for Establishing Reference Intervals of Human Urine Proteome for Monitoring Physiological and Pathological Changes.

Urine as a true non-invasive sampling source holds great potential for biomarker discovery. While approximately 2000 proteins can be detected by mass spectrometry in urine from healthy people, the amount of these proteins vary considerably. A systematic evaluation of a large number of samples is needed to determine the range of the variations. Current biomarker studies often measure limited number of urine samples in the discovery phase, which makes it difficult to determine whether proteins differentially expressed between control and disease groups represent actual difference, or are just physiological variations among the individuals, leads to failures in the validation phase with the increased sample numbers. Here, we report a streamlined workflow with capacity of measuring 8 urine proteomes per day at the coverage of >1500 proteins. With this workflow, we evaluated variations in 497 urine proteomes from 167 healthy donors, establishing reference intervals (RIs) that covered urine protein variations. We demonstrated that RIs could be used to monitor physiological changes by detecting transient outlier proteins. Furthermore, we provided a RIs-based algorithm for biomarker discovery and validation to screen for diseases such as cancer. This study provided a proof-of-principle workflow for the use of urine proteome for health monitoring and disease screening.

Acetylation-regulated interaction between p53 and SET reveals a widespread regulatory mode.

Although lysine acetylation is now recognized as a general protein modification for both histones and non-histone proteins, the mechanisms of acetylation-mediated actions are not completely understood. Acetylation of the C-terminal domain (CTD) of p53 (also known as TP53) was an early example of non-histone protein acetylation and its precise role remains unclear. Lysine acetylation often creates binding sites for bromodomain-containing 'reader' proteins. Here we use a proteomic screen to identify the oncoprotein SET as a major cellular factor whose binding with p53 is dependent on CTD acetylation status. SET profoundly inhibits p53 transcriptional activity in unstressed cells, but SET-mediated repression is abolished by stress-induced acetylation of p53 CTD. Moreover, loss of the interaction with SET activates p53, resulting in tumour regression in mouse xenograft models. Notably, the acidic domain of SET acts as a 'reader' for the unacetylated CTD of p53 and this mechanism of acetylation-dependent regulation is widespread in nature. For example, acetylation of p53 also modulates its interactions with similar acidic domains found in other p53 regulators including VPRBP (also known as DCAF1), DAXX and PELP1 (refs. 7, 8, 9), and computational analysis of the proteome has identified numerous proteins with the potential to serve as acidic domain readers and lysine-rich ligands. Unlike bromodomain readers, which preferentially bind the acetylated forms of their cognate ligands, the acidic domain readers specifically recognize the unacetylated forms of their ligands. Finally, the acetylation-dependent regulation of p53 was further validated in vivo by using a knock-in mouse model expressing an acetylation-mimicking form of p53. These results reveal that acidic-domain-containing factors act as a class of acetylation-dependent regulators by targeting p53 and, potentially, other proteins.

A divergent role of the SIRT1-TopBP1 axis in regulating metabolic checkpoint and DNA damage checkpoint.

DNA replication is executed only when cells have sufficient metabolic resources and undamaged DNA. Nutrient limitation and DNA damage cause a metabolic checkpoint and DNA damage checkpoint, respectively. Although SIRT1 activity is regulated by metabolic stress and DNA damage, its function in these stress-mediated checkpoints remains elusive. Here we report that the SIRT1-TopBP1 axis functions as a switch for both checkpoints. With glucose deprivation, SIRT1 is activated and deacetylates TopBP1, resulting in TopBP1-Treslin disassociation and DNA replication inhibition. Conversely, SIRT1 activity is inhibited under genotoxic stress, resulting in increased TopBP1 acetylation that is important for the TopBP1-Rad9 interaction and activation of the ATR-Chk1 pathway. Mechanistically, we showed that acetylation of TopBP1 changes the conformation of TopBP1, thereby facilitating its interaction with distinct partners in DNA replication and checkpoint activation. Taken together, our studies identify the SIRT1-TopBP1 axis as a key signaling mode in the regulation of the metabolic checkpoint and the DNA damage checkpoint.

Analysis of gene expression changes in Trichophyton rubrum after skin interaction.

Trichophyton rubrum, an anthropophilic and cosmopolitan fungus, is the most common agent of superficial mycoses. In this study, T. rubrum infection was modelled by adding human skin sections to a limited medium containing glucose and cDNA microarrays were used to monitor T. rubrum gene expression patterns on a global level. We observed that exposure to human skin resulted in upregulation of the expression levels of T. rubrum genes related to many cellular and biological processes, including transcription and translation, metabolism and secondary transport, the stress response, and signalling pathways. These results provide a reference set of T. rubrum genes whose expression patterns change upon infection and reveal previously unknown genes that most likely correspond to proteins that should be considered as virulence factor candidates and potential new drug targets for T. rubrum infection.

A comprehensive proteomic analysis of Mycobacterium bovis bacillus Calmette-Guérin using high resolution Fourier transform mass spectrometry.

Since 1921, Mycobacterium bovis bacillus Calmette-Guérin (BCG) has been recognized as an important vaccine to prevent tuberculosis worldwide. Nonetheless, a global analysis of BCG proteome has not been clearly investigated. In this study, we performed an in-depth proteomic analysis of BCG under an in vitro cultivation condition using SDS-PAGE and high resolution Fourier transform mass spectrometry. In total, 3434 proteins (35,259 unique peptides) including 512 transmembrane proteins were identified, covering ~87% of the predicted BCG proteome. Seven pseudogene protein products were also obtained and validated by RT-PCR at gene transcript level. Additionally, translational start sites of 832 proteins were confirmed and 186 were extended using N-terminus-derived peptides. The physicochemical characteristics of all identified proteins were determined. Some predominant proteins, including PE and PPE family proteins, lipoproteins, heat shock proteins, transport proteins and low molecular weight protein antigens, are discussed, which represent potential prominent antigens in the humoral and cellular immune response. This study represents the most comprehensive BCG proteome to date, which will likely facilitate the design of vaccination and immunodiagnostic strategies against TB.

A proteogenomic analysis of Shigella flexneri using 2D LC-MALDI TOF/TOF.

BACKGROUND: New strategies for high-throughput sequencing are constantly appearing, leading to a great increase in the number of completely sequenced genomes. Unfortunately, computational genome annotation is out of step with this progress. Thus, the accurate annotation of these genomes has become a bottleneck of knowledge acquisition. RESULTS: We exploited a proteogenomic approach to improve conventional genome annotation by integrating proteomic data with genomic information. Using Shigella flexneri 2a as a model, we identified total 823 proteins, including 187 hypothetical proteins. Among them, three annotated ORFs were extended upstream through comprehensive analysis against an in-house N-terminal extension database. Two genes, which could not be translated to their full length because of stop codon 'mutations' induced by genome sequencing errors, were revised and annotated as fully functional genes. Above all, seven new ORFs were discovered, which were not predicted in S. flexneri 2a str.301 by any other annotation approaches. The transcripts of four novel ORFs were confirmed by RT-PCR assay. Additionally, most of these novel ORFs were overlapping genes, some even nested within the coding region of other known genes. CONCLUSIONS: Our findings demonstrate that current Shigella genome annotation methods are not perfect and need to be improved. Apart from the validation of predicted genes at the protein level, the additional features of proteogenomic tools include revision of annotation errors and discovery of novel ORFs. The complementary dataset could provide more targets for those interested in Shigella to perform functional studies.

Global gene expression profiles for the growth phases of Trichophyton rubrum.

Trichophyton rubrum (T. rubrum) is a common superficial fungus. Molecular and genetic studies of T. rubrum are still limited. In this paper, we report the global analysis of gene expression profiles at different growth phases using cDNA microarray technology. A total of 2044 differentially expressed genes were obtained and clustered into three expression patterns. Our data confirmed previous results that many mRNAs were pre-stored in the conidia of T. rubrum. Transcriptional profiling and function analysis showed that some glycolytic enzymes share similar expression patterns and may be coregulated during the transition of growth phases. Some genes involved in small GTPase signaling pathways, and in cAMP-dependent and MAPK regulation pathways were induced in response to the growth dynamics of T. rubrum. Although the detailed biological roles of these T. rubrum genes are still unknown, our results suggest that these genes may be involved in regulation mechanisms in the life cycle of the fungus.

Combining blue native polyacrylamide gel electrophoresis with liquid chromatography tandem mass spectrometry as an effective strategy for analyzing potential membrane protein complexes of Mycobacterium bovis bacillus Calmette-Guérin.

BACKGROUND: Tuberculosis is an infectious bacterial disease in humans caused primarily by Mycobacterium tuberculosis, and infects one-third of the world's total population. Mycobacterium bovis bacillus Calmette-Guérin (BCG) vaccine has been widely used to prevent tuberculosis worldwide since 1921. Membrane proteins play important roles in various cellular processes, and the protein-protein interactions involved in these processes may provide further information about molecular organization and cellular pathways. However, membrane proteins are notoriously under-represented by traditional two-dimensional polyacrylamide gel electrophoresis (2-D PAGE) and little is known about mycobacterial membrane and membrane-associated protein complexes. Here we investigated M. bovis BCG by an alternative proteomic strategy coupling blue native PAGE to liquid chromatography tandem mass spectrometry (LC-MS/MS) to characterize potential protein-protein interactions in membrane fractions. RESULTS: Using this approach, we analyzed native molecular composition of protein complexes in BCG membrane fractions. As a result, 40 proteins (including 12 integral membrane proteins), which were organized in 9 different gel bands, were unambiguous identified. The proteins identified have been experimentally confirmed using 2-D SDS PAGE. We identified MmpL8 and four neighboring proteins that were involved in lipid transport complexes, and all subunits of ATP synthase complex in their monomeric states. Two phenolpthiocerol synthases and three arabinosyltransferases belonging to individual operons were obtained in different gel bands. Furthermore, two giant multifunctional enzymes, Pks7 and Pks8, and four mycobacterial Hsp family members were determined. Additionally, seven ribosomal proteins involved in polyribosome complex and two subunits of the succinate dehydrogenase complex were also found. Notablely, some proteins with high hydrophobicity or multiple transmembrane helixes were identified well in our work. CONCLUSIONS: In this study, we utilized LC-MS/MS in combination with blue native PAGE to characterize modular components of multiprotein complexes in BCG membrane fractions. The results demonstrated that the proteomic strategy was a reliable and reproducible tool for analysis of BCG multiprotein complexes. The identification in our study may provide some evidence for further study of BCG protein interaction.

Transcriptional profile of the Shigella flexneri response to an alkaloid: berberine.

Berberine, a natural isoquinoline alkaloid found in many medicinal herbs, has been shown to be active against a variety of microbial infections. To examine the potential effects of berberine on Shigella flexneri, a whole-genome DNA microarray was constructed and a transcriptome analysis of the cellular responses of S. flexneri when exposed to berberine chloride (BC) was performed. Our data revealed that BC upregulated a group of genes involved in DNA replication, repair and division. Intriguingly, the expression of many genes related to cell envelope biogenesis was increased. In addition, many genes involved in cell secretion, nucleotide metabolism, translation, fatty acid metabolism and the virulence system were also induced by the drug. However, more genes from the functional classes of carbohydrate metabolism, energy production and conversion as well as amino acid metabolism were significantly repressed than were induced. These results provide a comprehensive view of the changes in gene expression when S. flexneri was exposed to BC, and shed light on its complicated effects on this pathogen.

Expression dynamics of secreted protease genes in Trichophyton rubrum induced by key host's proteinaceous components.

Trichophyton rubrum is the most common agent of dermatophytosis, a disease that affects millions of individuals worldwide. Its molecular pathogenicity mechanisms are still not completely elucidated. It has been widely recognized that proteases secreted by T. rubrum are the key virulence factors during host infection. However, our knowledge about the expression of its secreted proteases in host infection is still obscure. This investigation provides the expression patterns and dynamics of secreted protease genes belonging to the subtilisins (SUB) and metalloproteases (MEP) gene families in T. rubrum. The data was obtained under simulated host infection conditions through relative quantification of real time PCR. Keratin, collagen, and elastin induced the expression of similar protease genes, and the expression patterns and dynamics of these protease genes in media containing human skin sections were different from those in media containing individual protein substrates. According to the expression dynamics of these protease genes, we conclude that Sub3, Sub4, and Mep4 may be the dominant proteases secreted by T. rubrum during host infection, and that these proteases could be good targets for new antifungal chemotherapy and molecular diagnostic markers. This work presents useful molecular details to further our understanding of the pathogenesis of dermatophytosis.

Recent dermatophyte divergence revealed by comparative and phylogenetic analysis of mitochondrial genomes.

BACKGROUND: Dermatophytes are fungi that cause superficial infections of the skin, hair, and nails. They are the most common agents of fungal infections worldwide. Dermatophytic fungi constitute three genera, Trichophyton, Epidermophyton, and Microsporum, and the evolutionary relationships between these genera are epidemiologically important. Mitochondria are considered to be of monophyletic origin and mitochondrial sequences offer many advantages for phylogenetic studies. However, only one complete dermatophyte mitochondrial genome (E. floccosum) has previously been determined. RESULTS: The complete mitochondrial DNA sequences of five dermatophyte species, T. rubrum (26,985 bp), T. mentagrophytes (24,297 bp), T. ajelloi (28,530 bp), M. canis (23,943 bp) and M. nanum (24,105 bp) were determined. These were compared to the E. floccosum sequence. Mitochondrial genomes of all 6 species were found to harbor the same set of genes arranged identical order indicating that these dermatophytes are closely related. Genome size differences were largely due to variable lengths of non-coding intergenic regions and the presence/absence of introns. Phylogenetic analyses based on complete mitochondrial genomes reveals that the divergence of the dermatophyte clade was later than of other groups of pathogenic fungi. CONCLUSION: This is the first systematic comparative genomic study on dermatophytes, a highly conserved and recently-diverged lineage of ascomycota fungi. The data reported here provide a basis for further exploration of interrelationships between dermatophytes and will contribute to the study of mitochondrial evolution in higher fungi.

Transcriptional profiles of Trichophyton rubrum in response to itraconazole.

Trichophyton rubrum is the predominant causative agent for superficial dermatomycosis. In order to understand how triazole antifungal agents interact with dermatophytes, the gene expression response of T. rubrum to itraconazole was studied by large-scale gene expression profiling. A total of 670 genes were found to be responsive to itraconazole, including 305 that were up-regulated and 365 down-regulated. Most genes involved in lipid metabolism and especially in ergosterol biosynthesis were up-regulated in response to itraconazole, including ERG6, ERG7, ERG11, ERG24, ERG25 and ERG26. In addition, transcription of some genes involved in cell stress response, drug efflux, and small molecule transport was also affected by itraconazole. Differential expression of selected genes was confirmed by real-time reverse transcription polymerase chain reaction (RT-PCR). This is the first microarray hybridization analysis of T. rubrum exposed to a triazole antifungal agent.

Proteomic profile of dormant Trichophyton rubrum conidia.

BACKGROUND: Trichophyton rubrum is the most common dermatophyte causing fungal skin infections in humans. Asexual sporulation is an important means of propagation for T. rubrum, and conidia produced by this way are thought to be the primary cause of human infections. Despite their importance in pathogenesis, the conidia of T. rubrum remain understudied. We intend to intensively investigate the proteome of dormant T. rubrum conidia to characterize its molecular and cellular features and to enhance the development of novel therapeutic strategies. RESULTS: The proteome of T. rubrum conidia was analyzed by combining shotgun proteomics with sample prefractionation and multiple enzyme digestion. In total, 1026 proteins were identified. All identified proteins were compared to those in the NCBI non-redundant protein database, the eukaryotic orthologous groups database, and the gene ontology database to obtain functional annotation information. Functional classification revealed that the identified proteins covered nearly all major biological processes. Some proteins were spore specific and related to the survival and dispersal of T. rubrum conidia, and many proteins were important to conidial germination and response to environmental conditions. CONCLUSION: Our results suggest that the proteome of T. rubrum conidia is considerably complex, and that the maintenance of conidial dormancy is an intricate and elaborate process. This data set provides the first global framework for the dormant T. rubrum conidia proteome and is a stepping stone on the way to further study of the molecular mechanisms of T. rubrum conidial germination and the maintenance of conidial dormancy.

Membrane subproteomic analysis of Mycobacterium bovis bacillus Calmette-Guérin.

Mycobacterium bovis bacillus Calmette-Guérin (BCG) vaccine has been known for a long time to prevent tuberculosis (TB) worldwide since 1921. Nonetheless, we know little about BCG membrane proteome. In the present study, we utilized alkaline incubation and Triton X-114-based methods to enrich BCG membrane proteins and subsequently digested them using proteolytic enzyme. The recovered peptides were further separated by 2-D LC and identified by ESI-MS/MS. As a result, total 474 proteins were identified, including 78 integral membrane proteins (IMPs). Notably, 18 BCG IMPs were described for the first time in mycobacterium. Further analysis of the 78 IMPs indicated that the theoretical molecular mass distribution of them ranged from 8.06 to 167.86 kDa and pI scores ranged from 4.40 to 11.60. Functional classification revealed that a large proportion of the identified IMPs (67.9%, 53 out of 78) were involved in cell wall and cell processes functional group. In conclusion, here we reported a comprehensive profile of the BCG membrane subproteome. The present investigation may allow the identification of some valuable vaccine and drug target candidates and thus provide basement for future designing of preventive, diagnostic, and therapeutic strategies against TB.

Transcriptional profile induced by furazolidone treatment of Shigella flexneri.

Shigella flexneri is a facultative intracellular pathogen responsible for endemic shigellosis especially in developing countries. Furazolidone, a nitrofuran derivative, is very effective against the infection with S. flexneri. To examine potential effects of furazolidone on this germ, a whole-genome DNA microarray was constructed and transcriptional profiles of the responses to furazolidone were determined. The expressing data revealed adaptive responses of S. flexneri to oxidative stress induced by furazolidone treatment. Iron metabolism was found to be disturbed by furazolidone through derepression of the iron uptake regulon. In addition, energy metabolism, amino acid metabolism, cofactors metabolism, and DNA repair system were also affected by the drug. These data establish a potential for furazolidone to enhance free radical reactions through reductive activation by oxygen-sensitive nitroreductase. Moreover, we provide evidence that furazolidone is able to cause metabolic dysfunction, which cannot always be attributed to oxidative stress, and interactions between reductive metabolites of furazolidone and S. flexneri should be considered.

Comparison between gene expression of conidia and germinating phase in Trichophyton rubrum.

Trichophyton rubrum is a dominating superficial dermatophyte, whose conidial germination is correlated to pathopoiesis and a highly important developmental process. To investigate the changes of physiology, biochemistry and cytology during the germination, we selected 3364 function identified ESTs from T. rubrum cDNA library to construct cDNA microarrays, and compared the gene expression levels of conidia and germinating phase. Data analysis indicated that 335 genes were up-regulated during the germination, which mainly encoded translated, modified proteins and structural proteins. The constituents of cell wall and cell membrane were synthetized abundantly, suggesting that they are the foundation of cell morphogenesis. The ingredients of the two-component signal transduction system were up-regulated, presuming that they were important for the conidial germination. Genes of various metabolic pathways were expressed prosperously, especially the genes that participated in glycolysis and oxidative phosphorylation were up-regulated on the whole, demonstrating that in the environment with sufficient oxygen and glucose, conidia obtained energy through aerobic respiration. This paper provides important clues which are helpful to understanding the changes in gene expression, signal conduction and metabolism characteristics during T. rubrum conidial germination, and possess significant meaning to the study of other superficial dermatophytes.

The use of global transcriptional analysis to reveal the biological and cellular events involved in distinct development phases of Trichophyton rubrum conidial germination.

BACKGROUND: Conidia are considered to be the primary cause of infections by Trichophyton rubrum. RESULTS: We have developed a cDNA microarray containing 10250 ESTs to monitor the transcriptional strategy of conidial germination. A total of 1561 genes that had their expression levels specially altered in the process were obtained and hierarchically clustered with respect to their expression profiles. By functional analysis, we provided a global view of an important biological system related to conidial germination, including characterization of the pattern of gene expression at sequential developmental phases, and changes of gene expression profiles corresponding to morphological transitions. We matched the EST sequences to GO terms in the Saccharomyces Genome Database (SGD). A number of homologues of Saccharomyces cerevisiae genes related to signalling pathways and some important cellular processes were found to be involved in T. rubrum germination. These genes and signalling pathways may play roles in distinct steps, such as activating conidial germination, maintenance of isotropic growth, establishment of cell polarity and morphological transitions. CONCLUSION: Our results may provide insights into molecular mechanisms of conidial germination at the cell level, and may enhance our understanding of regulation of gene expression related to the morphological construction of T. rubrum.

cDNA microarray analysis of the expression profiles of Trichophyton rubrum in response to novel synthetic fatty acid synthase inhibitor PHS11A.

Trichophyton rubrum (T. rubrum) is a major pathogen responsible for dermatophytosis. Because of potential relapse of disease with current antifungal therapy protocols, there is a need for additional and/or alternative antifungal agents for the treatment of disease caused by T. rubrum. We synthesized a potent fungal fatty acid synthase inhibitor, PHS11A, based on the structure of fungal fatty acid synthase. The antifungal activities of PHS11A were tested against 38 clinical isolates of T. rubrum and compared with those of ketoconazole and terbinafine, the MIC(50) and MIC(90) of PHS11A on the isolates were 2 and 4 microg/ml, respectively. We evaluated the transcriptional response of T. rubrum hyphae exposed to PHS11A using 11,232-spot cDNA microarrays. PHS11A exposure increased transcription of fatty acid synthases (FASs) genes FAS1 and FAS2. PHS11A also affected transcription of some genes involved in lipid metabolism, cAMP and MAPK pathways, and multidrug resistance. Quantitative real-time PCR was performed for selected genes to verify the microarray results.