• Animal Bioscience
• /
• v.36 no.9
• /
• pp.1336-1349
• /
• 2023

Objective: The study was conducted to screen differentially expressed miRNAs in sows at early pregnancy by high-throughput sequencing and explore its mechanism of action on embryo implantation. Methods: The blood serum of pregnant and non-pregnant Landrace×Yorkshire sows were collected 14 days after artificial insemination, and exosomal miRNAs were purified for high throughput miRNA sequencing. The expression patterns of 10 differentially expressed (DE) miRNAs were validated by quantitative reverse transcription-polymerase chain reaction (qRT-PCR). The qRT-PCR quantified the abundance of serum exosomal miR-192 in pregnant and control sows, and the diagnostic power was assessed by receiver operating characteristic (ROC) analysis. The target genes of DE miRNAs were predicted with bioinformatics software, and the functional and pathway enrichment analysis was performed on gene ontology and the Kyoto encyclopedia of genes and genomes terms. Furthermore, a luciferase reporter system was used to identify the target relation between miR-192 and integrin alpha 4 (ITGA4), a gene influencing embryo implantation in pigs. Finally, the expression levels of miRNAs and the target gene ITGA4 were analyzed by qRT-PCR, and western blot, with the proliferation of BeWo cells detected by cell counting kit-8 (CCK-8). Results: A total of 221 known miRNAs were detected in the libraries of the pregnant and non-pregnant sows, of which 55 were up-regulated and 67 were down-regulated in the pregnant individuals compared with the non-pregnant controls. From these, the expression patterns of 10 DE miRNAs were validated. The qRT-PCR analysis further confirmed a significantly higher expression of miR-192 in the serum exosomes extracted from pregnant sows, when compared to controls. The ROC analysis revealed that miR-192 provided excellent diagnostic accuracy for pregnancy (area under the ROC curve [AUC]=0.843; p>0.001). The dual-luciferase reporter assay indicated that miR-192 directly targeted ITGA4. The protein expression of ITGA4 was reduced in cells that overexpressed miR-192. Overexpression of miR-192 resulted in the decreased proliferation of BeWo cells and regulated the expression of cell cycle-related genes. Conclusion: Serum exosomal miR-192 could serve as a potential biomarker for early pregnancy in pigs. miR-192 targeted ITGA4 gene directly, and miR-192 can regulate cellular proliferation.

• Korean Journal of Clinical Laboratory Science
• /
• v.55 no.4
• /
• pp.235-243
• /
• 2023

As next-generation sequencing has been developed and used widely, RNA-sequencing (RNA-seq) has rapidly emerged as the first choice of tools to validate global transcriptome profiling. With the significant advances in RNA-seq, various types of RNA-seq have evolved in conjunction with the progress in bioinformatic tools. On the other hand, it is difficult to interpret the complex data underlying the biological meaning without a general understanding of the types of RNA-seq and bioinformatic approaches. In this regard, this paper discusses the two main sections of RNA-seq. First, two major variants of RNA-seq are described and compared with the standard RNA-seq. This provides insights into which RNA-seq method is most appropriate for their research. Second, the most widely used RNA-seq data analyses are discussed: (1) exploratory data analysis and (2) pathway enrichment analysis. This paper introduces the most widely used exploratory data analysis for RNA-seq, such as principal component analysis, heatmap, and volcano plot, which can provide the overall trends in the dataset. The pathway enrichment analysis section introduces three generations of pathway enrichment analysis and how they generate enriched pathways with the RNA-seq dataset.

• Korean Journal of Clinical Laboratory Science
• /
• v.56 no.1
• /
• pp.10-20
• /
• 2024

RNA-sequencing (RNA-seq) is a technique used for providing global patterns of transcriptomes in samples. However, it can only provide the average gene expression across cells and does not address the heterogeneity within the samples. The advances in single-cell RNA sequencing (scRNA-seq) technology have revolutionized our understanding of heterogeneity and the dynamics of gene expression at the single-cell level. For example, scRNA-seq allows us to identify the cell types in complex tissues, which can provide information regarding the alteration of the cell population by perturbations, such as genetic modification. Since its initial introduction, scRNA-seq has rapidly become popular, leading to the development of a huge number of bioinformatic tools. However, the analysis of the big dataset generated from scRNA-seq requires a general understanding of the preprocessing of the dataset and a variety of analytical techniques. Here, we present an overview of the workflow involved in analyzing the scRNA-seq dataset. First, we describe the preprocessing of the dataset, including quality control, normalization, and dimensionality reduction. Then, we introduce the downstream analysis provided with the most commonly used computational packages. This review aims to provide a workflow guideline for new researchers interested in this field.

• Korean Journal of Poultry Science
• /
• v.37 no.3
• /
• pp.275-282
• /
• 2010

Chicken is an important livestock as a valuable biomedical model as well as food for human, and there is a strong rationale for improving our understanding on metabolism and physiology of this organism. The first draft of chicken genome assembly was released in 2004, which enables elaboration on the linkage between genetic and metabolic traits of chicken. The objectives of this study were thus to reconstruct metabolic pathway of the chicken genome and to construct a chicken specific pathway genome database (PGDB). We developed a comprehensive genome database for chicken by integrating all the known annotations for chicken genes and proteins using a pipeline written in Perl. Based on the comprehensive genome annotations, metabolic pathways of the chicken genome were reconstructed using the PathoLogic algorithm in Pathway Tools software. We identified a total of 212 metabolic pathways, 2,709 enzymes, 71 transporters, 1,698 enzymatic reactions, 8 transport reactions, and 1,360 compounds in the current chicken genome build, Gallus_gallus-2.1. Comparative metabolic analysis with the human, mouse and cattle genomes revealed that core metabolic pathways are highly conserved in the chicken genome. It was indicated the quality of assembly and annotations of the chicken genome need to be improved and more researches are required for improving our understanding on function of genes and metabolic pathways of avian species. We conclude that the chicken PGDB is useful for studies on avian and chicken metabolism and provides a platform for comparative genomic and metabolic analysis of animal biology and biomedicine.

• Asian Pacific Journal of Cancer Prevention
• /
• v.15 no.12
• /
• pp.4773-4780
• /
• 2014

Background: To investigate the effects of small interference RNA (siRNA) targeting BCR/ABL mRNA on proliferation and apoptosis in the K562 human chronic myeloid leukemia (CML) cell line and to provide a theoretical rationale and experimental evidence for its potential clinical application for anti-CML treatment. Materials and Methods: The gene sequence for BCR/ABL mRNA was found from the GeneBank. The target gene site on the BCR/ABL mRNA were selected according to Max-Planck-Institute (MPI) and rational siRNA design rules, the secondary structure of the candidate targeted mRNA was predicted, the relevant thermodynamic parameters were analyzed, and the targeted gene sequences were compared with BLAST to eliminate any sequences with significant homology. Inhibition of proliferation was evaluated by MTT assay and colony-formation inhibiting test. Apoptosis was determined by flow cytometry (FCM) and the morphology of apoptotic cells was identified by Giemsa-Wright staining. Western blotting was used to analyze the expression of BCR/ABL fusion protein in K562 cells after siRNA treatment. Results: The mRNA local secondary structure calculated by RNA structure software, and the optimal design of specific siRNA were contributed by bioinformatics rules. Five sequences of BCR/ABL siRNAs were designed and synthesized in vitro. Three sequences, siRNA1384, siRNA1276 and siRNA1786, which showed the most effective inhibition of K562 cell growth, were identified among the five candidate siRNAs, with a cell proliferative inhibitory rate nearly 50% after exposure to 12.5nmol/L~50nmol/L siRNA1384 for 24,48 and 72 hours. The 50% inhibitory concentrations ($IC_{50}$) of siRNA1384, siRNA1276 and siRNA1786 for 24hours were 46.6 nmol/L, 59.3 nmol/L and 62.6 nmol/L, respectively, and 65.668 nmol/L, 76.6 nmol/L, 74.4 nmol/L for 72 hours. The colony-formation inhibiting test also indicated that, compared with control, cell growth of siRNA treated group was inhibited. FCM results showed that the rate of cell apoptosis increased 24 hours after transfecting siRNA. The results of annexinV/PI staining indicated that the rate of apoptosis imcreased (1.53%, 15.3%, 64.5%, 57.5% and 21.5%) following treamtne with siRNAs (siRNA34, siRNA372, siRNA1384, siRNA1276 and siRNA1786). Morphological analysis showed td typical morphologic changes of apoptosis such as shrunken, fragmentation nucleus as well as "apoptotic bodies" after K562 cell exposure to siRNA. Western blot analysis showed that BCR/ABL protein was reduced sharply after a single dose of 50nmol/L siRNA transfection. Conclusions: Proliferation of K562 cells was remarkbly inhibited by siRNAs (siRNA1384, siRNA1276 and siRNA1786) in a concentration-dependent manner in vitro, with effective induction of apoptosis at a concentration of 50 nmol/L. One anti-leukemia mechanism in K562 cells appeared that BCR/ABL targeted protein was highly down-regulated. The siRNAs (siRNA1384, siRNA1276 and siRNA1786) may prove valuable in the treatment of CML.

• Journal of KIISE:Software and Applications
• /
• v.33 no.2
• /
• pp.143-153
• /
• 2006

Microarray which enables us to obtain hundreds and thousands of expression of gene or genotype at once is an epoch-making technology in comparative analysis of genes. First of all, we have to measure the intensity of each gene in an microarray image from the experiment to gain the expression level of each gene. But it is difficult to analyze the microarray image in manual because it has a lot of genes. Meta-gridding method and various auto-gridding methods have been proposed for this, but thew still have some problems. For example, meta-gridding requires manual-work due to some variations in spite of experiment in same microarray, and auto-gridding nay not carried out fully or correctly when an image has a lot of noises or is lowly expressed. In this article, we propose Hierarchical Grid Alignment algorithm for new methodology combining meta-gridding method with auto-gridding method. In our methodology, we necd a meta-grid as an input, and then align it with the microarray image automatically. Experimental results show that the proposed method serves more robust and reliable gridding result than the previous methods. It is also possible for user to do more reliable batch analysis by using our algorithm.

• Interdisciplinary Bio Central
• /
• v.1 no.1
• /
• pp.3.1-3.6
• /
• 2009

Introduction: GTPases known as translation factor play a vital role as ribosomal subunit assembly chaperone. The bacterial Obg proteins ($Spo{\underline{0B}}$-associated ${\underline{G}}TP$-binding protein) belong to the subfamily of P-loop GTPase proteins and now it is considered as one of the new target for antibacterial drug. The majority of bacterial Obgs have been commonly found to be associated with ribosome, implying that these proteins may play a fundamental role in ribosome assembly or maturation. In addition, one of the experimental evidences suggested that Bacillus subtilis Obg (BsObg) protein binds to the L13 ribosomal protein (BsL13) which is known to be one of the early assembly proteins of the 50S ribosomal subunit in Escherichia coli. In order to investigate binding mode between the BsObg and the BsL13, protein-protein docking simulation was carried out after generating 3D structure of the BsL13 structure using homology modeling method. Materials and Methods: Homology model structure of BsL13 was generated using the EcL13 crystal structure as a template. Protein-protein docking of BsObg protein with ribosomal protein BsL13 was performed by DOT, a macro-molecular docking software, in order to predict a reasonable binding mode. The solvated energy minimization calculation of the docked conformation was carried out to refine the structure. Results and Discussion: The possible binding conformation of BsL13 along with activated Obg fold in BsObg was predicted by computational docking study. The final structure is obtained from the solvated energy minimization. From the analysis, three important H-bond interactions between the Obg fold and the L13 were detected: Obg:Tyr27-L13:Glu32, Obg:Asn76-L13:Glu139, and Obg:Ala136-L13:Glu142. The interaction between the BsObg and BsL13 structures were also analyzed by electrostatic potential calculations to examine the interface surfaces. From the results, the key residues for hydrogen bonding and hydrophobic interaction between the two proteins were predicted. Conclusion and Prospects: In this study, we have focused on the binding mode of the BsObg protein with the ribosomal BsL13 protein. The interaction between the activated Obg and target protein was investigated with protein-protein docking calculations. The binding pattern can be further used as a base for structure-based drug design to find a novel antibacterial drug.