BackgroundThe forest musk deer, Moschus berezovskii, is one of seven musk deer (Moschus spp.) and is distributed in Southwest China. Akin to other musk deer, the forest musk deer has been traditionally and is currently hunted for its musk (i.e., global perfume industry). Considerable hunting pressure and habitat loss have caused significant population declines. Consequently, the Chinese government commenced captive breeding programs for musk harvesting in the 1950s. However, the prevalence of fatal diseases is considerably restricting population increases. Disease severity and extent are exacerbated by inbreeding and genetic diversity declines in captive musk deer populations. It is essential that knowledge of captive and wild forest musk deer populations' immune system and genome be gained in order to improve their physical and genetic health. We have thus sequenced the whole genome of the forest musk deer, completed the genomic assembly and annotation, and performed preliminary bioinformatic analyses.FindingsA total of 407 Gb raw reads from whole-genome sequencing were generated using the Illumina HiSeq 4000 platform. The final genome assembly is around 2.72 Gb, with a contig N50 length of 22.6 kb and a scaffold N50 length of 2.85 Mb. We identified 24,352 genes and found that 42.05% of the genome is composed of repetitive elements. We also detected 1,236 olfactory receptor genes. The genome-wide phylogenetic tree indicated that the forest musk deer was within the order Artiodactyla, and it appeared as the sister clade of four members of Bovidae. In total, 576 genes were under positive selection in the forest musk deer lineage.ConclusionsWe provide the first genome sequence and gene annotation for the forest musk deer. The availability of these resources will be very useful for the conservation and captive breeding of this endangered and economically important species and for reconstructing the evolutionary history of the order Artiodactyla.

The Chinese forest musk deer ( Moschus berezovskii ; FMD) is an artiodactyl mammal and is both economically valuable and highly endangered. To investigate the genetic mechanisms of musk secretion and adaptive immunity in FMD, we compared its genome to nine other artiodactyl genomes. Comparative genomics demonstrated that eight positively selected genes (PSGs) in FMD were annotated in three KEGG pathways that were related to metabolic and synthetic activity of musk, similar to previous transcriptome studies. Functional enrichment analysis indicated that many PSGs were involved in the regulation of immune system processes, implying important reorganization of the immune system in FMD. FMD-specific missense mutations were found in two PSGs ( MHC class II antigen DRA and ADA ) that were classified as deleterious by PolyPhen-2, possibly contributing to immune adaptation to infectious diseases. Functional assessment showed that the FMD-specific mutation enhanced the ADA activity, which was likely to strengthen the immune defense against pathogenic invasion. Single nucleotide polymorphism-based inference showed the recent demographic trajectory for FMD. Our data and findings provide valuable genomic resources not only for studying the genetic mechanisms of musk secretion and adaptive immunity, but also for facilitating more effective management of the captive breeding programs for this endangered species.

Raptors are carnivorous birds including accipitrids (Accipitridae, Accipitriformes) and owls (Strigiformes), which are diurnal and nocturnal, respectively. To examine the evolutionary basis of adaptations to different light cycles and hunting behavior between accipitrids and owls, we de novo assembled besra (Accipiter virgatus, Accipitridae, Accipitriformes) and oriental scops owl (Otus sunia, Strigidae, Strigiformes) draft genomes. Comparative genomics demonstrated four PSGs (positively selected genes) (XRCC5, PRIMPOL, MDM2, and SIRT1) related to the response to ultraviolet (UV) radiation in accipitrids, and one PSG (ALCAM) associated with retina development in owls, which was consistent with their respective diurnal/nocturnal predatory lifestyles. We identified five accipitrid-specific and two owl-specific missense mutations and most of which were predicted to affect the protein function by PolyPhen-2. Genome comparison showed the diversification of raptor olfactory receptor repertoires, which may reflect an important role of olfaction in their predatory lifestyle. Comparison of TAS2R gene (i.e. linked to tasting bitterness) number in birds with different dietary lifestyles suggested that dietary toxins were a major selective force shaping the diversity of TAS2R repertoires. Fewer TAS2R genes in raptors reflected their carnivorous diet, since animal tissues are less likely to contain toxins than plant material. Our data and findings provide valuable genomic resources for studying the genetic mechanisms of raptors’ environmental adaptation, particularly olfaction, nocturnality and response to UV radiation.

No abstract

Since the release of the complete human genome, the priority of human genomic study has now been shifting towards closing gaps in ethnic diversity. Here, we present a fully phased and well-annotated diploid human genome from a Han Chinese male individual (CN1), in which the assemblies of both haploids achieve the telomere-to-telomere (T2T) level. Comparison of this diploid genome with the CHM13 haploid T2T genome revealed significant variations in the centromere. Outside the centromere, we discovered 11,413 structural variations, including numerous novel ones. We also detected thousands of CN1 alleles that have accumulated high substitution rates and a few that have been under positive selection in the East Asian population. Further, we found that CN1 outperforms CHM13 as a reference genome in mapping and variant calling for the East Asian population owing to the distinct structural variants of the two references. Comparison of SNP calling for a large cohort of 8869 Chinese genomes using CN1 and CHM13 as reference respectively showed that the reference bias profoundly impacts rare SNP calling, with nearly 2 million rare SNPs miss-called with different reference genomes. Finally, applying the CN1 as a reference, we discovered 5.80 Mb and 4.21 Mb putative introgression sequences from Neanderthal and Denisovan, respectively, including many East Asian specific ones undetected using CHM13 as the reference. Our analyses reveal the advances of using CN1 as a reference for population genomic studies and paleo-genomic studies. This complete genome will serve as an alternative reference for future genomic studies on the East Asian population.

The monal genus (Lophophorus) is a branch of Phasianidae and its species inhabit the high-altitude mountains of the Qinghai-Tibet Plateau. The Chinese monal, L. lhuysii, is a threatened endemic bird of China that possesses high-altitude adaptability, diversity of plumage color and potentially low reproductive life history. This is the first study to describe the monal genome using next generation sequencing technology. The Chinese monal genome size is 1.01 Gb, with 16,940 protein-coding genes. Gene annotation yielded 100.93 Mb (9.97%) repeat elements, 785 ncRNA, 5,465,549 bp (0.54%) SSR and 15,550 (92%) genes in public databases. Compared to other birds and mammals, the genome evolution analysis showed numerous expanded gene families and positive selected genes involved in high-altitude adaptation, especially related to the adaptation of low temperature and hypoxia. Consequently, this gene data can be used to investigate the molecular evolution of high-altitude adaptation in future bird research. Our first published genome of the genus Lophophorus will be integral for the study of monal population genetic diversity and conservation, genomic evolution and Galliformes species differentiation in the Qinghai-Tibetan Plateau.

All the genes of Paradoxornis gularis and Niltava davidi complete mitochondrial genomes encoded on the H-strand, apart from one PCG (ND6) and eight tRNAs (tRNA Gln , tRNA Ala , tRNA Asn , tRNA Cys , tRNA Tyr , tRNA Ser(UCN) , tRNA Pro , and tRNA Glu). Strikingly large tandem repeats were found in the control region (CR) of P. gularis mitogenome. The phylogenetic analysis demonstrated that P. gularis possessed a basal phylogenetic position within the genus Paradoxornis, and N. davidi was more closely related to Cyanoptila cyanomelana than other Muscicapidae species.