SPI: a tool for incorporating gene expression data into a four-dimensional database of Caenorhabditis elegans embryogenesis

MOTIVATION: A comprehensive gene expression database is essential for computer modeling and simulation of biological phenomena, including development. Development is a four-dimensional (4D; 3D structure and time course) phenomenon. We are constructing a 4D database of gene expression for the early embryogenesis of the nematode Caenorhabditis elegans. As a framework of the 4D database, we have constructed computer graphics (CG), into which we will incorporate the expression data of a number of genes at the subcellular level. However, the assignment of 3D distribution of gene products (protein, mRNA), of embryos at various developmental stages, is both difficult and tedious. We need to automate this process. For this purpose, we developed a new system, named SPI after superimposing fluorescent confocal microscopic data onto a CG framework. RESULTS: The scheme of this system comprises the following: (1) acquirement of serial sections (40 slices) of fluorescent confocal images of three colors (4',6'-diamino-2-phenylindole (DAPI) for nuclei, indodicarbocyanine (Cy-3) for the internal marker, which is a germline-specific protein POS-1 and indocarbocyanine (Cy-5) for the gene product to be examined); (2) identification of several features of the stained embryos, such as contour, developmental stage and position of the internal marker; (3) selection of CG images of the corresponding stage for template matching; (4) superimposition of serial sections onto the CG; (5) assignment of the position of superimposed gene products. The Snakes algorithm identified the embryo contour. The detection accuracy of embryo contours was 92.1% when applied to 2- to 28-cell-stage embryos. The accuracy of the developmental stage prediction method was 81.2% for 2- to 8-cell-stage embryos. We manually judged only the later stage embryos because the accuracy for embryos at the later stages was unsatisfactory due to experimental noise effects. Finally, our system chose the optimal CG and performed the superposition and assignment of gene product distribution. We established an initial 4D gene expression database with 56 maternal gene products. AVAILABILITY: This system is available at http://anti.lab.nig.ac.jp/spi/ and http://anti.lab.nig.ac.jp/4ddb/     

Mosaic gene expression in nuclear transfer-derived embryos and the production of cloned transgenic pigs from ear-derived fibroblasts

Genetically modified domestic animals have many potential applications ranging from basic research to production agriculture. One of the goals in transgenic animal production schemes is to reliably predict the expression pattern of the foreign gene. Establishing a method to screen genetically modified embryos for transgene expression before transfer to surrogates may improve the likelihood of producing offspring with the desired expression pattern. In order to determine how transgene expression may be regulated in the early embryo, we generated porcine embryos from two distinct genetically modified cell lines by using the nuclear transfer (NT) technique. Both cell lines expressed the enhanced green fluorescent protein (eGFP); the first was a fibroblast cell line derived from the skin of a newborn pig that expressed eGFP, whereas the second was a fetal derived fibroblast cell line into which the eGFP gene was introduced by a retroviral vector. The reconstructed embryos were activated by electrical pulses and cultured in NCSU23. Although the in vitro developmental ability of each group of NT embryos was not different, the eGFP expression pattern was different. All embryos produced from the transduced fetal cell line fluoresced, but only 26% of the embryos generated from the newborn cell line fluoresced, and among those that did express eGFP, more than half had a mosaic expression pattern. This was unexpected because the fetal cell line was not clonally selected, and each cell had potentially different sites of integration. Embryos generated from the newborn cell line were surgically transferred to five surrogate gilts. One gilt delivered four female piglets, all of which expressed eGFP, and all had microsatellites identical to the donor. Here we demonstrate that transgene expression in all the blastomeres of an NT embryo is not uniform. In addition, transgene expression in a genetically manipulated embryo may not be an accurate indicator of expression in the resulting offspring.     

利用斑马鱼胚胎快速鉴定真核质粒中目的基因表达的研究

目的建立利用斑马鱼胚胎快速鉴定真核质粒中目的基因表达的实验体系。方法选20枚斑马鱼受精卵,在显微镜下每隔1h记录胚胎的发育情况。另选250枚单细胞期斑马鱼胚胎,平均分成5组,一组胚胎作为对照,剩余4组分别向胚胎的单细胞内注射pEGFP-N1(真核表达质粒)、pCMV-DsRed-Express2(真核表达质粒)、pET28-GFP(原核表达质粒)、pET28-RFP(原核表达质粒)质粒,在不同时间点连续观察绿色荧光及红色荧光的表达情况。另选600枚单细胞期斑马鱼胚胎,平均分成3组,一组胚胎作为对照,一组向胚胎单细胞内注射pEGFP-N1质粒,另外一组向胚胎单细胞内注射pEGFP-N1-MUC1外源基因融合重组质粒,注射4h后在荧光显微镜下观察绿色荧光的表达情况,并用RT-PCR的方法检测目的基因MUC1mRNA的转录情况。结果注射pEGFP-N1、pCMV-DsRed-Express2真核表达质粒的胚胎,注射4h后分别观察到很强的绿色荧光及红色荧光;注射pET28-GFP、pET28-RFP原核表达质粒的胚胎,10h内都未观察到绿色荧光及红色荧光;注射pEGFP-N1-MUC1外源基因融合质粒,注射4h后同样...     

根据GO分类体系选择与实验条件相关的差异表达功能结点进行基因功能深化预测

利用有限个实验条件下的基因表达谱数据,只能对与实验条件相关的基因功能类进行有效预测,所以有必要限定可预测的基因功能类范围。据此,首先基于GeneOntology(GO)选择富集差异表达基因与实验条件相关的功能类。再通过支持向量机分类器,深化预测迄今只注释到实验条件相关功能类的父结点的基因是否属于该实验条件相关功能类。应用于一套酵母基因表达谱数据,结果显示,在剔除了高度不平衡的训练集合后,平均真阳性率(precision)与平均覆盖率(recall)都分别达到了71%与47%以上。     

A method for two-dimensional registration and construction of the two-dimensional atlas of gene expression patterns in situ

We apply the fast redundant dyadic wavelet transform to the spatial registration of two-dimensional gene expression patterns of 736 Drosophila melanogaster embryos. This method is superior to the Fourier transform or windowed Fourier transform because of its ability to reduce noise and is of high resolution. In registration of the dataset we use two cost functions based on computing the Euclidean or Mahalanobis distance. The algorithm shows a high level of accuracy. For early temporal classes the cost function based on Mahalanobis distance gives better results. We have reported a method for construction of an integrated dataset elsewhere. In this paper the method is extended to the two-dimensional case. The procedure for data assembly provides for the preservation of some aspects of the nuclear structure of a two-dimensional gene expression pattern. It is based on creating an averaged model that reproduces the spatial distribution of nuclei over the embryo image. The average concentrations of each protein in each averaged nucleus are computed from the series of embryos of the same age.     

Specific gene silencing in the pre-implantation stage mouse embryo by an siRNA expression vector system

Recently, small interfering RNAs (siRNAs) have become a powerful and widely used tool for the analysis of gene function in mammalian cells. Here we report that the microinjection of an siRNA expression vector into the nucleus is an efficient and powerful method of specific gene silencing in pre-implantation mouse embryos. We used this method to examine the expression of two genes EGFP and Oct4. Vectors encoding siRNAs targeted against EGFP or Oct4 were injected into the pronucleus or nucleus of zygotes, which were then cultured until the blastocyst stage. When the effects of RNAi were examined in blastocyst stage eggs, there was robust inhibition of the gene product in a concentration-dependent manner at both the mRNA and the protein level. The expression of other endogenous genes was not affected, showing the specificity of the vector-mediated RNAi. In addition, this method was effective for inhibiting maternally expressed mRNA. To demonstrate that RNAi of Oct4 induced a similar phenotype to that of Oct4-null embryos, the blastocysts were further cultured in ES medium. After the fourth day of culture, the embryos either had outgrown only a layer of trophoblast cells or showed developmental arrest at the blastocyst stage (>90%). Moreover, concomitant with Oct4 suppression at the blastocyst stage, we observed inhibition of Fgf4, a gene that is known to be induced downstream of Oct4 expression. Taken together, these results demonstrate that the use of siRNA expression vector is a powerful way to achieve gene silencing in the pre-implantation stage embryo.     

Gene expression of bovine embryos developing at the air-liquid interface on oviductal epithelial cells (ALI-BOEC)

We recently developed an air-liquid interface long-term culture of differentiated bovine oviductal epithelial cells (ALI-BOEC). This ex vivo oviduct epithelium is capable of supporting embryo development in co-culture up to the blastocyst stage without addition of embryo culture medium. However, blastocyst rates in co-culture were markedly lower than in conventional in vitro embryo production procedures. In the present study, we assessed target gene expression of ALI-BOEC derived embryos to test their similarity to embryos from conventional in vitro embryo culture. We screened previously published data from developing bovine embryos and selected 41 genes which are either differentially expressed during embryo development, or reflect differences between various in vitro culture conditions or in vitro and in vivo embryos. Target gene expression was measured in 8-cell embryos and blastocysts using a 48.48 Dynamic Array? on a Biomark HD instrument. For comparison with the ALI-BOEC system, we generated embryos by two different standard IVP protocols. The culture conditions lead to differential gene expression in both 8-cell embryos and blastocysts. Across the expression of all target genes the embryos developing on ALI-BOEC did not depart from conventional IVP embryos. These first results prove that gene expression in ALI-BOEC embryos is not largely aberrant. However, there was no clear indication for a more in vivo-like target gene expression of these embryos. This calls for further optimization of the ALI-BOEC system to increase its efficiency both quantitatively and qualitatively.     

pRRophetic: An R Package for Prediction of Clinical Chemotherapeutic Response from Tumor Gene Expression Levels

We recently described a methodology that reliably predicted chemotherapeutic response in multiple independent clinical trials. The method worked by building statistical models from gene expression and drug sensitivity data in a very large panel of cancer cell lines, then applying these models to gene expression data from primary tumor biopsies. Here, to facilitate the development and adoption of this methodology we have created an R package called pRRophetic. This also extends the previously described pipeline, allowing prediction of clinical drug response for many cancer drugs in a user-friendly R environment. We have developed several other important use cases; as an example, we have shown that prediction of bortezomib sensitivity in multiple myeloma may be improved by training models on a large set of neoplastic hematological cell lines. We have also shown that the package facilitates model development and prediction using several different classes of data.     

Integration of Steady-State and Temporal Gene Expression Data for the Inference of Gene Regulatory Networks

We develop a new regression algorithm, cMIKANA, for inference of gene regulatory networks from combinations of steady-state and time-series gene expression data. Using simulated gene expression datasets to assess the accuracy of reconstructing gene regulatory networks, we show that steady-state and time-series data sets can successfully be combined to identify gene regulatory interactions using the new algorithm. Inferring gene networks from combined data sets was found to be advantageous when using noisy measurements collected with either lower sampling rates or a limited number of experimental replicates. We illustrate our method by applying it to a microarray gene expression dataset from human umbilical vein endothelial cells (HUVECs) which combines time series data from treatment with growth factor TNF and steady state data from siRNA knockdown treatments. Our results suggest that the combination of steady-state and time-series datasets may provide better prediction of RNA-to-RNA interactions, and may also reveal biological features that cannot be identified from dynamic or steady state information alone. Finally, we consider the experimental design of genomics experiments for gene regulatory network inference and show that network inference can be improved by incorporating steady-state measurements with time-series data.     

A new method to efficiently produce transgenic embryos and mice from low-titer lentiviral vectors

Vector injection into the perivitelline space has emerged as the standard delivery method to transduce lentivirus to mammalian oocytes or one-cell embryos, but its application is limited by the need for high titers of lentivirus. Herein we developed a new method by using a Piezo impact micro-manipulator for injecting low titer of lentivirus into the subzonal space of two-cell embryos or the perivitelline space of one-cell embryos that were shrunk with a highly concentrated sucrose solution. The survival rate of embryos was greater than 98% using this micromanipulation strategy, which was increased compared to the normal one-cell embryo injection method. More than 90% of injected embryos were GFP positive after subzonal injection of a lentivirus vector carrying the GFP gene with titers of 2 × 10⁸ I.U./ml. Even when a low titer of lentivirus (2 × 10⁶ I.U./ml) was used, 53.26% and 40.85% transgenic embryos were obtained after two-cell embryonic injection and one-cell sucrose treated embryonic injection, respectively. The GFP-positive rates were also greater than in the conventional method of injecting one-cell embryos (25.39%). In addition, blastocysts from the two-cell embryo injection group displayed stronger GFP fluorescence than the one-cell embryo injection groups treated with or without the sucrose solution. Increased expression of GFP suggests that the embryos obtained from this injection method have higher exogenous gene expression levels compared to previous methods. Therefore, in contrast with the traditional injection method, we have demonstrated a simplified and efficient lentivirus-mediated gene transfer method based on a low-titer virus preparation.     

Predicting protein submitochondria locations by combining different descriptors into the general form of Chou’s pseudo amino acid composition

Knowledge of the submitochondria location of protein is integral to understanding its function and a necessity in the proteomics era. In this work, a new submitochondria data set is constructed, and an approach for predicting protein submitochondria locations is proposed by combining the amino acid composition, dipeptide composition, reduced physicochemical properties, gene ontology, evolutionary information, and pseudo-average chemical shift. The overall prediction accuracy is 93.57% for the submitochondria location and 97.79% for the three membrane protein types in the mitochondria inner membrane using the algorithm of the increment of diversity combined with the support vector machine. The performance of the pseudo-average chemical shift is excellent. For contrast, the method is also used to predict submitochondria locations in the data set constructed by Du and Li; an accuracy of 94.95% is obtained by our method, which is better than that of other existing methods.     

Gene manipulation of chick embryos in vitro, early chick culture, and long survival in transplanted eggs

We introduce a revolutionary gene transfer system in chick: transfect chick embryos at early developmental stage by electroporation in vitro, Early Chick (EC) culture, and transplant to the egg to let the embryo survive until E5.5. Referring to the fate map, we could target the tissues of transfection, or transfect large areas of the embryo. We could get tissue-specific expression of a transgene by tissue-specific promoter. This method is very convenient and rapid, but allows us to get stable expression of the transgene in combination with transposon system.     

Registration of the expression patterns of Drosophila segmentation genes by two independent methods

MOTIVATION: To construct an integrated map of Drosophila segmentation gene expression from partial data taken from individual embryos. RESULTS: Spline and wavelet based registration techniques were developed to register Drosophila segmentation gene expression data. As ground control points for registration we used the locations of extrema on gene expression patterns, represented in 1D. The registration method was characterized by unprecedented high accuracy. A method for constructing the integrated pattern of gene expression at cellular resolution was designed. These patterns were constructed for 9 segmentation genes belonging to gap and pair-rule classes.