Search

Abstract 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 ( $$\:{F}_{HOM}$$ , $$\:{F}_{L\&H}$$ , $$\:{F}_{VR1}$$ , $$\:{F}_{VR2}$$ , and $$\:{F}_{YAN}$$ ), as well as two segment-based estimators of runs of homozygosity (ROH, $$\:{F}_{ROH}$$ ) and homozygous-by-descent (HBD, $$\:{F}_{HBD}$$ ). 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 ( $$\:{F}_{HOM}$$ ), 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 $$\:{F}_{ROH}$$ and $$\:{F}_{HBD}$$ , 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.

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.

Abstract 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.

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.

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.

Long non-coding RNAs (lncRNAs) have been essentially involved in all kinds of important biological processes, such as cell differentiation and development. In this study, we comprehensively profiled genome-wide lncRNAs of bovine skeletal muscle in two Chinese cattle breeds using RNA-seq approach. The longissimus dorsi samples were collected from three Charolais (CH) and three Tongjiang (TJ) adult cattle, respectively. Approximately 80 million reads were obtained per library, and by which we assembled 31,009 unique transcripts. We successfully identified a total of 4155 different lncRNAs in six muscle samples, consisting of 3204 intronic lncRNAs, 737 intergenic lncRNAs, 113 sense lncRNAs and 101 antisense lncRNAs. Among them, 32 lncRNAs were differentially expressed (Fold Change >1.0, FDR < 0.05) between CH and TJ. Among these lncRNAs, 1744 identified bovine lncRNAs were located adjacent to 1146 protein-coding genes (

Abstract Parentage assignment is defined as the identification of the true parents of one focal offspring among a list of candidates and has been commonly used in zoological, ecological, and agricultural studies. Although likelihood‐based parentage assignment is the preferred method in most cases, it requires genotyping a predefined set of DNA markers and providing their population allele frequencies. In the present study, we proposed an alternative method of parentage assignment that does not depend on genotype data and prior information of allele frequencies. Our method employs the restriction site‐associated DNA sequencing (RAD‐seq) reads for clustering into the overlapped RAD loci among the compared individuals, following which the likelihood ratio of parentage assignment could be directly calculated using two parameters—the genome heterozygosity and error rate of sequencing reads. This method was validated on one simulated and two real data sets with the accurate assignment of true parents to focal offspring. However, our method could not provide a statistical confidence to conclude that the first ranked candidate is a true parent.

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.

Ketosis is a common metabolic disease in dairy cows during early lactation. However, information about the metabolomic and proteomic profiles associated with the incidence and progression of ketosis is still limited. In this study, an integrated metabolomics and proteomics approach was performed on blood serum sampled from cows diagnosed with clinical ketosis (case, ≥ 2.60 mmol/L plasma β-hydroxybutyrate; BHBA) and healthy controls (control, < 1.0 mmol/L BHBA). Samples were taken 2 weeks before parturition and 2 weeks after parturition from 19 animals (nine cases, 10 controls). All serum samples (n = 38) were subjected to Liquid Chromatography-Mass Spectrometry (LC-MS) based metabolomic analysis, and 20 samples underwent Data-Independent Acquisition (DIA) LC-MS based proteomic analysis. A total of 97 metabolites and 540 proteins were successfully identified, and multivariate analysis revealed significant differences in both metabolomic and proteomic profiles between cases and controls. We investigated clinical ketosis-associated metabolomic and proteomic changes using statistical analyses. Correlation analysis of statistically significant metabolites and proteins showed 78 strong correlations (correlation coefficient, R ≥ 0.7) between 38 metabolites and 25 proteins, which were then mapped to pathways using IMPaLA. Results showed that ketosis altered a wide range of metabolic pathways, such as metabolism, metabolism of proteins, gene expression and post-translational protein modification, vitamin metabolism, signaling, and disease related pathways. Findings presented here are relevant for identifying molecular targets for ketosis and biomarkers for ketosis detection during the transition period.

Abstract 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.

Abstract It is widely acknowledged that transcriptional diversity largely contributes to biological regulation in eukaryotes. Since the advent of second-generation sequencing technologies, a large number of RNA sequencing studies have considerably improved our understanding of transcriptome complexity. However, it still remains a huge challenge for obtaining full-length transcripts because of difficulties in the short read-based assembly. In the present study we employ PacBio single-molecule long-read sequencing technology for whole-transcriptome profiling in rabbit (Oryctolagus cuniculus). We totally obtain 36,186 high-confidence transcripts from 14,474 genic loci, among which more than 23% of genic loci and 66% of isoforms have not been annotated yet within the current reference genome. Furthermore, about 17% of transcripts are computationally revealed to be non-coding RNAs. Up to 24,797 alternative splicing (AS) and 11,184 alternative polyadenylation (APA) events are detected within this de novo constructed transcriptome, respectively. The results provide a comprehensive set of reference transcripts and hence contribute to the improved annotation of rabbit genome.

There are a few well-known indigenous breeds of Chinese rabbits in Sichuan and Fujian provinces, for which the genetic diversity and population structure have been poorly investigated. In the present study, we successfully employed the restriction-site-associated DNA sequencing (RAD-seq) approach to comprehensively discover genome-wide SNPs of 104 rabbits from four Chinese indigenous breeds: 30 Sichuan White, 34 Tianfu Black, 32 Fujian Yellow and eight Fujian Black. A total of 7,055,440 SNPs were initially obtained, from which 113,973 high-confidence SNPs (read depth ≥ 3, calling rate = 100% and biallelic SNPs) were selected to study the genetic diversity and population structure. The mean polymorphism information content (PIC) and nucleotide diversity (π) of each breed slightly varied with ranging from 0.2000 to 0.2281 and from 0.2678 to 0.2902, respectively. On the whole, Fujian Yellow rabbits showed the highest genetic diversity, which was followed by Tianfu Black and Sichuan White rabbits. The principal component analysis (PCA) revealed that the four breeds were clearly distinguishable. Our results first reveal the genetic differences among these four rabbit breeds in the Sichuan and Fujian provinces and also provide a high-confidence set of genome-wide SNPs for Chinese indigenous rabbits that could be employed for gene linkage and association analyses in the future.

Abstract Background Uncoupling proteins 2 (UCP2) plays an important role in energy regulation, previous studies suggested that UCP2 is an excellent candidate gene for human obesity and growth-related traits in cattle and chicks. The current study was designed to detect the genetic variation of UCP2 gene, and to explore the association between polymorphism of UCP2 gene and growth, carcass and meat quality traits in rabbits. Results A synonymous mutation in exon 1 and four variants in the first intron of the UCP2 gene were identified by using PCR-sequencing. The synonymous mutation c.72G>A was subsequently genotyped by MassArray system (Sequenom iPLEXassay) in 248 samples from three meat rabbit breeds (94 Ira rabbits, 83 Champagne rabbits, and 71 Tianfu black rabbits). Association analysis suggested that the individuals with AA and AG genotypes showed greater 70 d body weight (P

Growth, carcass characteristics and meat quality are the most important traits used in the rabbit industry. Identification of the candidate markers and genes significantly associated with these traits will be beneficial in rabbit breeding. In this study, we enrolled 465 rabbits, including 16 male Californian rabbits and 17 female Kangda5 line rabbits as the parental generation, along with their offspring (232 male and 200 female), in a genome-wide association study (GWAS) based on SLAF-seq technology. Bodyweight at 35, 42, 49, 56, 63 and 70 d was recorded for growth traits; and slaughter liveweight (84 d) and dressing out percentage were measured as carcass traits; and cooking loss and drip loss were measured as meat quality traits. A total of 5,223,720 SLAF markers were obtained by digesting the rabbit genome using RsaI + EcoRV-HF® restriction enzymes. After quality control, a subset of 317,503 annotated single-nucleotide polymorphisms (SNPs) was retained for subsequent analysis. A total of 28, 81 and 10 SNPs for growth, carcass and meat quality traits, respectively, were identified based on genome-wide significance (p < 3.16 × 10−7). Additionally, 16, 71 and 9 candidate genes were identified within 100 kb upstream or downstream of these SNPs. Further analysis is required to determine the biological roles of these candidate genes in determining rabbit growth, carcass traits and meat quality.

The domestic water buffalo is native to the Asian continent but through historical migrations and recent importations, nowadays has a worldwide distribution. The two types of water buffalo, i.e., river and swamp, display distinct morphological and behavioral traits, different karyotypes and also have different purposes and geographical distributions. River buffaloes from Pakistan, Iran, Turkey, Egypt, Romania, Bulgaria, Italy, Mozambique, Brazil and Colombia, and swamp buffaloes from China, Thailand, Philippines, Indonesia and Brazil were genotyped with a species-specific medium-density 90K SNP panel. We estimated the levels of molecular diversity and described population structure, which revealed historical relationships between populations and migration events. Three distinct gene pools were identified in pure river as well as in pure swamp buffalo populations. Genomic admixture was seen in the Philippines and in Brazil, resulting from importations of animals for breed improvement. Our results were largely consistent with previous archeological, historical and molecular-based evidence for two independent domestication events for river- and swamp-type buffaloes, which occurred in the Indo-Pakistani region and close to the China/Indochina border, respectively. Based on a geographical analysis of the distribution of diversity, our evidence also indicated that the water buffalo spread out of the domestication centers followed two major divergent migration directions: river buffaloes migrated west from the Indian sub-continent while swamp buffaloes migrated from northern Indochina via an east-south-eastern route. These data suggest that the current distribution of water buffalo diversity has been shaped by the combined effects of multiple migration events occurred at different stages of the post-domestication history of the species.

The polled phenotype of cattle is increasingly becoming favourable mainly because of the enhanced emphasis on animal welfare, for which the causative mutations have been reported during the past years. The Shuxuan cattle are a new synthetic breed by crossing the indigenous cattle with both Simmental and Holstein semen in Sichuan of Southwest China, in which about 15% of polled individuals have newly emerged. Because official record about POLLED genotypes for the historically imported sires is unavailable, we therefore genotyped the proposed POLLED variants of P202ID, P80kbID and P219ID among 48 polled and 16 horned Shuxuan cattle. It was first revealed that all three candidate mutations have been simultaneously introgressed into Shuxuan cattle, whereas the P202ID mutation is dominant. Furthermore, one polled animal still remains to carry none of the three candidate mutations, which suggests that further mutation(s) would also exist. Additionally, we sequenced mitochondrial DNA and found that Shuxuan cattle are composed of two matrilineal origins of Bos taurus (65.6%) and B. indicus (34.4%); and there is no origin-biased distribution of polled phenotype. In conclusion, our study first supports the recently reported novel candidate mutation of P219ID and detects simultaneous presences of all three known POLLED mutations within a cattle breed.

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 (pC 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

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.

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.

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

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

Although various computer tools have been elaborately developed to calculate a series of statistics in molecular population genetics for both small- and large-scale DNA data, there is no efficient and easy-to-use toolkit available yet for exclusively focusing on the steps of mathematical calculation. Here, we present PopSc, a bioinformatic toolkit for calculating 45 basic statistics in molecular population genetics, which could be categorized into three classes, including (i) genetic diversity of DNA sequences, (ii) statistical tests for neutral evolution, and (iii) measures of genetic differentiation among populations. In contrast to the existing computer tools, PopSc was designed to directly accept the intermediate metadata, such as allele frequencies, rather than the raw DNA sequences or genotyping results. PopSc is first implemented as the web-based calculator with user-friendly interface, which greatly facilitates the teaching of population genetics in class and also promotes the convenient and straightforward calculation of statistics in research. Additionally, we also provide the Python library and R package of PopSc, which can be flexibly integrated into other advanced bioinformatic packages of population genetics analysis.

In addition to the important role in maintaining calcium homeostasis, vitamin D3 (VD3) has also been suggested to participate in the regulation of innate and adaptive immune responses in human by enhancing antimicrobial peptides expression. In this study, we investigated the effect of supplementation of 0, 800, 1,600, 3,200 and 6,400 IU/kg of VD3 in the basal diet on the mRNA expression of antimicrobial peptides (AvBD-1 and chCATH-1) in chicken by using quantitative real-time RT-PCR. The supplementation of VD3 in the basal diet could significantly promote the mRNA expression of both AvBD-1 and chCATH-1 genes (P

This study aimed to investigate the possible involvement of Chemokine receptor 6 (CCR6) gene in genetic susceptibility to digestive disorders in rabbits. A total of 346 samples (161 cases vs 185 controls) were used for this case-control association study. One coding SNP (c.96C>T) spanning the CCR6 gene coding region was identified by direct sequencing and then genotyped using PCR-SSCP. The results revealed that individuals with allele "T" was significantly protective against digestive disorders (odds ratio=0.66, 95% confidence interval: 0.46-0.94, P<0.05), and TC/TT genotype can significantly reduce the incidence of digestive disorders (odds ratio=0.64, 95% CI: 0.41-0.99, P<0.05). we also investigated CCR6 gene mRNA expression level in different severity status of digestive disorders. The significantly higher mRNA expression level of sacculus rotundus in severe group was observed compared with mild and health group (P<0.05), in addition, the individuals with CC genotype had a higher mRNA expression than TC and TT genotype (P<0.05). [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Abstract  Tongjiang cattle are a local cattle population of Sichuan Province, China, numbering approximately half a million in 2005. They have long been grouped into the Bashan breed, although they have a unique breeding history and phenotypic characteristics, as well as a restricted geographic distribution. Morphologically, they can be divided into two groups based on the basic coat color (black and russet). In order to dissect the matrilineal components of Tongjiang cattle and to compare the body size traits of the two morphological groups, we measured five body size traits among 59 Tongjiang cattle samples and further sequenced the mtDNA D-loop sequence of 54 individuals. Among the 54 mtDNAs, 37 (68.5%) were Bos taurus types and 17 (31.5%) were Bos indicus types. Four known B. taurus haplogroups (T1–T4) and one B. indicus haplogroup (I1) were detected in these samples. Two body size traits differed significantly (P B. taurus haplogroups in one local Chinese cattle population. Our results suggest that the contribution of different matrilineal lineages to Chinese cattle might be more complex than we originally thought. [ABSTRACT FROM AUTHOR]

Many studies showed that long intergenic non-coding RNAs (lincRNAs) have multifaceted roles in various biological processes. However, roles and profiles of lincRNAs in clinical ketosis of dairy cows are still poorly understood. In this study, blood samples of both nine cows affected with clinical ketosis (i.e., case group) and ten healthy ones as control group at 2 weeks before as well as 2 weeks after parturition were collected for genome-wide analysis of lincRNA transcriptome. A total of 2154 novel intergenic transcripts were obtained from all the 38 RNA-Seq libraries, with 795 lincRNAs were identified. By comparing to protein coding transcripts, the characteristics reported in long non-coding RNAs were observed in our identified lincRNAs, such as shorter in transcript length, lower expression level, and fewer exons. Furthermore, we found thirteen differentially expressed lincRNAs (DELs) between case and control group. The cis- and trans-regulatory analysis obtained 124 target genes for these 13 DELs that were functionally involved in metal ion binding, cation binding, biological stress response, and energy metabolism. This study provides a comprehensive expression profile of lincRNAs associated with the progression and onset of clinical ketosis in transition dairy cows, which will facilitate further research on the functional roles of non-coding RNAs on cow ketosis biology. [ABSTRACT FROM AUTHOR]

Insulin receptor substrate (IRS) proteins play key roles in signal transduction in insulin and insulin-like growth factor signaling to control postnatal growth. The fat mass and obesity-associated (FTO) protein also play an essential role in postnatal growth. The aim of this study was to investigate the association between the IRS1 and FTO genes and the regulation of growth traits in rabbits. A total of nine synonymous SNPs were detected in the IRS1 coding sequence using direct sequencing, and the c.189G>T and c.2574G>A SNPs from two linkage disequilibrium blocks were further genotyped for association analysis in 216 New Zealand rabbits. The association results revealed that the TT genotype of c.189G>T and the AA genotype of c.2574G>A were significantly associated with higher body weight at 70 (BW70) and 84 (BW84) days of age and with higher average daily gain (P < 0.05). Linear-regression analysis revealed that the two-gene combination model of FTO c.499G>A and IRS1 c.2574G>A was associated with BW70 and BW84. The combination model of the GA genotype of FTO c.499G>A with the AA genotype of IRS1 c.2574G>A was associated with preferred values for BW70 and BW84. The performance values for the FTO c.499G>A genotypes after stratification with regard to the IRS1 c.189G>T genotypes revealed that the TT genotype of IRS1 c.189G>T reduced the FTO c.499G>A significance associated with BW70 and BW84. Together, our data indicated that the IRS1 gene was associated with growth traits in rabbits. The IRS1 and FTO combination model may be exploited to assist breeders in selecting rabbits with preferred body weight. [ABSTRACT FROM AUTHOR]