IAnimal: a cross-species omics knowledgebase for animals

The role of NUDT3 in lipid accumulation and its functional variants related to backfat thickness in pigs

NUDT3 is a leading candidate gene that strongly linked to pig fatness traits, however, its function in porcine adipocytes remains poorly understood. Here, the percentage of EdU+ cells was significantly reduced when NUDT3 was knocked down, as was the expression of cell cycle repressors. NUDT3 overexpression yielded the opposite outcome. Moreover, the knockdown of NUDT3 resulted in more lipid droplets in adipocytes, whereas its enforced expression had the reverse effect. In addition, exogenous expression of NUDT3 in adipose tissue significantly reduced fat expansion triggered by a high-fat diet in mice. At molecular level, integrative RIP-seq and RNA-seq analysis revealed that genes influenced by NUDT3 overexpression or knockdown were significantly enriched in the PI3K-AKT signaling pathway, and western blot confirmed that AKT phosphorylation was significantly increased by NUDT3 knockdown, while the phosphorylation levels of PI3K, AKT, and mTOR were significantly decreased by the enforced NUDT3 expression both ex vivo and in vivo. Notably, rs694899689 was identified as a potential genetic variant for modulates NUDT3 expression and impacting backfat thickness in pigs through analysis of multi-omics data, CRISPRi (CRISPR interference) and dual luciferase reporter assays. Overall, our work established NUDT3 as a novel negative regulator of adipogenesis and lipid deposition and revealed that rs694899689 might serve as a potential molecular marker for pig breeding.

AnimalGWASAtlas: Annotation and prioritization of GWAS loci and quantitative trait loci for animal complex traits

Genome-wide association study (GWAS) and quantitative trait locus (QTL) mapping methods provide valuable insights and opportunities for identifying functional gene underlying phenotype formation. However, the majority of GWAS risk loci and QTLs located in noncoding regions poses significant challenges in pinpointing the protein-coding genes associated with specific traits. Moreover, growing evidence suggests not all GWAS risk loci and QTLs are functional, emphasizing the critical need for prioritizing causal sites-a task of paramount importance for biologists. The accumulation of publicly available multiomics data provides an unprecedented opportunity to annotate and prioritize GWAS risk loci and QTLs. Therefore, we developed a comprehensive multiomics database encompassing four major agricultural species-pig, sheep, cattle, and chicken. This database integrates publicly accessible datasets, including 140 GWAS studies (covering 471 traits), 2625 QTL datasets (spanning 1235 traits), 86 Hi-C datasets (from eight cells/tissue types), 95 epigenomic datasets (from four cells/tissue types), and 769 transcription factor motifs. The database aims to link GWAS-QTL loci located in the noncoding regions to the target genes they regulate and prioritize functional and causal regulatory elements. Ultimately, it provides a valuable resource and potential validation targets for elucidating the genes and molecular pathways underlying economically important traits in agricultural animals.

Genome sequencing and assembly of feral chickens in the wild of Sulawesi, Indonesia

The feralization of domestic chicken makes the conservation and management of red jungle fowl (Gallus gallus) more complicated and challenging. We collected two Sulawesi feral chickens, located east of the Wallace Line, for whole-genome sequencing and de novo genome assembly. Phylogenetic and f4-statistics analyses indicated that the Sulawesi feralized domestic chickens (G. g. domesticus) received gene flow from G. g. gallus. We integrated ~45× ultra-long Oxford Nanopore Technology reads and ~28× PacBio HiFi reads to generate a de novo genome assembly of a female Sulawesi feral chicken (GGsula) with a contig N50 of 19.88 Mbp. We characterized structural variations in GGsula, and found some were related to nervous system. Our study provides the first genome assembly of feral chickens, which is a unique genomic resource to explore the process of chicken domestication and feralization.

Experimental reporting of fish transcriptomic responses in environmental toxicology and ecotoxicology

Due to its increasing affordability and efforts to understand transcriptional responses of organisms to biotic and abiotic stimuli, transcriptomics has become an important tool with significant impact on toxicological investigations and hazard and risk assessments, especially during development and application of new approach methodologies (NAMs). Data generated using transcriptomic methodologies have directly informed adverse outcome pathway frameworks, chemical and biological read across, and aided in the identification of points of departure. Using data reporting frameworks for transcriptomics data offers improved transparency and reproducibility of research and an opportunity to identify barriers to adoption of these NAMs, especially in environmental toxicology and ecotoxicology with aquatic models. Improved reporting also allows for reexamination of existing data, limiting needs for experiment replication and further reducing animal experimentation. Here, we use a standardized form of data reporting for omics-based studies, the Organisation for Economic Co-operation and Development omics reporting framework, which specifically reports on a list of parameters that should be included in transcriptomics studies used in a regulatory context. We focused specifically on fish studies using RNA- Sequencing (Seq)/microarray technologies within a toxicology context. Inconsistencies in reporting and methodologies among the experimental designs (toxicology vs. molecular characterization) were observed in addition to foundational differences in reporting of sample concentration or preparation or quality assessments, which can affect reproducibility and read across, confidence in results, and contribute substantially to understanding molecular mechanisms of toxicants and toxins. Our findings present an opportunity for improved research reporting. We also provide several recommendations as logical steps to reduce barriers to adoption of transcriptomics within environmental toxicology and ecotoxicology.

A genome-wide association study identified candidate regions and genes for commercial traits in a Landrace population

Backfat thickness (BFT) and feed conversion ratio (FCR) are important commercial traits in the pig industry. With the increasing demand for human health and meat production, identifying functional genomic regions and genes associated with these commercial traits is critical for enhancing production efficiency. In this research, we conducted a genome-wide association study (GWAS) on a Landrace population comprising 4,295 individuals with chip data for BFT and FCR. Our analysis revealed a total of 118 genome-wide significant signals located on chromosomes SSC1, SSC2, SSC7, SSC12, and SSC13, respectively. Furthermore, we identified 10 potential regions associated with the two traits and annotated the genes within these regions. In addition, enrichment analysis was also performed. Notably, candidate genes such as SHANK2, KCNQ1, and ABL1 were found to be associated with BFT, whereas NAP1L4, LSP1, and PPFIA1 genes were related to the FCR. Our findings provide valuable insights into the genetic architecture of these two traits and offer guidance for future pig breeding efforts.

iDog: a multi-omics resource for canids study

iDog (https://ngdc.cncb.ac.cn/idog/) is a comprehensive public resource for domestic dogs (Canis lupus familiaris) and wild canids, designed to integrate multi-omics data and provide data services for the worldwide canine research community. Notably, iDog 2.0 features a 15-fold increase in genomic samples, including 29.55 million single nucleotide polymorphisms (SNPs) and 16.54 million insertions/deletions (InDels) from 1929 modern samples and 29.09 million SNPs from 111 ancient Canis samples. Additionally, 43487 breed-specific SNPs and 530 disease/trait-associated variants have been identified and integrated. The platform also includes data from 141 BioProjects involving gene expression analyses and a single-cell transcriptome module containing data from 105 057 Beagle hippocampus cells. iDog 2.0 also includes an epignome module that evaluates DNA methylation patterns across 547 samples and chromatin accessibility across 87 samples for the analysis of gene expression regulation. Additionally, it provies phenotypic data for 897 dog diseases, 3207 genotype-to-phenotype (G2P) pairs, and 349 dog disease-associated genes, along with two newly constructed ontologies for breed and disease standardization. Finally, 13 new analytical tools have been added. Given these enhancements, the updated iDog 2.0 is an invaluable resource for the global cannie research community.

The effects of runs-of-homozygosity on pig domestication and breeding

BACKGROUND

Since their domestication, recent inbreeding together with intensive artificial selection and population bottlenecks have allowed the prevalence of deleterious mutations and the increase of runs-of-homozygosity (ROH) in domestic pigs. This makes pigs a good model to understand the genetic underpinnings of inbreeding depression.

RESULTS

Here we integrated a comprehensive dataset comprising 7239 domesticated pigs and wild boars genotyped by single nucleotide polymorphism (SNP) chips, along with phenotypic data encompassing growth, reproduction and disease-associated traits. Our study revealed differential ROH landscapes during domestication and artificial selection of Eurasian pigs. We observed associations between ROH burden and phenotypic traits such as body conformation and susceptibility to diseases like scrotal hernia. By examining associations of whole-genome and regional ROH burden with gene expression, we identified specific genes and pathways affected by inbreeding depression. Associations of regional ROH burden with gene expression also enabled the discovery of novel regulatory elements. Lastly, we inferred recessive lethal mutations by examining depletion of ROH in an inbred population with relatively small sample size, following by fine mapping with sequencing data.

CONCLUSIONS

These findings suggested that both phenotypic and genetic variations have been reshaped by inbreeding, and provided insights to the genetic mechanisms underlying inbreeding depression.

Genomic prediction accounting for dominance and epistatic genetic effects on litter size traits in Large White pigs

Litter size traits of sows are crucial for the economic benefits of the pig industry. Three phenotypic traits of 1,206 Large White pigs, the total number born (TNB), number born alive (NBA), and number of healthy piglets (NHP), were recorded. We evaluated a series of genomic best linear unbiased prediction models that sequentially added additive effects (model A), dominance effects (model A + D), and epistatic effects (model A + D + AA, model A + D + AA + AD, and model A + D + AA + AD + DD) using chip data and imputed whole-genome sequencing (WGS) data to estimate genetic parameters and predictive accuracy. The reproductive traits of sows showed low heritability in this study, with narrow heritability of the 3 traits ranging from 0.030 to 0.064, and broad heritability ranging from 0.125 to 0.145. The inclusion of nonadditive effects in the model improved the accuracy of genomic selection. In the chip data, compared with that of the A model, the A + D + AA + AD + DD model showed the greatest increase in accuracy for TNB, NBA, and NHP, with improvements of 1.78%, 1.67%, and 1.74%, respectively. Additionally, the accuracy of the imputed WGS data was greater compared to the chip data. For the TNB, NBA, and NHP traits, the predictive accuracy of the imputed WGS data improved by 3.26%, 7.72%, and 3.00%, respectively, compared with that of the chip data. In summary, these results suggest that nonadditive effects in genomic selection could improve prediction accuracy and should be considered in pig genomic evaluation procedures.

The 1000 Chinese Indigenous Pig Genomes Project provides insights into the genomic architecture of pigs

Pigs play a central role in human livelihoods in China, but a lack of systematic large-scale whole-genome sequencing of Chinese domestic pigs has hindered genetic studies. Here, we present the 1000 Chinese Indigenous Pig Genomes Project sequencing dataset, comprising 1011 indigenous individuals from 50 pig populations covering approximately two-thirds of China's administrative divisions. Based on the deep sequencing (~25.95×) of these pigs, we identify 63.62 million genomic variants, and provide a population-specific reference panel to improve the imputation performance of Chinese domestic pig populations. Using a combination of methods, we detect an ancient admixture event related to a human immigration climax in the 13th century, which may have contributed to the formation of southeast-central Chinese pig populations. Analyzing the haplotypes of the Y chromosome shows that the indigenous populations residing around the Taihu Lake Basin exhibit a unique haplotype. Furthermore, we find a 13 kb region in the THSD7A gene that may relate to high-altitude adaptation, and a 0.47 Mb region on chromosome 7 that is significantly associated with body size traits. These results highlight the value of our genomic resource in facilitating genomic architecture and complex traits studies in pigs.

Epigenomic features associated with body temperature stabilize tissues during cold exposure in cold-resistant pigs

Cold stress in low-temperature environments can trigger changes in gene expression, but epigenomics regulation of temperature stability in vital tissues, including the fat and diencephalon, is still unclear. Here, we explore the cold-induced changes in epigenomic features in the diencephalon and fat tissues of two cold-resistant Chinese pig breeds, Min and Enshi black (ES) pigs, utilizing H3K27ac CUT&Tag, RNA-seq, and selective signature analysis. Our results show significant alterations in H3K27ac modifications in the diencephalon of Min pigs and the fat of ES pigs after cold exposure. Dramatic changes in H3K27ac modifications in the diencephalon of Min pig are primarily associated with genes involved in energy metabolism and hormone regulation, whereas those in the fat of ES pig are primarily associated with immunity-related genes. Moreover, transcription factors PRDM1 and HSF1, which show evidence of selection, are enriched in genomic regions presenting cold-responsive alterations in H3K27ac modification in the Min pig diencephalon and ES pig fat, respectively. Our results indicate the diversity of epigenomic response mechanisms to cold exposure between Min and ES pigs, providing unique epigenetic resources for studies of low-temperature adaptation in large mammals.

An overview of recent technological developments in bovine genomics

Cattle are regarded as highly valuable animals because of their milk, beef, dung, fur, and ability to draft. The scientific community has tried a number of strategies to improve the genetic makeup of bovine germplasm. To ensure higher returns for the dairy and beef industries, researchers face their greatest challenge in improving commercially important traits. One of the biggest developments in the last few decades in the creation of instruments for cattle genetic improvement is the discovery of the genome. Breeding livestock is being revolutionized by genomic selection made possible by the availability of medium- and high-density single nucleotide polymorphism (SNP) arrays coupled with sophisticated statistical techniques. It is becoming easier to access high-dimensional genomic data in cattle. Continuously declining genotyping costs and an increase in services that use genomic data to increase return on investment have both made a significant contribution to this. The field of genomics has come a long way thanks to groundbreaking discoveries such as radiation-hybrid mapping, in situ hybridization, synteny analysis, somatic cell genetics, cytogenetic maps, molecular markers, association studies for quantitative trait loci, high-throughput SNP genotyping, whole-genome shotgun sequencing to whole-genome mapping, and genome editing. These advancements have had a significant positive impact on the field of cattle genomics. This manuscript aimed to review recent advances in genomic technologies for cattle breeding and future prospects in this field.

Research Progress and Applications of Bovine Genome in the Tribe Bovini

Various bovine species have been domesticated and bred for thousands of years, and they provide adequate animal-derived products, including meat, milk, and leather, to meet human requirements. Despite the review studies on economic traits in cattle, the genetic basis of traits has only been partially explained by phenotype and pedigree breeding methods, due to the complexity of genomic regulation during animal development and growth. With the advent of next-generation sequencing technology, genomics projects, such as the 1000 Bull Genomes Project, Functional Annotation of Animal Genomes project, and Bovine Pangenome Consortium, have advanced bovine genomic research. These large-scale genomics projects gave us a comprehensive concept, technology, and public resources. In this review, we summarize the genomics research progress of the main bovine species during the past decade, including cattle (Bos taurus), yak (Bos grunniens), water buffalo (Bubalus bubalis), zebu (Bos indicus), and gayal (Bos frontalis). We mainly discuss the development of genome sequencing and functional annotation, focusing on how genomic analysis reveals genetic variation and its impact on phenotypes in several bovine species.

Comprehensive analysis of the gene expression profile of the male and female BTBR mice with diabetic nephropathy

Comprehensive insight into the gender-based gene expression-related omics data in a rodent model of diabetic nephropathy (DN) is scarce. In the present study, the gender-based genes regulating different pathways involved in the progression of DN were explored through an unbiased RNA sequence of kidneys from BTBR mice with DN. We identified 17,739 and 17,981 genes in male and female DN mice; 1121 and 655 genes were expressed differentially (DEGs, differentially expressed genes) in male and female DN mice; both genders displayed only 195 DEGs. In the male DN mice, the number of upregulated genes was nearly the same as that of the down-regulated genes. In contrast, the number of upregulated genes was lesser than that of the down-regulated genes in the female DN mice, manifesting a remarkable gender disparity during the progression of DN in this animal model. Gene Ontology (GO) and KEGG-enriched results showed that most of these DEGs were related to the critical biological processes, including metabolic pathways, natural oxidation, bile secretion, and PPAR signaling; all are highly associated with DN. Notably, the DEGs significantly enriched for steroid hormone biosynthesis pathway were identified in both genders; the number of DEGs increased was 22 in male DN mice and 14 in female DN mice. Specifically, the Ugt1a10, Akr1c12, and Akr1c14 were upregulated in both genders. Interestingly, the Hsd11b1 gene was upregulated in female DN mice but downregulated in male DN mice. These results suggest that a significant gender-based variance in the gene expression occurs during the progression of DN and may be playing a role in the advancement of DN in the BTBR mouse model.

TCOD: an integrated resource for tropical crops

Tropical crops are vital for tropical agriculture, with resource scarcity, functional diversity and extensive market demand, providing considerable economic benefits for the world's tropical agriculture-producing countries. The rapid development of sequencing technology has promoted a milestone in tropical crop research, resulting in the generation of massive amount of data, which urgently needs an effective platform for data integration and sharing. However, the existing databases cannot fully satisfy researchers' requirements due to the relatively limited integration level and untimely update. Here, we present the Tropical Crop Omics Database (TCOD, https://ngdc.cncb.ac.cn/tcod), a comprehensive multi-omics data platform for tropical crops. TCOD integrates diverse omics data from 15 species, encompassing 34 chromosome-level de novo assemblies, 1 255 004 genes with functional annotations, 282 436 992 unique variants from 2048 WGS samples, 88 transcriptomic profiles from 1997 RNA-Seq samples and 13 381 germplasm items. Additionally, TCOD not only employs genes as a bridge to interconnect multi-omics data, enabling cross-species comparisons based on homology relationships, but also offers user-friendly online tools for efficient data mining and visualization. In short, TCOD integrates multi-species, multi-omics data and online tools, which will facilitate the research on genomic selective breeding and trait biology of tropical crops.

AGIDB: a versatile database for genotype imputation and variant decoding across species

The high cost of large-scale, high-coverage whole-genome sequencing has limited its application in genomics and genetics research. The common approach has been to impute whole-genome sequence variants obtained from a few individuals for a larger population of interest individually genotyped using SNP chip. An alternative involves low-coverage whole-genome sequencing (lcWGS) of all individuals in the larger population, followed by imputation to sequence resolution. To overcome limitations of processing lcWGS data and meeting specific genotype imputation requirements, we developed AGIDB (https://agidb.pro), a website comprising tools and database with an unprecedented sample size and comprehensive variant decoding for animals. AGIDB integrates whole-genome sequencing and chip data from 17 360 and 174 945 individuals, respectively, across 89 species to identify over one billion variants, totaling a massive 688.57 TB of processed data. AGIDB focuses on integrating multiple genotype imputation scenarios. It also provides user-friendly searching and data analysis modules that enable comprehensive annotation of genetic variants for specific populations. To meet a wide range of research requirements, AGIDB offers downloadable reference panels for each species in addition to its extensive dataset, variant decoding and utility tools. We hope that AGIDB will become a key foundational resource in genetics and breeding, providing robust support to researchers.

Epigenomic analysis of the myometrium during late implantation revealed regulatory elements in genes related to the cellular zinc homeostasis pathway in pigs

The myometrium, composed of the inner circular muscle (CM) and outer longitudinal muscle (LM), is crucial in establishing and maintaining early pregnancy. However, the molecular mechanisms involved are not well understood. In this study, we identified the transcriptomic features of the CM and LM collected from the mesometrial (M) and anti-mesometrial (AM) sides of the pig uterus on day 18 of pregnancy during the placentation initiation phase. Some genes in the cellular zinc ion level regulatory pathways (MT-1A, MT-1D, MT-2B, SLC30A2, and SLC39A2) were spatially and highly enriched in uterine CM at the mesometrial side. In addition, the histone modification profiles of H3K27ac and H3K4me3 in uterine CM and LM collected from the mesometrial side were characterized. Genomic regions associated with the expression of genes regulating the cellular zinc ion level were detected. Moreover, six highly linked variants in the H3K27ac-enriched region of the pig SLC30A2 gene were identified and found to be significantly associated with the total number born at the second parity (P < 0.05). In conclusion, the genes in the pathways of cellular zinc homeostasis and their regulatory elements identified have implications for pig reproduction trait improvement and warrant further investigations.