The Gerbera cDNA Microarray: A Tool for Large-Scale Identification of Genes Involved in Flower Development

During the past ten years, large-scale transcript analysis using microarrays has become a powerful tool to identify and predict functions for new genes. It allows simultaneous monitoring of the expression of thousands of genes and has become a routinely used tool in laboratories worldwide. Microarray analysis will, together with other functional genomics tools, take us closer to understanding the functions of all genes in genomes of living organisms. Flower development is a genetically regulated process which has mostly been studied in the traditional model species Arabidopsis thaliana, Antirrhinum majus and Petunia hybrida. The molecular mechanisms behind flower development in them are partly applicable in other plant systems. However, not all biological phenomena can be approached with just a few model systems. In order to understand and apply the knowledge to ecologically and economically important plants, other species also need to be studied. Sequencing of 17 000 ESTs from nine different cDNA libraries of the ornamental plant Gerbera hybrida made it possible to construct a cDNA microarray with 9000 probes. The probes of the microarray represent all different ESTs in the database. From the gerbera ESTs 20% were unique to gerbera while 373 were specific to the Asteraceae family of flowering plants. Gerbera has composite inflorescences with three different types of flowers that vary from each other morphologically. The marginal ray flowers are large, often pigmented and female, while the central disc flowers are smaller and more radially symmetrical perfect flowers. Intermediate trans flowers are similar to ray flowers but smaller in size. This feature together with the molecular tools applied to gerbera, make gerbera a unique system in comparison to the common model plants with only a single kind of flowers in their inflorescence. In the first part of this thesis, conditions for gerbera microarray analysis were optimised including experimental design, sample preparation and hybridization, as well as data analysis and verification. Moreover, in the first study, the flower and flower organ-specific genes were identified. After the reliability and reproducibility of the method were confirmed, the microarrays were utilized to investigate transcriptional differences between ray and disc flowers. This study revealed novel information about the morphological development as well as the transcriptional regulation of early stages of development in various flower types of gerbera. The most interesting finding was differential expression of MADS-box genes, suggesting the existence of flower type-specific regulatory complexes in the specification of different types of flowers. The gerbera microarray was further used to profile changes in expression during petal development. Gerbera ray flower petals are large, which makes them an ideal model to study organogenesis. Six different stages were compared and specifically analysed. Expression profiles of genes related to cell structure and growth implied that during stage two, cells divide, a process which is marked by expression of histones, cyclins and tubulins. Stage 4 was found to be a transition stage between cell division and expansion and by stage 6 cells had stopped division and instead underwent expansion. Interestingly, at the last analysed stage, stage 9, when cells did not grow any more, the highest number of upregulated genes was detected. The gerbera microarray is a fully-functioning tool for large-scale studies of flower development and correlation with real-time RT-PCR results show that it is also highly sensitive and reliable. Gene expression data presented here will be a source for gene expression mining or marker gene discovery in the future studies that will be performed in the Gerbera Laboratory. The publicly available data will also serve the plant research community world-wide.

BRCA1 and BRCA2 germline mutation analysis among Indian women from south India: identification of four novel mutations and high-frequency occurrence of 185delAG mutation

Mutations in the BRCA1 and BRCA2 genes profoundly increase the risk of developing breast and/or ovarian cancer among women. To explore the contribution of BRCA1 and BRCA2 mutations in the development of hereditary breast cancer among Indian women, we carried out mutation analysis of the BRCA1 and BRCA2 genes in 61 breast or ovarian cancer patients from south India with a positive family history of breast and/or ovarian cancer. Mutation analysis was carried out using conformation-sensitive gel electrophoresis (CSGE) followed by sequencing. Mutations were identified in 17 patients (28.0%); 15 (24.6%) had BRCA1 mutations and two (3.28%) had BRCA2 mutations. While no specific association between BRCA1 or BRCA2 mutations with cancer type was seen, mutations were more often seen in families with ovarian cancer. While 40% (4/10) and 30.8% (4/12) of families with ovarian or breast and ovarian cancer had mutations, only 23.1% (9/39) of families with breast cancer carried mutations in the BRCA1 and BRCA2 genes. In addition, while BRCA1 mutations were found in all age groups, BRCA2 mutations were found only in the age group of <= 40 years. Of the BRCA1 mutations, there were three novel mutations (295delCA; 4213T -> A; 5267T -> G) G) and three mutations that have been reported earlier. Interestingly, 185delAG, a BRCA1 mutation which occurs at a very high frequency in Ashkenazi Jews, was found at a frequency of 16.4% (10/61). There was one novel mutation (4866insT) and one reported mutation in BRCA2. Thus, our study emphasizes the importance of mutation screening in familial breast and/or ovarian cancers, and the potential implications of these findings in genetic counselling and preventive therapy.

Identification of IDUA and WNT16 phosphorylation-related non-synonymous polymorphisms for bone mineral density in meta-analyses of genome-wide association studies

Protein phosphorylation regulates a wide variety of cellular processes. Thus, we hypothesize that single-nucleotide polymorphisms (SNPs) that may modulate protein phosphorylation could affect osteoporosis risk. Based on a previous conventional genome-wide association (GWA) study, we conducted a three-stage meta-analysis targeting phosphorylation-related SNPs (phosSNPs) for femoral neck (FN)-bone mineral density (BMD), total hip (HIP)-BMD, and lumbar spine (LS)-BMD phenotypes. In stage 1, 9593 phosSNPs were meta-analyzed in 11,140 individuals of various ancestries. Genome-wide significance (GWS) and suggestive significance were defined by α = 5.21 × 10–6 (0.05/9593) and 1.00 × 10–4, respectively. In stage 2, nine stage 1–discovered phosSNPs (based on α = 1.00 × 10–4) were in silico meta-analyzed in Dutch, Korean, and Australian cohorts. In stage 3, four phosSNPs that replicated in stage 2 (based on α = 5.56 × 10–3, 0.05/9) were de novo genotyped in two independent cohorts. IDUA rs3755955 and rs6831280, and WNT16 rs2707466 were associated with BMD phenotypes in each respective stage, and in three stages combined, achieving GWS for both FN-BMD (p = 8.36 × 10–10, p = 5.26 × 10–10, and p = 3.01 × 10–10, respectively) and HIP-BMD (p = 3.26 × 10–6, p = 1.97 × 10–6, and p = 1.63 × 10–12, respectively). Although in vitro studies demonstrated no differences in expressions of wild-type and mutant forms of IDUA and WNT16B proteins, in silico analyses predicts that WNT16 rs2707466 directly abolishes a phosphorylation site, which could cause a deleterious effect on WNT16 protein, and that IDUA phosSNPs rs3755955 and rs6831280 could exert indirect effects on nearby phosphorylation sites. Further studies will be required to determine the detailed and specific molecular effects of these BMD-associated non-synonymous variants. © 2015 American Society for Bone and Mineral Research.

Identification and functional characterization of a cis-acting positive DNA element regulating CYP 2B1/B2 gene transcription in rat liver

A positive cis-acting DNA element in the near 5'-upstream region of the CYP2B1/B2 genes in rat liver was found to play an important role in the transcription of these genes. An oligonucleotide covering -69 to -98 nt mimicked the gel mobility shift pattern given by the fragment -179 to +29 nt, which was earlier found adequate to confer the regulatory features of this gene. Two major complexes were seen, of which the slower and faster moving complexes became intense under uninduced and Phenobarbitone-induced conditions respectively. Minigene cloned DNA plasmid covering -179 to +181 nt in pUC 19 and Bal 31 mutants derived from this parent were transcribed in whole nuclei and cell free transcription extracts and mutants containing only upto -75 nt of the upstream were poorly transcribed. Transcription extracts from phenobarbitone-injected rat liver nuclei were significantly more active than extracts from uninduced rats in transcribing the minigene constructs. Addition of the oligonucleotide (-69 to -98nt) specifically inhibited the transcription of the minigene construct (-179 to +181 nt) in the cell free transcription system. It is therefore, concluded that the region -69 to -98 nt acts as a positive cis-acting element in the transcription of the CYP2B1/B2 genes and in mediating the inductive effects of phenobarbitone.

Mapping of assembled epitopes with microgram quantities of antigen: Identification of an epitope at the receptor binding region of human follicle stimulating hormone

Identification of epitopes by modification studies has been reported by us recently. The method requires milligram quantities of antigen and since several proteins are not available in large quantities they are not amenable for such an investigation. One such protein is human follicle stimulating hormone (hFSH) whose mapping of epitopes is of importance in reproductive biology. Here we report a method that uses microgram quantities of hFSH to map a beta-specific epitope located at the receptor binding region. This identification has also been validated by the chemical modification method using heterologous antigen ovine follicle stimulating hormone (oFSH).

Golgi proteomics : Identification of a novel cartilage-specific Golgi protein GoPro49

The Golgi complex is a central organelle of the secretory pathway, responsible for a range of post-translational modifications, as well as for membrane traffic to the plasma membrane and to the endosomal-lysosomal pathway. In addition, this organelle has roles in cell migration, in the regulation of traffic, and as a mitotic check point. The structure of the Golgi complex is highly dynamic and able to respond to the amount of cargo being transported and the stage of the cell cycle. The Golgi proteome reflects the functions and structure of this organelle, and can be divided into three major groups: the Golgi resident proteins (e.g. modification enzymes), the Golgi matrix proteins (involved in structure and tethering events), and trafficking proteins (e.g. vesicle coat proteins and Rabs). The Golgi proteome has been studied on several occasions, from both rat liver and mammary gland Golgi membranes using proteomic approaches, but still little more than half of the estimated Golgi proteome is known. Nevertheless, methodological improvements and introduction of shotgun proteomics have increased the number of identified proteins, and especially the number of identified transmembrane proteins. Cartilage, even though not a typical tissue in which to study membrane traffic, secretes large amounts of extracellular matrix proteins that are extensively modified, especially by amino acid hydroxylation, glycosylation and sulfation. Furthermore, the cartilage ECM contains several, large oligomeric proteins (such as collagen II) that are difficult to assemble and transport. Indeed, cartilage has been shown to be susceptible to changes both in secretory pathway (e.g. the COPII coat assembly) and in post-translational modifications (e.g. heparan sulfate formation). Dental follicle, and the periodontal ligament (PDL) that it forms, are another type of connective tissue, and they have a role in anchoring teeth to bone. This anchorage is achieved by numerous matrix fibres that connect the bone matrix with the cementum. These tissues have in common the secretion of large matrix molecules. In this study the Golgi proteome was analysed from purified, stacked Golgi membranes isolated from rat liver. The identified, extensive proteome included a protein similar to Ab2-095, or Golgi protein 49kDa (GoPro49), which was shown to localise to the Golgi complex as an EGFP fusion protein. Surprisingly, in situ hybridisation showed the GoPro49 expression to be highly restricted to different mesenchymal tissues, especially in cartilage, and this expression pattern was clearly developmentally regulated. In addition to cartilage, GoPro49 was also expressed in the dental follicle, but was not observed in the mature PDL. Importantly, GoPro49 is the first specific marker for the dental follicle. Endogenous GoPro49 protein co-localised with β-COP in both chondrosarcoma and primary dental follicle cell lines. The COPI staining in these cells was highly dynamic, showing a number of tubules. This may reflect the type of secretory cargo they secrete. Currently GoPro49 is the only Golgi protein with such a restricted expression pattern.

Identification of the fungal catabolic D-galacturonate pathway

Pectin is a natural polymer consisting mainly of D-galacturonic acid monomers. Microorganisms living on decaying plant material can use D-galacturonic acid for growth. Although bacterial pathways for D-galacturonate catabolism had been described previously, no eukaryotic pathway for D-galacturonate catabolism was known at the beginning of this work. The aim of this work was to identify such a pathway. In this thesis the pathway for D-galacturonate catabolism was identified in the filamentous fungus Trichoderma reesei. The pathway consisted of four enzymes: NADPH-dependent D-galacturonate reductase (GAR1), L-galactonate dehydratase (LGD1), L-threo-3-deoxy-hexulosonate aldolase (LGA1) and NADPH-dependent glyceraldehyde reductase (GLD1). In this pathway D-galacturonate was converted to pyruvate and glycerol via L-galactonate, L-threo-3-deoxy-hexulosonate and L-glyceraldehyde. The enzyme activities of GAR1, LGD1 and LGA1 were present in crude mycelial extract only when T. reesei was grown on D-galacturonate. The activity of GLD1 was equally present on all the tested carbon sources. The corresponding genes were identified either by purifying and sequencing the enzyme or by expressing genes with homology to other similar enzymes in a heterologous host and testing the activities. The new genes that were identified were expressed in Saccharomyces cerevisiae and resulted in active enzymes. The GAR1, LGA1 and GLD1 were also produced in S. cerevisiae as active enzymes with a polyhistidine-tag, and purified and characterised. GAR1 and LGA1 catalysed reversible reactions, whereas only the forward reactions were observed for LGD1 and GLD1. When gar1, lgd1 or lga1 was deleted in T. reesei the deletion strain was unable to grow with D-galacturonate as the only carbon source, demonstrating that all the corresponding enzymes were essential for D-galacturonate catabolism and that no alternative D-galacturonate pathway exists in T. reesei. A challenge for biotechnology is to convert cheap raw materials to useful and more valuable products. Filamentous fungi are especially useful for the conversion of pectin, since they are efficient producers of pectinases. Identification of the fungal D-galacturonate pathway is of fundamental importance for the utilisation of pectin and its conversion to useful products.

Solution NMR studies of acetohydroxy acid synthase I: Identification of the sites of inter-subunit interactions using multidimensional NMR methods

The novel multidomain organization in the multimeric Escherichia coli AHAS I (ilvBN) enzyme has been dissected to generate polypeptide fragments. These fragments when cloned, expressed and purified reassemble in the presence of cofactors to yield a catalytically competent enzyme. Structural characterization of AHAS has been impeded due to the fact that the holoenzyme is prone to dissociation leading to heterogeneity in samples. Our approach has enabled the structural characterization using high-resolution nuclear magnetic resonance methods. Near complete sequence specific NMR assignments for backbone H-N, N-15, C-13 alpha and C-13(beta) atoms of the FAD binding domain of ilvB have been obtained on samples isotopically enriched in H-2, C-13 and N-15. The secondary structure determined on the basis of observed C-13(alpha) secondary chemical shifts and sequential NOEs indicates that the secondary structure of the FAD binding domain of E. coli AHAS large Subunit (ilvB) is similar to the structure of this domain in the catalytic subunit of yeast AHAS. Protein-protein interactions involving the regulatory subunit (ilvN) and the domains of the catalytic subunit (ilvB) were studied using circular dichroic and isotope edited solution nuclear magnetic resonance spectroscopic methods. Observed changes in circular dichroic spectra indicate that the regulatory subunit (ilvN) interacts with ilvB alpha and ilvB beta domains of the catalytic subunit and not with the ilvB gamma domain. NMR chemical shift mapping methods show that ilvN binds close to the FAD binding site in ilvB beta and proximal to the intrasubunit ilvB alpha/ilvB beta domain interface. The implication of this interaction on the role of the regulatory subunit oil the activity of the holoenzyme is discussed. NMR studies of the regulatory domains show that these domains are structured in solution. Preliminary evidence for the interaction of ilvN with the metabolic end product of the pathway, viz., valine is also presented.

A semi-analytical particle filter for identification of nonlinear oscillators

Particle filters find important applications in the problems of state and parameter estimations of dynamical systems of engineering interest. Since a typical filtering algorithm involves Monte Carlo simulations of the process equations, sample variance of the estimator is inversely proportional to the number of particles. The sample variance may be reduced if one uses a Rao-Blackwell marginalization of states and performs analytical computations as much as possible. In this work, we propose a semi-analytical particle filter, requiring no Rao-Blackwell marginalization, for state and parameter estimations of nonlinear dynamical systems with additively Gaussian process/observation noises. Through local linearizations of the nonlinear drift fields in the process/observation equations via explicit Ito-Taylor expansions, the given nonlinear system is transformed into an ensemble of locally linearized systems. Using the most recent observation, conditionally Gaussian posterior density functions of the linearized systems are analytically obtained through the Kalman filter. This information is further exploited within the particle filter algorithm for obtaining samples from the optimal posterior density of the states. The potential of the method in state/parameter estimations is demonstrated through numerical illustrations for a few nonlinear oscillators. The proposed filter is found to yield estimates with reduced sample variance and improved accuracy vis-a-vis results from a form of sequential importance sampling filter.

Identification of a novel FGFRL1 MicroRNA target site polymorphism for bone mineral density in meta-analyses of genome-wide association studies

MicroRNAs (miRNAs) are critical post-transcriptional regulators. Based on a previous genome-wide association (GWA) scan, we conducted a polymorphism in microRNAs' Target Sites (poly-miRTS)-centric multistage meta-analysis for lumbar spine (LS)-, total hip (HIP)-, and femoral neck (FN)-bone mineral density (BMD). In stage I, 41,102 poly-miRTSs were meta-analyzed in 7 cohorts with a genome-wide significance (GWS) α=0.05/41,102=1.22×10-6. By applying α=5×10-5 (suggestive significance), 11 poly-miRTSs were selected, with FGFRL1 rs4647940 and PRR5 rs3213550 as top signals for FN-BMD (P-value=7.67×10-6 and 1.58×10-5) in gender-combined sample. In stage II in silico replication (two cohorts), FGFRL1 rs4647940 was the only signal marginally replicated for FN-BMD (P-value=5.08×10-3) at α=0.10/11=9.09×10-3. PRR5 rs3213550 was also selected based on biological significance. In stage III de novo genotyping replication (two cohorts), FGFRL1 rs4647940 was the only signal significantly replicated for FN-BMD (P-value=7.55×10-6) at α=0.05/2=0.025 in gender-combined sample. Aggregating three stages, FGFRL1 rs4647940 was the single stage I-discovered and stages II- and III-replicated signal attaining GWS for FN-BMD (P-value=8.87×10-12). Dual-luciferase reporter assays demonstrated that FGFRL1 3' untranslated region harboring rs4647940 appears to be hsa-miR-140-5p's target site. In a zebrafish microinjection experiment, dre-miR-140-5p is shown to exert a dramatic impact on craniofacial skeleton formation. Taken together, we provided functional evidence for a novel FGFRL1 poly-miRTS rs4647940 in a previously known 4p16.3 locus, and experimental and clinical genetics studies have shown both FGFRL1 and hsa-miR-140-5p are important for bone formation. © The Author 2015. Published by Oxford University Press. All rights reserved.

Identification of alpha- and beta-Hydroxy Acid Containing Cyclodepsipeptides in Natural Peptide Mixtures Using Negative Ion Mass Spectrometry

Natural peptide libraries often contain cyclodepsipeptides containing alpha or beta hydroxy residues. Extracts of fungal hyphae of Isaria yield a microheterogenous cyclodepsipeptide mixture in which two classes of molecules can be identified by mass spectral fragmentation of negative ions. In the case of isaridins, which contain an alpha-hydroxy residue and a beta-amino acid residue, a characteristic product ion corresponding to a neutral loss of 72 Da is obtained. hi addition, neutral loss of water followed by a 72 Da loss is also observed. Two distinct modes of fragmentation rationalize the observed product ion distribution. The neutral loss of 72 Da has also been obtained for a roseotoxin component, which is also an alpha-hydroxy residue containing cyclodepsipeptide. In the case of isariins, which contain a beta-hydroxy acid residue, ring opening and subsequent loss of the terminal residue as an unsaturated ketene fragment, rationalizes the observed product ion formation. Fragmentation of negative ions provide characteristic neutral losses, which are diagnostic of the presence of alpha-hydroxy or beta-hydroxy residues.

Identification and Characterization of a Novel Deoxyhypusine Synthase in Leishmania donovani

Deoxyhypusine synthase, an NAD(+)-dependent enzyme, catalyzes the first step in the post-translational synthesis of an unusual amino acid, hypusine (N-epsilon-(4-amino-2-hydroxybutyl)lysine), in the eukaryotic initiation factor 5A precursor protein. Two putative deoxyhypusine synthase (DHS) sequences have been identified in the Leishmania donovani genome, which are present on chromosomes 20: DHSL20 (DHS-like gene from chromosome 20) and DHS34 (DHS from chromosome 34). Although both sequences exhibit an overall conservation of key residues, DHSL20 protein lacks a critical lysine residue, and the recombinant protein showed no DHS activity in vitro. However, DHS34 contains the critical lysine residue, and the recombinant DHS34 effectively catalyzed deoxyhypusine synthesis. Furthermore, in vivo labeling confirmed that hypusination of eukaryotic initiation factor 5A occurs in intact Leishmania parasites. Interestingly, the DHS34 is much longer, with 601 amino acids, compared with the human DHS enzyme (369 amino acids) and contains several unique insertions. To study the physiological role of DHS34 in Leishmania, gene deletion mutations were attempted via targeted gene replacement. However, chromosomal null mutants of DHS34 could only be obtained in the presence of a DHS34-containing episome. The present data provide evidence that DHS34 is essential for L. donovani and that structural differences in the human and leishmanial DHS enzyme may be exploited for designing selective inhibitors against the parasite.