Transcriptome-wide study revealed m6A regulation of embryonic muscle development in Dingan goose (Anser cygnoides orientalis).

BACKGROUND: The number of myofiber is determined during the embryonic stage and does not increase during the postnatal period for birds, including goose. Thus, muscle production of adult goose is pre-determined during embryogenesis. Previous studies show N6-methyladenosine (m6A) is an important regulator for skeletal muscle development of birds and miRNAs play as a co-regulator for the skeletal muscle development in birds. Herein, we sequenced m6A and miRNA transcriptomes to investigate the profiles of m6A and their potential mechanism of regulating breast muscle development in Dingan Goose. RESULTS: We selected embryonic 21th day (E21) and embryonic 30th day (E30) to investigate the roles of transcriptome-wide m6A modification combining with mRNAs and miRNAs in goose breast muscle development. In this study, m6A peaks were mainly enriched in coding sequence (CDS) and start codon and397 genes were identified as differentially methylated genes (DMGs). GO and KEGG analysis showed that DMGs were highly related to cellular and metabolic process and that most DMGs were enriched in muscle-related pathways including Wnt signaling pathway, mTOR signaling and FoxO signaling pathway. Interestingly, a negative correlation between m6A methylation level and mRNA abundance was found through the analysis of m6A-RNA and RNA-seq data. Besides, we found 26 muscle-related genes in 397 DMGs. We also detected 228 differentially expressed miRNAs (DEMs), and further found 329 genes shared by the target genes of DEMs and DMGs (m6A-miRNA-genes), suggesting a tightly relationship between DEMs and DMGs. Among the m6A-miRNA-genes, we found 10 genes are related to breast muscle development. We further picked out an m6A-miRNA-gene, PDK3, from the 10 genes to visualize it and the result showed differentially methylated peaks on the mRNA transcript consistent with our m6A-seq results. CONCLUSION: GO and KEGG of DMGs between E21 and E30 showed most DMGs were muscle-related. In total, 228 DEMs were found, and the majority of DMGs were overlapped with the targets of DEGs. The differentially methylated peaks along with an m6A-miRNA-gene, PDK3, showed the similar results with m6A-seq results. Taken together, the results presented here provide a reference for further investigation of embryonic skeletal muscle development mechanism in goose.

m6A and miRNA jointly regulate the development of breast muscles in duck embryonic stages.

N6-methyladenosine (m6A) is an abundant internal mRNA modification and plays a crucial regulatory role in animal growth and development. In recent years, m6A modification has been found to play a key role in skeletal muscles. However, whether m6A modification contributes to embryonic breast muscle development of Pekin ducks has not been explored. To explore the role of m6A in embryonic breast muscle development of ducks, we performed m6A sequencing and miRNA sequencing for the breast muscle of duck embryos on the 19th (E19) and 27th (E27) days. A total of 12,717 m6A peaks were identified at E19, representing a total of 7,438 gene transcripts. A total of 14,703 m6A peaks were identified, which overlapped with the transcripts of 7,753 genes at E27. Comparing E19 and E27, we identified 2,347 differential m6A peaks, which overlapped with 1,605 m6A-modified genes (MMGs). Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses revealed that MMGs were enriched in multiple muscle- or fat-related pathways, which was also revealed from our analysis of differentially expressed genes (DEGs). Conjoint analysis of m6A-seq and RNA-seq data showed that pathways related to ß-oxidation of fatty acids and skeletal muscle development were significantly enriched, suggesting that m6A modification is involved in the regulation of fat deposition and skeletal muscle development. There were 90 upregulated and 102 downregulated miRNAs identified between the E19 and E27 stages. Through overlapping analysis of genes shared by MMGs and DEGs and the targets of differentially expressed miRNAs (DEMs), we identified six m6A-mRNA-regulated miRNAs. Finally, we found that m6A modification can regulate fat deposition and skeletal muscle development. In conclusion, our results suggest that m6A modification is a key regulator for embryonic breast muscle development and fat deposition of ducks by affecting expressions of mRNAs and miRNAs. This is the first study to comprehensively characterize the m6A patterns in the duck transcriptome. These data provide a solid basis for future work aimed at determining the potential functional roles of m6A modification in adipose deposition and muscle growth.

Genetic characteristics of Jiaji Duck by whole genome re-sequencing.

Jiaji Duck (JJ) is a Muscovy duck species that possesses many superior characteristics, and it has become an important genetic resource in China. However, to date, its genetic characteristics and genetic relationship with other duck breeds have not been explored yet, which greatly limits the utilization of JJ. In the present study, we investigated the genome sequences of 15 individual ducks representing five different duck populations, including JJ, French Muscovy duck (FF), mallard (YD), hong duck (HD) and Beijing duck (BD). Moreover, we investigated the characteristics of JJ-specific single nucleotide polymorphisms (SNPs) and compared the genome sequences of JJ vs. YD and JJ vs. BD using integrated strategies, including mutation detection, selective screening, and Gene Ontology (GO) analysis. More than 40 Gb of clean data were obtained for each population (mean coverage of 13.46 Gb per individual). A total number of 22,481,367 SNPs and 4,156,829 small insertion-deletions (Indels) were identified for the five duck populations, which could be used as molecular markers in breeding and utilization of JJ. Moreover, we identified 1,447,932 JJ-specific SNPs, and found that genes covering at least one JJ-specific SNP mainly involved in protein phosphorylation and dephosphorylation, as well as DNA modification. Phylogenetic tree and principal components analysis (PCA) revealed that the genetic relationship of JJ was closest to FF, while it was farthest to BD. A total of 120 and 111 genes were identified as positive selection genes for JJ vs. BD and JJ vs. YD, respectively. GO and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses showed that the positive selection genes for JJ vs. BD ducks mainly involved in pigmentation, muscle contraction and stretch, gland secretion, and immunology, while the positive selection genes obtained from JJ vs. YD ducks mainly involved in embryo development, muscle contraction and stretch, and gland secretion. Taken together, our findings enabled us to better understand the characteristics of JJ and provided a molecular basis for the breeding and hybrid utilization of JJ in the future.

Genotyping and zoonotic potential of Enterocytozoon bieneusi in cattle farmed in Hainan Province, the southernmost region of China.

Enterocytozoon bieneusi is an intestinal pathogen that infects a wide range of species, including humans. Cattle constitute an important host for E. bieneusi; however, there is a scarcity of information on the prevalence and genotyping of E. bieneusi in cattle in the Hainan Province of China. In this study, PCR analysis of 314 fecal samples from cattle in six cities of Hainan was performed for genotype identification. The average prevalence of E. bieneusi in these animals was 9.9% (31/314), and ranged from 0.0% (0/12) to 20.5% (8/39). Five known genotypes - EbpC (n = 14), BEB4 (n = 12), J (n = 2), I (n = 1), and CHG5 (n = 1) - and a novel genotype: HNC-I (n = 1) - were identified. Genotypes EbpC and HNC-I were placed in zoonotic Group 1, and the remaining four genotypes (BEB4, J, I, and CHG5) were placed in Group 2. Since 93.5% of the genotypes found in the cattle (29/31) (EbpC, BEB4, J, and I) have previously been found in humans, these genotypes are probably involved in the transmission of microsporidiosis to humans.

TITLE: Génotypage et potentiel zoonotique d'Enterocytozoon bieneusi chez les bovins élevés dans la province de Hainan, la région la plus au sud de la Chine. ABSTRACT: Enterocytozoon bieneusi est un pathogène intestinal qui infecte un large éventail d'espèces, y compris les humains. Le bétail constitue un hôte important pour E. bieneusi, mais les informations sur la prévalence et le génotypage d'E. bieneusi chez les bovins de la province de Hainan en Chine sont rares. Dans cette étude, une analyse PCR de 314 échantillons fécaux provenant de bovins dans six villes de Hainan a été réalisée pour l'identification du génotype. La prévalence moyenne d'E. bieneusi chez ces animaux était de 9,9 % (31/314), et variait de 0,0 % (0/12) à 20,5 % (8/39). Cinq génotypes connus, EbpC (n = 14), BEB4 (n = 12), J (n = 2), I (n = 1) et CHG5 (n = 1), et un nouveau génotype, HNC-I (n = 1), ont été identifiés. Les génotypes EbpC et HNC-I sont placés dans le groupe zoonotique 1, et les quatre génotypes restants (BEB4, J, I et CHG5) sont placés dans le groupe 2. Puisque 93,5 % (29/31) (EbpC, BEB4, J et I) des génotypes trouvés chez les bovins ont déjà été trouvés chez l'homme, ces génotypes sont probablement impliqués dans la transmission de la microsporidiose à l'homme.