Analysis of 16S rRNA gene sequencing data for the taxonomic characterization of the vaginal and the fecal microbial communities in Hanwoo.

OBJECTIVE: The study of Hanwoo (Korean native cattle) has mainly been focused on meat quality and productivity. Recently the field of microbiome research has increased dramatically. However, the information on the microbiome in Hanwoo is still insufficient, especially relationship between vagina and feces. Therefore, the purpose of this study is to examine the microbial community characteristics by analyzing the 16S rRNA sequencing data of Hanwoo vagina and feces, as well as to confirm the difference and correlation between vaginal and fecal microorganisms. As a result, the goal is to investigate if fecal microbiome can be used to predict vaginal microbiome. METHODS: A total of 31 clinically healthy Hanwoo that delivered healthy calves more than once in Cheongju, South Korea were enrolled in this study. During the breeding season, we collected vaginal and fecal samples and sequenced the microbial 16S rRNA genes V3-V4 hypervariable regions from microbial DNA of samples. RESULTS: The results revealed that the phylum-level microorganisms with the largest relative distribution were Firmicutes, Actinobacteria, Bacteroidetes, and Proteobacteria in the vagina, and Firmicutes, Bacteroidetes, and Spirochaetes in the feces, respectively. In the analysis of alpha, beta diversity, and effect size measurements (LefSe), the results showed significant differences between the vaginal and fecal samples. We also identified the function of these differentially abundant microorganisms by functional annotation analyses. But there is no significant correlation between vaginal and fecal microbiome. CONCLUSION: There is a significant difference between vaginal and fecal microbiome, but no significant correlation. Therefore, it is difficult to interrelate vaginal microbiome as fecal microbiome in Hanwoo. In a further study, it will be necessary to identify the genetic relationship of the entire microorganism between vagina and feces through the whole metagenome sequencing analysis and meta-transcriptome analysis to figure out their relationship.

Integrative methylome and transcriptome analysis of porcine abdominal fat indicates changes in fat metabolism and immune responses during different development.

Fat is involved in synthesizing fatty acids (FAs), FA circulation, and lipid metabolism. Various genetic studies have been conducted on porcine fat but understanding the growth and specific adipose tissue is insufficient. The purpose of this study is to investigate the epigenetic difference in abdominal fat according to the growth of porcine. The samples were collected from the porcine abdominal fat of different developmental stages (10 and 26 weeks of age). Then, the samples were sequenced using MBD-seq and RNA-seq for profiling DNA methylation and RNA expression. In 26 weeks of age pigs, differentially methylated genes (DMGs) and differentially expressed genes (DEGs) were identified as 2,251 and 5,768, compared with 10 weeks of age pigs, respectively. Gene functional analysis was performed using GO and KEGG databases. In functional analysis results of DMGs and DEGs, immune responses such as chemokine signaling pathways, B cell receptor signaling pathways, and lipid metabolism terms such as PPAR signaling pathways and fatty acid degradation were identified. It is thought that there is an influence between DNA methylation and gene expression through changes in genes with similar functions. The effects of DNA methylation on gene expression were investigated using cis-regulation and trans-regulation analysis to integrate and interpret different molecular layers. In the cis-regulation analysis using 629 overlapping genes between DEGs and DMGs, immune response functions were identified, while in trans-regulation analysis through the TF-target gene network, the co-expression network of lipid metabolism-related functions was distinguished. Our research provides an understanding of the underlying mechanisms for epigenetic regulation in porcine abdominal fat with aging.

Fat is involved in the synthesis of new fatty acids (FAs), FA circulation, and lipid metabolism. Various genetic studies have been conducted on porcine fat but understanding the growth and specific adipose tissue is insufficient. The purpose of this study is to investigate the epigenetic difference in abdominal fat according to the growth of porcine. Modifications in DNA methylation and expression values were confirmed epigenetically with growth. Changed genes in each DNA and RNA showed identical trends in the function of immune response and lipid metabolism. The effects of DNA methylation on gene expression were investigated using cis-regulation (functional enrichment analysis of overlapping genes) and trans-regulation (transcription factor and target gene networking) analysis to integrate and interpret different molecular layers. Our research provides an understanding of the underlying mechanisms for epigenetic regulation in porcine abdominal fat with aging.

Insight into the potential candidate genes and signaling pathways involved in lymphoma disease in dogs using a comprehensive whole blood transcriptome analysis.

Lymphoma is one of the most prevalent hematological cancers, accounting for 15-20 % of new cancer diagnoses in dogs. Therefore, this study aims to explore the important genes and pathways involved in canine lymphoma progression and understand the underlying molecular mechanisms using RNA sequencing. In this study, RNAs acquired from seven pairs of lymphoma and non-lymphoma blood samples were sequenced from different breeds of dogs. Sequencing reads were preprocessed, aligned with the reference genome, assembled and expressions were estimated through bioinformatics approaches. At a false discovery rate (FDR) < 0.05 and fold change (FC) ≥ 1.5, a total of 625 differentially expressed genes (DEGs) were identified between lymphoma and non-lymphoma samples, including 347 up-regulated DEGs such as SLC38A11, SCN3A, ZIC5 etc. and 278 down-regulated DEGs such as LOC475937, CSMD1, KRT14 etc. GO enrichment analysis showed that these DEGs were highly enriched for molecular function of ATP binding and calcium ion binding, cellular process of focal adhesion, and biological process of immune response, and defense response to virus. Similarly, KEGG pathways analysis revealed 11 significantly enriched pathways such as ECM-receptor interaction, cell cycle, PI3K-Akt signaling pathway, ABC transporters etc. In the protein-protein interaction (PPI) network, CDK1 was found to be a top hub gene with highest degree of connectivity. Three modules selected from the PPI network showed that canine lymphoma was highly associated with cell cycle, ECM-receptor interaction, hypertrophic cardiomyopathy, dilated cardiomyopathy and RIG-I-like receptor signaling pathway. Overall, our findings highlighted new candidate therapeutic targets for further testing in canine lymphoma and facilitate the understanding of molecular mechanism of lymphoma's progression in dogs.

Comprehensive Transcriptomic Comparison between Porcine CD8- and CD8+ Gamma Delta T Cells Revealed Distinct Immune Phenotype.

We aimed to comprehensively understand the functional mechanisms of immunity, especially of the CD8+/- subsets of gamma delta (γδ) T cells, using an RNA-sequencing analysis. Herein, γδ T cells were obtained from bronchial lymph node tissues of 38-day-old (after weaning 10-day: D10) and 56-day-old (after weaning 28-day: D28) weaned pigs and sorted into CD8+ and CD8- groups. Differentially expressed genes (DEGs) were identified based on the CD8 groups at D10 and D28 time points. We confirmed 1699 DEGs between D10 CD8+ versus D10 CD8- groups and 1784 DEGs between D28 CD8+ versus D28 CD8- groups; 646 upregulated and 561 downregulated DEGs were common. The common upregulated DEGs were enriched in the cytokine-cytokine receptor interaction and T cell receptor (TCR) signaling pathway, and the common downregulated DEGs were enriched in the B cell receptor signaling pathway. Further, chemokine-related genes, interferon gamma, and CD40 ligand were involved in the cytokine-cytokine receptor interaction and TCR signaling pathway, which are associated with inter-regulation in immunity. We expect our results to form the basic data required for understanding the mechanisms of γδ T cells in pigs; however, further studies are required in order to reveal the dynamic changes in γδ T cells under pathogenic infections, such as those by viruses.

Identification of Monoallelically Expressed Genes Associated with Economic Traits in Hanwoo (Korean Native Cattle).

Hanwoo, an indigenous Korean cattle breed, has been genetically improved by selecting superior sires called Korean-proven bulls. However, cows still contribute half of the genetic stock of their offspring, and allelic-specific expressed genes have potential, as selective targets of cows, to enhance genetic gain. The aim of this study is to identify genes that have MAEs based on both the genome and transcriptome and to estimate their effects on breeding values (BVs) for economically important traits in Hanwoo. We generated resequencing data for the parents and RNA-sequencing data for the muscle, fat, and brain tissues of the offspring. A total of 3801 heterozygous single nucleotide polymorphisms (SNPs) in offspring were identified and they were located in 1569 genes. Only 14 genes showed MAE (seven expressing maternal alleles and seven expressing paternal alleles). Tissue-specific MAE was observed, and LANCL1 showed maternal allele expression across all tissues. MAE genes were enriched for the biological process of cell death and angiogenesis, which included ACKR3 and PDCL3 genes, whose SNPs were significantly associated with BVs of lean meat production-related traits, such as weight at 12 months of age, carcass weight, and loin eye area. In the current study, monoallelically expressed genes were identified in various adult tissues and these genes were associated with genetic capacity in Hanwoo.

Transcriptome analysis of pig macrophages expressing porcine reproductive and respiratory syndrome virus non-structural protein 1.

Porcine reproductive and respiratory syndrome virus (PRRSV) is a causative pathogen of PRRS, one of the most economically disastrous swine diseases. Non-structural protein 1 (NSP1) of PRRSV consists of NSP1α and NSP1ß which exhibit papain like cysteine protease activity. Recent evidence demonstrates that PRRSV NSP1 may be participated in modulating host immunity, but very few host proteins were discovered as targets for NSP1. In this study, we used RNA-seq to investigate the functional role of PRRSV NSP1 in porcine alveolar macrophages, 3D4/31 cells. Compared to empty vector (mock) transfectant, NSP1, NSP1α, and NSP1ß expressing 3D4/31 cells displayed a total of 60 genes, 63 genes, and 80 genes as differentially expressed genes (DEGs), respectively. Most of DEGs are involved in early inflammatory responses including interleukin (IL)-17 signaling pathway, chemokine signaling pathway, tumor necrosis factor (TNF)-α signaling pathway, and cell adhesion molecules. Interestingly, PRRSV NSP1 expression in 3D4/31 cells decreased mRNA transcripts of Fosb and Gdf15 known to be involved in host cell signaling or host cell protection during inflammation. Therefore, PRRSV NSP1 might block the signaling involved in host immune surveillance. Further study is required to define the mechanism on how PRRSV NSP1 protein represses mRNA transcripts of specific host genes.

Author Correction: Comprehensive genome and transcriptome analyses reveal genetic relationship, selection signature, and transcriptome landscape of small-sized Korean native Jeju horse.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Integrated time-serial transcriptome networks reveal common innate and tissue-specific adaptive immune responses to PRRSV infection.

Porcine reproductive and respiratory syndrome virus (PRRSV) infection is the most important viral disease causing severe economic losses in the swine industry. However, mechanisms underlying gene expression control in immunity-responsible tissues at different time points during PRRSV infection are poorly understood. We constructed an integrated gene co-expression network and identified tissue- and time-dependent biological mechanisms of PRRSV infection through bioinformatics analysis using three tissues (lungs, bronchial lymph nodes [BLNs], and tonsils) via RNA-Seq. Three groups with specific expression patterns (i.e., the 3-dpi, lung, and BLN groups) were discovered. The 3 dpi-specific group showed antiviral and innate-immune signalling similar to the case for influenza A infection. Moreover, we observed adaptive immune responses in the lung-specific group based on various cytokines, while the BLN-specific group showed down-regulated AMPK signalling related to viral replication. Our study may provide comprehensive insights into PRRSV infection, as well as useful information for vaccine development.

Population analysis of the Korean native duck using whole-genome sequencing data.

BACKGROUND: Advances in next-generation sequencing technologies have provided an opportunity to perform population-level comparative genomic analysis to discover unique genomic characteristics of domesticated animals. Duck is one of the most popular domesticated waterfowls, which is economically important as a source of meat, eggs, and feathers. The objective of this study is to perform population and functional analyses of Korean native duck, which has a distinct meat flavor and texture phenotype, using whole-genome sequencing data. To study the distinct genomic features of Korean native duck, we conducted population-level genomic analysis of 20 Korean native ducks together with 15 other duck breeds. RESULTS: A total of 15.56 million single nucleotide polymorphisms were detected in Korean native duck. Based on the unique existence of non-synonymous single nucleotide polymorphisms in Korean native duck, a total of 103 genes related to the unique genomic characteristics of Korean native duck were identified in comparison with 15 other duck breeds, and their functions were investigated. The nucleotide diversity and population structures among the used duck breeds were then compared, and their phylogenetic relationship was analyzed. Finally, highly differentiated genomic regions among Korean native duck and other duck breeds were identified, and functions of genes in those regions were examined. CONCLUSIONS: This is the first study to compare the population of Korean native duck with those of other duck breeds by using whole-genome sequencing data. Our findings can be used to expand our knowledge of genomic characteristics of Korean native duck, and broaden our understanding of duck breeds.

Comprehensive genome and transcriptome analyses reveal genetic relationship, selection signature, and transcriptome landscape of small-sized Korean native Jeju horse.

The Jeju horse, indigenous to the Jeju Island in Korea may have originated from Mongolian horses. Adaptations to the local harsh environment have conferred Jeju horse with unique traits such as small-sized body, stocky head, and shorter limbs. These characteristics have not been studied previously at the genomic level. Therefore, we sequenced and compared the genome of 41 horses belonging to 6 breeds. We identified numerous breed-specific non-synonymous SNPs and loss-of-function mutants. Demographic and admixture analyses showed that, though Jeju horse is genetically the closest to the Mongolian breeds, its genetic ancestry is independent of that of the Mongolian breeds. Genome wide selection signature analysis revealed that genes such as LCORL, MSTN, HMGA2, ZFAT, LASP1, PDK4, and ACTN2, were positively selected in the Jeju horse. RNAseq analysis showed that several of these genes were also differentially expressed in Jeju horse compared to Thoroughbred horse. Comparative muscle fiber analysis showed that, the type I muscle fibre content was substantially higher in Jeju horse compared to Thoroughbred horse. Our results provide insights about the selection of complex phenotypic traits in the small-sized Jeju horse and the novel SNPs identified will aid in designing high-density SNP chip for studying other native horse breeds.

Genome-Wide Analysis of Allele-Specific Expression Patterns in Seventeen Tissues of Korean Cattle (Hanwoo).

The functional hemizygosity could be caused by the MAE of a given gene and it can be one of the sources to affect the phenotypic variation in cattle. We aimed to identify MAE genes across the transcriptome in Korean cattle (Hanwoo). For three Hanwoo family trios, the transcriptome data of 17 tissues were generated in three offspring. Sixty-two MAE genes had a monoallelic expression in at least one tissue. Comparing genotypes among each family trio, the preferred alleles of 18 genes were identified (maternal expression, n = 9; paternal expression, n = 9). The MAE genes are involved in gene regulation, metabolic processes, and immune responses, and in particular, six genes encode transcription factors (FOXD2, FOXM1, HTATSF1, SCRT1, NKX6-2, and UBN1) with tissue-specific expression. In this study, we report genome-wide MAE genes in seventeen tissues of adult cattle. These results could help to elucidate epigenetic effects on phenotypic variation in Hanwoo.

Comparison of Bacterial Populations in the Ceca of Swine at Two Different Stages and their Functional Annotations.

The microbial composition in the cecum of pig influences host health, immunity, nutrient digestion, and feeding requirements significantly. Advancements in metagenome sequencing technologies such as 16S rRNAs have made it possible to explore cecum microbial population. In this study, we performed a comparative analysis of cecum microbiota of crossbred Korean native pigs at two different growth stages (stage L = 10 weeks, and stage LD = 26 weeks) using 16S rRNA sequencing technology. Our results revealed remarkable differences in microbial composition, α and ß diversity, and differential abundance between the two stages. Phylum composition analysis with respect to SILVA132 database showed Firmicutes to be present at 51.87% and 48.76% in stages L and LD, respectively. Similarly, Bacteroidetes were present at 37.28% and 45.98% in L and LD, respectively. The genera Prevotella, Anaerovibrio, Succinivibrio, Megasphaera were differentially enriched in stage L, whereas Clostridium, Terrisporobacter, Rikenellaceae were enriched in stage LD. Functional annotation of microbiome by level-three KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway analysis revealed that glycine, serine, threonine, valine, leucine, isoleucine arginine, proline, and tryptophan metabolism were differentially enriched in stage L, whereas alanine, aspartate, glutamate, cysteine, methionine, phenylalanine, tyrosine, and tryptophan biosynthesis metabolism were differentially enriched in stage LD. Through machine-learning approaches such as LEfSe (linear discriminant analysis effect size), random forest, and Pearson's correlation, we found pathways such as amino acid metabolism, transport systems, and genetic regulation of metabolism are commonly enriched in both stages. Our findings suggest that the bacterial compositions in cecum content of pigs are heavily involved in their nutrient digestion process. This study may help to meet the demand of human food and can play significant roles in medicinal application.

Transcriptome analysis to identify long non coding RNA (lncRNA) and characterize their functional role in back fat tissue of pig.

Long non coding RNAs (lncRNA) have been previously found to be involved in important cellular activities like epigenetics, implantation, cell growth etc. in pigs. However, comprehensive analysis of lncRNA in back fat tissues at different developmental stages in pigs is still lacking. In this study we conducted transcriptome analysis in the back fat tissue of a F1 crossbred Korean Native Pig (KNP) × Yorkshire Pig to identify lncRNA. We investigated their role in 16 pigs at two different growth stages; stage 1 (10 weeks, n = 8) and stage 2 (26 weeks, n = 8). After quality assessment of sequencing reads, we got a total of 1,641,165 assembled transcripts out of eight paired end read from each stage. Among them, 6808 lncRNA transcripts were identified by filtering on the basis of multiple parameters like read length ≥ 200 nucleotides, exon numbers ≥2, FPKM ≥0.5, coding potential score < 0 etc. PFAM and RFAM were used to filter out all possible protein coding genes and housekeeping RNAs respectively. A total of 103 lncRNAs and 1057 mRNAs were found to be differentially expressed (DE) between the two stages (|log2FC| > 2, q < 0.05). We also identified 306 genes located around 100 kb upstream and 234 genes downstream around these DE lncRNA transcripts. The expression of top eleven DE lncRNAs (COL4A6, LY7S, MYH2, OXCT1, SMPDL3A, TMEM182, TTC36, RFOOOO4, RFOOO15, RFOOO45, CADM2) had been validating by qRT-PCR. Pathway and GO terms analysis showed that, positive regulation of biosynthetic process, Wnt signaling pathway, cellular protein modification process, and positive regulation of nitrogen compound were differentially enriched. Our results suggested that, KEGG pathways such as protein digestion and absorption, Arrhythmogenic right ventricular cardiomyopathy (ARVC) to be significantly enriched in both DE lncRNAs as well as DE mRNAs and involved in back fat tissues development. It also suggests that, identified lncRNAs are involved in regulation of important adipose tissues development pathways.

Return to Play after Modified Broström Operation for Chronic Ankle Instability in Elite Athletes.

BACKGROUND: This study assessed the average time to return to training and official game participation after modified Broström operation (MBO) in elite athletes. METHODS: Sixty athletes diagnosed with lateral ankle instability underwent MBO from October 2011 to December 2013. Their average age was 19.3 years, and the average follow-up time was 28.8 months. We measured the time sequence of three phases of rehabilitation: start of personal training, start of team training, and start of the first official game after recovery. Patients were divided into an early return to play (RTP) group and late RTP group. The groups were compared to identify possible risk factors affecting the RTP time. RESULTS: The mean length of time to return to personal training was 1.9 months, return to team training was 2.9 months, and return to competitive play was 3.9 months. There were no significant differences of any variables including age, sex, body mass index, level of sports, grade of instability, presence of os subfibulare, and preoperative functional score between the early RTP and late RTP groups. CONCLUSIONS: The RTP was 83.3% at 4 months after lateral ankle ligament repair and 100% at 8 months postoperatively. The results provide reference data for orthopedic surgeons in evaluating surgical results and informing patients about expectations after surgery in terms of performance level and timing of return to sports.

Identification and characterization of a novel KG42 xylanase (GH10 family) isolated from the black goat rumen-derived metagenomic library.

This study was conducted to isolate and functionally characterize a novel xylan-degrading enzyme from the microbial metagenomes of black goat rumens. A novel gene, KG42, was isolated from one of the 17 xylan-degrading metagenomic fosmid clones obtained from black goat rumens. The KG42 gene, comprising a 1107 bp open reading frame, encodes a protein composed of 368 amino acids (41 kDa) with a glycosyl hydrolase family 10 (GH10) domain, consisting of a "salad-bowl" shaped tertiary structure (a typical 8-fold α/ß barrel (α/ß)8) and two catalytic residues. KG42 xylanase protein has at best 40% sequence identity with other homologous GH10 xylanase proteins. The enzyme displayed its optimum activity at pH 5.0 and 50 °C. The enzyme was thermally stable at pH and temperature ranges of 5.0-10.0 and 20-60 °C, respectively. Substrate specificity and hydrolytic patterns implied that the KG42 xylanase functions as an endo-ß-1,4-xylanase (EC 3.2.1.8). The KG42 xylanase was also used for the preparation of bifidogenic xylan hydrolysates, demonstrating its potential applications toward preparing prebiotic xylooligosaccharides.

Non-Coding Transcriptome Maps across Twenty Tissues of the Korean Black Chicken, Yeonsan Ogye.

Yeonsan Ogye is a rare Korean domestic chicken breed whose entire body, including feathers and skin, has a unique black coloring. Although some protein-coding genes related to this unique feature have been examined, non-coding elements have not been widely investigated. Thus, we evaluated coding and non-coding transcriptome expression and identified long non-coding RNAs functionally linked to protein-coding genes in Ogye. High-throughput RNA sequencing and DNA methylation sequencing were performed to profile the expression of 14,264 Ogye protein-coding and 6900 long non-coding RNA (lncRNA) genes and detect DNA methylation in 20 different tissues of an individual Ogye. Approximately 75% of Ogye lncRNAs and 45% of protein-coding genes showed tissue-specific expression. For some genes, tissue-specific expression levels were inversely correlated with DNA methylation levels in their promoters. Approximately 39% of tissue-specific lncRNAs displayed functional associations with proximal or distal protein-coding genes. Heat shock transcription factor 2-associated lncRNAs appeared to be functionally linked to protein-coding genes specifically expressed in black skin tissues, more syntenically conserved in mammals, and differentially expressed in black relative to in white tissues. Pending experimental validation, our findings increase the understanding of how the non-coding genome regulates unique phenotypes and can be used for future genomic breeding of chickens.

Metagenomic Mining and Functional Characterization of a Novel KG51 Bifunctional Cellulase/Hemicellulase from Black Goat Rumen.

A novel KG51 gene was isolated from a metagenomic library of Korean black goat rumen and its recombinant protein was characterized as a bifunctional enzyme (cellulase/hemicellulase). In silico sequence and domain analyses revealed that the KG51 gene encodes a novel carbohydrate-active enzyme that possesses a salad-bowl-like shaped glycosyl hydrolase family 5 (GH5) catalytic domain but, at best, 41% sequence identity with other homologous GH5 proteins. Enzymatic profiles (optimum pH values and temperatures, as well as pH and thermal stabilities) of the recombinant KG51 bifunctional enzyme were also determined. On the basis of the substrate specificity data, the KG51 enzyme exhibited relatively strong cellulase (endo-ß-1,4-glucanase [EC 3.2.1.4]) and hemicellulase (mannan endo-ß-1,4-mannosidase [EC 3.2.1.78] and endo-ß-1,4-xylanase [EC 3.2.1.8]) activities, but no exo-ß-1,4-glucanase (EC 3.2.1.74), exo-ß-1,4-glucan cellobiohydrolase (EC 3.2.1.91), and exo-1,4-ß-xylosidase (EC 3.2.1.37) activities. Finally, the potential industrial applicability of the KG51 enzyme was tested in the preparation of prebiotic konjac glucomannan hydrolysates.

Whole genome and transcriptome maps of the entirely black native Korean chicken breed Yeonsan Ogye.

Background: Yeonsan Ogye (YO), an indigenous Korean chicken breed (Gallus gallus domesticus), has entirely black external features and internal organs. In this study, the draft genome of YO was assembled using a hybrid de novo assembly method that takes advantage of high-depth Illumina short reads (376.6X) and low-depth Pacific Biosciences (PacBio) long reads (9.7X). Findings: The contig and scaffold NG50s of the hybrid de novo assembly were 362.3 Kbp and 16.8 Mbp, respectively. The completeness (97.6%) of the draft genome (Ogye_1.1) was evaluated with single-copy orthologous genes using Benchmarking Universal Single-Copy Orthologs and found to be comparable to the current chicken reference genome (galGal5; 97.4%; contigs were assembled with high-depth PacBio long reads (50X) and scaffolded with short reads) and superior to other avian genomes (92%-93%; assembled with short read-only or hybrid methods). Compared to galGal4 and galGal5, the draft genome included 551 structural variations including the fibromelanosis (FM) locus duplication, related to hyperpigmentation. To comprehensively reconstruct transcriptome maps, RNA sequencing and reduced representation bisulfite sequencing data were analyzed from 20 tissues, including 4 black tissues (skin, shank, comb, and fascia). The maps included 15,766 protein-coding and 6,900 long noncoding RNA genes, many of which were tissue-specifically expressed and displayed tissue-specific DNA methylation patterns in the promoter regions. Conclusions: We expect that the resulting genome sequence and transcriptome maps will be valuable resources for studying domestic chicken breeds, including black-skinned chickens, as well as for understanding genomic differences between breeds and the evolution of hyperpigmented chickens and functional elements related to hyperpigmentation.

Integrated transcriptomes throughout swine oestrous cycle reveal dynamic changes in reproductive tissues interacting networks.

Female fertility is a highly regulated process involving the synchronized activities of multiple tissues. The underlying genomic regulation of the tissue synchronization is poorly understood. To understand this better we investigated the transcriptomes of the porcine ovary, endometrium, and oviduct at days 0, 3, 6, 9, 12, 15, or 18 of the oestrous cycle. We analysed the transcriptome profiles of the individual tissues and focus on the bridging genes shared by two or more tissues. The three tissue-networks were connected forming a triangular shape. We identified 65 bridging genes with a high level of connectivity to all other genes in the network. The expression levels showed negative correlations between the ovary and the other two tissues, and low correlations between endometrium and oviduct. The main functional annotations involved biosynthesis of steroid hormones, cell-to-cell adhesion, and cell apoptosis, suggesting that regulation of steroid hormone synthesis and tissue viability are major regulatory mechanisms.