First de novo whole genome sequencing and assembly of the bar-headed goose

Whole-genome resequencing reveals melanin deposition candidate genes of Luning chicken

BACKGROUND

Melanin in the black-bone chicken's body is considered the material basis for its medicinal effects and is an economically important trait. Therefore, improving the melanin content is a crucial focus in the breeding process of black-bone chickens. Luning chickens are black-bone chickens, with black beaks, skin, and meat. To investigate the genetic diversity and molecular mechanisms of melanin deposition in Luning chickens, we conducted whole-genome resequencing to analyze their breeding history and identify candidate genes influencing their black phenotype, along with transcriptome sequencing of dorsal skin tissues of male Luning chickens.

RESULTS

Population structure analysis revealed that Luning chickens tend to cluster independently and are closely related to Tibetan chickens. Runs of homozygosity analysis suggested potential inbreeding in the Luning chicken and Tibetan chicken population. By combining genetic differentiation index (Fst) and nucleotide diversity (θπ) ratios, we pinpointed selected regions associated with melanin deposition. Gene annotation identified 540 genes with the highest Fst value in LOC101750371 and LOC121108313, located on the 68.24-68.58 Mb interval of chromosome Z. Combining genomic and transcriptomic data, we identified ATP5E, EDN3, and LOC101750371 as candidate genes influencing skin color traits in black-bone chickens.

CONCLUSIONS

This study characterized the evolutionary history of Luning chickens and preliminarily excavated candidate genes influencing the genetic mechanism of pigmentation in black-bone chickens, providing valuable insights for the study of animal melanin deposition.

Chromosome-level genome sequencing and multi-omics of the Hungarian White Goose (Anser anser domesticus) reveals novel miRNA-mRNA regulation mechanism of waterfowl feather follicle development

The Hungarian White Goose (Anser anser domesticus) is an excellent European goose breed, with high feather and meat production. Despite its importance in the poultry industry, no available genome assembly information has been published. This study aimed to present Chromosome-level and functional genome sequencing of the Hungarian White Goose. The results showed that the genome assembly has a total length of 1115.82 Mb, 39 pairs of chromosomes, 92.98% of the BUSCO index, and contig N50 and scaffold N50 were up to 2.32 Mb and 60.69 Mb, respectively. Annotation of the genome assembly revealed 19550 genes, 286 miRNAs, etc. We identified 235 expanded and 1,167 contracted gene families in this breed compared with the other 16 species. We performed a positive selection analysis between this breed and four species of Anatidae to uncover the genetic information underlying feather follicle development. Further, we detected the function of miR-199-x, miR-143-y, and miR-23-z on goose embryonic skin fibroblast. In summary, we have successfully generated a highly complete genome sequence of the Hungarian white goose, which will provide a great resource to improve our understanding of gene functions and enhance the studies on feather follicle development at the genomic level.

Genetic characterization of Tibetan pigs adapted to high altitude under natural selection based on a large whole-genome dataset

The Qinghai-Tibet Plateau is a valuable genetic resource pool, and the high-altitude adaptation of Tibetan pigs is a classic example of the adaptive evolution of domestic animals. Here, we report the presence of Darwinian positive selection signatures in Tibetan pigs (TBPs) using 348 genome-wide datasets (127 whole-genome sequence datasets (WGSs) and 221 whole-genome single-nucleotide polymorphism (SNP) chip datasets). We characterized a high-confidence list of genetic signatures related response to high-altitude adaptation in Tibetan pigs, including 4,598 candidate SNPs and 131 candidate genes. Functional annotation and enrichment analysis revealed that 131 candidate genes are related to multiple systems and organs in Tibetan pigs. Notably, eight of the top ten novel genes, RALB, NBEA, LIFR, CLEC17A, PRIM2, CDH7, GK5 and FAM83B, were highlighted and associated with improved adaptive heart functions in Tibetan pigs high-altitude adaptation. Moreover, genome-wide association analysis revealed that 29 SNPs were involved in 13 candidate genes associated with at least one adaptive trait. In particular, among the top ten candidate genes, CLEC17A is related to a reduction in hemoglobin (HGB) in Tibetan pigs. Overall, our study provides a robust SNP/gene list involving genetic adaptation for Tibetan pig high-altitude adaptation, and it will be a valuable resource for future Tibetan pig studies.

A de novo assembled genome of the Tibetan Partridge (Perdix hodgsoniae) and its high-altitude adaptation

The Tibetan Partridge (Perdix hodgsoniae) is an endemic species distributed in high-altitude areas of 3600-5600 m on the Qinghai-Tibet Plateau. To explore how the species is adapted to the high elevation environment, we assembled a draft genome based on both the Illumina and PacBio sequencing platforms with its population genetics and genomics analysis. In total, 134.74 Gb short reads and 30.81 Gb long reads raw data were generated. The 1.05-Gb assembled genome had a contig N50 of 4.56 Mb, with 91.94% complete BUSCOs. The 17 457 genes were annotated, and 11.35% of the genome was composed of repeat sequences. The phylogenetic tree showed that P. hodgsoniae was located at the basal position of the clade, including Golden Pheasant (Chrysolophus pictus), Common Pheasant (Phasianus colchicus), and Mikado Pheasant (Syrmaticus mikado). We found that 1014, 2595, and 2732 of the 6641 one-to-one orthologous genes were under positive selection in P. hodgsoniae, detected using PAML, BUSTED, and aBSREL programs, respectively, of which 965 genes were common under positive selection with 3 different programs. Several positively selected genes and immunity pathways relevant to high-altitude adaptation were detected. Gene family evolution showed that 99 gene families experienced significant expansion events, while 6 gene families were under contraction. The total number of olfactory receptor genes was relatively low in P. hodgsoniae. Genomic data provide an important resource for a further study on the evolutionary history of P. hodgsoniae, which provides a new insight into its high-altitude adaptation mechanisms.

Chromosome-level genome assembly of the bar-headed goose (Anser indicus)

Bar-headed geese (Anser indicus) are adaptable to plateau environments. In this study, we sequenced and assembled a high-quality chromosome-level genome of the bar-headed goose using PacBio long reads and Hi-C technique, and generated 115.73 Gb of Illumina short-reads and 95.89 Gb of PacBio long-reads. The assembled bar-headed goose genome, with a contig N50 of 5.734 Mb and a scaffold N50 of 65.77 Mb, is 1.129 Gb in length and includes 33 chromosomes and 451 fragments. BUSCO assessment yielded a completeness score of 94.4%. In total, 15,376 protein-coding genes were predicted, of which 94.95% had homologs in protein databases. We identified 78 positively selected genes (PSGs) in the bar-headed goose genome, which were mainly enriched in calcium ion and ATP-binding. This bar-headed goose genome will be an important resource for increasing our understanding regarding the genetic basis of adaptation to life at a high altitude.

Chromosome-Level Genome Assembly of the Common Chaffinch (Aves: Fringilla coelebs): A Valuable Resource for Evolutionary Biology

The common chaffinch, Fringilla coelebs, is one of the most common, widespread, and well-studied passerines in Europe, with a broad distribution encompassing Western Europe and parts of Asia, North Africa, and the Macaronesian archipelagos. We present a high-quality genome assembly of the common chaffinch generated using Illumina shotgun sequencing in combination with Chicago and Hi-C libraries. The final genome is a 994.87-Mb chromosome-level assembly, with 98% of the sequence data located in chromosome scaffolds and a N50 statistic of 69.73 Mb. Our genome assembly shows high completeness, with a complete BUSCO score of 93.9% using the avian data set. Around 7.8% of the genome contains interspersed repetitive elements. The structural annotation yielded 17,703 genes, 86.5% of which have a functional annotation, including 7,827 complete universal single-copy orthologs out of 8,338 genes represented in the BUSCO avian data set. This new annotated genome assembly will be a valuable resource as a reference for comparative and population genomic analyses of passerine, avian, and vertebrate evolution.