Improving the genome assembly of rabbits with long-read sequencing

Genomic insights into adaptation to bipedal saltation and desert-like habitats of jerboas

Jerboas is a lineage of small rodents displaying atypical mouse-like morphology with elongated strong hindlimbs and short forelimbs. They have evolved obligate bipedal saltation and acute senses, and been well-adapted to vast desert-like habitats. Using a newly sequenced chromosome-scale genome of the Mongolian five-toed jerboa (Orientallactaga sibirica), our comparative genomic analyses and in vitro functional assays showed that the genetic innovations in both protein-coding and non-coding regions played an important role in jerboa morphological and physiological adaptation. Jerboa-specific amino acid substitutions, and segment insertions/deletions (indels) in conserved non-coding elements (CNEs) were found in components of proteoglycan biosynthesis pathway (XYLT1 and CHSY1), which plays an important role in limb development. Meanwhile, we found specific evolutionary changes functionally associated with energy or water metabolism (e.g., specific amino acid substitutions in ND5 and indels in CNEs physically near ROR2) and senses (e.g., expansion of vomeronasal receptors and the FAM136A gene family) in jerboas. Further dual-luciferase reporter assay verified that some of the CNEs with jerboa-specific segment indels exerted a significantly different influence on luciferase activity, suggesting changes in their regulatory function in jerboas. Our results revealed the potential molecular mechanisms underlying jerboa adaptation since the divergence from the Eocene-Oligocene transition, and provided more resources and new insights to enhance our understanding of the molecular basis underlying the phenotypic diversity and the environmental adaptation of mammals.

Comparative Bioinformatic Analysis of the Proteomes of Rabbit and Human Sex Chromosomes

Studying proteins associated with sex chromosomes can provide insights into sex-specific proteins. Membrane proteins accessible through the cell surface may serve as excellent targets for diagnostic, therapeutic, or even technological purposes, such as sperm sexing technologies. In this context, proteins encoded by sex chromosomes have the potential to become targets for X- or Y-chromosome-bearing spermatozoa. Due to the limited availability of proteomic studies on rabbit spermatozoa and poorly annotated databases for rabbits compared to humans, a bioinformatic analysis of the available rabbit X chromosome proteome (RX), as well as the human X (HX) and Y (HY) chromosomes proteome, was conducted to identify potential targets that could be accessible from the cell surface and predict which of the potential targets identified in humans might also exist in rabbits. We identified 100, 211, and 3 proteins associated with the plasma membrane or cell surface for RX, HX, and HY, respectively, of which 61, 132, and 3 proteins exhibit potential as targets as they were predicted to be accessible from the cell surface. Cross-referencing the potential HX targets with the rabbit proteome revealed an additional 60 proteins with the potential to be RX targets, resulting in a total of 121 potential RX targets. In addition, at least 53 possible common HX and RX targets have been previously identified in human spermatozoa, emphasizing their potential as targets of X-chromosome-bearing spermatozoa. Further proteomic studies on rabbit sperm will be essential to identify and validate the usefulness of these proteins for application in rabbit sperm sorting techniques as targets of X-chromosome-bearing spermatozoa.

Differential gene expression and protein-protein interaction network profiling of sulfur mustard-exposed rabbit corneas employing RNA-seq data and bioinformatics tools

Sulfur mustard (SM) ocular exposure severely damages the cornea and causes vision impairment. At present, no specific therapy exists to mitigate SM-induced corneal injury and vision loss. This study performed transcriptome profiling of naïve, SM-damaged, and SM-undamaged rabbit corneas using RNA-seq analysis and bioinformatic tools to gain a better mechanistic understanding and develop SM-specific medical countermeasures. The mRNA profiles of rabbit corneas 4 weeks post SM vapor exposure were generated using Illumina-NextSeq deep sequencing (Gene Expression Omnibus accession # GSE127708). The RNA sequences of naïve (n = 4), SM-damaged (n = 5), and SM-undamaged (n = 5) corneas were subjected to differential expression (DE) analysis after quality control profiling with FastQC. DE analysis was performed using HISAT2, StringTie, and DESeq2. The log2(FC)±2 and adjusted p˂0.05 were chosen to identify the most relevant genes. A total of 5930 differentially expressed genes (DEGs) (upregulated: 3196, downregulated: 2734) were found in SM-damaged corneas compared to naïve corneas, whereas SM-undamaged corneas showed 1884 DEGs (upregulated: 1029, downregulated: 855) compared to naïve corneas. DE profiling of SM-damaged corneas to SM-undamaged corneas revealed 985 genes (upregulated: 308, downregulated: 677). The DE profiles were subsequently subjected to signaling pathway enrichment, and protein‒protein interactions (PPIs) were analyzed. Pathway enrichment was performed for the genes associated with cellular apoptosis, death, adhesion, migration, differentiation, proliferation, extracellular matrix, and tumor necrosis factor production. To identify novel targets, we narrowed the pathway analysis to upregulated and downregulated genes associated with cell proliferation and differentiation, and PPI networks were developed. Furthermore, protein targets associated with cell differentiation and proliferation that may play vital roles in corneal fibrosis and wound healing post SM injury were identified.

Temporal rate of postmortem DNA degradation in archived tissue samples: evidence from liver and muscle

Guidelines identifying best practices for harvesting tissues that lead to optimal DNA preservation are few but are important curatorial concerns for genetic resource collections. We conducted a temporal study to establish rate of DNA degradation of tissue samples extracted from field-caught museum specimens. Five individuals of Sigmodon hispidus were collected and their liver and muscle tissues were harvested. Each tissue type was sectioned into 15 subsamples, and each was preserved in liquid nitrogen at different time intervals (2, 4, 8, 16, and 32 min; 1, 2, 4, 8, and 16 h; and 1, 2, 4, 8, and 16 days) following death. DNA was extracted using an automated robotic instrument and molecular mass profiles were determined fluorometrically. Postmortem DNA degradation was continuous and dependent on time, but also was significantly affected by differences among individual cotton rats. DNA fragments of ≥10,000 base pairs in length were present in muscle samples across all time intervals, whereas DNA fragments of this size in liver samples were no longer present after 8–16 h postmortem. DNA molecular mass profiles showed that muscle samples retained 80% of their longest fragments (≥10,000 base pairs) until 1 day postmortem, whereas liver samples retained the same percentage only until 8 min after death. Although rates of decay were measured from samples in a laboratory (not field) setting, rates of decay presented here can guide field and museum workers in best practices. Results suggest that opportunistic samples, such as those from roadkill specimens, are more likely to be of use for a variety of molecular methods when muscle is preserved. Considerations of differences in rates of degradation may also guide selection of tissue types housed in genetic resource collections, especially under space-limited circumstances.

Combined analysis of differentially expressed lncRNAs and miRNAs in liver tissues of high-fat fed rabbits by transcriptome sequencing

High-fat diet could lead to a series of metabolic diseases, including obesity, and its mechanism is not clear. In this study, the rabbit individuals were fed with high-fat diet, the liver tissues were collected, high-throughput sequencing technology was used to reveal the expression of lncRNA and miRNA difference, and the molecular regulation mechanism of lncRNA-miRNA. A total of 24,615 DE lncRNAs and 52 DE miRNAs were identified, including 15 novel discovered DE miRNAs (5 upregulated and 10 downregulated). Furthermore, five miRNAs and three mRNAs were verified by qRT-PCR, and the results showed that the expression of the DE miRNAs and DE lncRNAs in the two groups was consistent with our sequencing results. GO and KEGG analyzed 7,57,139 target genes respectively, enriching the pathways related to lipid metabolism, including mucin O-glycan biosynthesis pathway, insulin resistance and glucagon signaling pathway. Moreover, 65 targeting relationships were obtained. Among them, LOC103348122/miR-450a-5p, LOC103350359/miR-450a-3p and LOC103350429/miR-148a-5p were proposed the first time. Significantly, LOC103348122/miR-450a-5p and LOC103350429/miR-148a-5p were related to lipid metabolism in the liver. This study is of great significance to the CeRNA regulatory network related to lipid metabolism in the liver of rabbits, and provides a basis for understanding hepatic steatosis in rabbits.