Search

With their high economic value and cultural significance, modern roses are one of the most important ornamental plants. Because of their complicated genetic background and tetraploid nature, the creation of high-density genetic maps of roses has been a challenge that has slowed the pace of molecular breeding for modern roses. The current construction of tetraploid genetic maps based on existing diploid rose genomes could lead to inaccurate marker information and genotyping results. Therefore, we generated the first high-quality tetraploid genome of Rosa chinensis ‘Yunzheng Xiawei.’ Utilizing Illumina, PacBio, and Hi-C sequencing technologies, we assembled a genome of 858.59 Mb with 14 pseudo-chromosomes. Mode of inheritance analysis using PolyOrigin indicated that modern roses show both quadrivalent and bivalent pairing. Based on this reference genome, high-density genetic maps were constructed using MSTmap with nearly saturated markers. Quantitative trait locus (QTL) analysis was conducted using WinQTLCart and R/qtl for flavonoids and carotenoids, and 11 QTL clusters were identified. By combining the genome annotation, phylogenetic analyses, and gene expression analyses, we were able to identify several key genes related to flavonoid and carotenoid biosynthesis. This study provides the basis for further genetic analyses of highly heterozygous tetraploid roses and could facilitate the progress of marker-assisted selection in modern roses.

Flavonoid glycosides are widespread in plants, and are of great interest owing to their diverse biological activities and effectiveness in preventing chronic diseases. Periploca forrestii, a renowned medicinal plant of the Apocynaceae family, contains diverse flavonoid glycosides and is clinically used to treat rheumatoid arthritis and traumatic injuries. However, the mechanisms underlying the biosynthesis of these flavonoid glycosides have not yet been elucidated. In this study, we used widely targeted metabolomics and full-length transcriptome sequencing to identify flavonoid diversity and biosynthetic genes in P. forrestii. A total of 120 flavonoid glycosides, including 21 C-, 96 O-, and 3 C/O-glycosides, were identified and annotated. Based on 24,123 full-length coding sequences, 99 uridine diphosphate sugar-utilizing glycosyltransferases (UGTs) were identified and classified into 14 groups. Biochemical assays revealed that four UGTs exhibited O-glycosyltransferase activity toward apigenin and luteolin. Among them, PfUGT74B4 and PfUGT92A8 were highly promiscuous and exhibited multisite O-glycosylation or consecutive glycosylation activities toward various flavonoid aglycones. These four glycosyltransferases may significantly contribute to the diversity of flavonoid glycosides in P. forrestii. Our findings provide a valuable genetic resource for further studies on P. forrestii and insights into the metabolic engineering of bioactive flavonoid glycosides.

Abstract Horse chestnut (Aesculus chinensis) is an important medicinal tree that contains various bioactive compounds, such as aescin, barrigenol-type triterpenoid saponins (BAT), and aesculin, a glycosylated coumarin. Herein, we report a 470.02 Mb genome assembly and characterize an Aesculus-specific whole-genome duplication event, which leads to the formation and duplication of two triterpenoid biosynthesis-related gene clusters (BGCs). We also show that AcOCS6, AcCYP716A278, AcCYP716A275, and AcCSL1 genes within these two BGCs along with a seed-specific expressed AcBAHD6 are responsible for the formation of aescin. Furthermore, we identify seven Aesculus-originated coumarin glycoside biosynthetic genes and achieve the de novo synthesis of aesculin in E. coli. Collinearity analysis shows that the collinear BGC segments can be traced back to early-diverging angiosperms, and the essential gene-encoding enzymes necessary for BAT biosynthesis are recruited before the splitting of Aesculus, Acer, and Xanthoceras. These findings provide insight on the evolution of gene clusters associated with medicinal tree metabolites.

Holarrhena pubescens Wall. ex G. Don, 1837 is an important medicinal plant belonging to the Holarrhena genus in the Apocynaceae family. In this study, the complete chloroplast (cp) genome sequence of H. pubescens was sequenced using the Illumina NovaSeq platform. The cp genome of H. pubescens was 160,108 bp in length with 37.21% overall GC content. The cp genome of H. pubescens containing a large single-copy region (LSC, 88,685 bp), a small single-copy region (SSC, 18,671 bp), and a pair of inverted repeat regions (SSC, 26,376 bp). The cp genome encoded 129 genes, including 84 protein-coding genes, 37 tRNA genes, and eight rRNA genes. Phylogenetic analysis based on complete protein coding genes sequences revealed that H. pubescens was closest to Beaumontia murtonii.

The basal eudicot Epimedium (Barrenwort) exhibits innovative floral morphology in the form of petal spurs filled with nectar and petaloid sepals. The B-class MADS-box genes APETALAS3 (AP3) and PISTILLA (PI) determine sepal and petal identity while the C-class gene AGAMOUS (AG) determines carpel and stamen identity in Arabidopsis. Complex histories of gene duplication resulted in subsequent subfunctionalization or neofunctionalization of paralogs. Here, a total of four B- and two C-class genes were successfully isolated from E. sagittatum. Phylogenetic analysis showed that EsAP3-1, EsAP3-2, and EsAP3-3 are part of the AP3 clade; EsPI is part of the PI group; and EsAG and EsAG11 clustered into the AG and AG11 groups, respectively. Quantitative real-time PCR was utilized to detect the expression patterns of these genes, and all B-class genes except for EsAP3-3 were found to be universally expressed. The transcribed EsAG and EsAG11 genes were confined to reproductive organs. In addition, yeast three and two hybrid assays were used to explore the status of protein complexes. EsAP3-1 was found to have broadly interactive partners, and EsPI can form a heterodimer with EsAP3-1 and EsAP3-2. Transgenic EsAG overexpression in wild-type Arabidopsis confirmed the conserved function determining carpel development.

Abstract Background Flavonoids are necessary for plant growth and resistance to adversity and stress. They are also an essential nutrient for human diet and health. Among the metabolites produced in Cannabis sativa (C. sativa), phytocannabinoids have undergone extensive research on their structures, biosynthesis, and biological activities. Besides the phytocannabinoids, C. sativa is also rich in terpenes, alkaloids, and flavonoids, although little research has been conducted in this area. Results In this study, we identified 11 classes of key enzyme-encoding genes, including 56 members involved in the flavonoid biosynthesis in C. sativa, from their physical characteristics to their expression patterns. We screened the potentially step-by-step enzymes catalyzing the precursor phenylalanine to the end flavonoids using a conjoin analysis of gene expression with metabolomics from different tissues and chemovars. Flavonol synthase (FLS), belonging to the 2-oxoglutarate-dependent dioxygenase (2-ODD) superfamily, catalyzes the dihydroflavonols to flavonols. In vitro recombinant protein activity analysis revealed that CsFLS2 and CsFLS3 had a dual function in converting naringenin (Nar) to dihydrokaempferol (DHK), as well as dihydroflavonols to flavonols with different substrate preferences. Meanwhile, we found that CsFLS2 produced apigenin (Api) in addition to DHK and kaempferol when Nar was used as the substrate, indicating that CsFLS2 has an evolutionary relationship with Cannabis flavone synthase I. Conclusions Our study identified key enzyme-encoding genes involved in the biosynthesis of flavonoids in C. sativa and highlighted the key CsFLS genes that generate flavonols and their diversified functions in C. sativa flavonoid production. This study paves the way for reconstructing the entire pathway for C. sativa’s flavonols and cannflavins production in heterologous systems or plant culture, and provides a theoretical foundation for discovering new cannabis-specific flavonoids.

Quinoa grains are gaining increasing popularity owing to their high nutritional merits. However, only limited information is available on the metabolic profiles of quinoa grains. In this study, we determined the metabolic profiles of black, red, and white quinoa grains via an ultraperformance liquid chromatography–electrospray ionization–tandem mass spectrometry (UPLC–ESI–MS/MS)-based metabolomics. A total of 689 metabolites were identified, among which 251, 182, and 317 metabolites displayed different accumulation patterns in the three comparison groups (Black vs Red, Black vs White, and Red vs White), respectively. In particular, flavonoid and phenolic acid contents displayed considerable differences, with 22 flavonoids, 5 phenolic acids, and 1 betacyanin being differentially accumulated among the three quinoa cultivars. Additionally, correlation analysis showed that flavonoids and phenolic acids could act as betanin co-pigments in quinoa grains. In conclusion, this study provides comprehensive insights into the adequate utilization and development of novel quinoa-based functional foods.

IntroductionFlower color is an ideal trait for studying the molecular basis for phenotypic variations in natural populations of species. Epimedium (Berberidaceae) species exhibit a wide range of flower colors resulting from the varied accumulation of anthocyanins and other pigments in their spur-like petals and petaloid sepals.MethodsIn this work, the anthocyanidins of eight different Epimedium species with different floral pigmentation phenotypes were analyzed using HPLC. Twelve genes involved in anthocyanin biosynthesis were cloned and sequenced, and their expression was quantified.ResultsThe expression levels of the catalytic enzyme genes DFR and ANS were significantly decreased in four species showing loss of floral pigmentation. Complementation of EsF3’H and EsDFR in corresponding Arabidopsis mutants together with overexpression of EsF3’5’H in wild type Arabidopsis analysis revealed that these genes were functional at the protein level, based on the accumulation of anthocyanin pigments.DiscussionThese results strongly suggest that transcriptional regulatory changes determine the loss of anthocyanins to be convergent in the floral tissue of Epimedium species.

The medicinal plant Cannabis sativa L. (C. sativa) accumulates plant cytotoxic but medicinally important cannabinoids in glandular trichomes and flowers of female plants. Although the major biosynthetic pathway of cannabinoids has been revealed, their transportation mechanism is still unknown. Multidrug and toxic compound extrusion proteins (MATEs) can transport plant metabolites, ions and phytohormones intra and inter-cellularly. MATEs could have the potential to translocate cannabinoids or their synthetic intermediates to cellular compartment, thus protecting them from unwanted modifications and cytotoxicity. In this study, we performed a genome-wide identification and expression analysis of Cannabis sativa MATEs (CsMATEs) and revealed 42 CsMATEs that were classified phylogenetically into four conserved subfamilies. Forty-two CsMATEs were unevenly distributed on 10 chromosomes, with 50% CsMATEs were physically adjacent to at least one another CsMATEs and 83% CsMATEs localized on plasma membrane. Tandem duplication is the major evolutionary driving force for CsMATEs expansion. Real-time quantitative PCR revealed CsMATE23, CsMATE28 and CsMATE34 mainly expressed in flower, whereas CsMATE17 and CsMATE27 showed strong transcription in root. Light responsive cis-acting element was most abundant in promoters of CsMATE23, CsMATE28 and CsMATE34. Finally, the contents of cannabinoids and corresponding biosynthetic intermediates as well as expressions of CsMATE28 and CsMATE34 were determined under UV-B treatment, among which strong correlation was found. Our results indicates that CsMATEs might involve in biosynthesis of cannabinoids and has the potential to be used in heterologous production of cannabinoids.

Tartary buckwheat sprouts have a high nutritional value and are gluten-free, and polyphenols are their main active constituents. However, information regarding the active constituents’ difference of Tartary buckwheat sprouts grown from seeds with different morphology, at different developmental stages and environments is limited. Here, we developed a LC–MS-based targeted metabolomics approach to analyze polyphenols (46 flavonoids and 6 anthraquinones) in 40 Tartary buckwheat sprouts varieties. Both flavonoids and anthraquinones contributed to significant differences in sprouts grown from seed with different color or shape. Twenty-seven differential compounds were all at a higher level in 3-day-old sprouts, and the fold change from 3-day-old to 8-day-old sprouts was 1.42–6.64. A total of 25 differential compounds were all significantly upregulated upon UV-B radiation, especially for epicatechin. This study is valuable not only for better breeding cultivars of Tartary buckwheat sprouts, but also assessing their metabolic quality.

Histone deacetylases (HDACs) play crucial roles nearly in all aspects of plant biology, including stress responses, development and growth, and regulation of secondary metabolite biosynthesis. The molecular functions of HDACs have been explored in depth in Arabidopsis thaliana, while little research has been reported in the medicinal plant Cannabis sativa L. Here, we excavated 14 CsHDAC genes of C. sativa L that were divided into three relatively conserved subfamilies, including RPD3/HDA1 (10 genes), SIR2 (2 genes), and HD2 (2 genes). Genes associated with the biosynthesis of bioactive constituents were identified by combining the distribution of cannabinoids with the expression pattern of HDAC genes in various organs. Using qRT-PCR and transcription group analysis, we verified the expression of candidate genes in different tissues. We found that the histone inhibitor Trichostatin A (TSA) affected the expression of key genes in the cannabinoid metabolism pathway and the accumulation of synthetic precursors, which indirectly indicates that histone inhibitor may regulate the synthesis of active substances in C. sativa L.

Roses are highly valuable within the flower industry. The metabolites of anthocyanins, flavonols, and carotenoids in rose petals are not only responsible for the various visible petal colors but also important bioactive compounds that are important for human health. In this study, we performed a QTL analysis on pigment contents to locate major loci that determine the flower color traits. An F1 population of tetraploid roses segregating for flower color was used to construct an ultra-high-density genetic linkage map using whole-genome resequencing technology to detect genome-wide SNPs. Previously developed SSR and SNP markers were also utilized to increase the marker density. Thus, a total of 9,259 markers were mapped onto seven linkage groups (LGs). The final length of the integrated map was 1285.11 cM, with an average distance of 0.14 cM between adjacent markers. The contents of anthocyanins, flavonols and carotenoids of the population were assayed to enable QTL analysis. Across the 33 components, 46 QTLs were detected, explaining 11.85–47.72% of the phenotypic variation. The mapped QTLs were physically clustered and primarily distributed on four linkage groups, namely LG2, LG4, LG6, and LG7. These results improve the basis for flower color marker-assisted breeding of tetraploid roses and guide the development of rose products.

Abstract The floral transition is a crucial developmental event, but little is known about the underlying regulatory networks in seasonally and continuously flowering roses. In this study, we compared the genetic basis of flowering in two rose species, Rosa chinensis ‘Old Blush’, which flowers continuously, and R. odorata var. gigantea, which blooms in early spring. Gene ontology (GO) terms related to methylation, light reaction, and starch metabolism were enriched in R. odorata var. gigantea and terms associated with sugar metabolism were enriched in R. chinensis ‘Old Blush’ during the floral transition. A MapMan analysis revealed that genes involved in hormone signaling mediate the floral transition in both taxa. Furthermore, differentially expressed genes (DEGs) involved in vernalization, photoperiod, gibberellin (GA), and starch metabolism pathways converged on integrators, e.g., LFY, AGL24, SOC1, CAL, and COLs, to regulate the floral transition in R. odorata var. gigantea, while DEGs related to photoperiod, sugar metabolism, and GA pathways, including COL16, LFY, AGL11, 6PGDH, GASA4, and BAM, modulated the floral transition in R. chinensis ‘Old Blush.’ Our analysis of the genes underlying the floral transition in roses with different patterns of flowering provides a basis for further functional studies.

Rose is one of the most valuable ornamental crops worldwide. In this study, the composition of hydrophilic and lipophilic pigments in petals of six rose cultivars at seven developing stages was investigated using high performance liquid chromatography and mass spectrometry. Four anthocyanins, 20 flavonols, and 10 carotenoids were detected in petals of tested cultivars. Major individual anthocyanin, flavonol, and carotenoid were cyanidin/pelargonidin 3,5-diglucoside, kaempferol 3-O-rhamnoside, and (9Z)-violaxanthin, respectively. Significant differences were observed in pigments content in petals of different rose cultivars. The yellow petals of YI and GC exhibited no to very small amounts of anthocyanins, moderate amount of total flavonols, and highest content of total carotenoids. Similarly, pink petals of PF, WQ, and YX showed average concentration of total anthocyanins, highest concentration of total flavonols, and small amount of carotenoids. Further, orange petals of CH showed highest content of total anthocyanins, lowest content of total flavonols, and average content of total carotenoids. Correlation analysis demonstrated that there were many pigments influencing petal colors. Moreover, multiple linear regression indicated that pelargonidin 3,5-diglucoside, total anthocyanins and (9Z)-violaxanthin were the major factors. In addition, this study showed that orange cultivar CH, pink cultivar PF and yellow cultivar YI can have great potential as a natural source for the extraction of pelargonidin 3-O-glucoside, kaempferol 3-O-rhamnoside, and (9Z)-violaxanthin, respectively. These investigations would contribute toward understanding the mechanism on the development of flower colors and provide a theoretical basis for the breeding of rose with specific color.

Rosa chinensis, which is a famous traditional flower in China, is a major ornamental plant worldwide. Long-term cultivation and breeding have resulted in considerable changes in the number of rose petals, while most wild Rosaceae plants have only one whorl consisting of five petals. The petals of double flowers reportedly originate from stamens, but the underlying molecular mechanism has not been fully characterized. In this study, we observed that the number of petals of R. chinensis ‘Old Blush’ flowers increased and decreased in response to low- and high-temperature treatments, respectively, similar to previous reports. We characterized these variations in further detail and found that the number of stamens exhibited the opposite trend. We cloned an APETALA2 homolog, RcAP2. A detailed analysis of gene structure and promoter cis-acting elements as well as RcAP2 temporospatial expression patterns and responses to temperature changes suggested that RcAP2 expression may be related to the number of petals from stamen origin. The overexpression of RcAP2 in Arabidopsis thaliana transgenic plants may induce the transformation of stamens to petals, thereby increasing the number of petals. Moreover, silencing RcAP2 in ‘Old Blush’ plants decreased the number of petals. Our results may be useful for clarifying the temperature-responsive mechanism involved in petaloid stamen production, which may be relevant for the breeding of new rose varieties with enhanced flower traits.

Abstract Background Chinese indigenous sheep are valuable resources with unique features and characteristics. They are distributed across regions with different climates in mainland China; however, few reports have analyzed the environmental adaptability of sheep based on their genome. We examined the variants and signatures of selection involved in adaptation to extreme humidity, altitude, and temperature conditions in 173 sheep genomes from 41 phenotypically and geographically representative Chinese indigenous sheep breeds to characterize the genetic basis underlying environmental adaptation in these populations. Results Based on the analysis of population structure, we inferred that Chinese indigenous sheep are divided into four groups: Kazakh (KAZ), Mongolian (MON), Tibetan (TIB), and Yunnan (YUN). We also detected a set of candidate genes that are relevant to adaptation to extreme environmental conditions, such as drought-prone regions (TBXT, TG, and HOXA1), high-altitude regions (DYSF, EPAS1, JAZF1, PDGFD, and NF1) and warm-temperature regions (TSHR, ABCD4, and TEX11). Among all these candidate genes, eight ABCD4, CNTN4, DOCK10, LOC105608545, LOC121816479, SEM3A, SVIL, and TSHR overlap between extreme environmental conditions. The TSHR gene shows a strong signature for positive selection in the warm-temperature group and harbors a single nucleotide polymorphism (SNP) missense mutation located between positions 90,600,001 and 90,650,001 on chromosome 7, which leads to a change in the protein structure of TSHR and influences its stability. Conclusions Analysis of the signatures of selection uncovered genes that are likely related to environmental adaptation and a SNP missense mutation in the TSHR gene that affects the protein structure and stability. It also provides information on the evolution of the phylogeographic structure of Chinese indigenous sheep populations. These results provide important genetic resources for future breeding studies and new perspectives on how animals can adapt to climate change.

AbstractContext di-(2-Ethylhexyl) phthalate (DEHP) has potential reproductive toxicity. Bu-Shen-Tian-Jing formulations (BSTJFs) are beneficial for female reproductive capacity. However, BSTJF2 has much lower cytotoxicity than BSTJF1.Objective To investigate the effects of BSTJFs on ovarian granulosa cells exposed to DEHP and determine the potential molecular mechanisms.Methods and materials Human granulosa-like tumor cell line (KGN) cells were divided into control, DEHP, BSTJF1 and BSTJF2 groups. The DEHP group were given 1 μM DEHP for 24 h. They were then given BSTJF1 at 200 μg/mL or BSTJF2 at 100 μg/mL for 24 h. The control group was treated with the same concentration of DMSO (0.1%). Oxidative stress and mitochondrial function were measured. The mRNA and protein expression levels of HDAC3 and HSP90AA were determined. Integrative network pharmacology analysis of BSTJF2 was also performed.Results DEHP (1 μM) significantly suppressed the proliferation of KGN cells by 17%, significantly increased ROS levels by 28% and MDA levels by 47%, significantly decreased MMP levels by 22% and mtDNA copy by 30%. DEHP significantly increased protein expression of HDAC3 by 21%and HSP90AA by 64%. All these changes were significantly reversed by BSTJFs. Integrative network pharmacology analysis revealed HSP90AA was a key target (degree = 8). Both RGFP966 and BSTJF2 significantly reversed the increased expression of HDAC3 and HSP90AA, attenuated oxidative stress, and mitochondrial damage which were induced by DEHP.Conclusion BSTJFs might have therapeutic potential on oxidative stress and mitochondrial damage through the HDAC3/HSP90AA pathway which encourages further clinical trials.

Due to the uncertain, diverse, and light-attenuating characteristics of the underwater environment, underwater images have low contrast and unclear problems. This paper proposes a histogram equalization algorithm based on optimized adaptive image quadruple segmentation and cropping (AQSCHE). Compared with the traditional histogram equalization underwater image enhancement algorithm, this algorithm introduces histogram quadruple segmentation and cropping technology. Using the exposure value and segmentation point calculation formula that optimizes the distribution range of the histogram, perform quadruple segmentation on the image to obtain a more refined histogram. The adaptive histogram clipping is realized by constructing the clipping parameter z to adjust the contrast and brightness of the image. The original image is enhanced by double equalization of the sub-histogram and the histogram of each channel. Finally, the simulation experiments verify the enhancement effect of the proposed algorithm AQSCHE on underwater images. The processed underwater image has higher contrast, is clearer and more natural in subjective evaluation, and has a better visual effect; in the image objective evaluation indicators, information entropy (Entropy), peak signal-to-noise ratio (PSNR), structural similarity index (SSIM) and universal color image quality evaluator (UCIQE), this algorithm also outperforms other common algorithms such as HE and CLAHE.

Due to the advantages of Mg on inclusions and carbides in H13 steel, Mg wire was added into industrial ladle after VD process to replace traditional Ca treatment in this paper, which is of great significance for industrial production of H13 steel. This paper focuses on the effects of Mg treatment and traditional Ca treatment on the inclusions, carbides, microstructure and mechanical properties in H13 steel. The results show that feeding Mg wire during EAF→ LF→ VD → CC steelmaking process can ensure the submerged nozzle not being blocked under the condition of continuous casting 15 heats of steel ladle, which means feeding Ca wire can be replaced with Mg wire from the aspect of steelmaking production. The primary carbides in Mg treatment steel are mainly precipitated with MgO or MgAl2O4 as the cores and small in size, while the primary carbides in Ca treatment steel are almost large-size carbides without core. The samples taken from the intermediate ladle show that the primary carbide types are all V-rich, Mo-rich and V–Mo-rich composite, and the size of primary carbides in Mg treatment steel is much smaller than that in Ca treatment steel. Adding Mg is helpful to improve the equiaxed crystal rate of continuous casting ingot. The annealed hot rolling samples show that Mg treatment could interrupt the M23C6 (M = Cr, Fe) chain carbides distributed along the grain boundary. The impact test shows that the dimple fracture surface become small and deep after adding Mg, which improves impact toughness significantly.

Compared to Chinese indigenous sheep, Western sheep have rapid growth rate, larger physique, and higher meat yield. These excellent Western sheep were introduced into China for crossbreeding to expedite the enhancement of production performance and mutton quality in local breeds. Here, we investigated population genetic structure and genome-wide selection signatures among the Chinese indigenous sheep and the introduced sheep based on whole-genome resequencing data. The PCA, N-J tree and ADMIXTURE results showed significant genetic difference between Chinese indigenous sheep and introduced sheep. The nucleotide diversity (π) and linkage disequilibrium (LD) decay results indicated that the genomic diversity of introduced breeds were lower. Then, Fst & π ratio, XP-EHH, and de-correlated composite of multiple signals (DCMS) methods were used to detect the selection signals. The results showed that we identified important candidate genes related to growth rate and body size in the introduced breeds. Selected genes with stronger selection signatures are associated with growth rate (CRADD), embryonic development (BVES, LIN28B, and WNT11), body size (HMGA2, MSRB3, and PTCH1), muscle development and fat metabolism (MSTN, PDE3A, LGALS12, GGPS1, and SAR1B), wool color (ASIP), and hair development (KRT71, KRT74, and IRF2BP2). Thus, these genes have the potential to serve as candidate genes for enhancing the growth traits of Chinese indigenous sheep. We also identified tail-length trait-related candidate genes (HOXB13, LIN28A, PAX3, and VEGFA) in Chinese long-tailed breeds. Among these genes, HOXB13 is the main candidate gene for sheep tail length phenotype. LIN28A, PAX3, and VEGFA are related to embryonic development and angiogenesis, so these genes may be candidate genes for sheep tail type traits. This study will serve as a foundation for further genetic improvement of Chinese indigenous sheep and as a reference for studies related to growth and development of sheep.

There is a genetic difference between Hu sheep (short/fat-tailed sheep) and Tibetan sheep (short/thin-tailed sheep) in tail type, because of fat metabolism. Previous studies have mainly focused directly on sheep tail fat, which is not the main organ of fat metabolism. The function of miRNAs in sheep liver fat metabolism has not been thoroughly elucidated. In this study, miRNA-Seq was used to identify miRNAs in the liver tissue of three Hu sheep (short/fat-tailed sheep) and three Tibetan sheep (short/thin-tailed sheep) to characterize the differences in fat metabolism of sheep. In our study, Hu sheep was in a control group, we identified 11 differentially expressed miRNAs (DE miRNAs), including six up-regulated miRNAs and five down-regulated miRNAs. Miranda and RNAhybrid were used to predict the target genes of DE miRNAs, obtaining 3,404 target genes. A total of 115 and 67 GO terms as well as 54 and 5 KEGG pathways were significantly (padj < 0.05) enriched for predicted 3,109 target genes of up-regulated and 295 target genes of down-regulated miRNAs, respectively. oar-miR-432 was one of the most up-regulated miRNAs between Hu sheep and Tibetan sheep. And SIRT1 is one of the potential target genes of oar-miR-432. Furthermore, functional validation using the dual-luciferase reporter assay indicated that the up-regulated miRNA; oar-miR-432 potentially targeted sirtuin 1 (SIRT1) expression. Then, the oar-miR-432 mimic transfected into preadipocytes resulted in inhibited expression of SIRT1. This is the first time reported that the expression of SIRT1 gene was regulated by oar-miR-432 in fat metabolism of sheep liver. These results could provide a meaningful theoretical basis for studying the fat metabolism of sheep.

Wool color is controlled by a variety of genes. Although the gene regulation of some wool colors has been studied in relative depth, there may still be unknown genetic variants and control genes for some colors or different breeds of wool that need to be identified and recognized by whole genome resequencing. Therefore, we used whole genome resequencing data to compare and analyze sheep populations of different breeds by population differentiation index and nucleotide diversity ratios (Fst and θπ ratio) as well as extended haplotype purity between populations (XP-EHH) to reveal selection signals related to wool coloration in sheep. Screening in the non-white wool color group (G1 vs. G2) yielded 365 candidate genes, among which PDE4B, GMDS, GATA1, RCOR1, MAPK4, SLC36A1, and PPP3CA were associated with the formation of non-white wool; an enrichment analysis of the candidate genes yielded 21 significant GO terms and 49 significant KEGG pathways (p < 0.05), among which 17 GO terms and 21 KEGG pathways were associated with the formation of non-white wool. Screening in the white wool color group (G2 vs. G1) yielded 214 candidate genes, including ABCD4, VSX2, ITCH, NNT, POLA1, IGF1R, HOXA10, and DAO, which were associated with the formation of white wool; an enrichment analysis of the candidate genes revealed 9 significant GO-enriched pathways and 19 significant KEGG pathways (p < 0.05), including 5 GO terms and 12 KEGG pathways associated with the formation of white wool. In addition to furthering our understanding of wool color genetics, this research is important for breeding purposes.

Wool fineness affects the quality of wool, and some studies have identified about forty candidate genes that affect sheep wool fineness, but these genes often reveal only a certain proportion of the variation in wool thickness. We further explore additional genes associated with the fineness of sheep wool. Whole-genome resequencing of eight sheep breeds was performed to reveal selection signals associated with wool fineness, including four coarse wool and four fine/semi-fine wool sheep breeds. Multiple methods to reveal selection signals (Fst and θπ Ratio and XP-EHH) were applied for sheep wool fineness traits. In total, 269 and 319 genes were annotated in the fine wool (F vs. C) group and the coarse wool (C vs. F) group, such as LGR4, PIK3CA, and SEMA3C and NFIB, OPHN1, and THADA. In F vs. C, 269 genes were enriched in 15 significant GO Terms (p < 0.05) and 38 significant KEGG Pathways (p < 0.05), such as protein localization to plasma membrane (GO: 0072659) and Inositol phosphate metabolism (oas 00562). In C vs. F, 319 genes were enriched in 21 GO Terms (p < 0.05) and 16 KEGG Pathways (p < 0.05), such as negative regulation of focal adhesion assembly (GO: 0051895) and Axon guidance (oas 04360). Our study has uncovered genomic information pertaining to significant traits in sheep and has identified valuable candidate genes. This will pave the way for subsequent investigations into related traits.

The fat tail is a unique characteristic of sheep that represents energy reserves and is a complex adaptative mechanism of fat-tailed sheep to environmental stress. MicroRNA plays a significant role as regulators at the posttranscriptional level, but no studies have explained the molecular mechanisms of miRNA which regulate fat deposition in sheep tails. In this study, mRNA and miRNA analysis examined tail fat tissue from three Hu fat-tailed and three Tibetan thin-tailed sheep. After aligning to the reference sequences, 2,108 differentially expressed genes and 105 differential expression miRNAs were identified, including 1,247 up- and 861 downregulated genes and 43 up- and 62 downregulated miRNAs. Among these differentially expressed miRNAs, oar-miR-432 was one of the most downregulated miRNAs between Hu sheep and Tibetan sheep, and 712 genes were predicted to be targeted by oar-miR-432, 80 of which overlapped with DEGs. The Gene Ontology analysis on these genes showed that BMP2, LEP, GRK5, BMP7, and RORC were enriched in fat cell differentiation terms. The genes for BMP2 targeted by oar-miR-432 were examined using dual-luciferase assay. The oar-miR-432 mimic transfected into preadipocytes resulted in increased expression of BMP2. The marker gene PPAR-γ of fat differentiation had a lower expression than the negative control on days 0, 2, and 4 after induced differentiation. The decrease in the number of lipids in the oar-miR-432 mimic group detected by oil red O stain was also less than that in the negative control. This is the first study to reveal the fat mechanisms by which oar-miR-432 inhibits fat differentiation and promotes the expression of BMP2 in sheep tails.