• 한국균학회소식:학술대회논문집
• /
• 2015.05a
• /
• pp.53-53
• /
• 2015

Agrobacterium tumefaciens-mediated transformation(ATMT) of Flammulina velutipes was used to produce a diverse number of transformants to discover the functions of gene that is vital for its variation color, spore pattern and cellulolytic activity. Futhermore, the transformant pool will be used as a good genetic resource for studying gene functions. Agrobacterium-mediated transformation was conducted in order to generate intentional mutants of F. velutipes strain KACC42777. Then Agrobacterium tumefaciens AGL-1 harboring pBGgHg was transformed into F. velutipes. This method is use to determine the functional gene of F. velutipes. Inverse PCR was used to insert T-DNA into the tagged chromosomal DNA segments and conducting sequence analysis of the F. velutipes. But this experiment had trouble in diverse morphological mutants because of dikaryotic nature of mushroom. It needed to make monokaryotic fruiting varients which introduced genes of compatible mating types. In this study, next generation sequencing data was generated from 28 strains of Flammulina velutipes with different phenotypes using Illumina Hiseq platform. Filtered short reads were initially aligned to the reference genome (KACC42780) to construct a SNP matrix. And then we built a phylogenetic tree based on the validated SNPs. The inferred tree represented that white- and brown- fruitbody forming strains were generally separated although three brown strains, 4103, 4028, and 4195, were grouped with white ones. This topological relationship was consistently reappeared even when we used randomly selected SNPs. Group I containing 4062, 4148, and 4195 strains and group II containing 4188, 4190, and 4194 strains formed early-divergent lineages with robust nodal supports, suggesting that they are independent groups from the members in main clades. To elucidate the distinction between white-fruitbody forming strains isolated from Korea and Japan, phylogenetic analysis was performed using their SNP data with group I members as outgroup. However, no significant genetic variation was noticed in this study. A total of 28 strains of Flammulina velutipes were analyzed to identify the genomic regions responsible for producing white-fruiting body. NGS data was yielded by using Illumina Hiseq platform. Short reads were filtered by quality score and read length were mapped on the reference genome (KACC42780). Between the white- and brown fruitbody forming strains. There is a high possibility that SNPs can be detected among the white strains as homozygous because white phenotype is recessive in F. velutipes. Thus, we constructed SNP matrix within 8 white strains. SNPs discovered between mono3 and mono19, the parental monokaryotic strains of 4210 strain (white), were excluded from the candidate. If the genotypes of SNPs detected between white and brown strains were identical with those in mono3 and mono19 strains, they were included in candidate as a priority. As a result, if more than 5 candidates SNPs were localized in single gene, we regarded as they are possibly related to the white color. In F. velutipes genome, chr01, chr04, chr07,chr11 regions were identified to be associated with white fruitbody forming. White and Brown Fruitbody strains can be used as an identification marker for F. veluipes. We can develop some molecular markers to identify colored strains and discriminate national white varieties against Japanese ones.

• Microbiology and Biotechnology Letters
• /
• v.46 no.3
• /
• pp.253-260
• /
• 2018

Beneficial microorganisms are widely used in the forestry, livestock, and, in particular, agricultural sectors to control soilborne diseases and promote plant growth. However, the industrial utilization of these microorganisms is very limited, mainly due to uncertainty concerning their ability to colonize and persist in soil. In this study, the survival of beneficial microorganisms in field soil microcosms was investigated for 13 days using quantitative PCR with B. subtilis group-specific primers. Bacterial community dynamics of the treated soils were analyzed using 16S ribosomal RNA (rRNA) gene amplicon sequencing on the Illumina MiSeq platform. The average 16S rRNA gene copy number per g dry soil of Bacillus spp. was $4.37{\times}10^6$ after treatment, which was 1,000 times higher than that of the control. The gene copy number was generally maintained for a week and was reduced thereafter, but remained 100 times higher than that of the control. Bacterial community analysis indicated that Acidobacteria ($26.3{\pm}0.9%$), Proteobacteria ($24.2{\pm}0.5%$), Chloroflexi ($11.1{\pm}0.4%$), and Actinobacteria ($9.7{\pm}2.5%$) were abundant phyla in both treated and non-treated soils. In the treated soils, the relative abundance of Actinobacteria was lower, whereas those of Bacteroidetes and Firmicutes were higher compared to the control. Differences in total relative abundances of operational taxonomic units belonging to several genera were observed between the treated and non-treated soils, suggesting that inoculation of soil with the Bacillus strains influenced the relative abundances of certain groups of bacteria and, therefore, the dynamics of resident bacterial communities. These changes in resident soil bacterial communities in response to inoculation of soil with beneficial Bacillus spp. provide important information for the use of beneficial microorganisms in soil for sustainable agriculture.

• Journal of Life Science
• /
• v.23 no.1
• /
• pp.116-124
• /
• 2013

The aim was to examine resistance exercise-related genes after 8 weeks of resistance training. Thirty-two male Sprague-Dawley rats were divided into four groups: 4 weeks sedentary (4 wks CON, n=8), 8 weeks sedentary (8 wks CON, n=8), 4 weeks exercise training (4 wks REG, n=8), and 8 weeks exercise training (8 wks REG, n=8). The rats were trained to climb a 1-m vertical incline (85-degree), with weights secured to their tails. They climbed 10 times, 3 days per week, for 8 consecutive weeks. Skeletal muscle was taken from the flexor halucis longus after the exercise training. After separating the total RNA, large-scale gene expression was investigated by beadarray (Illumina RatRef-12 Expression BeadChip) analysis, and qPCR was used to inspect the beadarray data and to analyze the RNA quantitatively. The detection p-value for the genes was p<0.01, the M-value {M=$log_2$(condition)-$log_2$(reference)} was >1.0, and the DiffScore was >20. In total, the expression of 30 genes significantly increased 4 weeks after the exercise training, and the expression of six genes decreased. At 8 weeks, the expression of five genes significantly increased and that of 12 decreased. Several genes are potentially involved in resistance exercise and muscle hypertrophy, including 1) regulation of cell growth (IGFBP1, PLA2G2A, OKL38); 2) myogenesis (CSRP3); 3) tissue regeneration and muscle development (MUSTN1, MYBPH); 4) hypertrophy (CYR61, ATF3, NR4A3); and 5) glucose metabolism (G6PC, PCK1). These results may help to explain previously reported physiological changes of the skeletal muscle and suggest new avenues for further investigation.

• Journal of Life Science
• /
• v.29 no.1
• /
• pp.1-8
• /
• 2019

In Hanwoo cattle (Korean native cattle), there is a scarcity of comparative analysis papers using highdepth sequencing and haplotype phasing, particularly a comparative analysis of the Truseq Synthetic Long-Read Haplotyping sequencing platform serviced by Illumina (TSLRH) versus the Chromium Genome Sequencing platform serviced by 10X Genomics (10XG). DNA was extracted from the sperm of a Hanwoo breeding bull (ID: TN1505D2184/27214) provided by Hanwoo research canter and used for the generation of sequence data from both the sequencing platforms. We then identified SNVs using an appropriate analysis pipeline tailored for each platform. The TSLRH and 10XG platforms generated a total of 355,208,304 and 1,632,772,004 reads, respectively, corresponding to a Q30 (%) of 89.04% and 88.60%, respectively, of which 351,992,768(99.09%) and 1,526,641,824(93.50%) were successfully mapped. For the TSLRH and 10XG platforms, the mean depth of the sequencing was 13.04X and 74.3X, the longest phase block was 1,982,706 bp and 1,480,081 bp, the N50 phase block was 57,637 bp and 114,394 bp, the total number of SNVs identified was 4,534,989 and 8,496,813, and the total phased rate was 72.29% and 87.67%, respectively. Moreover, for each chromosome, we identified unique and common SNVs using both sequencing platforms. The number of SNVs was directly proportional to the length of the chromosome. Based on our results, we recommend the use of the 10XG platform for haplotype phasing and SNV identification, as it generated a longer N50 phase block, in addition to a higher mean depth, total number of reads, total number of SNVs, and phase rate, than the TSLRH platform.

• IMMUNE NETWORK
• /
• v.11 no.5
• /
• pp.258-267
• /
• 2011

Background: Current management strategies attempt to diagnose rheumatoid arthritis (RA) at an early stage. Transcription profiling is applied in the search for biomarkers for detecting early-stage disease. Even though gene profiling has been reported using several animal models of RA, most studies were performed after the development of active arthritis, and conducted only on the peripheral blood and joint. Therefore, we investigated gene expression during the initial phase of collagen-induced arthritis (CIA) before the arthritic features developed in the thymus in addition to the peripheral blood and synovium. Methods: For gene expression analysis using cDNA microarray technology, samples of thymus, blood, and synovium were collected from CIA, rats immunized only with type II collagen (Cll), rats immunized only with adjuvant, and unimmunized rats on days 4 and 9 after the first immunization. Arrays were scanned with an Illumina bead array. Results: Of the 21,910 genes in the array, 1,243 genes were differentially expressed at least 2-fold change in various organs of CIA compared to controls. Among the 1,243 genes, 8 encode T-cell receptors (TCRs), including CD3${\zeta}$, CD3${\delta}$, CD3${\varepsilon}$, CD8${\alpha}$, and CD8${\beta}$ genes, which were down-regulated in CIA. The synovium was the organ in which the genes were differentially expressed between CIA and control group, and no difference were found in the thymus and blood. Further, we determined that the differential expression was affected by adjuvant more than Cll. The differential expression of genes as revealed by real-time RT-PCR, was in agreement with the microarray data. Conclusion: This study provides evidence that the genes encoding TCRs including CD3${\zeta}$, CD3${\delta}$, CD3${\varepsilon}$, CD8${\alpha}$, and CD8${\beta}$ genes were down-regulated during the initial phase of CIA in the synovium of CIA. In addition, adjuvant played a greater role in the down-regulation of the CD3 complex compared to CII. Therefore, the down-regulation of TCR gene expression occurred dominantly by adjuvant could be involved in the pathogenesis of the early stage at CIA.

• Asian-Australasian Journal of Animal Sciences
• /
• v.27 no.11
• /
• pp.1540-1547
• /
• 2014

The accuracy of genomic estimated breeding values (GEBV) was evaluated for sixteen meat quality traits in a Berkshire population (n = 1,191) that was collected from Dasan breeding farm, Namwon, Korea. The animals were genotyped with the Illumina porcine 62 K single nucleotide polymorphism (SNP) bead chips, in which a set of 36,605 SNPs were available after quality control tests. Two methods were applied to evaluate GEBV accuracies, i.e. genome based linear unbiased prediction method (GBLUP) and Bayes B, using ASREML 3.0 and Gensel 4.0 software, respectively. The traits composed different sets of training (both genotypes and phenotypes) and testing (genotypes only) data. Under the GBLUP model, the GEBV accuracies for the training data ranged from $0.42{\pm}0.08$ for collagen to $0.75{\pm}0.02$ for water holding capacity with an average of $0.65{\pm}0.04$ across all the traits. Under the Bayes B model, the GEBV accuracy ranged from $0.10{\pm}0.14$ for National Pork Producers Council (NPCC) marbling score to $0.76{\pm}0.04$ for drip loss, with an average of $0.49{\pm}0.10$. For the testing samples, the GEBV accuracy had an average of $0.46{\pm}0.10$ under the GBLUP model, ranging from $0.20{\pm}0.18$ for protein to $0.65{\pm}0.06$ for drip loss. Under the Bayes B model, the GEBV accuracy ranged from $0.04{\pm}0.09$ for NPCC marbling score to $0.72{\pm}0.05$ for drip loss with an average of $0.38{\pm}0.13$. The GEBV accuracy increased with the size of the training data and heritability. In general, the GEBV accuracies under the Bayes B model were lower than under the GBLUP model, especially when the training sample size was small. Our results suggest that a much greater training sample size is needed to get better GEBV accuracies for the testing samples.

• Asian-Australasian Journal of Animal Sciences
• /
• v.31 no.6
• /
• pp.784-790
• /
• 2018

Objective: The genetic effects of an individual on the phenotypes of its social partners, such as its pen mates, are known as social genetic effects. This study aims to identify the candidate genes for social (pen-mates') average daily gain (ADG) in pigs by using the genome-wide association approach. Methods: Social ADG (sADG) was the average ADG of unrelated pen-mates (strangers). We used the phenotype data (16,802 records) after correcting for batch (week), sex, pen, number of strangers (1 to 7 pigs) in the pen, full-sib rate (0% to 80%) within pen, and age at the end of the test. A total of 1,041 pigs from Landrace breeds were genotyped using the Illumina PorcineSNP60 v2 BeadChip panel, which comprised 61,565 single nucleotide polymorphism (SNP) markers. After quality control, 909 individuals and 39,837 markers remained for sADG in genome-wide association study. Results: We detected five new SNPs, all on chromosome 6, which have not been associated with social ADG or other growth traits to date. One SNP was inside the prostaglandin $F2{\alpha}$ receptor (PTGFR) gene, another SNP was located 22 kb upstream of gene interferon-induced protein 44 (IFI44), and the last three SNPs were between 161 kb and 191 kb upstream of the EGF latrophilin and seven transmembrane domain-containing protein 1 (ELTD1) gene. PTGFR, IFI44, and ELTD1 were never associated with social interaction and social genetic effects in any of the previous studies. Conclusion: The identification of several genomic regions, and candidate genes associated with social genetic effects reported here, could contribute to a better understanding of the genetic basis of interaction traits for ADG. In conclusion, we suggest that the PTGFR, IFI44, and ELTD1 may be used as a molecular marker for sADG, although their functional effect was not defined yet. Thus, it will be of interest to execute association studies in those genes.

• Korean Journal of Ichthyology
• /
• v.32 no.2
• /
• pp.39-48
• /
• 2020

The slender shinner (Pseudopungtungia tenuicorpa), a tiny freshwater fish of about 8 to 10 cm belonging to Cyprinidae, is an endangered species found only in the Han and Imjin Rivers on the Korean Peninsula. During the breeding season, this species spawns in nests of Coreoperca herzi, a predator of this species, or small crevices on rocks. This unique reproductive ecology can make this species more vulnerable to anthropogenic perturbance that can further limit the places to spawn. Here, mtDNA and microsatellite loci were analyzed to identify the genetic diversity and structure of slender shinners and further to provide the basic data necessary for the conservation planning of this species. A total of 28 polymorphic microsatellite markers were developed using Illumina paired-end sequencing, and 67 slender shinners collected from three localities in the Han River were genotyped using these loci. This species showed a remarkably high level of genetic diversity with mean expected heterozygosity of 0.914 and mean allele number per locus of 27.9, and no signature of drastic demographic decline was detected. As a result of our microsatellite analysis, the genetic structure between the two stems of the Han River, North Han and South Han, was prominent. Such a genetic structure was also evident in the sequence analysis of 14 haplotypes obtained from mtDNA control region. Although slender shinners are only found in very limited areas around the world, the genetic structure indicates that there is a block of gene flow among the populations, which should be reviewed in the future if management and restoration of this species is needed.

• Asian-Australasian Journal of Animal Sciences
• /
• v.30 no.9
• /
• pp.1215-1224
• /
• 2017

Objective: To assess the differences in ovarian transcriptomes in Shan Ma ducks between their peak and late stages of egg production, and to obtain new transcriptomic data of these egg-producing ducks. Methods: The Illumina HiSeq 2000 system was used for high throughput sequencing of ovarian transcriptomes from Shan Ma ducks at their peak or late stages of egg production. Results: Greater than 93% of the sequencing data had a base quality score (Q score) that was not less than 20 (Q20). From ducks at their peak stage of egg production, 42,782,676 reads were obtained, with 4,307,499,083 bp sequenced. From ducks at their late stage of egg production, 45,316,166 reads were obtained, with 4,562,063,363 bp sequenced. A comparison of the two datasets identified 2,002 differentially expressed genes, with 790 upregulated and 1,212 downregulated. Further analysis showed that 1,645 of the 2,002 differentially expressed genes were annotated in the non-redundant (NR) database, with 646 upregulated and 999 downregulated. Among the differentially expressed genes with annotations in the NR database, 696 genes were functionally annotated in the clusters of orthologous groups of proteins database, involving 25 functional categories. One thousand two hundred four of the differentially expressed genes with annotations in the NR database were functionally annotated in the gene ontology (GO) database, and could be divided into three domains and 56 categories. The three domains were cellular component, molecular function, and biological process. Among the genes identified in the GO database, 451 are involved in development and reproduction. Analysis of the differentially expressed genes with annotations in the NR database against the Kyoto encyclopedia of genes and genomes database revealed that 446 of the genes could be assigned to 175 metabolic pathways, of which the peroxisome proliferator-activated receptor signaling pathway, insulin signaling pathway, fructose and mannose metabolic pathways, gonadotropin releasing hormone signaling pathway and transforming growth factor beta signaling pathway were significantly enriched. Conclusion: The differences in ovarian transcriptomes in Shan Ma ducks between their peak and late stages of egg production were elucidated, which greatly enriched the ovarian transcriptomic information of egg-producing ducks.

• Asian-Australasian Journal of Animal Sciences
• /
• v.29 no.10
• /
• pp.1371-1382
• /
• 2016

MicroRNAs (miRNAs) are short non-coding RNAs that post-transcriptionally regulate expression of mRNAs in many biological pathways. Liver plays an important role in the feed efficiency of animals and high and low efficient cattle demonstrated different gene expression profiles by microarray. Here we report comprehensive miRNAs profiles by next-gen deep sequencing in Angus cattle divergently selected for residual feed intake (RFI) and identify miRNAs related to feed efficiency in beef cattle. Two microRNA libraries were constructed from pooled RNA extracted from livers of low and high RFI cattle, and sequenced by Illumina genome analyser. In total, 23,628,103 high quality short sequence reads were obtained and more than half of these reads were matched to the bovine genome (UMD 3.1). We identified 305 known bovine miRNAs. Bta-miR-143, bta-miR-30, bta-miR-122, bta-miR-378, and bta-let-7 were the top five most abundant miRNAs families expressed in liver, representing more than 63% of expressed miRNAs. We also identified 52 homologous miRNAs and 10 novel putative bovine-specific miRNAs, based on precursor sequence and the secondary structure and utilizing the miRBase (v. 21). We compared the miRNAs profile between high and low RFI animals and ranked the most differentially expressed bovine known miRNAs. Bovine miR-143 was the most abundant miRNA in the bovine liver and comprised 20% of total expressed mapped miRNAs. The most highly expressed miRNA in liver of mice and humans, miR-122, was the third most abundant in our cattle liver samples. We also identified 10 putative novel bovine-specific miRNA candidates. Differentially expressed miRNAs between high and low RFI cattle were identified with 18 miRNAs being up-regulated and 7 other miRNAs down-regulated in low RFI cattle. Our study has identified comprehensive miRNAs expressed in bovine liver. Some of the expressed miRNAs are novel in cattle. The differentially expressed miRNAs between high and low RFI give some insights into liver miRNAs regulating physiological pathways underlying variation in this measure of feed efficiency in bovines.