Evaluating genomic inbreeding of two Chinese yak (Bos grunniens) populations.

BACKGROUND: Yaks are a vital livestock in the Qinghai-Tibetan Plateau area for providing food products, maintaining sustainable ecosystems, and promoting cultural heritage. Because of uncontrolled mating, it is impossible to estimate inbreeding level of yak populations using the pedigree-based approaches. With the aims to accurately evaluate inbreeding level of two Chinese yak populations (Maiwa and Jiulong), we obtained genome-wide single nucleotide polymorphisms (SNPs) by DNA sequencing and calculated five SNP-by-SNP estimators ([Formula: see text], [Formula: see text], [Formula: see text], [Formula: see text], and [Formula: see text]), as well as two segment-based estimators of runs of homozygosity (ROH, [Formula: see text]) and homozygous-by-descent (HBD, [Formula: see text]). Functional implications were analyzed for the positional candidate genes located within the related genomic regions. RESULTS: A total of 151,675 and 190,955 high-quality SNPs were obtained from 71 Maiwa and 30 Jiulong yaks, respectively. Jiulong had greater genetic diversity than Maiwa in terms of allele frequency and nucleotide diversity. The two populations could be genetically distinguished by principal component analysis, with the mean differentiation index (Fst) of 0.0054. The greater genomic inbreeding levels of Maiwa yaks were consistently supported by all five SNP-by-SNP estimators. Based on simple proportion of homozygous SNPs ([Formula: see text]), a lower inbreeding level was indicated by three successfully sequenced old leather samples that may represent historical Maiwa yaks about five generations ago. There were 3304 ROH detected among all samples, with mean and median length of 1.97 Mb and 1.0 Mb, respectively. A total of 94 HBD segments were found among all samples, whereas 92 of them belonged to the shortest class with the mean length of 10.9 Kb. Based on the estimates of [Formula: see text] and [Formula: see text], however, there was no difference in inbreeding level between Maiwa and Jiulong yaks. Within the genomic regions with the significant Fst or enriched by ROH, we found several candidate genes and pathways that have been reported to be related to diverse production traits in farm animals. CONCLUSIONS: We successfully evaluated the genomic inbreeding level of two Chinese yak populations. Although different estimators resulted in inconsistent conclusions on their genomic inbreeding levels, our results may be helpful to implement the genetic conservation and utilization programs for the two yak populations.

Integrated analysis of microRNAs, circular RNAs, long non-coding RNAs, and mRNAs revealed competing endogenous RNA networks involved in brown adipose tissue whitening in rabbits.

BACKGROUND: The brown adipose tissue (BAT) is a target for treating obesity. BAT losses thermogenic capacity and gains a "white adipose tissue-like" phenotype ("BAT whitening") under thermoneutral environments, which is a potential factor causing a low curative effect in BAT-related obesity treatments. Circular RNAs (circRNAs) and long non-coding RNAs (lncRNAs) can act as competing endogenous RNAs (ceRNA) to mRNAs and function in various processes by sponging shared microRNAs (miRNAs). However, the roles of circRNA- and lncRNA-related ceRNA networks in regulating BAT whitening remain litter known. RESULTS: In this study, BATs were collected from rabbits at day0 (D0), D15, D85, and 2 years (Y2). MiRNA-seq was performed to investigate miRNA changes during BAT whitening. Then, a combined analysis of circRNA-seq and whole-transcriptome sequencing was used for circRNA assembly and quantification during BAT whitening. Our data showed that 1187 miRNAs and 6204 circRNAs were expressed in the samples, and many of which were identified as significantly changed during BAT whitening. Target prediction showed that D0-selective miRNAs were significantly enriched in the Ras, MAPK, and PI3K-Akt signaling pathways, and Y2-selective miRNAs were predicted to be involved in cell proliferation. The cyclization of several circRNAs could form novel response elements of key thermogenesis miRNAs at the back-splicing junction (BSJ) sites, and in combination with a dual-luciferase reporter assay confirmed the binding between the BSJ site of novel_circ_0013792 and ocu-miR-378-5p. CircRNAs and lncRNAs have high cooperativity in sponging miRNAs during BAT whitening. Both circRNA-miRNA-mRNA and lncRNA-miRNA-mRNA triple networks were significantly involved in immune response-associated biological processes. The D15-selective circRNA might promote BAT whitening by increasing the expression of IDH2. The Y2-selective circRNA-related ceRNA network and lncRNA-related ceRNA network might regulate the formation of the WAT-like phenotype of BAT via MAPK and Ras signaling pathways, respectively. CONCLUSIONS: Our work systematically revealed ceRNA networks during BAT whitening in rabbits and might provide new insight into BAT-based obesity treatments.

Genome-wide SNP discovery and genetic diversity evaluation of Liangshan cattle in China.

Liangshan cattle are a very small indigenous breed with adult weight of less than 300 Kg and have been mainly distributed in the Liangshan Yi Autonomous Prefecture of Southwestern Sichuan, China. Due to its long-term adaptation to local environments, Liangshan cattle is a valuable genetic resource and should be paid with more attentions. However, the genetic diversity of Liangshan cattle have not been specifically investigated yet, which would be required when designing the appropriate conservation and utilization programs. In this study, we successfully employed the restriction-site-associated DNA sequencing (RAD-seq) approach to explore a total of 84,854 genome-wide and high-confidence SNPs of Liangshan cattle. All these SNPs were evenly distributed through all chromosomes with an average of 98 SNPs per 1-Mb region. The nucleotide diversity, expected heterozygosity, polymorphism information content of Liangshan cattle were 0.227, 0.223 and 0.183, respectively. Furthermore, there was no obvious difference on the genetic diversity among the three studied geographical populations. In conclusion, we provided a list of SNPs that could be used in the follow-up studies for Liangshan cattle and revealed a relatively high genetic variation in this gene pool.

Genome-wide identification and characterization of long non-coding RNAs during differentiation of visceral preadipocytes in rabbit.

Emerging evidence suggests that long non-coding RNAs (lncRNAs) are critical regulators of diverse biological processes, including adipogenesis. Despite being considered an ideal animal model for studying adipogenesis, little is known about the roles of lncRNAs in the regulation of rabbit preadipocyte differentiation. In the present study, visceral preadipocytes isolated from newborn rabbits were cultured in vitro and induced for differentiation, and global lncRNA expression profiles of adipocytes collected at days 0, 3, and 9 of differentiation were analyzed by RNA-seq. A total of 2066 lncRNAs were identified from nine RNA-seq libraries. Compared to protein-coding transcripts, lncRNA transcripts exhibited characteristics of a longer length and lower expression level. Furthermore, 486 and 357 differentially expressed (DE) lncRNAs were identified when comparing day 3 vs. day 0 and day 9 vs. day 3, respectively. Target genes of DE lncRNAs were predicted by the cis-regulating approach. Prediction of functions revealed that DE lncRNAs when comparing day 3 vs. day 0 were involved in gene ontology (GO) terms of developmental growth, growth, developmental cell growth, and stem cell proliferation, and involved in Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways of PI3K-Akt signaling pathway, fatty acid biosynthesis, and the insulin signaling pathway. The DE lncRNAs when comparing day 9 vs. day 3 were involved in GO terms that associated with epigenetic modification and were involved in the KEGG pathway of cAMP signaling pathway. This study provides further insight into the regulatory function of lncRNAs in rabbit visceral adipose and facilitates a better understanding of different stages of preadipocyte differentiation.

Analysis of gene expression profiles at different stages during preadipocyte differentiation in rabbits.

Excessive accumulation of fat is harmful to human health. The preadipocyte differentiation is a critical process of fat development. Studying the expression profiles of genes related to preadipocyte differentiation contributes to understanding of the mechanism of fat accumulation. Despite being considered an ideal animal model for studying adipogenesis, little is known about the gene expression profiles at different stages during preadipocyte differentiation in rabbits. In the present study, rabbit preadipocytes were cultured in vitro and induced for differentiation, and gene expression profiles of adipocytes collected at days 0, 3, and 9 of differentiation were analyzed by RNA-seq. We identified 1352 differentially expressed genes (DEGs) when comparing day 3 with day 0 and identified 888 DEGs when comparing day 9 with day 3. Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed that the PPAR signaling pathway and PI3K-Akt signaling pathway were significantly enriched by the DEGs that up-regulated within the period of day 0 - day 3, and the GO terms and KEGG pathways that were associated with cell cycle were enriched by the DEGs that up-regulated within the period of day 3 - day 9. The DEGs that specifically up-regulated within the period of day 0 - day 3 might play roles in the cytoplasm, and the DEGs that specifically up-regulated within the period of day 3 - day 9 might act in the nucleus. The protein-protein interaction (PPI) network constructed by DEGs showed that hub node genes might modulate rabbit preadipocyte differentiation via regulating cell cycle.

Sequence and Evolutionary Features for the Alternatively Spliced Exons of Eukaryotic Genes.

Alternative splicing of pre-mRNAs is a crucial mechanism for maintaining protein diversity in eukaryotes without requiring a considerable increase of genes in the number. Due to rapid advances in high-throughput sequencing technologies and computational algorithms, it is anticipated that alternative splicing events will be more intensively studied to address different kinds of biological questions. The occurrences of alternative splicing mean that all exons could be classified to be either constitutively or alternatively spliced depending on whether they are virtually included into all mature mRNAs. From an evolutionary point of view, therefore, the alternatively spliced exons would have been associated with distinctive biological characteristics in comparison with constitutively spliced exons. In this paper, we first outline the representative types of alternative splicing events and exon classification, and then review sequence and evolutionary features for the alternatively spliced exons. The main purpose is to facilitate understanding of the biological implications of alternative splicing in eukaryotes. This knowledge is also helpful to establish computational approaches for predicting the splicing pattern of exons.

Genome-wide identification and characterization of long non-coding RNAs during postnatal development of rabbit adipose tissue.

BACKGROUND: The rabbit is widely used as an important experimental model for biomedical research, and shows low adipose tissue deposition during growth. Long non-coding RNAs (lncRNAs) are associated with adipose growth, but little is known about the function of lncRNAs in the rabbit adipose tissue. METHODS: Deep RNA-sequencing and comprehensive bioinformatics analyses were used to characterize the lncRNAs of rabbit visceral adipose tissue (VAT) at 35, 85 and 120 days after birth. Differentially expressed (DE) lncRNAs were identified at the three growth stages by DESeq. The cis and trans prediction ways predicted the target genes of the DE lncRNAs. To explore the function of lncRNAs, Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed on the candidate genes. RESULTS: A total of 991,157,544 clean reads were generated after RNA-Seq of the three growth stages, of which, 30,353 and 107 differentially expressed (DE) lncRNAs were identified. Compared to the protein-coding transcripts, the rabbit lncRNAs shared some characteristics such as shorter length and fewer exons. Cis and trans target gene prediction revealed, 43 and 64 DE lncRNAs respectively, corresponding to 72 and 20 protein-coding genes. GO enrichment and KEGG pathway analyses revealed that the candidate DE lncRNA target genes were involved in oxidative phosphorylation, glyoxylate and dicarboxylate metabolism, and other adipose growth-related pathways. Six DE lncRNAs were randomly selected and validated by q-PCR. CONCLUSIONS: This study is the first to profile the potentially functional lncRNAs in the adipose tissue growth in rabbits, and contributes to our understanding of mammalian adipogenesis.

Comparative Study on the Genetic Diversity of GHR Gene in Tibetan Cattle and Holstein Cows.

Due to the phenotype-based artificial selection in domestic cattle, the underlying functional genes may be indirectly selected and show decreasing diversity in theory. The growth hormone receptor (GHR) gene has been widely proposed to significantly associate with critical economic traits in cattle. In the present study, we comparatively studied the genetic diversity of GHR in Tibetan cattle (a traditional unselected breed, n = 93) and Chinese Holstein cow (the intensively selected breed, n = 94). The Tibetan yak (n = 38) was also included as an outgroup breed. A total of 21 variants were detected by sequencing 1279 bp genomic fragments encompassing the largest exon 9. Twelve haplotypes (H1∼H12) constructed by 15 coding SNPs were presented as a star-like network profile, in which haplotype H2 was located at the central position and almost occupied by Tibetan yaks. Furthermore, H2 was also identical to the formerly reported sequence specific to African cattle. Only haplotype H5 was simultaneously shared by all three breeds. Tibetan cattle showed higher nucleotide diversity (0.00215 ± 0.00015) and haplotype diversity (0.678 ± 0.026) than Holstein cow. Conclusively, we found Tibetan cattle have retained relatively high genetic variation of GHR. The predominant presence of African cattle specific H2 in the outgroup yak breed would highlight its ancestral relationship, which may be used as one informative molecular marker in the phylogenetic studies.

Expression of MyHC genes, composition of muscle fiber type and their association with intramuscular fat, tenderness in skeletal muscle of Simmental hybrids.

In adult bovine skeletal muscle, it expressed four isoforms of Myosin heavy chain (MyHC) gene, MyHC-I, MyHC-IIa, MyHC-IIb, and MyHC-IIx that are translated into different structural protein myofibrils, and then further form different types of muscle fiber. In the studies, our objective is to reveal the expression patterns of MyHC genes in longissimus dorsi (Ld), semitendinosus (Se) and soleus (Sol) of Simmental hybrids cattle, and their association with intramuscular fat (IMF) content and meat shearing force (MSF). The muscle tissue of Ld, Se and Sol were collected from 6, 12 and 36-month old Simmental hybrids respectively, then the expression levels of MyHCs were examined by real-time PCR, at the same time, IMF, MSF and muscle type were measured with chemical assessment, shearing force measurer and immunostaining respectively. Our results showed that t Ld, Se, and Sol expressed MyHC-I, MyHC-IIa and MyHC-IIx isoforms but not MyHC-IIb, furthermore MyHC-I, MyHC-IIa and MyHC-IIx had different expression patterns in different skeletal muscle. The expression of MyHC-I in Se and Sol, MyHC-IIa in Ld, Se, and Sol, and MyHC-IIx in Sol was decreased with increasing age. The highest expression of MyHC-I in Ld, and MyHC-IIx in Ld and Se was observed in 12-month-old animals. The percentage of type-IIa fiber approximately occupied 70-80 % among various muscle fiber of Ld, Se and Sol. The percentage of different type fiber was not related to IMF content and MSF, but the expression levels of MyHC-I and MyHC-IIa were negatively related to IMF content (r = -0.724, and -0.681, respectively) and MSF (r = -0.672, and -0.641, respectively). The expression level of MyHC-IIx was also negatively related to MSF (r = -0.655). In conclusion, MyHC gene might be considered as a negative factor in genetic selection of IMF content and MSF.

Effects of variants in proopiomelanocortin and neuropeptide y genes on growth, carcass, and meat quality traits in rabbits.

Appetite-related neuropeptides proopiomelanocortin (POMC) and Neuropeptide Y (NPY) are essential for regulating feeding behavior and energy homeostasis. The objective of this study was to evaluate the effects of variants in POMC and NPY genes on growth, carcass and meat quality traits in rabbits. A total of six SNPs were identified for POMC (n = 2) and NPY (n = 4) genes by direct sequencing. Three SNPs were subsequently genotyped by using MassArray system (Sequenom iPLEXassay) in 235 individuals, which belong to three meat rabbit breeds, including 93 Ira rabbits; 81 Champagne rabbits and 61 Tianfu black rabbits. The SNP c.112-12G>T was in intron-exon boundaries (intron 1) of POMC gene, and the association analysis showed that individuals with TT genotype had a greater 84 d body weight (BW84), eviscerated weight and semi-eviscerated weight than those with GT genotype (p<0.05); the TT individuals were also higher than those GG in the ripe meat ratio (RMR) (p<0.05). The g.1778G>C SNP, which was in complete linkage with other three SNPs (g.1491G>A, g.1525G>T and g.1530C>T) in intron 1 of NPY gene, was significantly correlated with eviscerated slaughter percentage and semi-eviscerated slaughter percentage in rabbits, and the individuals with CC genotype had a better performance than CG genotype (p<0.05). These findings would provide primary clues for the biological roles of POMC and NPY underlying the rabbit growth-related traits.

Parameter Estimation of Host Genomic and Gut Microbiota Contribution to Growth and Feed Efficiency Traits in Meat Rabbits.

Rabbits can efficiently utilize plant fibers that are indigestible to humans, and hence may contribute to the alleviation of feed-food competition. Therefore, it is economically and ecologically important to genetically improve the growth performance and feed efficiency of meat rabbits. In this study, we combined pedigree, genomic, and gut microbiota data to estimate genetic and microbial parameters for nine growth and feed efficiency traits of 739 New Zealand White rabbits, including body weight (BW) at 35 (BW35), 70 (BW70), and 84 (BW84) days of age, and average daily gain (ADG), feed conversion ratio (FCR), and residual feed intake (RFI) within two age intervals of 35-70 days (ADG70, FCR70, and RFI70) and 35-84 days (ADG84, FCR84, and RFI84). Based on single-step genomic best linear unbiased prediction, three BW traits and two ADG traits had the high estimates (±standard error, SE) of heritability, ranging from 0.44 ± 0.13 of BW35 to 0.66 ± 0.08 of BW70. Moderate heritabilities were observed for RFI70 (0.22 ± 0.07) and RFI84 (0.29 ± 0.07), whereas the estimates did not significantly deviate from zero for the two FCR traits. There was moderate positive genetic correlation (±SE) between BW70 and ADG70 (0.579 ± 0.086), but BW70 did not correlate with RFI70. Based on microbial best linear unbiased prediction, the estimates of microbiability did not significantly deviate from zero for any trait. Based on the combined use of genomic and gut microbiota data, the parameters obtained in this study could help us to implement efficient breeding schemes in meat rabbits.

Case-control study and mRNA expression analysis reveal the MyD88 gene is associated with digestive disorders in rabbit.

As in humans, significant associations between Toll-like receptor 4 (TLR4) and digestive disorders have been identified in rabbit and dog. However, as an essential adaptor downstream of TLR4, the genetic variation of myeloid differentiating factor 88 (MyD88) and its association with digestive disorders have remained unknown. In this study, we detected 10 single nucleotide polymorphisms (SNPs) in the entire genomic region of rabbit MyD88. The genetic variation in susceptibility to digestive disorders for the only coding SNP (synonymous c.699T>C) was studied in Yaan (183 cases and 142 controls) and Chengdu populations (145 cases and 140 controls). A case-control association study revealed that individuals with the C allele had significant protection against digestive disorders in the Yaan population (OR = 0.71; 95% CI, 0.51-0.99; P < 0.05), the Chengdu population (OR = 0.55; 95% CI, 0.39-0.78; P < 0.01) and for joint analysis (OR = 0.62; 95% CI, 0.49-0.79; P < 0.01). We also experimentally induced digestive disorders by feeding a fiber-deficient diet and found that increased susceptibility was significantly associated with higher MyD88 mRNA expression (P < 0.05). The lowest MyD88 mRNA expression was observed in individuals carrying the protective CC genotype. These results suggest that MyD88 is one of the most plausible candidate genes in relation to digestive disorders in rabbit. Further studies are required to explore the biological implications of MyD88 in the pathogenesis of digestive disorders.

Rapid Genotyping of MSTN Gene Polymorphism Using High-resolution Melting for Association Study in Rabbits.

The myostatin (MSTN) gene, as a negative regulator of skeletal muscle growth, has been proposed to be associated with production traits in farm animals. In the present study, a T/C variant at -125 bp (relative to ATG start codon) of 5'regulatory region of rabbit MSTN was identified by direct sequencing. Two hundred and twenty two rabbits, which were randomly sampled from 3 breeds (Ira rabbits, Champagne rabbits and Tianfu black rabbits), were genotyped by high-resolution melting (HRM). Comparing the genotyping results of 47 samples with direct sequencing, the HRM showed high sensitivity (0.96) and high specificity (0.98). In the three rabbit breeds, the allele C was the predominant allele. The polymorphic site showed high heterozygosity (He = 0.48) and high effective number of alleles (Ne = 1.91). The genetic diversity was reasonably informative (0.25

Polymorphism of NLRP3 Gene and Association with Susceptibility to Digestive Disorders in Rabbit.

NLR family pyrin domain containing 3 (NLRP3) is a key component of the inflammasome, whose assembly is a crucial part of the innate immune response. The aim of the present study was to evaluate the association between exon 3 polymorphisms of NLRP3 and the susceptibility to digestive disorders in rabbits. In total, five coding single-nucleotide polymorphisms (cSNPs) were identified; all of which are synonymous. Among them, c.456 C> and c.594 G> were further genotyped for association analysis based on case-control design (n =162 vs n =102). Meanwhile, growing rabbits were experimentally induced to digestive disorders by feeding a fiber-deficient diet, subsequently they were subjected to mRNA expression analysis. Association analysis revealed that haplotype H1 (the two cSNPs: GT) played a potential protective role against digestive disorders (p<0.001). The expression of NLRP3 in the group H1HX1 (H1HX1 is composed of H1H1, H1H3 and H1H4) was the lowest among four groups which were classified by different types of diplotypes. Those results suggested that the NLRP3 gene was significantly associated with susceptibility to digestive disorders in rabbit.

Identification and Association of SNPs in TBC1D1 Gene with Growth Traits in Two Rabbit Breeds.

The TBC1D1 plays a key role in body energy homeostasis by regulating the insulin-stimulated glucose uptake in skeletal muscle. The present study aimed to identify the association between genetic polymorphisms of TBC1D1 and body weight (BW) in rabbits. Among the total of 12 SNPs detected in all 20 exons, only one SNP was non-synonymous (c.214G>A. p.G72R) located in exon 1. c.214G>A was subsequently genotyped among 491 individuals from two rabbit breeds by the high-resolution melting method. Allele A was the predominant allele with frequencies of 0.7780 and 0.6678 in European white rabbit (EWR, n = 205) and New Zealand White rabbit (NZW, n = 286), respectively. The moderate polymorphism information content (0.250.05). Our results implied that the c.214G>A of TBC1D1 gene might be one of the candidate loci affecting the trait of 35 d BW in the rabbit.

Single Nucleotide Polymorphisms of NLRP12 Gene and Association with Non-specific Digestive Disorder in Rabbit.

The NLRP12 (NLR family, pyrin domain containing 12) serves as a suppressor factor in the inflammatory response and protects the host against inflammation-induced damage. In the present study, we aimed to study the polymorphisms of NLRP12 gene and its association with susceptibility to non-specific digestive disorder (NSDD) in rabbits. We re-sequenced the entire coding region of the rabbit NLRP12 gene and detected a total of 19 SNPs containing 14 synonymous and five non-synonymous variations. Among them, the coding SNP (c.1682A>G), which would carry a potential functional implication, was subsequently subjected to genotyping for case-control association study (272 cases and 267 controls). The results revealed that allele A was significantly protective against NSDD with an odds ratio value of 0.884 (95% confidence interval, 0.788 to 0.993; p = 0.038). We also experimentally induced NSDD in growing rabbits by feeding a fibre-deficient diet and subsequently investigated NLRP12 mRNA expression. The mRNA expression of NLRP12 in healthy status was significantly higher than that in severe NSDD (p = 0.0016). The highest expression was observed in individuals carrying the protective genotype AA (p = 0.0108). These results suggested that NLRP12 was significantly associated with the NSDD in rabbits. However, the precise molecular mechanism of NLRP12 involving in the development of rabbit NSDD requires further research.

Polymorphisms of mitochondrial ATPase 8/6 genes and association with milk production traits in Holstein cows.

The maternal effect has been widely proposed to affect the production traits in domestic animals. However, the sequence polymorphisms of mitochondrial DNA (mtDNA) and association with milk production traits in Holstein cows have remained unclear. In this study, we investigated the single nucleotide polymorphisms (SNPs) of mtDNA ATPase 8/6 genes and association with four milk production traits of interest in 303 Holstein cows. A total of 18 SNPs were detected among the 842 bp fragment of ATPase 8/6 genes, which determined six haplotypes of B. taurus (H1-H4) and B. indicus (H5-H6). The mixed model analysis revealed that there was significant association between haplotype and 305-day milk yield (MY). The highest MY was observed in haplotype H4. However, we did not detect statistically significant differences among haplotypes for the traits of milk fat (MF), milk protein (MP), and somatic cell count (SC). The overall haplotype diversity and nucleotide diversity of ATPase 8/6 genes were 0.563 ± 0.030 and 0.00609 ± 0.00043, respectively. The results suggested that mitochondrial ATPase 8/6 genes could be potentially used as molecular marker to genetically improve milk production in Holstein cows.

Favoring Expression of Yak Alleles in Interspecies F1 Hybrids of Cattle and Yak Under High-Altitude Environments.

Both cis- and trans-regulation could cause differential expression between the parental alleles in diploid species that might have broad biological implications. Due to the relatively distant genetic divergence between cattle and yak, as well as their differential adaptation to high-altitude environments, we investigated genome-wide allelic differential expression (ADE) in their F1 hybrids using Nanopore long-read RNA-seq technology. From adult F1 hybrids raised in high-altitude, ten lung and liver tissues were individually sequenced for producing 31.6 M full-length transcript sequences. Mapping against autosomal homologous regions between cattle and yak, we detected 17,744 and 14,542 protein-encoding genes expressed in lung and liver tissues, respectively. According to the parental assignments of transcript sequences, a total of 3,381 genes were detected to show ADE in at least one sample. There were 186 genes showing ubiquitous ADE in all the studied animals, and among them 135 and 37 genes had consistent higher expression of yak and cattle alleles, respectively. Functional analyses revealed that the genes with favoring expression of yak alleles have been involved in the biological progresses related with hypoxia adaptation and immune response. In contrast, the genes with favoring expression of cattle alleles have been enriched into different biological progresses, such as secretion of endocrine hormones and lipid metabolism. Our results would support unequal contribution of parental genes to environmental adaptation in the F1 hybrids of cattle and yak.

Dynamics of transcriptome and chromatin accessibility revealed sequential regulation of potential transcription factors during the brown adipose tissue whitening in rabbits.

Brown adipose tissue (BAT) represents a valuable target for treating obesity in humans. BAT losses of thermogenic capacity and gains a "white adipose tissue-like (WAT-like)" phenotype (BAT whitening) under thermoneutral environments, which could lead to potential low therapy responsiveness in BAT-based obesity treatments. However, the epigenetic mechanisms of BAT whitening remain largely unknown. In this study, BATs were collected from rabbits at day0 (D0), D15, D85, and 2 years (Y2). RNA-sequencing (RNA-seq) and the assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) were performed to investigate transcriptome and chromatin accessibility of BATs at the four whitening stages, respectively. Our data showed that many genes and chromatin accessible regions (refer to as "peaks") were identified as significantly changed during BAT whitening in rabbits. The BAT-selective genes downregulated while WAT-selective genes upregulated from D0 to Y2, and the de novo lipogenesis-related genes reached the highest expression levels at D85. Both the highly expressed genes and accessible regions in Y2 were significantly enriched in immune response-related signal pathways. Analysis of different relationships between peaks and their nearby genes found an increased proportion of the synchronous changes between chromatin accessibility and gene expression during BAT whitening. The synergistic changes between the chromatin accessibility of promoter and the gene expression were found in the key adipose genes. The upregulated genes which contained increased peaks were significantly enriched in the PI3K-Akt signaling pathway, steroid biosynthesis, TGF-beta signaling pathway, osteoclast differentiation, and dilated cardiomyopathy. Moreover, the footprinting analysis suggested that sequential regulation of potential transcription factors (TFs) mediated the loss of thermogenic phenotype and the gain of a WAT-like phenotype of BAT. In conclusion, our study provided the transcriptional and epigenetic frameworks for understanding BAT whitening in rabbits for the first time and might facilitate potential insights into BAT-based obesity treatments.

Genome-Wide Association Studies for Growth Curves in Meat Rabbits Through the Single-Step Nonlinear Mixed Model.

Growth is a complex trait with moderate to high heritability in livestock and must be described by the longitudinal data measured over multiple time points. Therefore, the used phenotype in genome-wide association studies (GWAS) of growth traits could be either the measures at the preselected time point or the fitted parameters of whole growth trajectory. A promising alternative approach was recently proposed that combined the fitting of growth curves and estimation of single-nucleotide polymorphism (SNP) effects into single-step nonlinear mixed model (NMM). In this study, we collected the body weights at 35, 42, 49, 56, 63, 70, and 84 days of age for 401 animals in a crossbred population of meat rabbits and compared five fitting models of growth curves (Logistic, Gompertz, Brody, Von Bertalanffy, and Richards). The logistic model was preferably selected and subjected to GWAS using the approach of single-step NMM, which was based on 87,704 genome-wide SNPs. A total of 45 significant SNPs distributed on five chromosomes were found to simultaneously affect the two growth parameters of mature weight (A) and maturity rate (K). However, no SNP was found to be independently associated with either A or K. Seven positional genes, including KCNIP4, GBA3, PPARGC1A, LDB2, SHISA3, GNA13, and FGF10, were suggested to be candidates affecting growth performances in meat rabbits. To the best of our knowledge, this is the first report of GWAS based on single-step NMM for longitudinal traits in rabbits, which also revealed the genetic architecture of growth traits that are helpful in implementing genome selection.