The dynamic process involving the selection and maturation of follicles is regulated and controlled by a highly synchronized and exquisitely timed cascade of gene expression. Studies have shown that long non-coding RNA (lncRNA) is essential for the normal maintenance of animal reproductive function and has an important regulatory function in ovarian development and hormone secretion. In this study, a total of 2076 lncRNAs (1362 known lncRNAs and 714 new lncRNAs) and 25,491 mRNAs were identified in libraries constructed from Duroc ovaries on days 0, 2 and 4 of follicle development. lncRNAs were shorter, had fewer exons, exhibited a shorter ORF (Open Reading Frame) length and lower expression levels, and were less conserved than mRNAs. Furthermore, 1694 transcripts (140 lncRNAs and 1554 mRNAs) were found to be differentially expressed in pairwise comparisons. A total of 6945 co-localized mRNAs were detected in cis in 2076 lncRNAs. The most enriched GO (Gene Ontology) terms were related to developmental processes. KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway analysis revealed that the differentially expressed lncRNAs targeted mRNAs, and the differentially expressed mRNAs were related to the TGF-β signaling pathway, the PI3K-Akt signaling pathway, the Retinol metabolic pathway and the Wnt signaling pathway. This study deepened our understanding of the genetic basis and molecular mechanisms of follicular development in pigs.

Lepidoptera (moths and butterflies) are widely distributed in the world, but high-level phylogeny in Lepidoptera remains uncertain. More mitochondrial genome (mitogenome) data can help to conduct comprehensive analysis and construct a robust phylogenetic tree. Here, we sequenced and annotated 17 complete moth mitogenomes and made comparative analysis with other moths. The gene order of trnM-trnI-trnQ in 17 moths was different from trnI-trnQ-trnM of ancestral insects. The number, type, and order of genes were consistent with reported moths. The length of newly sequenced complete mitogenomes ranged from 14,231 bp of Rhagastis albomarginatus to 15,756 bp of Numenes albofascia. These moth mitogenomes were typically with high A+T contents varied from 76.0% to 81.7% and exhibited negative GC skews. Among 13 protein coding genes (PCGs), some unusual initiations and terminations were found in part of newly sequenced moth mitogenomes. Three conserved gene-overlapping regions and one conserved intergenic region were detected among 17 mitogenomes. The phylogenetic relationship of major superfamilies in Macroheterocera was as follows: (Bombycoidea + Lasiocampoidea) + ((Drepanoidea + Geometroidea) + Noctuoidea)), which was different from previous studies. Moreover, the topology of Noctuoidea as (Notodontidae + (Erebidae + Noctuidae)) was supported by high Bayesian posterior probabilities (BPP = 1.0) and bootstrapping values (BSV = 100). This study greatly enriched the mitogenome database of moth and strengthened the high-level phylogenetic relationships of Lepidoptera.

Microplastics (MPs) are pollutants widely distributed in the aquatic environments and causing various degrees of aquatic toxicity to aquatic organisms, which has attracted global attention in recent years. Nano-selenium (NSe) has been shown to have the potential to mitigate the harmful impacts of toxic substances. However, there is currently no reported evidence regarding the protective influence of NSe against the adverse effects of MPs. The aim of this study is to determine whether NSe could ameliorate the polystyrene (PS)-MPs-induced injury in grass carp (Ctenopharyngodon idella). The individuals of grass carp were assigned into three groups: (1) the control group fed with basal diet, (2) the PS group fed with basal diet and exposed to PS-MPs, and (3) the NSe group fed with diet supplemented with NSe and exposed to PS-MPs. Our results indicated that NSe administration significantly alleviated the histological damage caused by the PS-MPs in the liver and intestine with lower goblet cell count and larger villus height in the intestine, and significantly lower damage score in the liver. Moreover, NSe mitigated PS-MPs-induced oxidative stress through restoring the activities of antioxidant enzymes (superoxide dismutase (SOD), catalase (CAT), and malondialdehyde (MDA)) except the intestinal CAT activity. Furthermore, NSe supplementation could help fish maintain lower transcriptional level of the immune-related genes (Toll-like receptor 4 (TLR4) and myeloid differentiation factor 88 (MyD88)), inflammation-related genes (major histocompatibility complex class II (MHC-II) and interleukin 8 (IL-8)) and antioxidant enzyme-related genes (nuclear factor (erythroid-derived 2)-like 2 (Nrf-2) and kelch-like ECH-associated protein 1 (Keap-1)) after PS-MPs exposure. Besides, NSe supplementation dramatically helped maintain the intestinal microbial composition, for example, the proportion of Proteobacteria in the grass carp intestine of the NSe group (41 %) was similar to that of the control group (34 %) while 85 % of the PS group. NSe also played a significant protective role in intestinal microbial diversity, effectively resisting the damage on intestinal microbial diversity due to PS-MPs exposure. PS-MPs reduced the beneficial bacteria and increased the pathogenic microorganism like Aeromonas, which was undeniable signs of intestinal dysbiosis. Functional analysis indicated that PS-MPs affected intestinal microbiota functions like inhibition of metabolism, while NSe could significantly alleviate the damage. Our findings suggested that NSe could ameliorate PS-MPs-induced injury, which could contribute to the better understanding of the ecotoxicological effects of MPs on fish and help develop relevant mitigation strategies.

Moniezia expansa (Cyclophyllidea: Anoplocephalidae) is a large species of tapeworm that occurs in sheep and cattle and inhabits the small intestine, causing diarrhea and weight declines, leading to stockbreeding losses. Interestingly, the body fat percentage of M. expansa, which lacks the ability to synthesize fatty acids, is as high as 78% (dry weight) and all of the proglottids of M. expansa exhibit a dynamic developmental process from top to bottom. The aim of this paper is to identify the molecular basis of this high body fat percentage, the dynamic expression of developmental genes and their expression regulation patterns. From 12 different proglottids (four sections: scolex and neck, immature, mature and gravid with three replicates), 13,874 transcripts and 680 differentially expressed genes (DEGs) were obtained. The gene expression patterns of the scolex and neck and immature proglottids were very similar, while those of the mature and gravid proglottids differed greatly. In addition, 13 lipid transport-related proteins were found in the DEGs, and the expression levels showed an increasing trend in the four proglottid types. Furthermore, it was shown that 33 homeobox genes, 9 of which were DEGs, had the highest expression in the scolex and neck section. The functional enrichment results of the DEGs were predominantly indicative of development-related processes, and there were also some signal transduction and metabolism results. The most striking result was the finding of Wnt signaling pathways, which appeared multiple times. Furthermore, the weighted gene co-expression networks were divided into 12 modules, of which the brown module was enriched with many development-related genes. We hypothesize that M. expansa uses lipid transport-associated proteins to transport lipids from the host gut to obtain energy to facilitate its high fecundity. In addition, homeobox genes and Wnt signaling pathways play a core role in development and regeneration. The results promote research on the cell differentiation involved in the continuous growth and extension of body structures.

The giant panda (Ailuropoda melanoleuca) is the epitome of a flagship species for wildlife conservation and also an ideal model of adaptive evolution. As an obligate bamboo feeder, the giant panda relies on the olfaction for food recognition. The number of olfactory receptor (OR) genes and the rate of pseudogenes are the main factors affecting the olfactory ability of animals. In this study, we used the chromosome-level genome of the giant panda to identify OR genes and compared the genome sequences of OR genes with five other Ursidae species (spectacled bear (Tremarctos ornatus), American black bear (Ursus americanus), brown bear (Ursus arctos), polar bear (Ursus maritimus) and Asian black bear (Ursus thibetanus)). The giant panda had 639 OR genes, including 408 functional genes, 94 partial OR genes and 137 pseudogenes. Among them, 222 OR genes were detected and distributed on 18 chromosomes, and chromosome 8 had the most OR genes. A total of 448, 617, 582, 521 and 792 OR genes were identified in the spectacled bear, American black bear, brown bear, polar bear and Asian black bear, respectively. Clustering analysis based on the OR protein sequences of the six species showed that the OR genes distributed in 69 families and 438 subfamilies based on sequence similarity, and the six mammals shared 72 OR gene subfamilies, while the giant panda had 31 unique OR gene subfamilies (containing 35 genes). Among the 35 genes, there are 10 genes clustered into 8 clusters with 10 known human OR genes (OR8J3, OR51I1, OR10AC1, OR1S2, OR1S1, OR51S1, OR4M1, OR4M2, OR51T1 and OR5W2). However, the kind of odor molecules can be recognized by the 10 known human OR genes separately, which needs further research. The phylogenetic tree showed that 345 (about 84.56%) functional OR genes were clustered as Class-II, while only 63 (about 15.44%) functional OR genes were clustered as Class-I, which required further and more in-depth research. The potential odor specificity of some giant panda OR genes was identified through the similarity to human protein sequences. Sequences similar to OR2B1, OR10G3, OR11H6 and OR11H7P were giant panda-specific lacking, which may be related to the transformation and specialization from carnivore to herbivore of the giant panda. Since our reference to flavoring agents comes from human research, the possible flavoring agents from giant panda-specific OR genes need further investigation. Moreover, the conserved motifs of OR genes were highly conserved in Ursidae species. This systematic study of OR genes in the giant panda will provide a solid foundation for further research on the olfactory function and variation of the giant panda.

The schizothoracine fishes, widely distributed in the Qinghai-Tibetan Plateau and its adjacent areas, are considered as ideal models for investigation of high-altitude adaptation. Schizophygopsis are one group of the highly specialized schizothoracine fishes, and the genetic basis for their high-altitude adaptation is poorly understood. In this study, we performed comparative genomics analyses to investigate the potential genetic mechanisms for high-altitude adaptation of Schizopygopsis malacanthus and S. pylzovi based on the chromosome-level genomes. Functional enrichment analysis revealed that many expanded gene families in Schizopygopsis were associated with immune response while many contracted gene families were functionally associated with olfaction. Among the 123 positively selected genes (PSGs), angpt2a was detected in HIF-1 signaling pathway and possibly related to the hypoxia adaptation of Schizopygopsis. Furthermore, two PSGs cox15 and ndufb10 were distributed in thermogenesis, and there was a Schizopygopsis-specific missense mutation in cox15 (Gln115Glu), which possibly contributed to the cold temperature adaptation of the Schizopygopsis. KEGG enrichment of the PSGs revealed three significant pathways including metabolic pathways, cell cycle, and homologous recombination and GO enrichment analysis of the PSGs revealed several categories associated with DNA repair, cellular response to DNA damage stimulus, and metabolic process. Chromosome-scale characterization of olfactory receptor repertoires indicated that Schizopygopsis had the least number of OR genes, and the OR gene contraction was possibly caused by the limited food variety and the environmental factors such as lower air pressure, lower humidity, and lower temperature. Our study will help expand our understanding of the potential adaptive mechanism of Schizopygopsis to cope with the high-altitude conditions.