The widely observed RNA-DNA differences (RDDs) have been found to be due to nucleotide alteration by RNA editing. Canonical RNA editing (i.e., A-to-I and C-to-U editing) mediated by the adenosine deaminases acting on RNA (ADAR) family and apolipoprotein B mRNA editing catalytic polypeptide-like (APOBEC) family during the transcriptional process is considered common and essential for the development of an individual. To date, an increasing number of RNA editing sites have been reported in human, rodents, and some farm animals; however, genome-wide detection of RNA editing events in sheep has not been reported. The aim of this study was to identify RNA editing events in sheep by comparing the RNA-seq and DNA-seq data from three biological replicates of the kidney and spleen tissues. A total of 607 and 994 common edited sites within the three biological replicates were identified in the ovine kidney and spleen, respectively. Many of the RDDs were specific to an individual. The RNA editing-related genes identified in the present study might be evolved for specific biological functions in sheep, such as structural constituent of the cytoskeleton and microtubule-based processes. Furthermore, the edited sites found in the ovine BLCAP and NEIL1 genes are in line with those in previous reports on the porcine and human homologs, suggesting the existence of evolutionarily conserved RNA editing sites and they may play an important role in the structure and function of genes. Our study is the first to investigate RNA editing events in sheep. We screened out 607 and 994 RNA editing sites in three biological replicates of the ovine kidney and spleen and annotated 164 and 247 genes in the kidney and spleen, respectively. The gene function and conservation analysis of these RNA editing-related genes suggest that RNA editing is associated with important gene function in sheep. The putative functionally important RNA editing sites reported in the present study will help future studies on the relationship between these edited sites and the genetic traits in sheep.

Black bone sheep was first discovered in Yunnan province of China in 1970, with unique black pigmentation on the body and internal organs. Endothelin 3 (EDN3) has been known as a key gene causing hyperpigmentation in black bone chicken, the Silky fowl. In this study, EDN3 was employed as a candidate gene for regulating black color pigmentation. First, EDN3 was cloned from sheep to obtain the full-length cDNA by using the rapid amplification of cDNA ends (RACE). Genomic EDN3 was screened and a total of thirty predicted single nucleotide polymorphisms (SNPs) were genotyped for allele and genotype frequency analysis in a case-control study involving two black bone sheep populations. Genomic copy number analysis of EDN3 in sheep was conducted to measure the variation in copy number. EDN3 expression levels were observed among the groups in adult liver, lymph node, and kidney tissues, as well as embryo kidney samples. Also, among the tissues of black bone and non-black bone sheep. The size of the full-length cDNA was 1,578 bp, which included 426 bp of 5′-untranslated region (5′-UTR), an open reading frame (ORF) of 639 bp encoding a protein of 212 amino acids, and a 3′-UTR of 513 bp. Genotype and allele frequencies of all the discovered SNPs were found insignificantly different in black bone and non-black bone sheep (P > 0.05). Genomic copy number analysis of EDN3 in sheep revealed no significant difference between the two sheep groups. No significant variations were found in the adult liver and kidney embryo samples. However, the expression in lymph node and kidney tissue was significantly higher in black bone sheep than that in non-black bone sheep (P < 0.05). Significant variations in the EDN3 expression levels were observed among the tissues of non-black bone sheep. The findings of the present study indicate that unlike in Silky chickens, EDN3 is not responsible for hyperpigmentation but may play a key functional role in immune and excretory systems of black bone sheep.

The eggshell blueness is an interesting object for chicken genetic studies and blue-shelled chicken industry, especially after the discovery of the causative mutation of chicken blue eggshell. In the present study, genome wide association study (GWAS) was conducted in Chinese Dongxiang blue-shelled chicken underlying four traits of blue eggshell pigments: quantity of biliverdin (QB), quantity of protoporphyrin (QP), quantity of total pigment (QT), and color density trait (CD). A total of 139 individuals were randomly collected for GWAS. We detected two SNPs in genome-wise significance and 35 in suggestive significance, 24 out of the 37 SNP were located either within intron/exon or near 15 genes in a range of ~1.17 Mb on GGA21. For further confirmation of the identified SNP loci by GWAS, the follow-up replication studies were performed in two populations. A total of 146 individuals of the second generation derived from the former GWAS population, as well as 280 individuals from an alternative independent population were employed for genotyping by MALDI-TOF MS in a genotype-phenotype association study. Eighteen SNPs evenly distributed on the GGA21 significant region were successfully genotyped in the two populations, of which 4 and 6 SNP loci were shown significantly associated with QB, QT and QP in the two repeat populations, respectively. Further, the SNPs were narrowed down to a region of ~ 653.819 Kb on GGA21 that harbors five candidate genes: AJAP1, TNFRSF9, C1ORF174, CAMTA1, and CEP104. Shell gland of chickens laying dark and light blue eggshell was chosen for detection of mRNA expression of the five candidate genes. The results showed differential expression levels of these genes in the two groups. The specific function of these genes has not yet been defined clearly in chickens and further in-depth studies are needed to explore the new functional role in chicken eggshell blueness.

Heat-stress resilience is vital for poultry in tropical/subtropical regions where high temperatures impair productivity. Ake chickens, as the only naked-neck chicken breed in China, exhibit robust resistance to heat stress, but this breed lacks clarity in its genetic origins. This study utilized the next-generation sequencing data from 22 chicken breeds to conduct phylogenetic and population analyses. Gene flow analysis revealed a gene migration event from Iranian naked-neck chickens and Indian local breeds to Ake chickens, and population separation estimates suggested that the naked-neck gene was introduced to China around 500–600 years ago. NJ-tree, PCA, and population structure analyses showed that Ake chickens cluster with Yunnan native breeds, which diverged only 100–200 years ago. A selective sweep in the candidate region on chromosome 3 (97.0–97.37 Mb) showed elevated genetic differentiation (F_ST) and reduced nucleotide diversity (π) compared to the genome-wide average, indicating rapid fixation of the trait under natural/artificial selection. Demographic reconstruction indicated that the current effective size of Ake chickens is stable at 2000–3000 individuals. These findings deepen our understanding of Ake chicken evolution and provide valuable insights for conservation and the development of heat-stress-resistant poultry breeds.