Integrated mRNA-miRNA transcriptome profiling of blood immune responses potentially related to pulmonary fibrosis in forest musk deer.

Background: Forest musk deer (FMD, Moschus Berezovskii) is a critically endangered species world-widely, the death of which can be caused by pulmonary disease in the farm. Pulmonary fibrosis (PF) was a huge threat to the health and survival of captive FMD. MicroRNAs (miRNAs) and messenger RNAs (mRNAs) have been involved in the regulation of immune genes and disease development. However, the regulatory profiles of mRNAs and miRNAs involved in immune regulation of FMD are unclear. Methods: In this study, mRNA-seq and miRNA-seq in blood were performed to constructed coexpression regulatory networks between PF and healthy groups of FMD. The hub immune- and apoptosis-related genes in the PF blood of FMD were explored through Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis. Further, protein-protein interaction (PPI) network of immune-associated and apoptosis-associated key signaling pathways were constructed based on mRNA-miRNA in the PF blood of the FMD. Immune hub DEGs and immune hub DEmiRNAs were selected for experimental verification using RT-qPCR. Results: A total of 2744 differentially expressed genes (DEGs) and 356 differentially expressed miRNAs (DEmiRNAs) were identified in the PF blood group compared to the healthy blood group. Among them, 42 DEmiRNAs were negatively correlated with 20 immune DEGs from a total of 57 correlations. The DEGs were significantly associated with pathways related to CD molecules, immune disease, immune system, cytokine receptors, T cell receptor signaling pathway, Th1 and Th2 cell differentiation, cytokine-cytokine receptor interaction, intestinal immune network for IgA production, and NOD-like receptor signaling pathway. There were 240 immune-related DEGs, in which 186 immune-related DEGs were up-regulated and 54 immune-related DEGs were down-regulated. In the protein-protein interaction (PPI) analysis of immune-related signaling pathway, TYK2, TLR2, TLR4, IL18, CSF1, CXCL13, LCK, ITGB2, PIK3CB, HCK, CD40, CD86, CCL3, CCR7, IL2RA, TLR3, and IL4R were identified as the hub immune genes. The mRNA-miRNA coregulation analysis showed that let-7d, miR-324-3p, miR-760, miR-185, miR-149, miR-149-5p, and miR-1842-5p are key miRNAs that target DEGs involved in immune disease, immune system and immunoregulation. Conclusion: The development and occurrence of PF were significantly influenced by the immune-related and apoptosis-related genes present in PF blood. mRNAs and miRNAs associated with the development and occurrence of PF in the FMD.

Identification of Potential miRNA-mRNA Regulatory Network Associated with Growth and Development of Hair Follicles in Forest Musk Deer.

In this study, sRNA libraries and mRNA libraries of HFs of FMD were constructed and sequenced using an Illumina HiSeq 2500, and the expression profiles of miRNAs and genes in the HFs of FMD were obtained at the anagen and catagen stages. In total, 565 differentially expressed unigenes (DEGs) were identified, 90 of which were upregulated and 475 of which were downregulated. In the BP category of GO enrichment, the DEGs were enriched in the processes related to HF development and differentiation, including the hair cycle regulation and processes, HF development, skin epidermis development, regulation of HF development, skin development, the Wnt signaling pathway, and the BMP signaling pathway. Through KEGG analysis it was found that DEGs were significantly enriched in pathways associated with HF development and growth. A total of 186 differentially expressed miRNAs (DEmiRNAs) were screened (p < 0.05) in the HFs of FMD at the anagen stage vs. the catagen stage, 33 of which were upregulated and 153 of which were downregulated. Through DEmiRNA-mRNA association analysis, we found DEmiRNAs and target genes that mainly play regulatory roles in HF development and growth. The enrichment analysis of DEmiRNA target genes revealed similarities with the enrichment results of DEGs associated with HF development. Notably, both sets of genes were enriched in key pathways such as the Notch signaling pathway, melanogenesis, the cAMP signaling pathway, and cGMP-PKG. To validate our findings, we selected 11 DEGs and 11 DEmiRNAs for experimental verification using RT-qPCR. The results of the experimental validation were consistent with the RNA-Seq results.

[Structure Characteristics and Driving Variables of Epilithic Algae Community in Lhasa River Basin of Qinghai-Tibet Plateau].

In order to explore the characteristics and driving factors of the epilithic algae community in the middle and lower reaches of the Lhasa River, epilithic algae was collected and identified in September 2019, the species composition and spatial distribution were analyzed, and the key environmental factors affecting the epilithic algae community were identified through redundancy analysis. The results indicated that 31 genera of epilithic algae belonging to six phyla were identified, and the average cell density of the epilithic algae was 1.92×106 cells·m-2. The number of species and cell density of the diatom phylum were the largest at each sampling point. The species number and cell density of the epilithic algae were significantly different between the main stream and tributaries. The main stream and tributaries contained relatively small differences in the dominant algae genera, which were Gomphonema, Fragilaria, Cymbella, and Planktolyngbya in the main stream and Gomphonema, Fragilaria, Cymbella, and Oscillatoria in the tributaries. The redundancy analysis revealed that temperature, pH, and dissolved oxygen were the main driving factors affecting the community structure of the dominant genus in the main stream, while HCO3- and flow velocity controlled the community structure of the dominant genus in the tributaries. This study provides basic data and the theoretical basis for the conservation of aquatic ecosystems and water quality management in the middle and lower reaches of the Lhasa River.

Distribution patterns of microsatellites and development of its marker in different genomic regions of forest musk deer genome based on high throughput sequencing.

Forest musk deer (Moschus berezovskii, FMD) is an endangered artiodactyl species, male FMD produce musk. We have sequenced the whole genome of FMD, completed the genomic assembly and annotation, and performed bioinformatic analyses. Our results showed that microsatellites (SSRs) displayed nonrandomly distribution in genomic regions, and SSR abundances were much higher in the intronic and intergenic regions compared to other genomic regions. Tri- and hexanucleotide perfect (P) SSRs predominated in coding regions (CDSs), whereas, tetra- and pentanucleotide P-SSRs were less abundant. Trifold P-SSRs had more GC-contents in the 5'-untranslated regions (5'UTRs) and CDSs than other genomic regions, whereas mononucleotide P-SSRs had the least GC-contents. The repeat copy numbers (RCN) of the same mono- to hexanucleotide P-SSRs had different distributions in different genomic regions. The RCN of trinucleotide P-SSRs had increased significantly in the CDSs compared to the transposable elements (TEs), intronic and intergenic regions. The analysis of coefficient of variability (CV) of P-SSRs showed that the RCN of mononucleotide P-SSRs had relative higher variation in different genomic regions, followed by the CV pattern of RCN: dinucleotide P-SSRs > trinucleotide P-SSRs > tetranucleotide P-SSRs > pentanucleotide P-SSRs > hexanucleotide P-SSRs. The CV variations of RCN of the same mono- to hexanucleotide P-SSRs were relative higher in the intron and intergenic regions, followed by that in the TEs, and the relative lower was in the 5'UTR, CDSs and 3'UTRs. 58 novel polymorphic SSR loci were detected based on genotyping DNA from 36 captive FMD and 22 SSR markers finally showed polymorphism, stability, and repetition.

Distribution patterns and variation analysis of simple sequence repeats in different genomic regions of bovid genomes.

As the first examination of distribution, guanine-cytosine (GC) pattern, and variation analysis of microsatellites (SSRs) in different genomic regions of six bovid species, SSRs displayed nonrandomly distribution in different regions. SSR abundances are much higher in the introns, transposable elements (TEs), and intergenic regions compared to the 3'-untranslated regions (3'UTRs), 5'UTRs and coding regions. Trinucleotide perfect SSRs (P-SSRs) were the most frequent in the coding regions, whereas, mononucleotide P-SSRs were the most in the introns, 3'UTRs, TEs, and intergenic regions. Trifold P-SSRs had more GC-contents in the 5'UTRs and coding regions than that in the introns, 3'UTRs, TEs, and intergenic regions, whereas mononucleotide P-SSRs had the least GC-contents in all genomic regions. The repeat copy numbers (RCN) of the same mono- to hexanucleotide P-SSRs showed significantly different distributions in different regions (P < 0.01). Except for the coding regions, mononucleotide P-SSRs had the most RCNs, followed by the pattern: di- > tri- > tetra- > penta- > hexanucleotide P-SSRs in the same regions. The analysis of coefficient of variability (CV) of SSRs showed that the CV variations of RCN of the same mono- to hexanucleotide SSRs were relative higher in the intronic and intergenic regions, followed by the CV variation of RCN in the TEs, and the relative lower was in the 5'UTRs, 3'UTRs, and coding regions. Wide SSR analysis of different genomic regions has helped to reveal biological significances of their distributions.

Distinct patterns of simple sequence repeats and GC distribution in intragenic and intergenic regions of primate genomes.

As the first systematic examination of simple sequence repeats (SSRs) and guanine-cytosine (GC) distribution in intragenic and intergenic regions of ten primates, our study showed that SSRs and GC displayed nonrandom distribution for both intragenic and intergenic regions, suggesting that they have potential roles in transcriptional or translational regulation. Our results suggest that the majority of SSRs are distributed in non-coding regions, such as the introns, TEs, and intergenic regions. In these primates, trinucleotide perfect (P) SSRs were the most abundant repeats type in the 5'UTRs and CDSs, whereas, mononucleotide P-SSRs were the most in the intron, 3'UTRs, TEs, and intergenic regions. The GC-contents varied greatly among different intragenic and intergenic regions: 5'UTRs > CDSs > 3'UTRs > TEs > introns > intergenic regions, and high GC-content was frequently distributed in exon-rich regions. Our results also showed that in the same intragenic and intergenic regions, the distribution of GC-contents were great similarity in the different primates. Tri- and hexanucleotide P-SSRs had the most GC-contents in the 5'UTRs and CDSs, whereas mononucleotide P-SSRs had the least GC-contents in the six genomic regions of these primates. The most frequent motifs for different length varied obviously with the different genomic regions.

[Correlation between distribution of rhizospheric microorganisms and contents of steroidal saponins of Paris polyphylla var. yunnanensis].

In this paper, the varying pattern of the amount of rhizospheric microorganisms, including bacteria, actinomycetes and fungus, was observed during the cultivation of Paris polyphylla var. yunnanensis. And the correlations between number of rhizospheric microorganisms and the quality of P. polyphylla var. yunnanensis were also studied. The results showed that the rhizospheric microorganism source of P. polyphylla var. yunnanensis was rich. The distribution of rhizospheric microorganisms (soil bacteria, fungus, actinomycetes, potassium-solubilizing bacteria, inorganic phosphorus-solubilizing bacteria, organic phosphorus-solubilizing bacteria) collected from different origin places existed significant difference (P < 0.05). The varying pattern for the amount of rhizospheric microorganisms was showed as following: the amount of bacteria > the amount of actinomycetes > the amount of fungus. The medicinal quality of P. polyphylla var. yunnanensis was influenced by their habits, and the increase of cultivation years caused the obvious decrease of the quality of P. polyphylla var. yunnanensis. Therefore, the increase of cultivation years will cause the variation of the soil micro-ecology flora, and decrease the nutrient absorption and the utilization of P. polyphylla var. yunnanensis, which will make the decrease of the medical quality of P. polyphylla var. yunnanensis.

Genome-Wide Survey and Analysis of Microsatellite Sequences in Bovid Species.

Microsatellites or simple sequence repeats (SSRs) have become the most popular source of genetic markers, which are ubiquitously distributed in many eukaryotic and prokaryotic genomes. This is the first study examining and comparing SSRs in completely sequenced genomes of the Bovidae. We analyzed and compared the number of SSRs, relative abundance, relative density, guanine-cytosine (GC) content and proportion of SSRs in six taxonomically different bovid species: Bos taurus, Bubalus bubalis, Bos mutus, Ovis aries, Capra hircus, and Pantholops hodgsonii. Our analysis revealed that, based on our search criteria, the total number of perfect SSRs found ranged from 663,079 to 806,907 and covered from 0.44% to 0.48% of the bovid genomes. Relative abundance and density of SSRs in these Bovinae genomes were non-significantly correlated with genome size (Pearson, r < 0.420, p > 0.05). Perfect mononucleotide SSRs were the most abundant, followed by the pattern: perfect di- > tri- > penta- > tetra- > hexanucleotide SSRs. Generally, the number of SSRs, relative abundance, and relative density of SSRs decreased as the motif repeat length increased in each species of Bovidae. The most GC-content was in trinucleotide SSRs and the least was in the mononucleotide SSRs in the six bovid genomes. The GC-contents of tri- and pentanucleotide SSRs showed a great deal of similarity among different chromosomes of B. taurus, O. aries, and C. hircus. SSR number of all chromosomes in the B. taurus, O.aries, and C. hircus is closely positively correlated with chromosome sequence size (Pearson, r > 0.980, p < 0.01) and significantly negatively correlated with GC-content (Pearson, r < -0.638, p < 0.01). Relative abundance and density of SSRs in all chromosomes of the three species were significantly negatively correlated with GC-content (Pearson, r < -0.333, P < 0.05) but not significantly correlated with chromosome sequence size (Pearson, r < -0.185, P > 0.05). Relative abundances of the same nucleotide SSR type showed great similarity among different chromosomes of B. taurus, O. aries, and C. hircus.

[Effects of mycorrhizal colonization and medicine quality of Paris polyphylla var. yunnanensis inoculated by different foreign AM fungi species].

After 28 foreign species of AM fungi were inoculated in sterilized soil, the effects of the AM mycorrhizal colonization and the medicine quality of Paris polyphylla var. yunnanensis were observed by combination of inoculation test in pot at room temperature and instrumental analysis. The results showed that, compared with control group (CK), the inoculation of foreign AM fungi in the soil influenced the spore density, mycorrhizal infection rate, and colonization intensity of AM fungi in root system of P. polyphylla var. yunnanensis. The inoculation of foreign AM fungi enhanced the mycorrhiza viability of P. polyphylla var. yunnanensis by increasing the activity of succinic dehydrogenase (SDH) and alkaline phosphatase (ALP) in intraradical hyphae. The content of single steroid saponin in rhizome of P. polyphylla var. yunnanensis showed variation after P. polyphylla var. yunnanensis was inoculated by different foreign species of AM fungi, which was beneficial for increasing the medicine quality; however, the kinds of steroid saponin showed no difference. In a degree, there was a selectivity of symbiosis between P. polyphylla var. yunnanensis and foreign AM fungi. And we found that the Claroideoglomus claroideum and Racocetra coralloidea were best foreign AM fungi species for cultivating P. polyphylla var. yunnanensis under field condition.

The reproductive performance of female Forest musk deer (Moschus berezovskii) in captivity.

Reproductive performance of 750 adult female Forest musk deer (FMD, Moschus berezovskii) was monitored (from 2005 to 2009) on two Forest musk deer farms in Sichuan province, China. The mean (± SEM) lengths of the estrous cycle, pregnancy, and consecutive fawning intervals were 17.5 ± 0.8 d (n = 64), 182.9 ± 0.8 d (n =190), and 363.2 ± 2.4 d (n = 120), respectively. Copulation occurred from the end of October to early March, with the majority (84.1 ± 1.9%) during November and December (P < 0.05). The fawning rate was 68.7 ± 2.2% in yearling hinds and 73.3 ± 1.3% in adult hinds. Fawning occurred from the end of April to early September, with the majority (88. 9 ± 2.2%) during May and June (P < 0.005), including a peak in May (52.0 ± 1.0%). Yearling hinds fawned later (7.3 ± 1.8 d) than adult hinds. Rates of premature delivery and dystocia were 2.7 ± 0.4 and 2.0 ± 0.3%, respectively. There was a 1.07:1 female-to-male ratio at birth (P > 0.05), with 67.7% twins, but only two cases of triplets and one set of quadruplets among 537 fawns. Birth weight averaged 518.7 ± 10.3 g, ranging from 304 to 775 g (n = 240), with male fawns slightly heavier than females (521.1 ± 14.8 vs 516.9 ± 14.3 g). Postnatal and pre-weaning mortality rates were 25.0 ± 1.2 and 17.9 ± 3.1%, respectively. Postnatal mortality and urinary calculus were much more common in male fawns (P < 0.01). The mortality rate due to suppuration and dyspepsia were different (P < 0.001) between postnatal and pre-weaning periods. Due to higher mortality of male fawns, the adult sex ratio was significantly female-biased. In conclusion, captive FMD had similar seasonal patterns of copulation and fawning as wild FMD and other musk deer species.

[Behavior ethogram and PAE coding system of Cervus nippon sichuanicus].

A monthly 5-day periodic observation at 06:00-18:00 from March to November 2007 was conducted to record the behavioral processes, contents, and results, and the surrounding habitats of Sichuan sika deer (Cervus nippon sichuanicus) in Donglie, Chonger, and Reer villages of Tiebu Natural Reserve of Sichuan Province. The behavioral ethogram, vigilance behaviors ethogram and its PAE (posture, act, and environment) coding system of the Sichuan sika deer were established, which filled the gap of the PAE coding of ungulates vigilance behaviors. A total of 11 kinds of postures, 83 acts, and 136 behaviors were recorded and distinguished, with the relative frequency of each behavior in relation to gender, age, and season described. Compared with other ungulates, the behavioral repertoire of Sichuan sika deer was mostly similar to that of other cervid animals.