Skeletal muscle development plays an important role in muscle quality and yield, which decides the economic value of livestock. Long non-coding RNAs (lncRNAs) have been reported to be associated with skeletal muscle development. However, little is revealed about the function of lncRNAs in rabbits’ muscle development. LncRNAs and mRNAs in two rabbit breeds (ZIKA rabbits (ZKR) and Qixin rabbits (QXR)) with different growth rates at three developmental stages (0 day, 35 days, and 84 days after birth) were researched by transcriptome sequencing. Differentially expressed lncRNAs and mRNAs were identified for two rabbit breeds at the same stages by DESeq package. Co-expression correlation analysis of differentially expressed lncRNAs and mRNAs were performed to construct lncRNA–mRNA pairs. To explore the function of lncRNAs, Gene Ontology (GO) analysis of co-expression mRNAs in lncRNA–mRNA pairs were performed. In three comparisons, there were 128, 109, and 115 differentially expressed lncRNAs, respectively. LncRNAs TCONS_00013557 and XR_518424.2 differentially expressed in the two rabbit breeds might play important roles in skeletal muscle development, for their co-expressed mRNAs were significantly enriched in skeletal muscle development related GO terms. This study provides potentially functional lncRNAs in skeletal muscle development of two rabbit breeds and might be beneficial to the production of rabbits.

The domestic rabbit (Oryctolagus cuniculus f. domesticus) is a very important variety in biomedical research and agricultural animal breeding. Due to the different geographical areas in which rabbit breeds originated, and the long history of domestication/artificial breeding, rabbits have experienced strong selection pressure, which has shaped many traits of most rabbit varieties, such as color and weight. An efficient genome-wide single-nucleotide polymorphism (SNP) detection strategy is genotyping-by-sequencing (GBS), which has been widely used in many organisms. This study attempted to explore bi-allelic SNPs associated with fur color and weight-related traits using GBS in five rabbit breeds. The data consisted of a total 831,035 SNPs in 150 individuals from Californian rabbits (CF), German Zika rabbits (ZK), Qixing rabbits (QX), Sichuan grey rabbits (SG), and Sichuan white rabbits (SW). In addition, these five breeds of rabbits were obviously independent populations, with high genetic differentiation among breeds and low genetic diversity within breeds. A total of 32,144 SNP sites were identified by selective sweep among the different varieties. The genes that carried SNP loci in these selected regions were related to important traits (fur color and weight) and signal pathways, such as the MAPK/ERK signaling pathway and the Hippo signaling pathway. In addition, genes related to fur color and weight were identified, such as ASIPs, MITFs and KITs, ADCY3s, YAPs, FASs, and ACSL5s, and they had more SNP sites. The research offers the foundation for further exploration of molecular genetic markers of SNPs that are related to traits.

Rabbit hemorrhagic disease (RHD) is an acute fatal disease caused by the rabbit hemorrhagic disease virus (RHDV). Since the first outbreaks of type 2 RHDV (RHDV2) in April 2020 in China, the persistence of this virus in the rabbit population has caused substantial economic losses in rabbit husbandry. Previous failures in preventing RHDV2 prompted us to further investigate the immune mechanisms underlying the virus’s pathogenicity, particularly concerning the spleen, a vital component of the mononuclear phagocyte system (MPS). For this, a previous RHDV2 isolate, CHN/SC2020, was utilized to challenge naive adult rabbits. Then, the splenic transcriptome was determined by RNA-Seq. This study showed that the infected adult rabbits had 3148 differentially expressed genes (DEGs), which were associated with disease, signal transduction, cellular processes, and cytokine signaling categories. Of these, 100 upregulated DEGs were involved in inflammatory factors such as IL1α, IL-6, and IL-8. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that these DEGs were significantly enriched in the cytokine–cytokine receptor interaction signaling pathway, which may play a vital role in CHN/SC2020 infection. At the same time, proinflammatory cytokines and chemokines were significantly increased in the spleen at the late stages of infection. These findings suggested that RHDV2 (CHN/SC2020) might induce dysregulation of the cytokine network and compromise splenic immunity against viral infection, which expanded our understanding of RHDV2 pathogenicity.