Genetic Basis of Phenotypic Differences Between Chinese Yunling Black Goats and Nubian Goats Revealed by Allele-Specific Expression in Their F1 Hybrids

Telomere-to-telomere genome assembly of a male goat reveals variants associated with cashmere traits

A complete goat (Capra hircus) reference genome enhances analyses of genetic variation, thus providing insights into domestication and selection in goats and related species. Here, we assemble a telomere-to-telomere (T2T) gap-free genome (2.86 Gb) from a cashmere goat (T2T-goat1.0), including a Y chromosome of 20.96 Mb. With a base accuracy of >99.999%, T2T-goat1.0 corrects numerous genome-wide structural and base errors in previous assemblies and adds 288.5 Mb of previously unresolved regions and 446 newly assembled genes to the reference genome. We sequence the genomes of five representative goat breeds for PacBio reads, and use T2T-goat1.0 as a reference to identify a total of 63,417 structural variations (SVs) with up to 4711 (7.42%) in the previously unresolved regions. T2T-goat1.0 was applied in population analyses of global wild and domestic goats, which revealed 32,419 SVs and 25,397,794 SNPs, including 870 SVs and 545,026 SNPs in the previously unresolved regions. Also, our analyses reveal a set of selective variants and genes associated with domestication (e.g., NKG2D and ABCC4) and cashmere traits (e.g., ABCC4 and ASIP).

Comprehensive analyses of 435 goat transcriptomes provides insight into male reproduction

A systematic analysis of genes related to reproduction is crucial for obtaining a comprehensive understanding of the molecular mechanisms that underlie male reproductive traits in mammals. Here, we utilized 435 goat transcriptome datasets to unveil the testicular tissue-specific genes (TSGs), allele-specific expression (ASE) genes and their uncharacterized transcriptional features related to male goat reproduction. Results showed a total of 1790 TSGs were identified in goat testis, which was the most among all tissues. GO enrichment analyses suggested that testicular TSGs were mainly involved in spermatogenesis, multicellular organism development, spermatid development, and flagellated sperm motility. Subsequently, a total of 95 highly conserved TSGs (HCTSGs), 508 middle conserved TSGs (MCTSGs) and 42 no conserved TSGs (NCTSGs) were identified in goat testis. GO enrichment analyses suggested that the HCTSGs and MCTSGs has a more important association with male reproduction than NCTSGs. Additionally, we identified 644 ASE genes, including 88 tissue-specific ASE (TS-ASE) genes (e.g., FSIP2, TDRD9). GO enrichment analyses indicated that both ASE genes and TS-ASE genes were associated with goat male reproduction. Overall, this study revealed an extensive gene set involved in the regulation of male goat reproduction and their dynamic transcription patterns. Data reported here provide valuable insights for a further improvement of the economic benefits of goats as well as future treatments for male infertility.

Diverse WGBS profiles of longissimus dorsi muscle in Hainan black goats and hybrid goats

BACKGROUND

Goat products have played a crucial role in meeting the dietary demands of people since the Neolithic era, giving rise to a multitude of goat breeds globally with varying characteristics and meat qualities. The primary objective of this study is to pinpoint the pivotal genes and their functions responsible for regulating muscle fiber growth in the longissimus dorsi muscle (LDM) through DNA methylation modifications in Hainan black goats and hybrid goats.

METHODS

Whole-genome bisulfite sequencing (WGBS) was employed to scrutinize the impact of methylation on LDM growth. This was accomplished by comparing methylation differences, gene expression, and their associations with growth-related traits.

RESULTS

In this study, we identified a total of 3,269 genes from differentially methylated regions (DMR), and detected 189 differentially expressed genes (DEGs) through RNA-seq analysis. Hypo DMR genes were primarily enriched in KEGG terms associated with muscle development, such as MAPK and PI3K-Akt signaling pathways. We selected 11 hub genes from the network that intersected the gene sets within DMR and DEGs, and nine genes exhibited significant correlation with one or more of the three LDM growth traits, namely area, height, and weight of loin eye muscle. Particularly, PRKG1 demonstrated a negative correlation with all three traits. The top five most crucial genes played vital roles in muscle fiber growth: FOXO3 safeguarded the myofiber's immune environment, FOXO6 was involved in myotube development and differentiation, and PRKG1 facilitated vasodilatation to release more glucose. This, in turn, accelerated the transfer of glucose from blood vessels to myofibers, regulated by ADCY5 and AKT2, ultimately ensuring glycogen storage and energy provision in muscle fibers.

CONCLUSION

This study delved into the diverse methylation modifications affecting critical genes, which collectively contribute to the maintenance of glycogen storage around myofibers, ultimately supporting muscle fiber growth.

Comparative and Functional Analysis of miRNAs and mRNAs Involved in Muscle Fiber Hypertrophy of Juvenile and Adult Goats

MicroRNAs (miRNAs) are small non-coding RNAs that post-transcriptionally regulate several pathway intermediates and affect the skeletal muscle development in mice, pigs, sheep, and cattle. However, to date, only a small number of miRNAs have been reported in the muscle development of goats. In this report, the longissimus dorsi transcripts of one- and ten-month-old goats were analyzed by sequencing RNAs and miRNAs. The results showed that the ten-month-old Longlin goats had 327 up- and 419 down-regulated differentially expressed genes (DEGs) compared with the one-month-old. In addition, 20 co-up-regulated and 55 co-down-regulated miRNAs involved in the muscle fiber hypertrophy of goats were identified in ten-month-old Longlin and Nubian goats compared with one-month-old. Five miRNA-mRNA pairs (chi-let-7b-3p-MIRLET7A, chi-miR193b-3p-MMP14, chi-miR-355-5p-DGAT2, novel_128-LOC102178119, novel_140-SOD3) involved in the goat skeletal muscle development were identified by miRNA-mRNA negative correlation network analysis. Our results provided new insight into the functional roles of goat muscle-associated miRNAs, allowing a deeper understanding of the transformation of miRNA roles during mammalian muscle development.

Genome-Wide Population Structure and Selection Signatures of Yunling Goat Based on RAD-seq

Simple Summary

Goats are important domestic animals that provide meat, milk, fur, and other products for humans. The demand for these products has increased in recent years. Disease resistance among goat breeds is different, but the genetic basis of the differences in resistance to diseases is still unclear and needs to be further studied. In this study, many genes and pathways related to immunity and diseases were identified to be under positive selection between Yunling and Nubian goats using RAD-seq technology. This study on the selection signatures of Yunling goats provides the scientific basis and technical support for the breeding of domestic goats for disease resistance, which has important social and economic significance.

Abstract

Animal diseases impose a huge burden on the countries where diseases are endemic. Conventional control strategies of vaccines and veterinary drugs are to control diseases from a pharmaceutical perspective. Another alternative approach is using pre-existing genetic disease resistance or tolerance. We know that the Yunling goat is an excellent local breed from Yunnan, southwestern China, which has characteristics of strong disease resistance and remarkable adaptability. However, genetic information about the selection signatures of Yunling goats is limited. We reasoned that the genes underlying the observed difference in disease resistance might be identified by investigating selection signatures between two different goat breeds. Herein, we selected the Nubian goat as the reference group to perform the population structure and selection signature analysis by using RAD-seq technology. The results showed that two goat breeds were divided into two clusters, but there also existed gene flow. We used Fst (F-statistics) and π (pi/θπ) methods to carry out selection signature analysis. Eight selected regions and 91 candidate genes were identified, in which some genes such as DOK2, TIMM17A, MAVS, and DOCK8 related to disease and immunity and some genes such as SPEFI, CDC25B, and MIR103 were associated with reproduction. Four GO (Gene Ontology) terms (GO:0010591, GO:001601, GO:0038023, and GO:0017166) were associated with cell migration, signal transduction, and immune responses. The KEGG (Kyoto Encyclopedia of Genes and Genomes) signaling pathways were mainly associated with immune responses, inflammatory responses, and stress reactions. This study preliminarily revealed the genetic basis of strong disease resistance and adaptability of Yunling goats. It provides a theoretical basis for the subsequent genetic breeding of disease resistance of goats.

Genetic Basis of Dorper Sheep (Ovis aries) Revealed by Long-Read De Novo Genome Assembly

Dorper sheep (Ovis aries) (DPS), developed in the 1930s by crossing Dorset Horn and Blackhead Persian sheep in South Africa, is a world-famous composite breed for mutton production. The genetic basis underlying this breed is yet to be elucidated. Here, we report the sequencing and assembly of a highly contiguous Dorper sheep genome via integration of Oxford Nanopore Technology (ONT) sequencing and Hi-C (chromatin conformation capture) approaches. The assembled genome was around 2.64 Gb with a contig N50 of 73.33 Mb and 140 contigs in total. More than 99.5% of the assembled sequences could be anchored to 27 chromosomes and they were annotated with 20,450 protein-coding genes. Allele-specific expression (ASE) genes of Dorper sheep were revealed through ASE analysis and they were involved in the immune system, lipid metabolism, and environmental adaptation. A total of 5,701 and 456 allelic sites were observed in the SNP and indels loci identified from relevant whole-genome resequencing data. These allelic SNP and INDEL sites were annotated in 1,002 and 294 genes, respectively. Moreover, we calculated the number of variant sites and related genes derived from the maternal and paternal ancestors, revealing the genetic basis of outstanding phenotypic performance of Dorper sheep. In conclusion, this study reports the first reference genome of Dorper sheep and reveals its genetic basis through ASE. This study also provides a pipeline for mining genetic information of composite breeds, which has an implication for future hybrid-breeding practices.

Patterns of Gene Content and Co-occurrence Constrain the Evolutionary Path toward Animal Association in Candidate Phyla Radiation Bacteria

Candidate Phyla Radiation (CPR) bacteria are small, likely episymbiotic organisms found across Earth's ecosystems. Despite their prevalence, the distribution of CPR lineages across habitats and the genomic signatures of transitions among these habitats remain unclear. Here, we expand the genome inventory for Absconditabacteria (SR1), Gracilibacteria, and Saccharibacteria (TM7), CPR bacteria known to occur in both animal-associated and environmental microbiomes, and investigate variation in gene content with habitat of origin. By overlaying phylogeny with habitat information, we show that bacteria from these three lineages have undergone multiple transitions from environmental habitats into animal microbiomes. Based on co-occurrence analyses of hundreds of metagenomes, we extend the prior suggestion that certain Saccharibacteria have broad bacterial host ranges and constrain possible host relationships for Absconditabacteria and Gracilibacteria. Full-proteome analyses show that animal-associated Saccharibacteria have smaller gene repertoires than their environmental counterparts and are enriched in numerous protein families, including those likely functioning in amino acid metabolism, phage defense, and detoxification of peroxide. In contrast, some freshwater Saccharibacteria encode a putative rhodopsin. For protein families exhibiting the clearest patterns of differential habitat distribution, we compared protein and species phylogenies to estimate the incidence of lateral gene transfer and genomic loss occurring over the species tree. These analyses suggest that habitat transitions were likely not accompanied by large transfer or loss events but rather were associated with continuous proteome remodeling. Thus, we speculate that CPR habitat transitions were driven largely by availability of suitable host taxa and were reinforced by acquisition and loss of some capacities. IMPORTANCE Studying the genetic differences between related microorganisms from different environment types can indicate factors associated with their movement among habitats. This is particularly interesting for bacteria from the Candidate Phyla Radiation because their minimal metabolic capabilities require associations with microbial hosts. We found that shifts of Absconditabacteria, Gracilibacteria, and Saccharibacteria between environmental ecosystems and mammalian mouths/guts probably did not involve major episodes of gene gain and loss; rather, gradual genomic change likely followed habitat migration. The results inform our understanding of how little-known microorganisms establish in the human microbiota where they may ultimately impact health.

Elimination of Reference Mapping Bias Reveals Robust Immune Related Allele-Specific Expression in Crossbred Sheep

Pervasive allelic variation at both gene and single nucleotide level (SNV) between individuals is commonly associated with complex traits in humans and animals. Allele-specific expression (ASE) analysis, using RNA-Seq, can provide a detailed annotation of allelic imbalance and infer the existence of cis-acting transcriptional regulation. However, variant detection in RNA-Seq data is compromised by biased mapping of reads to the reference DNA sequence. In this manuscript, we describe an unbiased standardized computational pipeline for allele-specific expression analysis using RNA-Seq data, which we have adapted and developed using tools available under open license. The analysis pipeline we present is designed to minimize reference bias while providing accurate profiling of allele-specific expression across tissues and cell types. Using this methodology, we were able to profile pervasive allelic imbalance across tissues and cell types, at both the gene and SNV level, in Texel×Scottish Blackface sheep, using the sheep gene expression atlas data set. ASE profiles were pervasive in each sheep and across all tissue types investigated. However, ASE profiles shared across tissues were limited, and instead, they tended to be highly tissue-specific. These tissue-specific ASE profiles may underlie the expression of economically important traits and could be utilized as weighted SNVs, for example, to improve the accuracy of genomic selection in breeding programs for sheep. An additional benefit of the pipeline is that it does not require parental genotypes and can therefore be applied to other RNA-Seq data sets for livestock, including those available on the Functional Annotation of Animal Genomes (FAANG) data portal. This study is the first global characterization of moderate to extreme ASE in tissues and cell types from sheep. We have applied a robust methodology for ASE profiling to provide both a novel analysis of the multi-dimensional sheep gene expression atlas data set and a foundation for identifying the regulatory and expressed elements of the genome that are driving complex traits in livestock.

Discovery of Genomic Characteristics and Selection Signatures in Korean Indigenous Goats Through Comparison of 10 Goat Breeds

Indigenous breeds develop their own genomic characteristics by adapting to local environments or cultures over long periods of time. Most of them are not particularly productive in commercial terms, but they have abilities to survive in harsh environments or tolerate to specific diseases. Their adaptive characteristics play an important role as genetic materials for improving commercial breeds. As a step toward this goal, we analyzed the genome of Korean indigenous goats within 10 goat breeds. We collected 136 goat individuals by sequencing 46 new goats and employing 90 publicly available goats. Our whole-genome data was comprised of three indigenous breeds (Korean indigenous goat, Iranian indigenous goat, and Moroccan indigenous goat; n = 29, 18, 20), six commercial breeds (Saanen, Boer, Anglo-Nubian, British Alpine, Alpine, and Korean crossbred; n = 16, 11, 5, 5, 2, 13), and their ancestral species (Capra aegagrus; n = 17). We identified that the Iranian indigenous goat and the Moroccan indigenous goat have relatively similar genomic characteristics within a large category of genomic diversity but found that the Korean indigenous goat has unique genomic characteristics distinguished from the other nine breeds. Through population analysis, we confirmed that these characteristics have resulted from a near-isolated environment with strong genetic drift. The Korean indigenous goat experienced a severe genetic bottleneck upon entering the Korean Peninsula about 2,000 years ago, and has subsequently rarely experienced genetic interactions with other goat breeds. From selection analysis and gene-set enrichment analysis, we revealed selection signals for Salmonella infection and cardiomyopathy in the genome of the Korean indigenous goat. These adaptive characteristics were further identified with genomic-based evidence. We uncovered genomic regions of selective sweeps in the LBP and BPI genes (Salmonella infection) and the TTN and ITGB6 genes (cardiomyopathy), among several candidate genes. Our research presents unique genomic characteristics and distinctive selection signals of the Korean indigenous goat based on the extensive comparison. Although the adaptive traits require further validation through biological experiments, our findings are expected to provide a direction for future biodiversity conservation strategies and to contribute another option to genomic-based breeding programmes for improving the viability of Capra hircus.