Search

A comprehensive transcriptomic survey of pigs can provide a mechanistic understanding of tissue specialization processes underlying economically valuable traits and accelerate their use as a biomedical model. Here we characterize four transcript types (lncRNAs, TUCPs, miRNAs, and circRNAs) and protein-coding genes in 31 adult pig tissues and two cell lines. We uncover the transcriptomic variability among 47 skeletal muscles, and six adipose depots linked to their different origins, metabolism, cell composition, physical activity, and mitochondrial pathways. We perform comparative analysis of the transcriptomes of seven tissues from pigs and nine other vertebrates to reveal that evolutionary divergence in transcription potentially contributes to lineage-specific biology. Long-range promoter–enhancer interaction analysis in subcutaneous adipose tissues across species suggests evolutionarily stable transcription patterns likely attributable to redundant enhancers buffering gene expression patterns against perturbations, thereby conferring robustness during speciation. This study can facilitate adoption of the pig as a biomedical model for human biology and disease and uncovers the molecular bases of valuable traits., A comprehensive transcriptomic survey of the pig could enable mechanistic understanding of tissue specialization and accelerate its use as a biomedical model. Here the authors characterize four distinct transcript types in 31 adult pig tissues to dissect their distinct structural and transcriptional features and uncover transcriptomic variability related to tissue physiology.

Enterocytozoon bieneusi is an obligate intracellular parasite and the most common pathogen of microsporidiosis in humans and animals. In this study, a total of 198 fecal samples were collected from 11 captive populations of Rhesus macaque in Chinese zoos, to investigate the prevalence and analyze the zoonotic potential of E. bieneusi by genotype of the Internal Transcribed Spacer (ITS) gene on the rRNA of E. bieneusi via nested PCR. Results showed that the average infection rate of E. bieneusi in the 11 populations was 13.6%, and the highest infection rate was 56.5% in the population of Xinjiang Tianshan Zoo. Seven genotypes were identified including 2 known genotypes (D and CM1) and 5 novel genotypes (Mul1, Mul2, Mul3, Mul4 and Mul5). Phylogenetic analysis revealed that the novel genotypes Mul2, Mul3, Mul4 and Mul5 belonged to Group 1 showed the zoonotic potential．These findings extend the distribution of E. bieneusi genotypes and provide baseline data for controlling E. bieneusi infection.

The Tibetan chicken is a native Chinese breed that lives at high elevations and has adapted to the extreme environmental conditions of the Tibetan Plateau. However, its hypoxic adaptation at the gene expression level is unclear. Here, we sequenced nine lung transcriptomes of the Tibetan chicken at three developmental stages (5 and 42 weeks and 4.5 years). A total of 1.02 billion clean reads were obtained. We identified 16 012 mRNAs and 6898 lncRNAs. The expression of mRNA showed that nine samples were significantly divided into three clusters, with higher correlation and closer relationship between the 5 and 42 week groups. We identified 399 differentially expressed genes (DEGs) between the 5 and 42 week groups, 3532 DEGs between the 5 week and 4.5 year groups, and 3909 DEGs between the 42 week and 4.5 year groups. The up-regulated DEGs in the 5 week group, compared with 42 week and 4.5 year groups, were enriched in GO terms associated with growth and development, whereas the up-regulated DEGs in the 4.5 year group were mainly enriched in many metabolic-related categories. Moreover, the enrichment results with up-regulated DEGs in the 5 and/or 42 week groups, compared with the 4.5 year group, were associated with hypoxic adaptation, such as oxygen transport, oxygen binding and oxygen carrier activity, and calcium signaling pathway. In addition, we identified 978 high-correlation lncRNA and protein-coding gene pairs, and 524 significant neighboring protein-coding genes were also DEGs. Our results provide new insights into gene expression of lung tissue in Tibetan chickens during the aging process.

Inhibin A regulates follicular development, and its expression level is related to physiological activities, such as the recruitment, selection, and predominance during follicular development. Therefore, examining inhibin A and its regulatory effects on the reproductive performance of poultry is crucial. In this study, we measured the mRNA and protein abundances of INHA and INHBA in the chicken reproductive system and determined the hormone secretion and apoptosis of follicular granulosa cells (GCs) after being treated with inhibin A protein, and flow cytometry was performed to analyze GC apoptosis in INHA-specific small RNA interference (siRNA). We detected that INHA and INHBA were mainly expressed in chicken follicles. The highest INHA mRNA abundance was found in the fifth largest preovulatory follicle (F5) (P

The genetic structure and evolutionary relationship of Dulong chicken with other native Chinese species remained unclear. In this study, the mitochondrial control region was analyzed in total of 343 samples comprising 59 from Dulong chicken and 284 from 8 other Chinese local breeds revealed 51 mutation sites that defined 42 haplotypes. The maximum genetic variation was observed between the Shimian caoke and Pengxian yellow chickens. Phylogenetic analysis revealed that these local chickens mainly scatter in two southwestern clades. Dulong chickens have close relationship with other native chicken. Finding of this study suggests a single matrilineal lineage of indigenous Dulong chickens.

Four and a half LIM domain protein 1 (FHL1) belongs to the FHL protein family and is predominantly expressed in skeletal and cardiac muscle. FHL1 acts as a scaffold during sarcomere assembly and plays a vital role in muscle growth and development. Autophagy is key to skeletal muscle development and regeneration, with its dysfunction associated with a range of muscular pathologies and disorders. In this study, we constructed FHL1-silenced or FHL1-overexpressed myoblasts to investigate its role in autophagy during the differentiation of chicken myoblasts into myotubules. Our data showed that FHL1 contributes to myoblast differentiation as measured through MyoG, MyoD, Myh3, and Mb mRNA expression, MyoG and MyHC protein expression and the morphological characteristics of myoblasts. The results showed that FHL1 silencing inhibited the expression of ATG5 and ATG7, meanwhile, immunofluorescence and immunoprecipitation showed that FHL1 and LC3 interacted to regulate the correct formation of autophagosomes. FHL1 inhibition increased cleaved caspase-3 and PARP abundance and promoted myoblast apoptosis. Furthermore, FHL1 rescued skeletal muscle atrophy through regulating the expression of Atrogin-1 and MuRF1. Taken together, these data suggested that FHL1 regulates chicken myoblast differentiation through its interaction with LC3.

FOXO3, which encodes the transcription factor forkhead box O-3 (FoxO3), is a member of the FOXO subfamily of the forkhead box (FOX) family. FOXO3 can be negatively regulated by its phosphorylation by the PI3K/Akt signaling pathway and ultimately drives apoptosis when activated. In mammalian ovaries, the FOXO3 protein regulates atresia and follicle growth by promoting apoptosis of ovarian granulosa cells. Nonetheless, the specific effects of the FOXO3 protein on granulosa apoptosis of avian ovaries have not been elucidated. Therefore, we studiedFOXO3expression in follicles with different organization and at all hierarchical levels of chicken follicles. Via an immunofluorescence assay, the chicken follicular theca at all hierarchical levels were found to be strongly stained with an anti-FOXO3 antibody. In chicken primary ovarian granulosa cells, mRNA levels of proapoptotic factors BNIP3 and BCL2L11 decreased in the absence of FOXO3, and so did PARP-1 and cleaved caspase 3 protein levels. After treatment with a recombinant FOXO3 protein, PARP-1 and caspase 3 protein levels increased, along with mRNA levels of Bnip3 and BCL2L11 (significantly,pFOXO3was downregulated in chicken granulosa cells when different estradiol or FSH concentrations were applied. In conclusion,FOXO3is expressed in chicken reproductive tissues, including follicles and ovarian granulosa cells, and promotes apoptosis of chicken ovarian granulosa cells.

Skeletal muscle satellite cell growth and development is a complicated process driven by multiple genes. The PDZ and LIM domain 5 (PDLIM5) gene has been proven to function in C2C12 myoblast differentiation and is involved in the regulation of skeletal muscle development. The role of PDLIM5 in chicken skeletal muscle satellite cells, however, is unclear. In this study, in order to determine whether the PDLIM5 gene has a function in chicken skeletal muscle satellite cells, we examined the changes in proliferation and differentiation of chicken skeletal muscle satellite cells (SMSCs) after interfering and overexpressing PDLIM5 in cells. In addition, the molecular pathways of the PDLIM5 gene regulating SMSC proliferation and differentiation were analyzed by transcriptome sequencing. Our results show that PDLIM5 can promote the proliferation and differentiation of SMSCs, furthermore, through transcriptome sequencing, it can be found that the differential genes are enriched in the MAPK signaling pathway after knocking down PDLIM5. Finally, it was verified that PDLIM5 played an active role in the proliferation and differentiation of chicken SMSCs by activating the p38-MAPK signaling pathway. These results indicate that PDLIM5 may be involved in the growth and development of chicken skeletal muscle.

Chicken atrophic ovaries have decreased volume and are indicative of ovarian failure, presence of a tumor, or interrupted ovarian blood supply. Ovarian tumor is accompanied by an increase in follicular atresia, granulosa cell (GC) apoptosis, and autophagy. In a previous study, we found using high throughput sequencing that miR-204 is highly expressed in chicken atrophic ovaries. Thus, in the present study, we further investigated its function in GC apoptosis and autophagy. We found that overexpression of miR-204 reduced mRNA and protein levels of proliferation-related genes and increased apoptosis-related genes. Cell counting kit-8 (CCK-8), 5-ethynyl-2-deoxyuridine (EdU), and flow cytometry assays revealed that miR-204 inhibited GC proliferation and promoted apoptosis. Furthermore, we confirmed with reporter gene assays that Forkhead box K2 (FOXK2) was directly targeted by miR-204. FOXK2, as a downstream regulator of phosphoinositide 3-kinase (PI3K)/AKT/mammalian target of rapamycin (mTOR) signal pathways, promoted GC proliferation and inhibited apoptosis. Subsequently, we observed that miR-204 was involved in GC autophagy by targeting Transient Receptor Potential Melastatin 3 (TRPM3). The luciferase activities of the two binding sites of TRPM3 were decreased in response to treatment with a miR-204 mimic, and the autophagic flux was increased after miR-204 inhibition. However, overexpression of miR-204 had opposite results in autophagosomes and autolysosomes. miR-204 inhibits GC autophagy by suppressing the protein expression of TRPM3/AMP-activated protein kinase (AMPK)/ULK signaling pathway components. Inhibition of miR-204 enhanced autophagy by accumulating and degrading the protein levels of LC3-II (Microtubule Associated Protein Light Chain 3B) and p62 (Protein of 62 kDa), respectively, whereas miR-204 overexpression was associated with contrary results. Immunofluorescence staining showed that there was a significant reduction in the fluorescent intensity of LC3B, whereas p62 protein was increased after TRPM3 silencing. Collectively, our results indicate that miR-204 is highly expressed in chicken atrophic ovaries, which promotes GC apoptosis via repressing FOXK2 through the PI3K/AKT/mTOR pathway and inhibits autophagy by impeding the TRPM3/AMPK/ULK pathway.

Circular RNAs (circRNAs) are recognized as functional non-coding transcripts; however, emerging evidence has revealed that some synthetic circRNAs generate functional peptides or proteins. Additionally, the diverse biological functions of circRNAs include acting as miRNA-binding sponges, RNA-binding protein regulators, and protein translation templates. Previously, we found that circular RNA circFAM188B is a stable circular RNA and differentially expressed between broiler chickens and layers during embryonic skeletal muscle development. In this study, we found that circFAM188B exhibited a unique pattern of sharply decreased expression from embryonic day 10 (E10) to Day 35 (D35) after hatching. Our experimental results showed that circFAM188B promotes the proliferation, but inhibits the differentiation of chicken skeletal muscle satellite cells (SMSCs). Bioinformatic analysis revealed circFAM188B contain an opening reading frame (ORF) which translate into circFAM188B-103aa, internal ribosome entry site (IRES) analysis further confirmed the coding potential of circFAM188B. In addition, western blot assay detected a flag tagged circFAM188B-103aa, and several peptides of circFAM188B-103aa were detected by LC-MS/MS analysis. We further verified that the role of circFAM188B-103aa in chicken myogenesis is consistent with that of its parent transcript circFAM188B, which facilitates proliferation, but represses differentiation of chicken SMSC. Taken together, these results suggested that a novel protein circFAM188B-103aa encoded by circFAM188B that promotes the proliferation but inhibits the differentiation of chicken SMSCs.

In order to breed new birds with strong disease resistance, it is necessary to first understand the mechanism of avian antiviral response. Interferon regulatory factor 7 (IRF7) is not only a member of type I interferons (IFNs) regulatory factor (IRFs) family, but also a major regulator of the IFN response in mammals. However, whether IRF7 is involved in the host innate immune response remains unclear in poultry, due to the absence of IRF3. Here, we first observed by HE stains that with the increase of the time of ALV-J challenge, the thymus was obviously loose and swollen, the arrangement of liver cell was disordered, and the bursa of fabricius formed vacuolated. Real-time PCR detection showed that the expression level of IRF7 gene and related immune genes in ALV-J group was significantly higher than that in control group (P

INHA, the gene encoding the inhibin alpha subunit, was involved in folliculogenesis in mammals, but no study was reported for its working pathway in birds. Here we hypothesize that gene polymorphism in INHA 3′UTR might influence miRNAs binding efficiency and further affect the function of this gene. Thus, we investigated the association between the 3′UTR single-nucleotide polymorphisms (SNPs) in INHA and the laying performance in chickens and further explore their possible molecular cascades in granulosa cells (GC). Five SNPs were detected in Tianfu green-shell layers and g. 22,178,975 G > A was significantly associated with total egg numbers at the age of 300 days (EN, n = 286). Birds carrying the AA genotype laid more EN than those with GG (P

As bioinformatic approaches have been developed, it has been demonstrated that microRNAs (miRNAs) are involved in the formation of muscles and play important roles in regulation of muscle cell proliferation and differentiation. Previously, it has been demonstrated that miR-148a-3p is one of the most abundant miRNAs in chicken skeletal muscle. Here, we build on that work and demonstrate that miR-148a-3p is important in the control of differentiation of chicken skeletal muscle satellite cells (SMSCs). Elevated expression of miR-148a-3p significantly promoted differentiation and inhibited apoptosis of SMSCs but did not affect proliferation. Furthermore, it was observed that the mesenchyme homeobox 2 (Meox2) is a target gene of miR-148a-3p and that miR-148a-3p can down-regulate expression of Meox2, which promote differentiation of SMSCs and suppress apoptosis. Furthermore, miR-148a-3p overexpression encouraged activation of the PI3K/AKT signaling pathway, which could be recovered by overexpression of Meox2. Overall, these findings suggest that microRNA-148a-3p is a potent promoter of myogenesis via direct targeting of Meox2 and increase of the PI3K/AKT signaling pathway in chicken SMSCs.

The hypothalamic-pituitary-ovarian (HPO) axis regulates the breeding process cycle of laying hens. However, the key regulatory genes of the HPO axis and pathways that drive chicken egg laying performance remain elusive. A total of 856 Chinese Luhua chicken was raised and the highest two hundred and the lowest two hundred chicken egg production were considered as high egg production (HEP) and low egg production (LEP) according to the total egg number at 300 days of age, respectively. RNA-seq sequencing (RNA-Seq) was conducted to explore the chicken transcriptome from the hypothalamus, pituitary gland and ovary tissue of 6 Chinese Luhua chicken with 3 high and low-rate egg production. In total, 76.09 Gb RNA-seq sequences were generated from 15 libraries with an average of 5.07 Gb for each library. Further analysis showed that 414, 356 and 10 differentially expressed genes (DEGs) were identified in pituitary gland, ovary and hypothalamus between HEP and LEP chickens, respectively. In pituitary gland, DEGs were involve in regulation of cellular glucose homeostasis, Ras protein signal transduction, negative regulation of hormone secretion. In Ovary DEGs were mainly involved in embryonic organ development, regulation of canonical Wnt signaling, response to peptide hormone. Our study identified DEGs that regulate mTOR signaling pathway, Jak-STAT signaling pathway, Tryptophan metabolism and PI3K-Akt signaling pathways at HPO-axis in laying hens. These important data contribute to improve our understanding of reproductive biology of chicken and isolating effective molecular markers that can be used for genetic selection in Chinese domestic Luhua chicken.

Noncoding RNAs, especially microRNAs (miRNAs), have been reported to play important roles during skeletal muscle development and regeneration. Our previous sequencing data revealed that miR-99a-5p is one of the most abundant miRNAs in chicken breast muscle. The purpose of this study was to reveal the regulatory mechanism of miR-99a-5p in the proliferation and differentiation of chicken skeletal muscle satellite cells (SMSCs). Through the investigation of cell proliferation activity, cell cycle progression, and 5-ethynyl-29-deoxyuridine (EdU) assay, we found that miR-99a-5p can significantly promote the proliferation of SMSCs. Moreover, we found that miR-99a-5p can inhibit myotube formation by decreasing the expression of muscle cell differentiation marker genes. After miR-99a-5p target gene scanning, we confirmed that miR-99a-5p directly targets the 3&prime, untranslated region (UTR) of myotubularin-related protein 3 (MTMR3) and regulates its expression level during chicken SMSC proliferation and differentiation. We also explored the role of MTMR3 in muscle development and found that its knockdown significantly facilitates the proliferation but represses the differentiation of SMSCs, which is opposite to the effects of miR-99a-5p. Overall, we demonstrated that miR-99a-5p regulates the proliferation and differentiation of SMSCs by targeting MTMR3.

MicroRNAs (miRNAs) are evolutionarily conserved, small noncoding RNAs that play critical post-transcriptional regulatory roles in skeletal muscle development. Chicken is an optimal model to study skeletal muscle formation because its developmental anatomy is similar to that of mammals. In this study, we identified potential miRNAs in the breast muscle of broilers and layers at embryonic day 10 (E10), E13, E16, and E19. We detected 1836 miRNAs, 233 of which were differentially expressed between broilers and layers. In particular, miRNA-200a-3p was significantly more highly expressed in broilers than layers at three time points. In vitro experiments showed that miR-200a-3p accelerated differentiation and proliferation of chicken skeletal muscle satellite cells (SMSCs) and inhibited SMSCs apoptosis. The transforming growth factor 2 (TGF-&beta, 2) was identified as a target gene of miR-200a-3p, and which turned out to inhibit differentiation and proliferation, and promote apoptosis of SMSCs. Exogenous TGF-&beta, 2 increased the abundances of phosphorylated SMAD2 and SMAD3 proteins, and a miR-200a-3p mimic weakened this effect. The TGF&beta, 2 inhibitor treatment reduced the promotional and inhibitory effects of miR-200a-3p on SMSC differentiation and apoptosis, respectively. Our results indicate that miRNAs are abundantly expressed during embryonic skeletal muscle development, and that miR-200a-3p promotes SMSC development by targeting TGF-&beta, 2 and regulating the TGF&beta, 2/SMAD signaling pathway.

Immunoglobulin superfamily containing leucine-rich repeat (Islr) contains an Ig-like domain, an LRR motif, and a transmembrane domain and is highly expressed in various chicken tissues. Although Islr has known roles in muscle regeneration, its role in the regulation of muscle atrophy has not been studied. In this study, we constructed Islr-silenced or Islr-overexpressed myoblasts to investigate its role during the differentiation of myoblasts into myotubes. The results showed that Islr was highly expressed in chicken skeletal muscle tissue and regulated myoblast differentiation, but not proliferation. Islr regulated the expression of atrophy-related genes including atrogin-1 and MuRF-1, and could rescue dexamethasone-induced atrophy in myoblasts and myotubes. Western blot analysis indicated that Islr participates in myoblast atrophy through IGF/PI3K/AKT-FOXO signaling. Meanwhile, the expression of caspase-8 and caspase-9 increased in Islr-silenced groups, indicating its role in cell viability. Taken together, these data suggested that Islr plays an important role in myoblasts differentiation, and which can alleviate skeletal muscle atrophy and prevents muscle cell apoptosis via IGF/PI3K/AKT-FOXO signaling pathway.

MicroRNA (miRNA) is known to be an important regulator of muscle growth and development. The regulation of microRNA on the skeletal muscle phenotype of animals is mainly achieved by regulating the proliferation and differentiation of myoblasts. In this study, we sequenced a total of 60 samples from 15 developing stages of the pectoral muscle and five other tissues at 300 days of Tibetan chicken. We characterized the expression patterns of miRNAs across muscle developmental stages, and found that the chicken growth and development stage was divided into early-embryonic and late-embryonic as well as postnatal stages. We identified 81 and 21 DE-miRNAs by comparing the miRNA profiles of pectoral muscle of three broad periods and different tissues, respectively; and 271 miRNAs showed time-course patterns. Their potential targets were predicted and used for functional enrichment to understand their regulatory functions. Significantly, GgmiRNA-454 is a time-dependent and tissue-differential expression miRNA. In order to elucidate the role of gga-miRNA-454 in the differentiation of myoblasts, we cultured chicken myoblasts in vitro. The results show that although gga-miRNA-454-3p initiates increase and thereafter decrease during the chicken myoblasts differentiation, it had no effect on primary myoblasts proliferation. Furthermore, we confirm that gga-miRNA-454 inhibits myoblast differentiation by targeting the myotube-associated protein SBF2.

We characterized 26 wild fruit flies comparative population genomics from six different altitude and latitude locations by whole genome resequencing. Genetic diversity was relatively higher in Ganzi and Chongqing populations. We also found 13 genes showing selection signature between different altitude flies and variants related to hypoxia and temperature stimulus, were preferentially selected during the flies evolution. One of the most striking selective sweeps found in all high altitude flies occurred in the region harboring Hsp70Aa and Hsp70Ab on chromosome 3R. Interestingly, these two genes are involved in GO terms including response to hypoxia, unfolded protein, temperature stimulus, heat, oxygen levels. Mutation in HPH gene, a candidate gene in the hypoxia inducible factor pathway, might contributes to hypoxic high-altitude adaptation. Intriguingly, some of the selected genes, primarily utilized in humans, were involved in the response to hypoxia, which could imply a conserved molecular mechanisms underlying high-altitude adaptation between insects and humans.

T-2 toxin is a trichothecene mycotoxin produced by fungi which are known to contaminate cereals, especially in wheat and corn. T-2 toxin is known to cause a range of toxic effects in humans and animals, including immunosuppression and carcinogenesis. Although the effects of T-2 toxin on condition of chickens’ spleens have been reported, there has been no systematic study of damage to the spleen of broiler chickens exposed to T-2 toxin. The purpose of the present study was to assess the effects of T-2 toxin on pathology, rates of apoptosis, oxidative stress, and T-lymphocyte subsets in the spleen of broiler chickens. One hundred and twenty male broiler chickens were randomly assigned to one of four groups (30 birds per group), fed 0 mg/kg (control), 0.5 mg/kg, 1 mg/kg, or 2 mg/kg T-2 toxin, respectively. After 21 days, chickens exposed to T-2 toxin demonstrated decreased relative weight and size of the spleen, increased percentage of apoptotic splenocytes, and evident lesions. Concentrations of reactive oxygen species and MDA content increased in splenocytes during T-2 toxin treatments, whereas activities of SOD, CAT, and GSH-PX decreased. The ratio of CD4+/CD8+ T cells also decreased as the dose of T-2 toxin increased. Overall, these results suggest that T-2 toxin causes oxidative stress, leading to increased rates of splenocyte apoptosis and might impair the splenic immune function of broiler chickens.

Animal gut microbiota begins to colonize after birth and is functionally indispensable for maintaining the health of the host. It has been reported that gender and age influence the composition of the intestinal microbiome. However, the effects of gender and age on the intestinal microorganism of forest musk deer (FMD) remain unclear. The aim of this study was to establish the relationship between the structure and composition of fecal microbiota of male and female forest musk deer with age. Here, Illumina Miseq 300PE sequencing platform targeting 16S rRNA V3–V4 hypervariable region applied to define the fecal microbiota of male and female FMD with two age groups, juvenile (age 1–2 years) and adult (age 4–10 years). Alpha diversity index did not show significant difference in bacterial diversity between the males and females or among age groups. The intestinal microbiota of FMD was dominated by three phyla, the Firmicutes, Proteobacteria and Bacteroidetes regardless of gender and different ages. Higher proportions of Proteobacteria were found in adult male and juvenile female individuals. The composition of Bacteroidetes was stable with the gender and age of FMD. Interestingly, the relative abundance of genera Clostridiales and Bacteroidales were higher in the juvenile FMD. Conversely, proportions of Pseudomonas and Lachnospiraceae were abundant in the adult FMD. Higher proportions of Ruminococcaceae, Dore, and 5-7N15 were found in the juvenile male groups. They may reflect the different immune resistance of male and female individuals at different stages of development. This study explored the fecal microbiota composition of forest musk deer in relation to gender and age, which may provide an effective strategy for developing intestinal microecological preparations and potential musk deer breeding.

Broodiness in laying hens results in atrophy of the ovary and consequently decreases productivity. However, the regulatory mechanisms that drive ovary development remain elusive. Thus, we collected atrophic ovaries (AO) from 380-day-old broody chickens (BC) and normal ovaries (NO) from even-aged egg-laying hens (EH) for RNA sequencing. We identified 3,480 protein-coding transcripts that were differentially expressed (DE), including 1,719 that were down-regulated and 1,761 that were up-regulated in AO. There were 959 lncRNA transcripts that were DE, including 56 that were down-regulated and 903 that were up-regulated. Among the116 miRNAs that were DE, 79 were down-regulated and 37 were up-regulated in AO. Numerous DE protein-coding transcripts and target genes for miRNAs/lncRNAs were significantly enriched in reproductive processes, cell proliferation, and apoptosis pathways. A miRNA-intersection gene-pathway network was constructed by considering target relationships and correlation of the expression levels between ovary development-related genes and miRNAs. We also constructed a competing endogenous RNA (ceRNA) network by integrating competing relationships between protein-coding genes and lncRNA transcripts, and identified several lncRNA transcripts predicted to regulate the CASP6, CYP1B1, GADD45, MMP2, and SMAS2 genes. In conclusion, we discovered protein-coding genes, miRNAs, and lncRNA transcripts that are candidate regulators of ovary development in broody chickens.

The goral (Naemorhedus caudatus), a species of the Caprinae, is classified as globally vulnerable on the IUCN Red List of Threatened Species because of its sharply declining population size, which is attributable to habitat destruction. In this study, we determined and described the complete mitochondrial genome of N. c. griseus, a goral subspecies from the western Sichuan of China; the mitochondrial genome of this goral is 16,557 bp in length (GenBank accession No: MG865962) and contains 13 protein-coding genes (PCGs), two rRNA genes, 23 tRNA genes, and one control region (D-loop). Ten PCGs start with ATG, but ND2, ND3, and ND5 use ATA as the start codon. The overall base composition of the mitochondrial DNA is 34.6% A, 28.6% T, 22.4% C, and 14.4% G. The G+C content is 36.8%. Phylogenetic analysis revealed that N. c. griseus clustered into a clade with other Naemorhedus species with high bootstrap support. This work provides additional molecular information for studying the conservation genetics and evolutionary relationships of this endangered species.

Knowledge about the impact of different geographical environments on rhesus macaque gut microbiota is limited. In this study, we compared the characteristics of gut microbiota in six different Chinese rhesus macaque populations, including Hainan, Nanning, Guizhou, Xichang, Jianchuan and Tibet. Through the composition analysis of operational taxonomic units (OTUs), we found that there were significant differences in the abundance of core overlapping OTUs in the six Chinese groups. Specifically, the Tibet population exhibited the highest gut microbial diversity and the most unique OTUs. Statistically significant differences in the composition of gut microbiota among the six groups at phylum and family level were evident. Specifically, Tibet had higher abundances of Firmicutes and lower abundances of Bacteroidetes than the other geographical groups, and the higher abundance of Firmicutes in the Tibetan group was mainly caused by a significant increase in the family Ruminococcaceae and Christensenellaceae. Phylogenetic investigation of communities by reconstruction of unobserved state analysis showed that the enrichment ratio for environmental information processing and organismal systems was the highest in the Tibet population. Additionally, our results suggested that in the adaptation process of rhesus macaques to different geographical environments, the abundance of the core common flora of the intestinal microbes had undergone varying degree of change and produced new and unique flora, both of which helped to reshape the gut microbiota of rhesus macaques. In particular, this change was more obvious for animals in the high-altitude environments.

The white-lipped deer (Cervus albirostris) is a rare and endangered species endemic to the Qinghai–Tibet Plateau in China. In this study, we determined the complete mitochondrial genome sequence of the white-lipped deer. The total length of the mitogenome is 16473 bp. It has the typical structure of a vertebrate mitochondrial genome including 13 protein-coding genes(PCGs), 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes and a control region (D-loop). tRNASer lacks the stem of the dihydrouridine (DHU) arm and fails to form a typical cloverleaf secondary structure. A phylogenetic tree based on the 13 PCG sequences of nine species closely related to white-lipped deer shows the classification status of the species within the family Cervidae. These results provide new molecular biology information to further understand the diversity of the white-lipped deer, which will help to protect this population.

RNF11 (RING finger protein 11) has been indicated to be related to numerous immunological diseases. As a neo-regulator in innate immunity, RNF11 was proved to inhibit antiviral signaling and interferon-β (IFN-β) production in human cell line. Chicken antiviral signaling pathway is different with other species. However, there is little study about the role of RNF11 in chicken antiviral signaling pathway and inflammatory response. This study was designed to identify the function of RNF11 in chicken antiviral signaling pathway and inflammatory response. We overexpressed and interfered RNF11 in chicken embryo fibroblasts (CEFs) and then used quantitative reverse transcription polymerase chain reaction (qRT-PCR) to investigate the poly (I:C)-induced production of IFN-β and inflammatory response factors IL2 and IL10. The results indicated that RNF11 extremely significantly decreased the production of IFN-β and highly significantly increased the messenger RNA (mRNA) levels of IL2 and IL10, but depletion of RNF11 obtained inverse results. Thus, we concluded that RNF11 not only can inhibit antiviral signaling but also plays an important role in inflammatory response in chicken.

Background Species living at high altitude are subject to strong selective pressures due to inhospitable environments (e.g., hypoxia, low temperature, high solar radiation, and lack of biological production), making these species valuable models for comparative analyses of local adaptation. Studies that have examined high-altitude adaptation have identified a vast array of rapidly evolving genes that characterize the dramatic phenotypic changes in high-altitude animals. However, how high-altitude environment shapes gene expression programs remains largely unknown. Findings We generated a total of 910 Gb of high-quality RNA-seq data for 180 samples derived from 6 tissues of 5 agriculturally important high-altitude vertebrates (Tibetan chicken, Tibetan pig, Tibetan sheep, Tibetan goat, and yak) and their cross-fertile relatives living in geographically neighboring low-altitude regions. Of these, ∼75% reads could be aligned to their respective reference genomes, and on average ∼60% of annotated protein coding genes in each organism showed FPKM expression values greater than 0.5. We observed a general concordance in topological relationships between the nucleotide alignments and gene expression–based trees. Tissue and species accounted for markedly more variance than altitude based on either the expression or the alternative splicing patterns. Cross-species clustering analyses showed a tissue-dominated pattern of gene expression and a species-dominated pattern for alternative splicing. We also identified numerous differentially expressed genes that could potentially be involved in phenotypic divergence shaped by high-altitude adaptation. Conclusions These data serve as a valuable resource for examining the convergence and divergence of gene expression changes between species as they adapt or acclimatize to high-altitude environments.

Skeletal muscle is the most abundant tissue in the human and animal body, loss of its function can lead to muscle aging and various myogenic diseases. The skeletal muscle development is a complex and tightly regulated process, which is driven by a variety of many factors, signaling pathways and regulatory mechanisms. Plectin (Plec), a cytolinker protein, is ubiquitously expressed in various tissues such as skin, muscle, plasma membrane, and most types of cells. Although known isoforms of Plec is well-characterized in muscle dystrophy, very little is known on the function of Plec in the skeletal muscle development. Here, we found that Plec plays a vital role in promoting C2C12 myoblasts differentiation and proliferation, but inhibits their apoptosis. Also, Plec regulates the expression of atrophy-related genes (atrogin-1 and muRF-1) to rescue muscle atrophy. Furthermore, we have demonstrated that Plec binds to Dishevelled-2 (Dvl-2) and forms a protein complex, which is then activate the canonical Wnt signaling. We also observed that Plec resists ubiquitination by stabilizing Dvl-2 and reduces the level of LC3-labeled Dvl-2 and antagonizes the autophagy system. In conclusion, our findings suggest that Plec regulates canonical Wnt signaling mediated skeletal development by stabilizing Dvl-2 and downregulating the cellular autophagic degradation system.

The Sichuan deer (Cervus elaphus macneilli) is a rare animal native to western Sichuan and southeastern Tibet, China. In this study, the entire Sichuan deer mitochondrial genome was sequenced; it was 16350 bp long (GenBank accession No: KX449334) and contained 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes, and one control region (D-loop). Ten PCGs started with ATG, but ND2, ND3, and ND5 used ATA as the start codon. PCG homology between C.e.macneilli and four other red deer subspecies ranged from 93.9 to 96%, which indicates a clear difference among subspecies. Phylogenetic analysis revealed that the Sichuan deer clustered into a clade with other red deer subspecies with high bootstrap support. This work provides additional molecular information for studying red deer conservation genetics and evolutionary relationships.

Forest musk deer (FMD) is an endangered species, and its population size has dropped dramatically. In this study, we determined the complete mitochondrial genome of a FMD. The genome was 16,353 bp long, and contained 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, two ribosomal RNA genes, and one control region. Most PCGs were distributed on the heavy strand except ND6 and eight tRNA genes, which were encoded on the light strand. The nucleotide composition was asymmetric, with an overall A+T content of 62.1%. Comparison with one previously described mitochondrial genome revealed 94.4% sequence homology and 1128 nucleotide mutation sites, which represents a substantial difference. A phylogenetic tree based on Cytb genes of eight closely related musk deer species showed that our sample clustered with two FMD subspecies from Yunnan, China. These results provide novel molecular information that can potentially be used for genetic diversity conservation of this species.

Tibetan chicken (TC) is a unique chicken breed that has inhabited the Tibetan Plateau for thousands of years and has adapted to the harsh environment of an altitude of over greater than 4000 m. The changes in blood characteristics that occur in high-altitude rearing of TCs are unknown. To obtain a deeper insight into the blood characteristics of chickens in low and high altitudes, we focussed on four chicken populations to compare the variations among them. In the present study, we found significant differences in white blood cell (WBC) count, red blood cell (RBC) count and haematocrit (HCT) values between TCs and other chicken populations. When TCs were reared at low altitude, there was a significant (p

GPR133 (G protein-coupled receptor 133) plays significant roles in various physiological processes. Alternatively splicing (AS) variants of GPR133 in many species have been predicted in multiple databases, but there is no available information about the AS events of chicken GPR133 (cGPR133). In the present study, two chicken GPR133 variants, cGPR133-va and cGPR133-vb, were identified by a combination of reverse transcription PCR (RT-PCR) and rapid amplification of cDNA 5'-ends (5' RACE). Sequence analysis shows that cGPR133-va and cGPR133-vb are resulting from different AS modules and their sequences are predicted to encode two distinct putative proteins, respectively. Quantitative real-time PCR (qRT-PCR) analysis revealed that cGPR133-va and cGPR133-vb are widely expressed in different tissues, while exhibiting specific expression profile. Altogether, our results first demonstrate the existence of novel cGPR133 variants and illustrate its transcriptional diversity and their widespread distribution, which provides a foundation for the further research of GPR133.

Background: Since the domestication of the red jungle fowls (Gallus gallus; dating back to ∼10 000 B.P.) in Asia, domestic chickens (Gallus gallus domesticus) have been subjected to the combined effects of natural selection and human-driven artificial selection; this has resulted in marked phenotypic diversity in a number of traits, including behavior, body composition, egg production, and skin color. Population genomic variations through diversifying selection have not been fully investigated. Findings: The whole genomes of 78 domestic chickens were sequenced to an average of 18-fold coverage for each bird. By combining this data with publicly available genomes of five wild red jungle fowls and eight Xishuangbanna game fowls, we conducted a comprehensive comparative genomics analysis of 91 chickens from 17 populations. After aligning ∼21.30 gigabases (Gb) of high-quality data from each individual to the reference chicken genome, we identified ∼6.44 million (M) single nucleotide polymorphisms (SNPs) for each population. These SNPs included 1.10 M novel SNPs in 17 populations that were absent in the current chicken dbSNP (Build 145) entries. Conclusions: The current data is important for population genetics and further studies in chickens and will serve as a valuable resource for investigating diversifying selection and candidate genes for selective breeding in chickens.

Trachypithecus poliocephalus : (Primates: Cercopithecidae), classified as critically endangered (IUCN) and in China was listed National First Class Protected Animal. In this study, the complete mitochondrial genome sequence of T. poliocephalus has been reported for the first time in China. The complete mitochondrial genome of T. poliocephalus was sequenced and it contained 13 protein-coding genes, 2 rRNA genes (12S rRNA and 16S rRNA), 22 tRNA genes, and one control region (D-loop).

We newly obtained the near complete mitochondrial genome of Oreolalax schmidti (18,481bp) by using polymerase chain reaction (PCR) in this study. It includes 13 protein-coding genes (PCG), two ribosomal RNA (rRNA) genes, and 23 transfer RNA (tRNA) genes. The phylogenetic tree indicates that the Oreolalax schmidti is closely related to the O. lichuanensis. This report will help the further studies of Oreolalax species classification and source protection.

The complete mitochondrial genome of the Oreolalax major (17,786 bp long) was obtained in this study. It includes 13 protein-coding genes (PCG), two ribosomal RNA (rRNA) genes, and 23 transfer RNA (tRNA) genes (GenBank accession number MN803320). The phylogenetic tree indicates that the O. major is closely related to the O. xiangchengensis.

Ferret-badger (Melogale moschata) is a widely distributed and rarely studied species. Therefore, we sequenced its complete mitochondrial genome, which is circular in shape with 16,511 bp length and have a structure typical to mammals. There is 98.67% similarity when compared with the reported sequence (HM106328). It contained 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes, and 1 control regions (D-loop). This study can provide novel molecular information and genetic diversity of ferret-badger.

Red-and-white giant flying squirrel (Petaurista alborufus) is a widely distributed but rarely studied species. Here, we sequenced the complete mitochondrial genome about it, which is a circular genome of 16,511 bp and have a typical structure with mammals. The base composition shows higher composition of A\T. There are 61 variations and eight base deletions between two sequences when compared with the reported sequence (JQ743657). In this study, we obtain a new mitochondrial genome sequence of red-and-white giant flying squirrel, which provide novel molecular information that could prove fundamental to our understanding of diversity research and phylogenetic structure.

The complete mitochondrial of Cynops cyanurus cyanurus was determined by using Sanger sequencing and the mitogenome size of C. c. cyanurus was 16,465 bp. It contained 13 protein-coding genes, 22 tRNAs, 2 rRNAs, 1 control region (D-loop) and a non-coding region (NC). In addition, the phylogenetic tree shows that C. c. cyanurus has the closest relationship with C. c. chuxiongensis. Hoping this study can improve our understanding of the species evolutionary relationship of Cynops.

Due to the multiple causes, the population of roe deer has declined significantly. In this study, we analyzed the complete mitogenome of Capreolus pygargus bedfordi, whose genome was 16,357 bp long. There were 13 protein-coding genes (PCG), two ribosomal RNA genes (12S rRNA and 16S rRNA), 22 transfer RNA genes, and one control region. Nine PCGs started with ATG, while NAD2, NAD3, and NAD5 genes commenced with ATA, and ND4L began with GTG. ND6 and eight tRNA genes were encoded on the L-strand. These results provide newer molecular information, which contribute to its molecular and phylogenetic studies, and genetic diversity conservation.

// Binlong Chen 1, * , Diyan Li 1, * , Miao Li 1 , Sichen Li 1 , Kenan Peng 1 , Xian Shi 1 , Lanyun Zhou 1 , Pu Zhang 1 , Zhongxian Xu 1 , Huadong Yin 1 , Yan Wang 1 , Xiaoling Zhao 1 , Qing Zhu 1 1 Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, Sichuan Agricultural University, Chengdu, P.R.China, 611130 * These authors contribute equally to this work Correspondence to: Qing Zhu, email: zhuqing@sicau.edu.cn Keywords: aflatoxin B 2 , mitochondria-mediated apoptosis, PI3K/Akt/mTOR-mediated autophagy, hepatocytes, broiler Received: May 18, 2016 Accepted: October 26, 2016 Published: November 15, 2016 ABSTRACT Aflatoxins have been shown to induce hepatotoxicity in animal models, but the effects of aflatoxin B 2 (AFB 2 ) on broiler hepatocytes is unclear. This study aimed to investigate the effects of AFB 2 on apoptosis and autophagy to provide an experimental basis for understanding the mechanism of aflatoxin-induced hepatotoxicity. One hundred-twenty Cobb500 broilers were allocated to four groups and exposed to 0 mg/kg, 0.2 mg/kg, 0.4 mg/kg, and 0.8 mg/kg of AFB2 per day for 21 d. AFB 2 exerted potent proapoptotic and proautophagic effects on hepatocytes, with increased numbers of apoptotic and autophagic hepatocytes. Poly ADP-ribose polymerase (PARP) was cleaved and caspase-3 was activated in experimental groups, showing that the apoptosis of hepatocytes was triggered by AFB 2 . Increased levels of the autophagy factors Beclin-1 and LC3-II/LC3-I, as well as down-regulation of p62, a marker of autophagic flux, provided additional evidence for AFB 2 -triggered autophagy. AFB2 induced mitochondria-mediated apoptosis via the production of reactive oxygen species (ROS) and promotion of the translocation of Bax and cytochrome c (cyt c) between mitochondria and the cytosol, triggering the formation of apoptosomes. AFB2 also inhibited the phosphoinositide 3-kinase/Akt/mammalian target of rapamycin (PI3K/Akt/mTOR) pathway by activating PI3K, Akt, and mTOR and inhibiting their phosphorylation, contributing to the proautophagic activity of AFB 2 . These findings provide new insights into the mechanisms involved in AFB 2 -induced hepatotoxicity in broilers.

In this study, the complete mitochondrial genome sequence of Prionodon pardicolor pardicolor has been reported for the first time in China. The total length of the mitogenome was 16,689 bp. It contained the typical structure, including 13 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes and 1 control region (D-loop). tRNASer lost the stem of DHC loop and failed to form a typical cloverleaf secondary structure. Mitochondrial control region of P. p. pardicolor was 1242 bp long, laying between the tRNAPro and tRNAPhe genes. Based on the 13 PCG sequences, homology and single-nucleotide mutation loci were performed between P. p. pardicolor and P. pardicolor. As a result, the homology between two species was 99.6% on average, and Nucleotide mutation loci was 41. Also, in view of their mitochondrial control region comparison, P. p. pardicolor drop 28 bp nucleotides.

As the rate-limiting enzyme of the mitochondrial respiratory chain, cytochrome c oxidase (COX) plays a crucial role in biological metabolism. "Living fossil" giant panda (Ailuropoda melanoleuca) is well-known for its special bamboo diet. In an effort to explore functional variation of COX1 in the energy metabolism behind giant panda's low-energy bamboo diet, we looked at genetic variation of COX1 gene in giant panda, and tested for its selection effect. In 1545 base pairs of the gene from 15 samples, 9 positions were variable and 1 mutation leaded to an amino acid sequence change. COX1 gene produces six haplotypes, nucleotide (pi), haplotype diversity (Hd). In addition, the average number of nucleotide differences (k) is 0.001629±0.001036, 0.8083±0.0694 and 2.517, respectively. Also, dN/dS ratio is significantly below 1. These results indicated that giant panda had a low population genetic diversity, and an obvious purifying selection of the COX1 gene which reduces synthesis of ATP determines giant panda's low-energy bamboo diet. Phylogenetic trees based on the COX1 gene were constructed to demonstrate that giant panda is the sister group of other Ursidae.

Cynops cyanurus chuxiongensis is distributed in Yunnan plateau. In this paper, we sequenced and determined the complete mitochondrial genome of C. c. chuxiongensis. The assembled mitogenome is 16,4...

In this study, we sequenced the complete mitochondrial genome of Argali sheep (Ovis ammon). The total length was 16,612 bp, which contained 13 protein-coding genes (PCGs), 22 tRNA genes, 2 rRNA genes, and a control region (D-loop). Eight tRNA genes with one PCG (ND6) encoded on the L-strand, others were encoded on the H-strand. The neighbor-joining analysis shows that Argali sheep has a close relationship with the same genus species of Ovis aries and Ovis orientalis. This study has provided new data for the phylogeny of Argali sheep.

Bitterness is an important taste sensation for chickens, which provides useful sensory information for acquisition and selection of diet, and warns them against ingestion of potentially harmful and noxious substances in nature. Bitter taste receptors (T2Rs) mediate the recognition of bitter compounds belonging to a family of proteins known as G-protein coupled receptors. The aim of this study was to identify and evaluate the expression ofT2R7in chicken tongue tissue and constructcT2R7-1andcT2R7-2-expressing HEK-293T cells to access the expression ofPLCβ2andITPR3after exposure with different concentrations of the bitter compounds. Using real-time PCR, we show that the relative expression level ofT2R7mRNA in 5, 1, 0.1, and 10−3 mM of camphor and erythromycin solutions and 5 mM of chlorpheniramine maleate solutions was significantly higher than that in 50 mM KCL solutions. We confirmed that the bitter taste receptor T2R7 and downstream signaling effectors are sensitive to different concentrations of bitter compounds. Moreover, T2R7-1 (corresponding to the unique haplotype of the Tibetan chicken) had higher sensitivity to bitter compounds compared with that of T2R7-2 (corresponding to the unique haplotype of the Jiuyuan black-chicken). These results provide great significance of taste response on dietary intake to improve chicken feeding efficiency in poultry production and have certain reference value for future taste research in other bird species.

Muscle LIM protein (MLP/CSRP3/CRP3) is a microtubule-associated protein preferentially expressed in cardiac and skeletal muscle and has a central role during muscle development and for architectural maintenance of muscle cells. LIM-domain proteins act as both modulators and downstream targets of TGF-β signaling, which is well documented to negatively regulate differentiation of myogenic precursor cells or myoblasts. Herein, we determined whether CSRP3 regulates chicken satellite cell proliferation and differentiation in vitro, and examined its mechanism of action by focusing on the TGF-β signaling pathway. Interference of CSRP3 mRNA expression had no effect on the proliferation of satellite cells, but significantly inhibited satellite cell differentiation into myotubes at 24, 48, and 72 h after initiation of differentiation. However, CSRP3 overexpression did not affect the proliferation or differentiation of satellite cells. Moreover, knockdown of CSRP3 caused up-regulation of TGF-β and Smad3 mRNA and protein levels. The phosphorylation of Smad3 in CSRP3-knockdown cells was greater than that in wild-type cells at 24, 48, and 72 h after initiation of differentiation. Collectively, knockdown of CSRP3 suppressed chicken satellite cell differentiation by regulating Smad3 phosphorylation in the TGF-β signaling pathway. Our results indicate that CSRP3 might play an important role in promoting satellite cell differentiation in chicken.

The Chinese forest musk deer (Moschus berezovskii) is an economically important species distributed throughout southwest China and northern Vietnam. Occurrence and development of disease are aggravated by inbreeding and genetic diversity declines in captive musk deer populations. Deep transcriptomics investigation may provide a promising way to improve genetic health of captive and wild FMD population. MicroRNAs (miRNAs), which regulate gene expression by targeting and suppressing of mRNAs, play an important role in physiology and organism development control. In this study, RNA-seq technology was adopted to characterize the miRNA transcriptome signature among six tissues (heart, liver, spleen, lung, kidney, and muscle) in Chinese forest musk deer at two years of age. Deep sequencing generated a total of 103,261,451 (~87.87%) good quality small RNA reads; of them 6,622,520 were unique across all six tissues. A total of 2890 miRNAs were identified, among them 1129 were found to be expressed in all tissues. Moreover, coexpression of 20 miRNAs (>2000RPM) in all six tissues and top five highly expressed miRNAs in each tissue implied the crucial and particular function of them in FMD physiological processes. Our findings of forest musk deer miRNAs supplement the database of transcriptome information for this species and conduce to our understanding of forest musk deer biology.

Inhibin α (INHA) is a candidate gene controlling ovulation in poultry. As the functional center of inhibin, INHA is a molecular marker associated with egg-laying performance. The objective of the current study was to analyze the expression differences of INHA in reproductive system and single nucleotide polymorphisms (SNPs) associations with reproductive traits in chickens. A total of 260 LuHua chickens (barred-feather chicken) were adopted. Twelve SNPs were detected in INHA gene. Among the exonic SNPs, three (g. 22177991A>G, g. 22178249G>C, and g. 22178414G>A) were missense mutations, resulting in the amino acid substitutions Val→Ala, Ala→Gly, and Ala→Gly, respectively. Four SNPs in the 3＇ untranslated region of INHA were predicted to either disturb or create microRNA-target interactions. Five SNPs (g. 22176870T>C, g. 22177100T>C, g. 22177149T>C, g. 22177991A>G, and g. 22178975G>A) were significantly associated with the number of eggs at 300 d of age (EN) (P < 0.05). Birds carrying GA genotype exhibited more EN than those with AA genotype (P < 0.01). In addition, quantitative real-time PCR revealed that INHA is mainly expressed in follicles on d 300 in chickens. Firstly, INHA expression increased and then decreased. The highest INHA mRNA abundance was found in the fifth largest preovulatory follicle (F5) (P < 0.01). In the prehierarchical follicles, INHA mRNA expression increased dramatically in small yellow follicles (SYF) (P < 0.01). Western blotting analysis showed that the INHA protein expression profile in the follicle was similar to its mRNA counterpart with greater expression in F5 and SYF follicles and lowest expression in F1 follicles (P < 0.05). These results suggest that INHA is a potential candidate gene improving reproductive traits in chickens.

Summary The genetic footprints of adaptations to naturally occurring tropical stress along with domestication are poorly reported in chickens. Here, by conducting population genomic analyses of 67 chickens inhabiting distinct climates, we found signals of gene flow from Tibetan chickens to Sri Lankan and Saudi Arabian breeds and identified 12 positively selected genes that are likely involved in genetic adaptations to both tropical desert and tropical monsoon island climates. Notably, in tropical desert climate, advantageous alleles of TLR7 and ZC3HAV1, which could inhibit replication of viruses in cells, suggest immune adaptation to the defense against zoonotic diseases in chickens. Furthermore, comparative genomic analysis showed that four genes (OC90, PLA2G12B, GPR17 and TNFRSF11A) involved in arachidonic acid metabolism have undergone convergent adaptation to tropical desert climate between birds and mammals. Our study offers insights into the genetic mechanisms of adaptations to tropical climates in birds and other animals and provides practical value for breeding design and medical research on avian viruses., Graphical Abstract, Highlights • Gene flow was detected from Tibetan chickens to Sri Lankan and Saudi Arabian chickens • Twelve genes showed chicken adaptations to both tropical desert and monsoon climates • Advantageous alleles of TLR7 and ZC3HAV1 showed immune adaptation in chickens • Convergent adaptation to tropical desert climate was found between birds and mammals, genomics; genetics; biological sciences; zoology.