The First High-quality Reference Genome of Sika Deer Provides Insights into High-tannin Adaptation.

Sika deer are known to prefer oak leaves, which are rich in tannins and toxic to most mammals; however, the genetic mechanisms underlying their unique ability to adapt to living in the jungle are still unclear. In identifying the mechanism responsible for the tolerance of a highly toxic diet, we have made a major advancement by explaining the genome of sika deer. We generated the first high-quality, chromosome-level genome assembly of sika deer and measured the correlation between tannin intake and RNA expression in 15 tissues through 180 experiments. Comparative genome analyses showed that the UGT and CYP gene families are functionally involved in the adaptation of sika deer to high-tannin food, especially the expansion of the UGT family 2 subfamily B of UGT genes. The first chromosome-level assembly and genetic characterization of the tolerance to a highly toxic diet suggest that the sika deer genome may serve as an essential resource for understanding evolutionary events and tannin adaptation. Our study provides a paradigm of comparative expressive genomics that can be applied to the study of unique biological features in non-model animals.

Genome-wide study on genetic diversity and phylogeny of five species in the genus Cervus.

BACKGROUND: Previous investigations of phylogeny in Cervus recovered many clades without whole genomic support. METHODS: In this study, the genetic diversity and phylogeny of 5 species (21 subspecies/populations from C. unicolor, C. albirostris, C. nippon, C. elaphus and C. eldii) in the genus Cervus were analyzed using reduced-representation genome sequencing. RESULTS: A total of 197,543 SNPs were identified with an average sequencing depth of 16 x. A total of 21 SNP matrices for each subspecies/population and 1 matrix for individual analysis were constructed, respectively. Nucleotide diversity and heterozygosity analysis showed that all 21 subspecies/populations had different degrees of genetic diversity. C. eldii, C. unicolor and C. albirostris showed relatively high expected and observed heterozygosity, while observed heterozygosity in C. nippon was the lowest, indicating there was a certain degree of inbreeding rate in these subspecies/populations. Phylogenetic ML tree of all Cervus based on the 21 SNP matrices showed 5 robustly supported clades that clearly separate C. eldii, C. unicolor, C. albirostris, C. elaphus and C. nippon. Within C. elaphus clade, 4 subclades were well differentiated and statistically highly supported: C. elaphus (New Zealand), C. e. yarkandensis, C. c. canadensis and the other grouping the rest of C. canadensis from China. In the C. nippon clade, 2 well-distinct subclades corresponding to C. n. aplodontus and other C. nippon populations were separated. Phylogenetic reconstruction indicated that the first evolutionary event of the genus Cervus occurred approximately 7.4 millions of years ago. The split between C. elaphus and C. nippon could be estimated at around 3.6 millions of years ago. Phylogenetic ML tree of all samples based on individual SNP matrices, together with geographic distribution, have shown that there were 3 major subclades of C. elaphus and C. canadensis in China, namely C. e. yarkandensis (distributed in Tarim Basin), C. c. macneilli/C. c. kansuensis/C. c. alashanicus (distributed in middle west of China), and C. c. songaricus/C. c. sibiricus (distributed in northwest of China). Among them, C. e. yarkandensis was molecularly the most primitive subclade, with a differentiation dating back to 0.8-2.2 Myr ago. D statistical analysis showed that there was high probability of interspecific gene exchange between C. albirostris and C. eldii, C. albirostris and C. unicolor, C. nippon and C. unicolor, and there might be 2 migration events among 5 species in the genus Cervus. CONCLUSIONS: Our results provided new insight to the genetic diversity and phylogeny of Cervus deer. In view of the current status of these populations, their conservation category will need to be reassessed.

[Effect of anti-androgen and progesterone on development of regrowth antlers in 6 year old sika deer].

To investigate the effects of anti-androgen drugs and melengestrol acetate (MGA) on development of regrowth antlers in 6 year old sika deer, twenty healthysika deerwith similar body weight and antler weightwere randomly divided into five groups by using single factor test design: flutamide (n=4), bicalutamide (n=4), progesterone acetate (CPA, n=4), melengestrol acetate (MGA, n=4), control(n=4). All deer were fed with same diets and were housed outside together in an opened fence of 15 m×30 m with free access to water and feed. Treatment groups were injected subcutaneously sustained-release agents of the four drugs respectively when two-branched antlers were harvested. The control group had no special treatment. In the experiment period of 60 d, blood sampleswere collected for 4 times for each deer. The concentration of testosterone in plasma was tested and analyzed to compare the changes between different groups. Development of regrowth antlers was observed. At the end of the experiment, regrowth antlers were weighted and analyzed. The resultsshowed that the weights of regrowth antlers in treatment groups were significantly greater than those from control group and the weight gain (as compared with the control group) was 100.50%, 64.46%, 87.16% and 117.46% respectively in flutamide group, bicalutamide group, progesterone acetate group and melengestrol acetate group. For plasma testosterone concentration, it was not significantly different in the early stage (in the first 35 d), but at the end of the experimen, it was significantly higher than that of earlier stage (P<0.01) in various groups. Testosterone concentration of flutamide treated group was significantly lower than that of the other groups (P<0.01), while the level inbicalutamide and MGA treated groups was significantly higher than that in other groups (P<0.01). The results showed that both anti-androgen drugs and MGA treatment promoted the development of regrowth antlers and increased the weight of regrowth antlers, where the effect was most significant by MGA treatment. From the morphological observation of the antlers, it was found that anti-androgen and MGA treatments prolonged the growth period of regrowth antlers through delaying the ossification of antlers. However, plasma testosterone concentration was not affected by the treatments.

Histone deacetylase inhibitor stimulates E2 and P4 secretion in sika deer ovarian granulosa cells at a moderate dose.

The histone deacetylase inhibitor (HDACi) and tumor suppressor play an important role in genome reorganization and epigenetic regulation. In this study, granulosa cells (GCs) isolated from sika deer ovaries were cultured and treated with different concentrations of trichostatin A (TSA) for 48 h. It was found that TSA inhibited GCs proliferation and induced GCs apoptosis by upregulating expression of BAX, meanwhile, downregulating expression of GLUT3, GLUT8, BCL-XL. In addition, TSA caused cell cycle arrest at the G1 and G2/M phase accompanied by reducing expression of Cyclin D2 and CDK4. TSA pretreatment increased DNMT3a, DNMT1, HDAC1, and HAT1 expression, and attenuated them when TAS higher than 50 nM. The protein levels of H3K9ac and H4K8ac in GCs were increased at 48 h after TSA treatment. TSA stimulated the secretion of estradiol and progesterone at a moderate dose. Our data suggest that TSA is important as a regulator of steroid hormone synthesis in granulosa cells during follicular development in the sika deer ovary.

Screening weight related genes of velvet antlers by whole genome re-sequencing.

The velvet antler is a special organ that has important biological significance for deer, and its growth is a complicated biological metabolism process. Growing evidence suggests that genetics factors play essential roles in the weight of velvet antlers. In this study, we investigated five sika deer (Cervus nippon) populations under the same feeding condition, and screened genetic variations in the 100 samples (including 50 heavy and 50 light velvet antler weight samples) by whole genome re-sequencing. The results showed that 94 genetic variations were related to the velvet antler weight, among which two single nucleotide polymorphism (SNP) sites were located on the exon regions of OAS2 and ALYREF/THOC4, respectively. Furthermore, ALYREF/THOC4 is highly expressed in the velvet antler. The biological functions of these genetic variations were highly related to the growth and development of deer velvet antlers. Collectively, we screened genes related to the velvet antler weight in sika deer populations by whole genome re-sequencing and identified 94 sites as candidate genetic variations related to the velvet antler weight. We hope that it will contribute to further mechanistic studies of velvet antler development and weight variations.

Transcriptome analysis of sika deer in China.

Sika deer is of great commercial value because their antlers are used in tonics and alternative medicine and their meat is healthy and delicious. The goal of this study was to generate transcript sequences from sika deer for functional genomic analyses and to identify the transcripts that demonstrate tissue-specific, age-dependent differential expression patterns. These sequences could enhance our understanding of the molecular mechanisms underlying sika deer growth and development. In the present study, we performed de novo transcriptome assembly and profiling analysis across ten tissue types and four developmental stages (juvenile, adolescent, adult, and aged) of sika deer, using Illumina paired-end tag (PET) sequencing technology. A total of 1,752,253 contigs with an average length of 799 bp were generated, from which 1,348,618 unigenes with an average length of 590 bp were defined. Approximately 33.2 % of these (447,931 unigenes) were then annotated in public protein databases. Many sika deer tissue-specific, age-dependent unigenes were identified. The testes have the largest number of tissue-enriched unigenes, and some of them were prone to develop new functions for other tissues. Additionally, our transcriptome revealed that the juvenile-adolescent transition was the most complex and important stage of the sika deer life cycle. The present work represents the first multiple tissue transcriptome analysis of sika deer across four developmental stages. The generated data not only provide a functional genomics resource for future biological research on sika deer but also guide the selection and manipulation of genes controlling growth and development.

Bacterial community composition and fermentation patterns in the rumen of sika deer (Cervus nippon) fed three different diets.

Sika deer (Cervus nippon) rely on microorganisms living in the rumen to convert plant materials into chemical compounds, such as volatile fatty acids (VFAs), but how the rumen bacterial community is affected by different forages and adapt to altered diets remains poorly understood. The present study used 454-pyrosequencing of bacterial 16S ribosomal RNA (rRNA) genes to examine the relationship between rumen bacterial diversity and metabolic phenotypes using three sika deer in a 3 × 3 latin square design. Three sika deer were fed oak leaves (OL), corn stover (CS), or corn silage (CI), respectively. After a 7-day feeding period, when compared to the CS and CI groups, the OL group had a lower proportion of Prevotella spp. and a higher proportion of unclassified bacteria belonging to the families Succinivibrionaceae and Paraprevotellaceae (P<0.05). Meanwhile, the concentration of isobutyrate was significantly lower (P<0.05) in the OL group than in the CS and CI groups. There was no significant change of dominant bacterial genera in the OL group after 28 days of feeding. Conversely, total volatile fatty acids (TVFAs) showed an increase after 28 days of feeding, mainly due to the increasing of acetate, propionate, and valerate (P<0.05). The interplay between bacteria and metabolism in the OL group differed from that in the CS and CI groups, especially for the interaction of TVFAs and acetate/propionate. Overall, the current study suggested that Prevotella spp. played critical roles in the fermentation of feed in the rumen of sika deer. However, the differences in interplay patterns between rumen bacterial community composition and metabolic phenotypes were altered in the native and domesticated diets indicating the changed fermentation patterns in the rumen of sika deer.

Immunotoxicity of aluminum.

Aluminum (Al) is present in the daily life of all humans. With the incidence of Al contamination increased in recent years, the toxicity of Al on the immune function has attracted more attention. Even with this increased attention, the mechanism of Al immunotoxicity still remains unclear. The mechanism of Al immunotoxicity reviewed herein focused on the effects of Al on the splenic trace elements, the status of α-naphthyl acetate esterase (ANAE) cells, cytokines, complement and immunoglobulins, as well as macrophages. The studies in the literature showed that Al decreased splenic iron (Fe) and zinc (Zn) levels, but the effects of Al on splenic copper (Cu) level was ambiguous and controversial. Al exposure inhibited levels of ANAE(+) cells, the production of interleukin (IL)-2 and the functions of macrophages. With respect to other key cytokines, studies showed that Al suppressed the production of tumor necrosis factor (TNF)-α in vitro; effects of Al on TNF-α formation in vivo were less overt. Al exposure reduced complement 3 (C3) level, but effects of Al exposure on complement 4 (C4) level were not as clear-cut. Lastly, the effects of Al exposure on the IgG, IgM and IgA levels were conflicting. Taken in totality, the results of several studies in the literature demonstrated that Al could impart adverse effects on immune function.