Transcriptional atlas analysis from multiple tissues reveals the expression specificity patterns in beef cattle

Genetic parameters and genome-wide association studies including the X chromosome for various reproduction and semen quality traits in Nellore cattle

BACKGROUND

The profitability of the beef industry is directly influenced by the fertility rate and reproductive performance of both males and females, which can be improved through selective breeding. When performing genomic analyses, genetic markers located on the X chromosome have been commonly ignored despite the X chromosome being one of the largest chromosomes in the cattle genome. Therefore, the primary objectives of this study were to: (1) estimate variance components and genetic parameters for eighteen male and five female fertility and reproductive traits in Nellore cattle including X chromosome markers in the analyses; and (2) perform genome-wide association studies and functional genomic analyses to better understand the genetic background of male and female fertility and reproductive performance traits in Nellore cattle.

RESULTS

The percentage of the total direct heritability (h2total) explained by the X chromosome markers (h2x) ranged from 3 to 32% (average: 16.4%) and from 9 to 67% (average: 25.61%) for female reproductive performance and male fertility traits, respectively. Among the traits related to breeding soundness evaluation, the overall bull and semen evaluation and semen quality traits accounted for the highest proportion of h2x relative to h2total with an average of 39.5% and 38.75%, respectively. The total number of significant genomic markers per trait ranged from 7 (seminal vesicle width) to 43 (total major defects). The number of significant markers located on the X chromosome ranged from zero to five. A total of 683, 252, 694, 382, 61, and 77 genes overlapped with the genomic regions identified for traits related to female reproductive performance, semen quality, semen morphology, semen defects, overall bulls' fertility evaluation, and overall semen evaluation traits, respectively. The key candidate genes located on the X chromosome are PRR32, STK26, TMSB4X, TLR7, PRPS2, SMS, SMARCA1, UTP14A, and BCORL1. The main gene ontology terms identified are "Oocyte Meiosis", "Progesterone Mediated Oocyte Maturation", "Thermogenesis", "Sperm Flagellum", and "Innate Immune Response".

CONCLUSIONS

Our findings indicate the key role of genes located on the X chromosome on the phenotypic variability of male and female reproduction and fertility traits in Nellore cattle. Breeding programs aiming to improve these traits should consider adding the information from X chromosome markers in their genomic analyses.

Transcriptionally distinct B cell profiles in systemic immune tissues and peritoneal cavity of Atlantic salmon (Salmo salar) infected with salmonid alphavirus subtype 3

Teleost B cells producing neutralizing antibodies contribute to protection against salmonid alphavirus (SAV) infection, the etiological agent of pancreas disease, thereby reducing mortality and disease severity. Our previous studies show differences in B cell responses between the systemic immune tissues (head kidney (HK) and spleen) and the peritoneal cavity (PerC) after intraperitoneal SAV3 infection in Atlantic salmon (Salmo salar) where the response in PerC dominates at the late time points. By employing the same infection model, we aimed to further characterize these B cells. Immunophenotyping of teleost B cells is challenging due to limited availability of markers; however, RNA-seq opens an opportunity to explore differences in transcriptomic responses of these cells. Our analysis identified 334, 259 and 613 differentially expressed genes (DEGs) in Atlantic salmon IgM+IgD+ B cells from HK, spleen, and PerC, respectively, at 6 weeks post SAV3 infection. Of these, only 34 were common to all the three immune sites. Additionally, out of the top 100 genes with the highest fold change in expression, only four genes were common across B cells from the three sites. Functional enrichment analyses of DEGs using KEGG and GO databases demonstrated differences in enriched innate immune signaling and the cytokine-cytokine interaction pathways in B cells across the sites, with varying numbers of genes involved. Overall, these findings show the presence of transcriptionally distinct B cell subsets with innate immune functions in HK, spleen and PerC of SAV3-infected Atlantic salmon. Further, our data provide new insights into the immunoregulatory role of fish B cells through the differential expression of various cytokine ligands and receptors and will be a useful resource for further studies into B cell immune compartments.

Expression Profiles of Housekeeping Genes and Tissue-Specific Genes in Different Tissues of Chinese Sturgeon (Acipenser sinensis)

The Chinese sturgeon (Acipenser sinensis) is an ancient, complex autooctoploid fish species that is currently facing conservation challenges throughout its distribution. To comprehensively characterize the expression profiles of genes and their associated biological functions across different tissues, we performed a transcriptome-scale gene expression analysis, focusing on housekeeping genes (HKGs), tissue-specific genes (TSGs), and co-expressed gene modules in various tissues. We collected eleven tissues to establish a transcriptomic repository, including data from Pacific Biosciences isoform sequencing (PacBio Iso-seq) and RNA sequencing (RNA-seq), and then obtained 25,434 full-length transcripts, with lengths from 307 to 9515 bp and an N50 of 3195 bp. Additionally, 20,887 transcripts were effectively identified and classified as known homologous genes. We also identified 787 HKGs, and the number of TSGs varied from 25 in the liver to 2073 in the brain. TSG functions were mainly enriched in certain signaling pathways involved in specific physiological processes, such as voltage-gated potassium channel activity, nervous system development, glial cell differentiation in the brain, and leukocyte transendothelial migration in the spleen and pronephros. Meanwhile, HKGs were highly enriched in some pathways involved in ribosome biogenesis, proteasome core complex, spliceosome activation, elongation factor activity, and translation initiation factor activity, which have been strongly implicated in fundamental biological tissue functions. We also predicted five modules, with eight hub genes in the brown module, most of which (such as rps3a, rps7, rps23, rpl11, rpl17, rpl27, and rpl28) were linked to ribosome biogenesis. Our results offer insights into ribosomal proteins that are indispensable in ribosome biogenesis and protein synthesis, which are crucial in various cell developmental processes and neural development of Chinese sturgeon. Overall, these findings will not only advance the understanding of fundamental biological functions in Chinese sturgeon but also supply a valuable genetic resource for characterizing this extremely important species.

Gene coexpression network analysis reveals the genes and pathways in pectoralis major muscle and liver associated with wooden breast in broilers

Wooden breast (WB) is a myopathy mainly affecting pectoralis major (PM) muscle in modern commercial broiler chickens, causing enormous economic losses in the poultry industry. Recent studies have observed hepatic and PM muscle injury in broilers affected by WB, but the relationships between WB and the 2 tissues are mostly unclear. In the current study, the RNA-seq raw data of PM muscle and liver were downloaded from GSE144000, and we constructed the gene coexpression networks of PM muscle and liver to explore the relationships between WB and the 2 tissues using the weighted gene coexpression network analysis (WGCNA) method. Six and 2 gene coexpression modules were significantly correlated with WB in the PM muscle and liver networks, respectively. TGF-beta signaling, Toll-like receptor signaling and mTOR signaling pathways were significantly enriched in the genes within the 6 gene modules of PM muscle network. Meanwhile, mTOR signaling pathway was significantly enriched in the genes within the 2 gene modules of liver network. In the consensus gene coexpression network across the 2 tissues, salmon module (r = -0.5 and p = 0.05) was significantly negatively correlated with WB, in which Toll-like receptor signaling, apoptosis, and autophagy pathways were significantly enriched. The genes related with the 3 pathways, myeloid differentiation primary response 88 (MYD88), interferon regulatory factor 7 (IRF7), mitogen-activated protein kinase 14 (MAPK14), FBJ murine osteosarcoma viral oncogene homolog (FOS), jun proto-oncogene (JUN), caspase-10, unc-51 like autophagy activating kinase 2 (ULK2) and serine/threonine kinase 11 (LKB1), were identified in salmon module. In this current study, we found that the signaling pathways related with cell inflammation, apoptosis and autophagy might influence WB across 2 tissues in broilers.

Expression patterns of housekeeping genes and tissue-specific genes in black goats across multiple tissues

Black goats are a significant meat breed in southern China. To investigate the expression patterns and biological functions of genes in various tissues of black goats, we analyzed housekeeping genes (HKGs), tissue-specific genes (TSGs), and hub genes (HUBGs) across 23 tissues. Additionally, we analyzed HKGs in 13 tissues under different feeding conditions. We identified 2968 HKGs, including six important ones. Interestingly, HKGs in grazing black goats demonstrated higher and more stable expression levels. We discovered a total of 9912 TSGs, including 134 newly identified ones. The number of TSGs for mRNA and lncRNA were nearly equal, with 127 mRNA TSGs expressed solely in one tissue. Additionally, the predicted functions of tissue-specific long non-coding RNAs (lncRNAs) targeting mRNAs corresponded with the physiological functions of the tissues.Weighted gene co-expression network analysis (WGCNA) constructed 30 modules, where the dark grey module consists almost entirely of HKGs, and TSGs are located in modules most correlated with their respective tissues. Additionally, we identified 289 HUBGs, which are involved in regulating the physiological functions of their respective tissues. Overall, these identified HKGs, TSGs, and HUBGs provide a foundation for studying the molecular mechanisms affecting the genetic and biological processes of complex traits in black goats.

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.

Transcriptomic profiling of gastrointestinal tracts in dairy cattle during lactation reveals molecular adaptations for milk synthesis

During lactation, dairy cattle's digestive tract requires significant adaptations to meet the increased nutrient demands for milk production. As we attempt to improve milk-related traits through selective pressure, it is crucial to understand the biological functions of the epithelia of the rumen, small intestine, and colonic tissues in response to changes in physiological state driven by changes in nutrient demands for milk synthesis. In this study, we obtained a total of 108 transcriptome profiles from three tissues (epithelia of the colon, duodenum, and rumen) of five Holstein cows, spanning eight time points from the early, mid, late lactation periods to the dry period. On average 97.06% of reads were successfully mapped to the reference genome assembly ARS-UCD1.2. We analyzed 27,607 gene expression patterns at multiple periods, enabling direct comparisons within and among tissues during different lactation stages, including early and peak lactation. We identified 1645, 813, and 2187 stage-specific genes in the colon, duodenum, and rumen, respectively, which were enriched for common or specific biological functions among different tissues. Time series analysis categorized the expressed genes within each tissue into four clusters. Furthermore, when the three tissues were analyzed collectively, 36 clusters of similarly expressed genes were identified. By integrating other comprehensive approaches such as gene co-expression analyses, functional enrichment, and cell type deconvolution, we gained profound insights into cattle lactation, revealing tissue-specific characteristics of the gastrointestinal tract and shedding light on the intricate molecular adaptations involved in nutrient absorption, immune regulation, and cellular processes for milk synthesis during lactation.

A review of the role of epigenetic studies for intramuscular fat deposition in beef cattle

Intramuscular fat (IMF) deposition profoundly influences meat quality and economic value in beef cattle production. Meanwhile, contemporary developments in epigenetics have opened new outlooks for understanding the molecular basics of IMF regulation, and it has become a key area of research for world scholars. Therefore, the aim of this paper was to provide insight and synthesis into the intricate relationship between epigenetic mechanisms and IMF deposition in beef cattle. The methodology involves a thorough analysis of existing literature, including pertinent books, academic journals, and online resources, to provide a comprehensive overview of the role of epigenetic studies in IMF deposition in beef cattle. This review summarizes the contemporary studies in epigenetic mechanisms in IMF regulation, high-resolution epigenomic mapping, single-cell epigenomics, multi-omics integration, epigenome editing approaches, longitudinal studies in cattle growth, environmental epigenetics, machine learning in epigenetics, ethical and regulatory considerations, and translation to industry practices from perspectives of IMF deposition in beef cattle. Moreover, this paper highlights DNA methylation, histone modifications, acetylation, phosphorylation, ubiquitylation, non-coding RNAs, DNA hydroxymethylation, epigenetic readers, writers, and erasers, chromatin immunoprecipitation followed by sequencing, whole genome bisulfite sequencing, epigenome-wide association studies, and their profound impact on the expression of crucial genes governing adipogenesis and lipid metabolism. Nutrition and stress also have significant influences on epigenetic modifications and IMF deposition. The key findings underscore the pivotal role of epigenetic studies in understanding and enhancing IMF deposition in beef cattle, with implications for precision livestock farming and ethical livestock management. In conclusion, this review highlights the crucial significance of epigenetic pathways and environmental factors in affecting IMF deposition in beef cattle, providing insightful information for improving the economics and meat quality of cattle production.

De novo transcriptome assembly database for 100 tissues from each of seven species of domestic herbivore

Domesticated herbivores are an important agricultural resource that play a critical role in global food security, particularly as they can adapt to varied environments, including marginal lands. An understanding of the molecular basis of their biology would contribute to better management and sustainable production. Thus, we conducted transcriptome sequencing of 100 to 105 tissues from two females of each of seven species of herbivore (cattle, sheep, goats, sika deer, horses, donkeys, and rabbits) including two breeds of sheep. The quality of raw and trimmed reads was assessed in terms of base quality, GC content, duplication sequence rate, overrepresented k-mers, and quality score distribution with FastQC. The high-quality filtered RNA-seq raw reads were deposited in a public database which provides approximately 54 billion high-quality paired-end sequencing reads in total, with an average mapping rate of ~93.92%. Transcriptome databases represent valuable resources that can be used to study patterns of gene expression, and pathways that are related to key biological processes, including important economic traits in herbivores.

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.

Pig-eRNAdb: a comprehensive enhancer and eRNA dataset of pigs

Enhancers and the enhancer RNAs (eRNAs) have been strongly implicated in regulations of transcriptions. Based the multi-omics data (ATAC-seq, ChIP-seq and RNA-seq) from public databases, Pig-eRNAdb is a dataset that comprehensively integrates enhancers and eRNAs for pigs using the machine learning strategy, which incorporates 82,399 enhancers and 37,803 eRNAs from 607 samples across 15 tissues of pigs. This user-friendly dataset covers a comprehensive depth of enhancers and eRNAs annotation for pigs. The coordinates of enhancers and the expression patterns of eRNAs are downloadable. Besides, thousands of regulators on eRNAs, the target genes of eRNAs, the tissue-specific eRNAs, and the housekeeping eRNAs are also accessible as well as the sequence similarity of eRNAs with humans. Moreover, the tissue-specific eRNA-trait associations encompass 652 traits are also provided. It will crucially facilitate investigations on enhancers and eRNAs with Pig-eRNAdb as a reference dataset in pigs.

Large-scale genomic and transcriptomic analyses elucidate the genetic basis of high meat yield in chickens

INTRODUCTION

Investigating the genetic markers and genomic signatures related to chicken meat production by combing multi-omics methods could provide new insights into modern chicken breeding technology systems.

OBJECT

Chicken is one of the most efficient and environmentally friendly livestock, especially the fast-growing white-feathered chicken (broiler), which is well known for high meat yield, but the underlying genetic basis is poorly understood.

METHOD

We generated whole-genome resequencing of three purebred broilers (n = 748) and six local breeds/lines (n = 114), and sequencing data of twelve chicken breeds (n = 199) were obtained from the NCBI database. Additionally, transcriptome sequencing of six tissues from two chicken breeds (n = 129) at two developmental stages was performed. A genome-wide association study combined with cis-eQTL mapping and the Mendelian randomization was applied.

RESULT

We identified > 17 million high-quality SNPs, of which 21.74% were newly identified, based on 21 chicken breeds/lines. A total of 163 protein-coding genes underwent positive selection in purebred broilers, and 83 genes were differentially expressed between purebred broilers and local chickens. Notably, muscle development was proven to be the major difference between purebred broilers and local chickens, or ancestors, based on genomic and transcriptomic evidence from multiple tissues and stages. The MYH1 gene family showed the top selection signatures and muscle-specific expression in purebred broilers. Furthermore, we found that the causal gene SOX6 influenced breast muscle yield and also related to myopathy occurrences. A refined haplotype was provided, which had a significant effect on SOX6 expression and phenotypic changes.

CONCLUSION

Our study provides a comprehensive atlas comprising the typical genomic variants and transcriptional characteristics for muscle development and suggests a new regulatory target (SOX6-MYH1s axis) for breast muscle yield and myopathy, which could aid in the development of genome-scale selective breeding aimed at high meat yield in broiler chickens.

Characterization and analysis of transcriptomes of multiple tissues from estrus and diestrus in pigs

A comprehensive understanding of the complex regulatory mechanisms governing estrus and ovulation across multiple tissues in mammals is imperative to improve the reproductive performance of livestock and mitigate ovulation-related disorders in humans. To comprehensively elucidate the regulatory landscape, we analyzed the transcriptome of protein-coding genes and long intergenic non-coding RNAs (lincRNAs) in 58 samples (including the hypothalamus, pituitary, ovary, vagina, and vulva) derived from European Large White gilts and Chinese Mi gilts during estrus and diestrus. We constructed an intricate regulatory network encompassing 358 hub genes across the five examined tissues. Furthermore, our investigation identified 85 differentially expressed lincRNAs that are predicted to target 230 genes associated with critical functions including behavior, receptors, and apoptosis. Importantly, we found that vital components of estrus and ovulation events involve "Apoptosis" pathway in the hypothalamus, "Autophagy" in the ovary, as well as "Hypoxia" and "Angiogenesis" in the vagina and vulva. We have identified several differentially expressed transcription factors (TFs), such as SPI1 and HES2, which regulate these pathways. SPI1 may suppress transcription in the autophagy pathway, promoting apoptosis and inhibiting the proliferation of ovarian granulosa cells. Our study provides the most comprehensive transcriptional profiling information related to estrus and ovulation events.

Transcriptomic analysis reveals diverse expression patterns underlying the fiber diameter of oxidative and glycolytic skeletal muscles in steers

Skeletal muscles consist of heterogeneous fibers with various contractile and metabolic properties that affect meat quality. The size of muscle fibers contributes to muscle mass and myopathy. Thus, improved understanding of the expression patterns underlying fiber size might open possibilities to change them using genetic methods. The aim of this study was to reveal transcriptomic landscapes of one oxidative (Psoas major) and three glycolytic (Longissimus lumborum, Triceps brachii, and Semimembranosus) muscles. Principal component analysis (PCA) showed significant differences in gene expression among the four muscles. Specifically, 2777 differentially expressed genes (DEGs) were detected between six pairwise comparisons of the four muscles. Weighted gene co-expression network analysis (WGCNA) identified six modules, which were significantly associated with muscle fiber diameter. We also identified 23 candidate genes, and enrichment analysis showed that biosynthesis of amino acids (bta01230), sarcomere (GO:0030017), and regulation of actin cytoskeleton (bta04810) overlapped in DEGs and WGCNA. Nineteen of these genes (e.g., EEF1A2, FARSB, and PINK1) have been reported to promote or inhibit muscle growth and development. Our findings contribute to the understanding of fiber size differences among oxidative and glycolytic muscles, which may provide a basis for breeding to improve meat yield.

Comprehensive Analysis of CircRNA Expression Profiles in Multiple Tissues of Pigs

Circular RNAs (circRNAs) are a class of non-coding RNAs with diverse functions, and previous studies have reported that circRNAs are involved in the growth and development of pigs. However, studies about porcine circRNAs over the past few years have focused on a limited number of tissues. Based on 215 publicly available RNA sequencing (RNA-seq) samples, we conducted a comprehensive analysis of circRNAs in nine pig tissues, namely, the gallbladder, heart, liver, longissimus dorsi, lung, ovary, pituitary, skeletal muscle, and spleen. Here, we identified a total of 82,528 circRNAs and discovered 3818 novel circRNAs that were not reported in the CircAtlas database. Moreover, we obtained 492 housekeeping circRNAs and 3489 tissue-specific circRNAs. The housekeeping circRNAs were enriched in signaling pathways regulating basic biological tissue activities, such as chromatin remodeling, nuclear-transcribed mRNA catabolic process, and protein methylation. The tissue-specific circRNAs were enriched in signaling pathways related to tissue-specific functions, such as muscle system process in skeletal muscle, cilium organization in pituitary, and cortical cytoskeleton in ovary. Through weighted gene co-expression network analysis, we identified 14 modules comprising 1377 hub circRNAs. Additionally, we explored circRNA-miRNA-mRNA networks to elucidate the interaction relationships between tissue-specific circRNAs and tissue-specific genes. Furthermore, our conservation analysis revealed that 19.29% of circRNAs in pigs shared homologous positions with their counterparts in humans. In summary, this extensive profiling of housekeeping, tissue-specific, and co-expressed circRNAs provides valuable insights into understanding the molecular mechanisms of pig transcriptional expression, ultimately deepening our understanding of genetic and biological processes.

Comparative analyses of dynamic transcriptome profiles highlight key response genes and dominant isoforms for muscle development and growth in chicken

BACKGROUND

Modern breeding strategies have resulted in significant differences in muscle mass between indigenous chicken and specialized broiler. However, the molecular regulatory mechanisms that underlie these differences remain elusive. The aim of this study was to identify key genes and regulatory mechanisms underlying differences in breast muscle development between indigenous chicken and specialized broiler.

RESULTS

Two time-series RNA-sequencing profiles of breast muscles were generated from commercial Arbor Acres (AA) broiler (fast-growing) and Chinese indigenous Lushi blue-shelled-egg (LS) chicken (slow-growing) at embryonic days 10, 14, and 18, and post-hatching day 1 and weeks 1, 3, and 5. Principal component analysis of the transcriptome profiles showed that the top four principal components accounted for more than 80% of the total variance in each breed. The developmental axes between the AA and LS chicken overlapped at the embryonic stages but gradually separated at the adult stages. Integrative investigation of differentially-expressed transcripts contained in the top four principal components identified 44 genes that formed a molecular network associated with differences in breast muscle mass between the two breeds. In addition, alternative splicing analysis revealed that genes with multiple isoforms always had one dominant transcript that exhibited a significantly higher expression level than the others. Among the 44 genes, the TNFRSF6B gene, a mediator of signal transduction pathways and cell proliferation, harbored two alternative splicing isoforms, TNFRSF6B-X1 and TNFRSF6B-X2. TNFRSF6B-X1 was the dominant isoform in both breeds before the age of one week. A switching event of the dominant isoform occurred at one week of age, resulting in TNFRSF6B-X2 being the dominant isoform in AA broiler, whereas TNFRSF6B-X1 remained the dominant isoform in LS chicken. Gain-of-function assays demonstrated that both isoforms promoted the proliferation of chicken primary myoblasts, but only TNFRSF6B-X2 augmented the differentiation and intracellular protein content of chicken primary myoblasts.

CONCLUSIONS

For the first time, we identified several key genes and dominant isoforms that may be responsible for differences in muscle mass between slow-growing indigenous chicken and fast-growing commercial broiler. These findings provide new insights into the regulatory mechanisms underlying breast muscle development in chicken.

Genome-wide scan for runs of homozygosity in South American Camelids

BACKGROUND

Alpaca (Vicugna pacos), llama (Lama glama), vicugna (Vicugna vicugna) and guanaco (Lama guanicoe), are the camelid species distributed over the Andean high-altitude grasslands, the Altiplano, and the Patagonian arid steppes. Despite the wide interest on these animals, most of the loci under selection are still unknown. Using whole-genome sequencing (WGS) data we investigated the occurrence and the distribution of Runs Of Homozygosity (ROHs) across the South American Camelids (SACs) genome to identify the genetic relationship between the four species and the potential signatures of selection.

RESULTS

A total of 37 WGS samples covering the four species was included in the final analysis. The multi-dimensional scaling approach showed a clear separation between the four species; however, admixture analysis suggested a strong genetic introgression from vicugna and llama to alpaca. Conversely, very low genetic admixture of the guanaco with the other SACs was found. The four species did not show significant differences in the number, length of ROHs (100-500 kb) and genomic inbreeding values. Longer ROHs (> 500 kb) were found almost exclusively in alpaca. Seven overlapping ROHs were shared by alpacas, encompassing nine loci (FGF5, LOC107034918, PRDM8, ANTXR2, LOC102534792, BSN, LOC116284892, DAG1 and RIC8B) while nine overlapping ROHs were found in llama with twenty-five loci annotated (ERC2, FZD9, BAZ1B, BCL7B, LOC116284208, TBL2, MLXIPL, PHF20, TRNAD-AUC, LOC116284365, RBM39, ARFGEF2, DCAF5, EXD2, HSPB11, LRRC42, LDLRAD1, TMEM59, LOC107033213, TCEANC2, LOC102545169, LOC116278408, SMIM15, NDUFAF2 and RCOR1). Four overlapping ROHs, with three annotated loci (DLG1, KAT6B and PDE4D) and three overlapping ROHs, with seven annotated genes (ATP6V1E1, BCL2L13, LOC116276952, BID, KAT6B, LOC116282667 and LOC107034552), were detected for vicugna and guanaco, respectively.

CONCLUSIONS

The signatures of selection revealed genomic areas potentially selected for production traits as well as for natural adaptation to harsh environment. Alpaca and llama hint a selection driven by environment as well as by farming purpose while vicugna and guanaco showed selection signals for adaptation to harsh environment. Interesting, signatures of selection on KAT6B gene were identified for both vicugna and guanaco, suggesting a positive effect on wild populations fitness. Such information may be of interest to further ecological and animal production studies.

The eQTL colocalization and transcriptome-wide association study identify potentially causal genes responsible for economic traits in Simmental beef cattle

BACKGROUND

A detailed understanding of genetic variants that affect beef merit helps maximize the efficiency of breeding for improved production merit in beef cattle. To prioritize the putative variants and genes, we ran a comprehensive genome-wide association studies (GWAS) analysis for 21 agronomic traits using imputed whole-genome variants in Simmental beef cattle. Then, we applied expression quantitative trait loci (eQTL) mapping between the genotype variants and transcriptome of three tissues (longissimus dorsi muscle, backfat, and liver) in 120 cattle.

RESULTS

We identified 1,580 association signals for 21 beef agronomic traits using GWAS. We then illuminated 854,498 cis-eQTLs for 6,017 genes and 46,970 trans-eQTLs for 1,903 genes in three tissues and built a synergistic network by integrating transcriptomics with agronomic traits. These cis-eQTLs were preferentially close to the transcription start site and enriched in functional regulatory regions. We observed an average of 43.5% improvement in cis-eQTL discovery using multi-tissue eQTL mapping. Fine-mapping analysis revealed that 111, 192, and 194 variants were most likely to be causative to regulate gene expression in backfat, liver, and muscle, respectively. The transcriptome-wide association studies identified 722 genes significantly associated with 11 agronomic traits. Via the colocalization and Mendelian randomization analyses, we found that eQTLs of several genes were associated with the GWAS signals of agronomic traits in three tissues, which included genes, such as NADSYN1, NDUFS3, LTF and KIFC2 in liver, GRAMD1C, TMTC2 and ZNF613 in backfat, as well as TIGAR, NDUFS3 and L3HYPDH in muscle that could serve as the candidate genes for economic traits.

CONCLUSIONS

The extensive atlas of GWAS, eQTL, fine-mapping, and transcriptome-wide association studies aid in the suggestion of potentially functional variants and genes in cattle agronomic traits and will be an invaluable source for genomics and breeding in beef cattle.

A Multi-Tissue Gene Expression Atlas of Water Buffalo (Bubalus bubalis) Reveals Transcriptome Conservation between Buffalo and Cattle

We generated 73 transcriptomic data of water buffalo, which were integrated with publicly available data in this species, yielding a large dataset of 355 samples representing 20 major tissue categories. We established a multi-tissue gene expression atlas of water buffalo. Furthermore, by comparing them with 4866 cattle transcriptomic data from the cattle genotype-tissue expression atlas (CattleGTEx), we found that the transcriptomes of the two species exhibited conservation in their overall gene expression patterns, tissue-specific gene expression and house-keeping gene expression. We further identified conserved and divergent expression genes between the two species, with the largest number of differentially expressed genes found in the skin, which may be related to structural and functional differences in the skin of the two species. This work provides a source of functional annotation of the buffalo genome and lays the foundations for future genetic and evolutionary studies in water buffalo.

Expression Profile of Housekeeping Genes and Tissue-Specific Genes in Multiple Tissues of Pigs

Pigs have become an ideal model system for human disease research and development and an important farm animal that provides a valuable source of nutrition. To profile the all-sided gene expression and their biological functions across multiple tissues, we conducted a comprehensive analysis of gene expression on a large scale around the side of housekeeping genes (HKGs), tissue specific genes (TSGs), and the co-expressed genes in 14 various tissues. In this study, we identified 2351 HKGs and 3018 TSGs across tissues, among which 4 HKGs (COX1, UBB, OAZ1/NPFF) exhibited low variation and high expression levels, and 31 particular TSGs (e.g., PDC, FKBP6, STAT2, and COL1A1) were exclusively expressed in several tissues, including endocrine brain, ovaries, livers, backfat, jejunum, kidneys, lungs, and longissimus dorsi muscles. We also obtained 17 modules with 230 hub genes (HUBGs) by weighted gene co-expression network analysis. On the other hand, HKGs functions were enriched in the signaling pathways of the ribosome, spliceosome, thermogenesis, oxidative phosphorylation, and nucleocytoplasmic transport, which have been highly suggested to involve in the basic biological tissue activities. While TSGs were highly enriched in the signaling pathways that were involved in specific physiological processes, such as the ovarian steroidogenesis pathway in ovaries and the renin-angiotensin system pathway in kidneys. Collectively, these stable, specifical, and co-expressed genes provided useful information for the investigation of the molecular mechanism for an understanding of the genetic and biological processes of complex traits in pigs.