The role of cathepsin K in lung development and newborn chronic lung injury.

Chronic lung diseases, specifically bronchopulmonary dysplasia (BPD), are still causing mortality and morbidity amongst newborn infants. High protease activity has been suggested to have a deleterious role in oxygen-induced lung injuries. Cathepsin K (CatK) is a potent protease found in fetal lungs, degrading collagen and elastin. We hypothesized that CatK may be an important modulator of chronic lung injury in newborn infants and neonatal mice. First we measured CatK protein levels in repeated tracheal aspirate fluid samples from 13 intubated preterm infants during the first two weeks of life. The amount of CatK at 9-13 days was low in infants developing chronic lung disease. Consequently, we studied CatK mRNA expression in oxygen-exposed wild-type (WT) rats at postnatal day (PN) 14 and found decreased pulmonary mRNA expression of CatK in whole lung samples. Thereafter we demonstrated that CatK deficiency modifies lung development by accelerating the thinning of alveolar walls in newborn mice. In hyperoxia-exposed newborn mice CatK deficiency resulted in increased number of pulmonary foam cells, macrophages and amount of reduced glutathione in lung homogenates indicating intensified pulmonary oxidative stress and worse pulmonary outcome due to CatK deficiency. Conversely, transgenic overexpression of CatK caused slight enlargement of distal airspaces with increased alveolar chord length in room air in neonatal mice. While hyperoxic exposure inhibited alveolarization and resulted in enlarged airspaces in wild-type mice, these changes were significantly milder in CatK overexpressing mice at PN7. Finally, we showed that the expression of macrophage scavenger receptor 2 (MSR2) mRNA was down-regulated in oxygen-exposed CatK-deficient mice analyzed by microarray analysis. Our results demonstrate that CatK seems to participate in normal lung development and its expression is altered during pulmonary injury. In the presence of pulmonary risk factors, like high oxygen exposure, low amount of CatK may contribute to aggravated lung injury while sustained or slightly elevated amount of CatK may even protect the newborn lungs from excessive injury. Besides collagen degrading and antifibrotic function of CatK in the lungs, it is obvious that CatK may affect macrophage activity and modify oxidative stress response. In conclusion, pulmonary proteases, specifically CatK, have distinct roles in lung homeostasis and injury development, and although suggested, broad range inhibition of proteases may not be beneficial in newborn lung injury.

Algorithmic Techniques in Gene Expression Processing. From Imputation to Visualization

The amount of biological data has grown exponentially in recent decades. Modern biotechnologies, such as microarrays and next-generation sequencing, are capable to produce massive amounts of biomedical data in a single experiment. As the amount of the data is rapidly growing there is an urgent need for reliable computational methods for analyzing and visualizing it. This thesis addresses this need by studying how to efficiently and reliably analyze and visualize high-dimensional data, especially that obtained from gene expression microarray experiments. First, we will study the ways to improve the quality of microarray data by replacing (imputing) the missing data entries with the estimated values for these entries. Missing value imputation is a method which is commonly used to make the original incomplete data complete, thus making it easier to be analyzed with statistical and computational methods. Our novel approach was to use curated external biological information as a guide for the missing value imputation. Secondly, we studied the effect of missing value imputation on the downstream data analysis methods like clustering. We compared multiple recent imputation algorithms against 8 publicly available microarray data sets. It was observed that the missing value imputation indeed is a rational way to improve the quality of biological data. The research revealed differences between the clustering results obtained with different imputation methods. On most data sets, the simple and fast k-NN imputation was good enough, but there were also needs for more advanced imputation methods, such as Bayesian Principal Component Algorithm (BPCA). Finally, we studied the visualization of biological network data. Biological interaction networks are examples of the outcome of multiple biological experiments such as using the gene microarray techniques. Such networks are typically very large and highly connected, thus there is a need for fast algorithms for producing visually pleasant layouts. A computationally efficient way to produce layouts of large biological interaction networks was developed. The algorithm uses multilevel optimization within the regular force directed graph layout algorithm.

The biological effects of high-pressure gas on the yeast transcriptome

The aim of the present study was to examine the feasibility of DNA microarray technology in an attempt to construct an evaluation system for determining gas toxicity using high-pressure conditions, as it is well known that pressure increases the concentration of a gas. As a first step, we used yeast (Saccharomyces cerevisiae) as the indicator organism and analyzed the mRNA expression profiles after exposure of yeast cells to nitrogen gas. Nitrogen gas was selected as a negative control since this gas has low toxicity. Yeast DNA microarray analysis revealed induction of genes whose products were localized to the membranes, and of genes that are involved in or contribute to energy production. Furthermore, we found that nitrogen gas significantly affected the transport system in the cells. Interestingly, nitrogen gas also resulted in induction of cold-shock responsive genes. These results suggest the possibility of applying yeast DNA microarray to gas bioassays up to 40 MPa. We therefore think that "bioassays" are ideal for use in environmental control and protection studies.

Relationship between Daily Habits, Oral Microbiota and Caries with Special Reference to Xylitol

Caries is a plaque-associated multifactorial chronic disease. Oral hygiene habits, sugar, and oral micobiota interactions are important for caries to occur. Xylitol has been shown to reduce caries mainly due to its effects on mutans streptococci (MS). The purpose of this study was to evaluate the relationship of daily oral health habits and bacterial level on the caries occurrence and to study the effect of xylitol on the composition of oral microflora. A total of 192, 10-12 years old, male school children had been screened for salivary MS. Healthy subjects with high MS counts participated in two parallel double-blinded, randomised, controlled trials. In the first 5-week trial, subjects were assigned into xylitol (n=35) and sorbitol gum (n=38) groups. At baseline, children were examined using International Caries Detection and Assessment System (ICDAS) criteria and interviewed for oral health habits. In the second 4-week trial, subjects were assigned into xylitol (n=25) and saccharine mouthrinse (n=25) groups. In the end of both interventions, saliva samples were collected. The samples were analysed for changes in MS counts and changes in the composition of the oral microbiota assessed by the Human Oral Microbe Identification Microarray (HOMIM). Relationships between daily habits, bacterial levels and caries were evaluated. Daily use of sweets and soft drinks were the habits significantly associated with caries severity measured by ICDAS Caries Index (CI), while toothbrushing was the only habit associated with the low caries severity. Abiotrophia defectiva and Actinomyces meyeri/ A. odontolyticus were significantly higher in caries-affected children while Shuttleworthia satelles was significantly higher in caries-free children. Xylitol showed significant reduction in salivary levels of MS in both trials. No significant effects on other members of the microbiota were found when evaluated by HOMIM. In conclusion, other members of oral microbiota than MS may be associated with caries occurrence or absence. The use of xylitol had significant effect on MS with no effects on the other members of the salivary microbiota.

Early Life Intestinal Microbiota in Health and in Atopic Eczema

Decrease in microbial contacts in affluent societies is considered to lie behind the rise in allergic and other chronic inflammatory diseases during the last decades. Indeed, deviations in the intestinal microbiota composition and diversity have been associated with several diseases, such as atopic eczema. However, there is no consensus yet on what would constitute a beneficial or harmful microbiota. The aim of this thesis was to study the microbiota development in healthy infants and to characterize intestinal microbiota signatures associated with disease status and severity in infants with atopic eczema. The methodological aim was to compare and optimize methods for DNA extraction from fecal samples to be used in high-throughput microbiota analyses. It was confirmed that the most critical step in successful microbial DNA extraction from fecal samples is the mechanical cell lysis procedure. Based on this finding, an efficient semi-automated extraction process was developed that can be scaled for use in high-throughput platforms such as phylogenetic microarray used in this series of studies. By analyzing a longitudinal motherchild cohort for 3 years it was observed that the microbiota development is a gradual process, where some bacterial groups reach the degree of adult-type pattern earlier than others. During the breast-feeding period, the microbiota appeared to be relatively simple, while major diversification was found to start during the weaning process. By the age of 3 years, the child’s microbiota composition started to resemble that of an adult, but the bacterial diversity has still not reached the full diversity, indicating that the microbiota maturation extends beyond this age. In addition, at three years of age, the child’s microbiota was more similar to mother’s microbiota than to microbiota of nonrelated women.In infants with atopic eczema, a high total microbiota diversity and abundance of butyrate-producing bacteria was found to correlate with mild symptoms at 6 months. At 18 months, infants with mild eczema had significantly higher microbiota diversity and aberrant microbiota composition when compared to healthy controls at the same age. In conclusion, the comprehensive phylogenetic microarray analysis of early life microbiota shows the synergetic effect of vertical transmission and shared environment on the intestinal microbiota development. By the age of three years, the compositional development of intestinal microbiota is close to adult level, but the microbiota diversification continues beyond this age. In addition, specific microbiota signatures are associated with the existence and severity of atopic eczema and intestinal microbiota seems to have a role in alleviating the symptoms of this disease.

Rapid, reliable and inexpensive quality assessment of biotinylated cRNA

The interpretation of oligonucleotide array experiments depends on the quality of the target cRNA used. cRNA target quality is assessed by quantitative analysis of the representation of 5' and 3' sequences of control genes using commercially available Test arrays. The Test array provides an economically priced means of determining the quality of labeled target prior to analysis on whole genome expression arrays. This manuscript validates the use of a duplex RT-PCR assay as a faster (6 h) and less expensive (microarrays from Affymetrix. Ten cRNA samples with a ß-actin 3'/5' ratio >6 were chosen and classified as degraded cRNAs, and 31 samples with a ß-actin 3'/5' ratio <6 were selected as good quality cRNAs. Blinded samples were then used for the RT-PCR assay. After gel electrophoresis, optical densities of the amplified 3' and 5' fragments of ß-actin were measured and the 3'/5' ratio was calculated. There was a strong correlation (r² = 0.6802) between the array and the RT-PCR ß-actin 3'/5' ratios. Moreover, the RT-PCR 3'/5' ratio was significantly different (P < 0.0001) between undegraded (mean ± SD, 0.34 ± 0.09) and degraded (1.71 ± 0.83) samples. None of the other parameters analyzed, such as i) the starting amount of RNA, ii) RNA quality assessed using the Bioanalyzer Chip technology, or iii) the concentration and OD260/OD280 ratio of the purified biotinylated cRNA, correlated with cRNA quality.

Gene expression profiling of clinical stages II and III breast cancer

Clinical stage (CS) is an established indicator of breast cancer outcome. In the present study, a cDNA microarray platform containing 692 genes was used to identify molecular differences between CSII and CSIII disease. Tumor samples were collected from patients with CSII or CSIII breast cancer, and normal breast tissue was collected from women without invasive cancer. Seventy-eight genes were deregulated in CSIII tumors and 22 in CSII tumors when compared to normal tissue, and 20 of them were differentially expressed in both CSII and CSIII tumors. In addition, 58 genes were specifically altered in CSIII and expression of 6 of them was tested by real time RT-PCR in another cohort of patients with CSII or CSIII breast cancer and in women without cancer. Among these genes, MAX, KRT15 and S100A14, but not APOBEC3G or KRT19, were differentially expressed on both CSIII and CSII tumors as compared to normal tissue. Increased HMOX1 levels were detected only in CSIII tumors and may represent a molecular marker of this stage. A clear difference in gene expression pattern occurs at the normal-to-cancer transition; however, most of the differentially expressed genes are deregulated in tumors of both CS (II and III) compared to normal breast tissue.

Dengue virus-induced regulation of the host cell translational machinery

Dengue virus (DV)-induced changes in the host cell protein synthesis machinery are not well understood. We investigated the transcriptional changes related to initiation of protein synthesis. The human hepatoma cell line, HepG2, was infected with DV serotype 2 for 1 h at a multiplicity of infection of one. RNA was extracted after 6, 24 and 48 h. Microarray results showed that 36.5% of the translation factors related to initiation of protein synthesis had significant differential expression (Z-score ≥ ±2.0). Confirmation was obtained by quantitative real-time reverse transcription-PCR. Of the genes involved in the activation of mRNA for cap-dependent translation (eIF4 factors), eIF4A, eIF4G1 and eIF4B were up-regulated while the negative regulator of translation eIF4E-BP3 was down-regulated. This activation was transient since at 24 h post-infection levels were not significantly different from control cells. However, at 48 h post-infection, eIF4A, eIF4E, eIF4G1, eIF4G3, eIF4B, and eIF4E-BP3 were down-regulated, suggesting that cap-dependent translation could be inhibited during the progression of infection. To test this hypothesis, phosphorylation of p70S6K and 4E-BP1, which induce cap-dependent protein synthesis, was assayed. Both proteins remained phosphorylated when assayed at 6 h after infection, while infection induced dephosphorylation of p70S6K and 4E-BP1 at 24 and 48 h of infection, respectively. Taken together, these results provide biological evidence suggesting that in HepG2 cells DV sustains activation of the cap-dependent machinery at early stages of infection, but progression of infection switches protein synthesis to a cap-independent process.

Cardiac gene expression and systemic cytokine profile are complementary in a murine model of post-ischemic heart failure

After myocardial infarction (MI), activation of the immune system and inflammatory mechanisms, among others, can lead to ventricular remodeling and heart failure (HF). The interaction between these systemic alterations and corresponding changes in the heart has not been extensively examined in the setting of chronic ischemia. The main purpose of this study was to investigate alterations in cardiac gene and systemic cytokine profile in mice with post-ischemic HF. Plasma was tested for IgM and IgG anti-heart reactive repertoire and inflammatory cytokines. Heart samples were assayed for gene expression by analyzing hybridization to AECOM 32k mouse microarrays. Ischemic HF significantly increased the levels of total serum IgM (by 5.2-fold) and total IgG (by 3.6-fold) associated with a relatively high content of anti-heart specificity. A comparable increase was observed in the levels of circulating pro-inflammatory cytokines such as IL-1β (3.8X) and TNF-α (6.0X). IFN-γ was also increased by 3.1-fold in the MI group. However, IL-4 and IL-10 were not significantly different between the MI and sham-operated groups. Chemokines such as MCP-1 and IL-8 were 1.4- and 13-fold increased, respectively, in the plasma of infarcted mice. We identified 2079 well annotated unigenes that were significantly regulated by post-ischemic HF. Complement activation and immune response were among the most up-regulated processes. Interestingly, 21 of the 101 quantified unigenes involved in the inflammatory response were significantly up-regulated and none were down-regulated. These data indicate that post-ischemic heart remodeling is accompanied by immune-mediated mechanisms that act both systemically and locally.

Genomic alterations detected by comparative genomic hybridization in ovarian endometriomas

Endometriosis is a complex and multifactorial disease. Chromosomal imbalance screening in endometriotic tissue can be used to detect hot-spot regions in the search for a possible genetic marker for endometriosis. The objective of the present study was to detect chromosomal imbalances by comparative genomic hybridization (CGH) in ectopic tissue samples from ovarian endometriomas and eutopic tissue from the same patients. We evaluated 10 ovarian endometriotic tissues and 10 eutopic endometrial tissues by metaphase CGH. CGH was prepared with normal and test DNA enzymatically digested, ligated to adaptors and amplified by PCR. A second PCR was performed for DNA labeling. Equal amounts of both normal and test-labeled DNA were hybridized in human normal metaphases. The Isis FISH Imaging System V 5.0 software was used for chromosome analysis. In both eutopic and ectopic groups, 4/10 samples presented chromosomal alterations, mainly chromosomal gains. CGH identified 11q12.3-q13.1, 17p11.1-p12, 17q25.3-qter, and 19p as critical regions. Genomic imbalances in 11q, 17p, 17q, and 19p were detected in normal eutopic and/or ectopic endometrium from women with ovarian endometriosis. These regions contain genes such as POLR2G, MXRA7 and UBA52 involved in biological processes that may lead to the establishment and maintenance of endometriotic implants. This genomic imbalance may affect genes in which dysregulation impacts both eutopic and ectopic endometrium.

Gene trio signatures as molecular markers to predict response to doxorubicin cyclophosphamide neoadjuvant chemotherapy in breast cancerpatients

In breast cancer patients submitted to neoadjuvant chemotherapy (4 cycles of doxorubicin and cyclophosphamide, AC), expression of groups of three genes (gene trio signatures) could distinguish responsive from non-responsive tumors, as demonstrated by cDNA microarray profiling in a previous study by our group. In the current study, we determined if the expression of the same genes would retain the predictive strength, when analyzed by a more accessible technique (real-time RT-PCR). We evaluated 28 samples already analyzed by cDNA microarray, as a technical validation procedure, and 14 tumors, as an independent biological validation set. All patients received neoadjuvant chemotherapy (4 AC). Among five trio combinations previously identified, defined by nine genes individually investigated (BZRP, CLPTM1,MTSS1, NOTCH1, NUP210, PRSS11, RPL37A, SMYD2, and XLHSRF-1), the most accurate were established by RPL37A, XLHSRF-1based trios, with NOTCH1 or NUP210. Both trios correctly separated 86% of tumors (87% sensitivity and 80% specificity for predicting response), according to their response to chemotherapy (82% in a leave-one-out cross-validation method). Using the pre-established features obtained by linear discriminant analysis, 71% samples from the biological validation set were also correctly classified by both trios (72% sensitivity; 66% specificity). Furthermore, we explored other gene combinations to achieve a higher accuracy in the technical validation group (as a training set). A new trio, MTSS1, RPL37 and SMYD2, correctly classified 93% of samples from the technical validation group (95% sensitivity and 80% specificity; 86% accuracy by the cross-validation method) and 79% from the biological validation group (72% sensitivity and 100% specificity). Therefore, the combined expression of MTSS1, RPL37 and SMYD2, as evaluated by real-time RT-PCR, is a potential candidate to predict response to neoadjuvant doxorubicin and cyclophosphamide in breast cancer patients.

Transcriptional regulation of bidirectional gene pairs by 17-β-estradiol in MCF-7 breast cancer cells

Using cDNA microarray analysis, we previously identified a set of differentially expressed genes in primary breast tumors based on the status of estrogen and progesterone receptors. In the present study, we performed an integrated computer-assisted and manual search of potential estrogen response element (ERE) binding sites in the promoter region of these genes to characterize their potential to be regulated by estrogen receptors (ER). Publicly available databases were used to annotate the position of these genes in the genome and to extract a 5’flanking region 2 kb upstream to 2 kb downstream of the transcription start site for transcription binding site analysis. The search for EREs and other binding sites was performed using several publicly available programs. Overall, approximately 40% of the genes analyzed were potentially able to be regulated by estrogen via ER. In addition, 17% of these genes are located very close to other genes organized in a head-to-head orientation with less than 1.0 kb between their transcript units, sharing a bidirectional promoter, and could be classified as bidirectional gene pairs. Using quantitative real-time PCR, we further investigated the effects of 17β-estradiol and antiestrogens on the expression of the bidirectional gene pairs in MCF-7 breast cancer cells. Our results showed that some of these gene pairs, such as TXNDC9/EIF5B, GALNS/TRAPPC2L, and SERINC1/PKIB, are modulated by 17β-estradiol via ER in MCF-7 breast cancer cells. Here, we also characterize the promoter region of potential ER-regulated genes and provide new information on the transcriptional regulation of bidirectional gene pairs.

Identification of microRNAs and mRNAs associated with multidrug resistance of human laryngeal cancer Hep-2 cells

Multidrug resistance (MDR) poses a serious impediment to the success of chemotherapy for laryngeal cancer. To identify microRNAs and mRNAs associated with MDR of human laryngeal cancer Hep-2 cells, we developed a multidrug-resistant human laryngeal cancer subline, designated Hep-2/v, by exposing Hep-2 cells to stepwise increasing concentrations of vincristine (0.02-0.96'µM). Microarray assays were performed to compare the microRNA and mRNA expression profiles of Hep-2 and Hep-2/v cells. Compared to Hep-2 cells, Hep-2/v cells were more resistant to chemotherapy drugs (∼45-fold more resistant to vincristine, 5.1-fold more resistant to cisplatin, and 5.6-fold more resistant to 5-fluorouracil) and had a longer doubling time (42.33±1.76 vs 28.75±1.12'h, P<0.05), higher percentage of cells in G0/G1 phase (80.98±0.52 vs69.14±0.89, P<0.05), increased efflux of rhodamine 123 (95.97±0.56 vs 12.40±0.44%, P<0.01), and up-regulated MDR1 expression. A total of 7 microRNAs and 605 mRNAs were differentially expressed between the two cell types. Of the differentially expressed mRNAs identified, regulator of G-protein signaling 10, high-temperature requirement protein A1, and nuclear protein 1 were found to be the putative targets of the differentially expressed microRNAs identified. These findings may open a new avenue for clarifying the mechanisms responsible for MDR in laryngeal cancer.

Immunohistochemical detection of virus through its nuclear cytopathic effect in idiopathic interstitial pneumonia other than acute exacerbation

Idiopathic interstitial pneumonias include complex diseases that have a strong interaction between genetic makeup and environmental factors. However, in many cases, no infectious agent can be demonstrated, and these clinical diseases rapidly progress to death. Theoretically, idiopathic interstitial pneumonias could be caused by the Epstein-Barr virus, cytomegalovirus, adenovirus, hepatitis C virus, respiratory syncytial virus, and herpesvirus, which may be present in such small amounts or such configuration that routine histopathological analysis or viral culture techniques cannot detect them. To test the hypothesis that immunohistochemistry provides more accurate results than the mere histological demonstration of viral inclusions, this method was applied to 37 open lung biopsies obtained from patients with idiopathic interstitial pneumonias. As a result, immunohistochemistry detected measles virus and cytomegalovirus in diffuse alveolar damage-related histological patterns of acute exacerbation of idiopathic pulmonary fibrosis and nonspecific interstitial pneumonia in 38 and 10% of the cases, respectively. Alveolar epithelium infection by cytomegalovirus was observed in 25% of organizing pneumonia patterns. These findings were coincident with nuclear cytopathic effects but without demonstration of cytomegalovirus inclusions. These data indicate that diffuse alveolar damage-related cytomegalovirus or measles virus infections enhance lung injury, and a direct involvement of these viruses in diffuse alveolar damage-related histological patterns is likely. Immunohistochemistry was more sensitive than the histological demonstration of cytomegalovirus or measles virus inclusions. We concluded that all patients with diffuse alveolar damage-related histological patterns should be investigated for cytomegalovirus and measles virus using sensitive immunohistochemistry in conjunction with routine procedures.

Whole transcriptome organisation in the dehydrated supraoptic nucleus

The supraoptic nucleus (SON) is part of the central osmotic circuitry that synthesises the hormone vasopressin (Avp) and transports it to terminals in the posterior lobe of the pituitary. Following osmotic stress such as dehydration, this tissue undergoes morphological, electrical and transcriptional changes to facilitate the appropriate regulation and release of Avp into the circulation where it conserves water at the level of the kidney. Here, the organisation of the whole transcriptome following dehydration is modelled to fit Zipf's law, a natural power law that holds true for all natural languages, that states if the frequency of word usage is plotted against its rank, then the log linear regression of this is -1. We have applied this model to our previously published euhydrated and dehydrated SON data to observe this trend and how it changes following dehydration. In accordance with other studies, our whole transcriptome data fit well with this model in the euhydrated SON microarrays, but interestingly, fit better in the dehydrated arrays. This trend was observed in a subset of differentially regulated genes and also following network reconstruction using a third-party database that mines public data. We make use of language as a metaphor that helps us philosophise about the role of the whole transcriptome in providing a suitable environment for the delivery of Avp following a survival threat like dehydration.