The network regulation mechanism of the effects of heat stress on the production performance and egg quality of Jinding duck was analyzed by miRNA‒mRNA.

To explore the differential regulation mechanism of heat stress on the egg production performance and egg quality of Jinding ducks, 200 Jinding ducks (360-day-old) in good health and with similar body weights and a normal appetite were selected and randomly divided into a control (normal temperature [NT]) group (20°C-25°C) and a heat stress (HS) group (32°C-36°C), with 4 replicates in each group and 25 ducks in each replicate. The pretrial period was 1 wk, and the formal trial period was 4 wk. At the end of the 4th wk, 12 duck eggs were collected from each replicate to determine egg quality. Pituitary and ovarian tissues of Jinding ducks were collected, transcriptome sequencing was performed to screen differentially expressed miRNAs and mRNAs related to high temperature and heat stress, and a competitive endogenous RNA regulatory network was constructed. The sequencing data were verified by qRT‒PCR method. The following results were obtained: (1) Compared with the NT group, the HS group had a significantly lower laying rate, total egg weight, average egg weight, total feed intake, and feed intake per duck (P < 0.01), an extremely significantly higher feed-to-egg ratio (P < 0.01), and a higher mortality rate. (2) Compared with the NT group, the HS group had an extremely significantly lower egg weight, egg yolk weight, eggshell weight, and eggshell strength (P < 0.01) and an extremely significantly lower yolk ratio and eggshell thickness (P < 0.01, P < 0.05); however, there was no significant difference in the egg shape index, Haugh unit or protein height (P > 0.05). (3) A total of 1,974 and 1,202 genes were identified in the pituitary and ovary, respectively, and there were 5 significantly differentially expressed miRNAs. The differentially expressed genes were involved in the arginine and proline metabolism pathways, ether lipid metabolism pathway, and drug metabolism-cytochrome P450 pathway, which are speculated to be related to the egg production performance of Jingding ducks under high-temperature heat stress. (4) Novel_221 may target the PRPS1 gene to participate in egg production performance; novel_168 and novel_289 may target PIGW; novel_289 may target Q3MUY2; and novel_289 and novel_208 may target PIGN or genes that may be related to high-temperature heat stress. (5) In pituitary tissue, upregulated novel_141 (center of the network) formed a regulatory network with HSPB1 and HSP30A, and downregulated novel_366 (center of the network) formed a regulatory network with the JIP1 gene. In ovarian tissue, downregulated novel_289 (center of the network) formed a regulatory network with the ZSWM7, ABI3, and K1C23 genes, novel_221 formed a regulatory network with the IGF1, BCL7B, SMC6, APOA4, and FARP2 genes, and upregulated novel_40 formed a regulatory network with the HA1FF10 gene. In summary, heat stress affects the production performance and egg quality of Jinding ducks by regulating the secretion of endocrine-related hormones and the release of neurotransmitters as well as the expression of miRNAs and mRNAs in pituitary and ovarian tissues. The miRNA‒mRNA regulatory network provides a theoretical basis for the molecular mechanism that regulates the stress response in pituitary and ovarian tissues, egg quality, and production performance under heat stress.

Analysis of circRNA-miRNA-mRNA regulatory network of embryonic gonadal development in Mulard duck.

The aim of the study was to explore the regulatory mechanism of differences in embryonic gonadal development between intergeneric distance hybrid offspring Mulard ducks and parent ducks. The morphological differences gonadal tissues of Muscovy ducks, Pekin ducks and Mulard ducks at 12.5-day embryonic age were observed by sectioning and hematoxylin-eosin (HE) staining. Then followed by transcriptome sequencing to screen for gonadal development-related differentially expressed circRNAs and mRNAs to construct a competitive endogenous RNA (ceRNA) regulatory network. Finally, qRT-PCR and luciferase reporter system were used to verify the sequencing data and targeting relationship of ceRNA pairs. The results showed that the seminiferous tubule lumen of Mulard ducks was not obvious, while there were obvious seminiferous tubules and tubular structures in testis of Pekin ducks and Muscovy ducks, with number and shape indicating maturity. There were 18 upregulated circRNAs and 16 downregulated circRNAs in Mulard ducks and Pekin ducks, respectively, and 39 upregulated circRNAs and 1 downregulated circRNA in Mulard ducks and Muscovy ducks, respectively. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis found that genes involves in dorso-ventral axis formation, for example, neurogenic locus notch homolog protein 1 (NOTCH1), were significantly enriched (P < 0.05). The novel_circ_0002265-gga-miR-122-5p-PAFAH1B2 regulatory network was constructed. The qRT-PCR results showed that the sequencing results were reliable. The dual-luciferase reporter assay showed that gga-miR-122-5p exists binding site of circ_0002265 and PAFAH1B2, indicating circ_0002265-gga-miR-122-5p-PAFAH1B2 targeting relationship. In summary, the embryonic gonadal development of intergeneric hybrid Mulard ducks may be regulated by differentially expressed circRNAs and genes, such as novel_circ_0000519, novel_circ_0003537, NOTCH1, FGFR2, PAFAH1B1, and PAFAH1B2, among which circ_0002265-gga-miR-122-5p-PAFAH1B2 may participate in the targeted regulation of gonadal development in Mulard ducks. The findings of this study are helpful for analyzing the mechanism of embryonic gonadal development differences in avians.

A novel splice variant of goat CPT1a gene and their diverse mRNA expression profiles.

The Alternative splicing (AS) of Carnitine palmitoyltransferase 1a (CPT1a) and their expression profiles had never been illuminated in goats until now. Herein, a novel splice transcript in the CPT1a gene that is predicted to result in the skipping of exons 6-19 (CPT1a-sv1) has been isolated in addition to the full-length transcript in goats. The result of RT-PCR showed that CPT1a-sv1 is 606 bp in length and consists of 6 exons. A novel exon 6 was consisted of partial exon 5 and partial exon 19, compared to that in CPT1a. RT-qPCR analysis showed that the expression patterns of CPT1a and CPT1a-sv1 are spatially different. In both kid and adult goats, the CPT1a transcript is strongly expressed in the liver, spleen, lung, kidney, and brain tissues. However, CPT1a-sv1 has a strong tissue-specific expression pattern, with moderate RNA levels in the liver and brain of kids, while highly expressed in the liver and minimally expressed in the brain of adults. We observed two transcripts to be involved in brain development. These findings improve our understanding of the function of the CPT1a gene in goats and provide information on the molecular mechanism of AS events.

Ovarian transcriptomic analysis of Shan Ma ducks at peak and late stages of egg production.

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.

[Dynamic changes of ecological footprint and ecological capacity in Fujian Province].

The analysis on the dynamic changes of ecological footprint and ecological capacity in Fujian Province showed that in 1999-2003, the ecological footprint per capita in the Province increased from 1.428 hm2 to 1.658 hm2, while the ecological capacity per capita decreased from 0.683 hm2 to 0.607 hm2, with an increased ecological deficit year after year. The contradiction between the ecological footprint and ecological capacity pricked up gradually, and the ecological environment was at risk. There existed a severe imbalance in the supply and demand of ecological footprint per capita. The main body of the demands was grassland and fossil fuel, accouting for 55.74% - 63.43% of the total, while their supply only occupied 0.77% - 0.82% and next to nothing of the ecological capacity per capita, respectively. As a whole, the ecological footprint per ten thousand yuan GDP declined in the five years, indicating that the resources use efficiency in the Province was improved gradually. Based on the analysis of the present situation of the economic development and resources distribution in the Province, the strategies on reducing ecological deficit were put forward.