STAT6, PBX2, and PBRM1 Emerge as Predicted Regulators of 452 Differentially Expressed Genes Associated With Puberty in Brahman Heifers

Identification of quantitative trait loci and associated candidate genes for pregnancy success in Angus-Brahman crossbred heifers

BACKGROUND

In beef cattle, more than 50% of the energy input to produce a unit of beef is consumed by the female that produced the calf. Development of genomic tools to identify females with high genetic merit for reproductive function could increase the profitability and sustainability of beef production.

RESULTS

Genome-wide association studies (GWAS) were performed using a single-step genomic best linear unbiased prediction approach on pregnancy outcome traits from a population of Angus-Brahman crossbred heifers. Furthermore, a validation GWAS was performed using data from another farm. Heifers were genotyped with the Bovine GGP F250 array that contains 221,077 SNPs. In the discovery population, heifers were bred in winter breeding seasons involving a single round of timed artificial insemination (AI) followed by natural mating for 3 months. Two phenotypes were analyzed: pregnancy outcome to first-service AI (PAI; n = 1,481) and pregnancy status at the end of the breeding season (PEBS; n = 1,725). The heritability was estimated as 0.149 and 0.122 for PAI and PEBS, respectively. In the PAI model, one quantitative trait locus (QTL), located between 52.3 and 52.5 Mb on BTA7, explained about 3% of the genetic variation, in a region containing a cluster of γ-protocadherin genes and SLC25A2. Other QTLs explaining between 0.5% and 1% of the genetic variation were found on BTA12 and 25. In the PEBS model, a large QTL on BTA7 was synonymous with the QTL for PAI, with minor QTLs located on BTA5, 9, 10, 11, 19, and 20. The validation population for pregnancy status at the end of the breeding season were Angus-Brahman crossbred heifers bred by natural mating. In concordance with the discovery population, the large QTL on BTA7 and QTLs on BTA10 and 12 were identified.

CONCLUSIONS

In summary, QTLs and candidate SNPs identified were associated with pregnancy outcomes in beef heifers, including a large QTL associated with a group of protocadherin genes. Confirmation of these associations with larger populations could lead to the development of genomic predictions of reproductive function in beef cattle. Moreover, additional research is warranted to study the function of candidate genes associated with QTLs.

Investigation of lncRNA in Bos taurus Mammary Tissue during Dry and Lactation Periods

This study aims to collect RNA-Seq data from Bos taurus samples representing dry and lactating mammary tissue, identify lncRNA transcripts, and analyze findings for their features and functional annotation. This allows for connections to be drawn between lncRNA and the lactation process. RNA-Seq data from 103 samples of Bos taurus mammary tissue were gathered from publicly available databases (60 dry, 43 lactating). The samples were filtered to reveal 214 dry mammary lncRNA transcripts and 517 lactating mammary lncRNA transcripts. The lncRNAs met common lncRNA characteristics such as shorter length, fewer exons, lower expression levels, and less sequence conservation when compared to the genome. Interestingly, several lncRNAs showed sequence similarity to genes associated with strong hair keratin intermediate filaments. Human breast cancer research has associated strong hair keratin filaments with mammary tissue cellular resilience. The lncRNAs were also associated with several genes/proteins that linked to pregnancy using expression correlation and gene ontology. Such findings indicate that there are crucial relationships between the lncRNAs found in mammary tissue and the development of the tissue, to meet both the animal's needs and our own production needs; these relationships should be further investigated to ensure that we continue to breed the most resilient, efficient dairy cattle.

Unrevealing functional candidate genes for bovine fertility through RNA sequencing meta-analysis and regulatory elements networks of co-expressed genes and lncRNAs

The high genetic heterogeneity and environmental effects of subfertility in livestock species make the elucidation of the genetic mechanisms associated with reproductive efficiency a difficult task. Network and co-expression network meta-analyses were applied alongside genetic variant calling and long non-coding RNA (lncRNA) characterization to identify functionally relevant target genes and regulatory subnetworks associated with fertility in dairy cattle. In total, 505 lncRNAs (441 previously annotated in the bovine reference genome ARS-UCD 1.2 and 64 novel lncRNAs) were identified. Seven differentially expressed genes between high-fertile (HF) and sub-fertile (SF) Holstein cows were identified in the network meta-analysis (CA5A, ENSBTAG00000051149, ENSBTAG00000003272, DEFB7, DIO2, TRPV3, and COL4A4). Additionally, seven functional candidate differentially co-expressed (DcoExp) modules with a differential regulatory pattern (|z-score|>2) were identified between HF and SF cows. The functional candidate genes and DcoExp modules identified were associated with fertility relevant processes such as the regulation of embryonic implantation and proliferation, interaction and molecule transfer between the fetus and the cow, and the immune system. These results help to better understand the genetic mechanisms associated with reproductive efficiency in dairy cattle through the identification of potential biomarkers and genetic variants associated with differentially expressed regulatory gene and lncRNAs regulatory element networks.

Proteomic Analysis of Hypothalamus and Pituitary Gland in Pre and Postpubertal Brahman Heifers

The hypothalamus and the pituitary gland are directly involved in the complex systemic changes that drive the onset of puberty in cattle. Here, we applied integrated bioinformatics to elucidate the critical proteins underlying puberty and uncover potential molecular mechanisms from the hypothalamus and pituitary gland of prepubertal (n = 6) and postpubertal (n = 6) cattle. Proteomic analysis in the hypothalamus and pituitary gland revealed 275 and 186 differentially abundant (DA) proteins, respectively (adjusted p-value < 0.01). The proteome profiles found herein were integrated with previously acquired transcriptome profiles. These transcriptomic studies used the same tissues harvested from the same heifers at pre- and post-puberty. This comparison detected a small number of matched transcripts and protein changes at puberty in each tissue, suggesting the need for multiple omics analyses for interpreting complex biological systems. In the hypothalamus, upregulated DA proteins at post-puberty were enriched in pathways related to puberty, including GnRH, calcium and oxytocin signalling pathways, whereas downregulated proteins were observed in the estrogen signalling pathway, axon guidance and GABAergic synapse. Additionally, this study revealed that ribosomal pathway proteins in the pituitary were involved in the pubertal development of mammals. The reported molecules and derived protein-protein networks are a starting point for future experimental approaches that might dissect with more detail the role of each molecule to provide new insights into the mechanisms of puberty onset in cattle.

Leveraging transcriptome and epigenome landscapes to infer regulatory networks during the onset of sexual maturation

Background

Despite sexual development being ubiquitous to vertebrates, the molecular mechanisms underpinning this fundamental transition remain largely undocumented in many organisms. We designed a time course experiment that successfully sampled the period when Atlantic salmon commence their trajectory towards sexual maturation.

Results

Through deep RNA sequencing, we discovered key genes and pathways associated with maturation in the pituitary-ovarian axis. Analyzing DNA methylomes revealed a bias towards hypermethylation in ovary that implicated maturation-related genes. Co-analysis of DNA methylome and gene expression changes revealed chromatin remodeling genes and key transcription factors were both significantly hypermethylated and upregulated in the ovary during the onset of maturation. We also observed changes in chromatin state landscapes that were strongly correlated with fundamental remodeling of gene expression in liver. Finally, a multiomic integrated analysis revealed regulatory networks and identified hub genes including TRIM25 gene (encoding the estrogen-responsive finger protein) as a putative key regulator in the pituitary that underwent a 60-fold change in connectivity during the transition to maturation.

Conclusion

The study successfully documented transcriptome and epigenome changes that involved key genes and pathways acting in the pituitary – ovarian axis. Using a Systems Biology approach, we identified hub genes and their associated networks deemed crucial for onset of maturation. The results provide a comprehensive view of the spatiotemporal changes involved in a complex trait and opens the door to future efforts aiming to manipulate puberty in an economically important aquaculture species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08514-8.

Identifying pleiotropic variants and candidate genes for fertility and reproduction traits in Holstein cattle via association studies based on imputed whole-genome sequence genotypes

Background

Genetic progress for fertility and reproduction traits in dairy cattle has been limited due to the low heritability of most indicator traits. Moreover, most of the quantitative trait loci (QTL) and candidate genes associated with these traits remain unknown. In this study, we used 5.6 million imputed DNA sequence variants (single nucleotide polymorphisms, SNPs) for genome-wide association studies (GWAS) of 18 fertility and reproduction traits in Holstein cattle. Aiming to identify pleiotropic variants and increase detection power, multiple-trait analyses were performed using a method to efficiently combine the estimated SNP effects of single-trait GWAS based on a chi-square statistic.

Results

There were 87, 72, and 84 significant SNPs identified for heifer, cow, and sire traits, respectively, which showed a wide and distinct distribution across the genome, suggesting that they have relatively distinct polygenic nature. The biological functions of immune response and fatty acid metabolism were significantly enriched for the 184 and 124 positional candidate genes identified for heifer and cow traits, respectively. No known biological function was significantly enriched for the 147 positional candidate genes found for sire traits. The most important chromosomes that had three or more significant QTL identified are BTA22 and BTA23 for heifer traits, BTA8 and BTA17 for cow traits, and BTA4, BTA7, BTA17, BTA22, BTA25, and BTA28 for sire traits. Several novel and biologically important positional candidate genes were strongly suggested for heifer (SOD2, WTAP, DLEC1, PFKFB4, TRIM27, HECW1, DNAH17, and ADAM3A), cow (ANXA1, PCSK5, SPESP1, and JMJD1C), and sire (ELMO1, CFAP70, SOX30, DGCR8, SEPTIN14, PAPOLB, JMJD1C, and NELL2) traits.

Conclusions

These findings contribute to better understand the underlying biological mechanisms of fertility and reproduction traits measured in heifers, cows, and sires, which may contribute to improve genomic evaluation for these traits in dairy cattle.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08555-z.

Extensive Variation in Gene Expression is Revealed in 13 Fertility-Related Genes Using RNA-Seq, ISO-Seq, and CAGE-Seq From Brahman Cattle

Fertility is a key driver of economic profitability in cattle production. A number of studies have identified genes associated with fertility using genome wide association studies and differential gene expression analysis; however, the genes themselves are poorly characterized in cattle. Here, we selected 13 genes from the literature which have previously been shown to have strong evidence for an association with fertility in Brahman cattle (Bos taurus indicus) or closely related breeds. We examine the expression variation of the 13 genes that are associated with cattle fertility using RNA-seq, CAGE-seq, and ISO-seq data from 11 different tissue samples from an adult Brahman cow and a Brahman fetus. Tissues examined include blood, liver, lung, kidney, muscle, spleen, ovary, and uterus from the cow and liver and lung from the fetus. The analysis revealed several novel isoforms, including seven from SERPINA7. The use of three expression characterization methodologies (5′ cap selected ISO-seq, CAGE-seq, and RNA-seq) allowed the identification of isoforms that varied in their length of 5′ and 3′ untranslated regions, variation otherwise undetectable (collapsed as degraded RNA) in generic isoform identification pipelines. The combinations of different sequencing technologies allowed us to overcome the limitations of relatively low sequence depth in the ISO-seq data. The lower sequence depth of the ISO-seq data was also reflected in the lack of observed expression of some genes that were observed in the CAGE-seq and RNA-seq data from the same tissue. We identified allele specific expression that was tissue-specific in AR, IGF1, SOX9, STAT3, and TAF9B. Finally, we characterized an exon of TAF9B as partially nested within the neighboring gene phosphoglycerate kinase 1. As this study only examined two animals, even more transcriptional variation may be present in a genetically diverse population. This analysis reveals the large amount of transcriptional variation within mammalian fertility genes and illuminates the fact that the transcriptional landscape cannot be fully characterized using a single technology alone.

RGD v2.0: a major update of the ruminant functional and evolutionary genomics database

Ruminant Genome Database (RGD; http://animal.nwsuaf.edu.cn/RGD) provides visualization and analysis tools for ruminant comparative genomics and functional annotations. As more high-quality ruminant genome assemblies have become available, we have redesigned the user interface, integrated and expanded multi-omics data, and developed novel features to improve the database. The new version, RGD v2.0, houses 78 ruminant genomes; 110-species synteny alignments for major livestock (including cattle, sheep, goat) and wild ungulates; 21 012 orthologous gene clusters with Gene Ontology and pathway annotation; ∼8 600 000 conserved elements; and ∼1 000 000 cis-regulatory elements by utilizing 1053 epigenomic data sets. The transcriptome data in RGD v2.0 has nearly doubled, currently with 1936 RNA-seq data sets, and 155 174 phenotypic data sets have been newly added. New and updated features include: (i) The UCSC Genome Browser, BLAT, BLAST and Table Browser tools were updated for six available ruminant livestock species. (ii) The LiftOver tool was newly introduced into our browser to allow coordinate conversion between different ruminant assemblies. And (iii) tissue specificity index, tau, was calculated to facilitate batch screening of specifically expressed genes. The enhanced genome annotations and improved functionality in RGD v2.0 will be useful for study of genome evolution, environmental adaption, livestock breeding and biomedicine.

Chromatin accessibility and regulatory vocabulary across indicine cattle tissues

BACKGROUND

Spatiotemporal changes in the chromatin accessibility landscape are essential to cell differentiation, development, health, and disease. The quest of identifying regulatory elements in open chromatin regions across different tissues and developmental stages is led by large international collaborative efforts mostly focusing on model organisms, such as ENCODE. Recently, the Functional Annotation of Animal Genomes (FAANG) has been established to unravel the regulatory elements in non-model organisms, including cattle. Now, we can transition from prediction to validation by experimentally identifying the regulatory elements in tropical indicine cattle. The identification of regulatory elements, their annotation and comparison with the taurine counterpart, holds high promise to link regulatory regions to adaptability traits and improve animal productivity and welfare.

RESULTS

We generate open chromatin profiles for liver, muscle, and hypothalamus of indicine cattle through ATAC-seq. Using robust methods for motif discovery, motif enrichment and transcription factor binding sites, we identify potential master regulators of the epigenomic profile in these three tissues, namely HNF4, MEF2, and SOX factors, respectively. Integration with transcriptomic data allows us to confirm some of their target genes. Finally, by comparing our results with Bos taurus data we identify potential indicine-specific open chromatin regions and overlaps with indicine selective sweeps.

CONCLUSIONS

Our findings provide insights into the identification and analysis of regulatory elements in non-model organisms, the evolution of regulatory elements within two cattle subspecies as well as having an immediate impact on the animal genetics community in particular for a relevant productive species such as tropical cattle.

Genome-wide association study to identify genomic regions and positional candidate genes associated with male fertility in beef cattle

Fertility plays a key role in the success of calf production, but there is evidence that reproductive efficiency in beef cattle has decreased during the past half-century worldwide. Therefore, identifying animals with superior fertility could significantly impact cow-calf production efficiency. The objective of this research was to identify candidate regions affecting bull fertility in beef cattle and positional candidate genes annotated within these regions. A GWAS using a weighted single-step genomic BLUP approach was performed on 265 crossbred beef bulls to identify markers associated with scrotal circumference (SC) and sperm motility (SM). Eight windows containing 32 positional candidate genes and five windows containing 28 positional candidate genes explained more than 1% of the genetic variance for SC and SM, respectively. These windows were selected to perform gene annotation, QTL enrichment, and functional analyses. Functional candidate gene prioritization analysis revealed 14 prioritized candidate genes for SC of which MAP3K1 and VIP were previously found to play roles in male fertility. A different set of 14 prioritized genes were identified for SM and five were previously identified as regulators of male fertility (SOD2, TCP1, PACRG, SPEF2, PRLR). Significant enrichment results were identified for fertility and body conformation QTLs within the candidate windows. Gene ontology enrichment analysis including biological processes, molecular functions, and cellular components revealed significant GO terms associated with male fertility. The identification of these regions contributes to a better understanding of fertility associated traits and facilitates the discovery of positional candidate genes for future investigation of causal mutations and their implications.

Gene regulation could be attributed to TCF3 and other key transcription factors in the muscle of pubertal heifers

Puberty is a whole-body event, driven by the hypothalamic integration of peripheral signals such as leptin or IGF-1. In the process of puberty, reproductive development is simultaneous to growth, including muscle growth. To enhance our understanding of muscle function related to puberty, we performed transcriptome analyses of muscle samples from six pre- and six post-pubertal Brahman heifers (Bos indicus). Our aims were to perform differential expression analyses and co-expression analyses to derive a regulatory gene network associate with puberty. As a result, we identified 431 differentially expressed (DEx) transcripts (genes and non-coding RNAs) when comparing pre- to post-pubertal average gene expression. The DEx transcripts were compared with all expressed transcripts in our samples (over 14,000 transcripts) for functional enrichment analyses. The DEx transcripts were associated with "extracellular region," "inflammatory response" and "hormone activity" (adjusted p < .05). Inflammatory response for muscle regeneration is a necessary aspect of muscle growth, which is accelerated during puberty. The term "hormone activity" may signal genes that respond to progesterone signalling in the muscle, as the presence of this hormone is an important difference between pre- and post-pubertal heifers in our experimental design. The DEx transcript with the highest average expression difference was a mitochondrial gene, ENSBTAG00000043574 that might be another important link between energy metabolism and puberty. In the derived co-expression gene network, we identified six hub genes: CDC5L, MYC, TCF3, RUNX2, ATF2 and CREB1. In the same network, 48 key regulators of DEx transcripts were identified, using a regulatory impact factor metric. The hub gene TCF3 was also a key regulator. The majority of the key regulators (22 genes) are members of the zinc finger family, which has been implicated in bovine puberty in other tissues. In conclusion, we described how puberty may affect muscle gene expression in cattle.

Genome-Wide Co-Expression Distributions as a Metric to Prioritize Genes of Functional Importance

Genome-wide gene expression analysis are routinely used to gain a systems-level understanding of complex processes, including network connectivity. Network connectivity tends to be built on a small subset of extremely high co-expression signals that are deemed significant, but this overlooks the vast majority of pairwise signals. Here, we developed a computational pipeline to assign to every gene its pair-wise genome-wide co-expression distribution to one of 8 template distributions shapes varying between unimodal, bimodal, skewed, or symmetrical, representing different proportions of positive and negative correlations. We then used a hypergeometric test to determine if specific genes (regulators versus non-regulators) and properties (differentially expressed or not) are associated with a particular distribution shape. We applied our methodology to five publicly available RNA sequencing (RNA-seq) datasets from four organisms in different physiological conditions and tissues. Our results suggest that genes can be assigned consistently to pre-defined distribution shapes, regarding the enrichment of differential expression and regulatory genes, in situations involving contrasting phenotypes, time-series, or physiological baseline data. There is indeed a striking additional biological signal present in the genome-wide distribution of co-expression values which would be overlooked by currently adopted approaches. Our method can be applied to extract further information from transcriptomic data and help uncover the molecular mechanisms involved in the regulation of complex biological process and phenotypes.

Weighted Gene Correlation Network Meta-Analysis Reveals Functional Candidate Genes Associated with High- and Sub-Fertile Reproductive Performance in Beef Cattle

Improved reproductive efficiency could lead to economic benefits for the beef industry, once the intensive selection pressure has led to a decreased fertility. However, several factors limit our understanding of fertility traits, including genetic differences between populations and statistical limitations. In the present study, the RNA-sequencing data from uterine samples of high-fertile (HF) and sub-fertile (SF) animals was integrated using co-expression network meta-analysis, weighted gene correlation network analysis, identification of upstream regulators, variant calling, and network topology approaches. Using this pipeline, top hub-genes harboring fixed variants (HF × SF) were identified in differentially co-expressed gene modules (DcoExp). The functional prioritization analysis identified the genes with highest potential to be key-regulators of the DcoExp modules between HF and SF animals. Consequently, 32 functional candidate genes (10 upstream regulators and 22 top hub-genes of DcoExp modules) were identified. These genes were associated with the regulation of relevant biological processes for fertility, such as embryonic development, germ cell proliferation, and ovarian hormone regulation. Additionally, 100 candidate variants (single nucleotide polymorphisms (SNPs) and insertions and deletions (INDELs)) were identified within those genes. In the long-term, the results obtained here may help to reduce the frequency of subfertility in beef herds, reducing the associated economic losses caused by this condition.

Identification and Annotation of Potential Function of Regulatory Antisense Long Non-Coding RNAs Related to Feed Efficiency in Bos taurus Bulls

Long non-coding RNAs (lncRNAs) can influence transcriptional and translational processes in mammalian cells and are associated with various developmental, physiological and phenotypic conditions. However, they remain poorly understood and annotated in livestock species. We combined phenotypic, metabolomics and liver transcriptomic data of bulls divergent for residual feed intake (RFI) and fat accretion. Based on a project-specific transcriptome annotation for the bovine reference genome ARS-UCD.1.2 and multiple-tissue total RNA sequencing data, we predicted 3590 loci to be lncRNAs. To identify lncRNAs with potential regulatory influence on phenotype and gene expression, we applied the regulatory impact factor algorithm on a functionally prioritized set of loci (n = 4666). Applying the algorithm of partial correlation and information theory, significant and independent pairwise correlations were calculated and co-expression networks were established, including plasma metabolites correlated with lncRNAs. The network hub lncRNAs were assessed for potential cis-actions and subjected to biological pathway enrichment analyses. Our results reveal a prevalence of antisense lncRNAs positively correlated with adjacent protein-coding genes and suggest their participation in mitochondrial function, acute phase response signalling, TCA-cycle, fatty acid β-oxidation and presumably gluconeogenesis. These antisense lncRNAs indicate a stabilizing function for their cis-correlated genes and a putative regulatory role in gene expression.

Enhanced identification of significant regulators of gene expression

Background

Diseases like cancer will lead to changes in gene expression, and it is relevant to identify key regulatory genes that can be linked directly to these changes. This can be done by computing a Regulatory Impact Factor (RIF) score for relevant regulators. However, this computation is based on estimating correlated patterns of gene expression, often Pearson correlation, and an assumption about a set of specific regulators, normally transcription factors. This study explores alternative measures of correlation, using the Fisher and Sobolev metrics, and an extended set of regulators, including epigenetic regulators and long non-coding RNAs (lncRNAs). Data on prostate cancer have been used to explore the effect of these modifications.

Results

A tool for computation of RIF scores with alternative correlation measures and extended sets of regulators was developed and tested on gene expression data for prostate cancer. The study showed that the Fisher and Sobolev metrics lead to improved identification of well-documented regulators of gene expression in prostate cancer, and the sets of identified key regulators showed improved overlap with previously defined gene sets of relevance to cancer. The extended set of regulators lead to identification of several interesting candidates for further studies, including lncRNAs. Several key processes were identified as important, including spindle assembly and the epithelial-mesenchymal transition (EMT).

Conclusions

The study has shown that using alternative metrics of correlation can improve the performance of tools based on correlation of gene expression in genomic data. The Fisher and Sobolev metrics should be considered also in other correlation-based applications.

Genomic measures of inbreeding coefficients and genome-wide scan for runs of homozygosity islands in Iranian river buffalo, Bubalus bubalis

Background

Consecutive homozygous fragments of a genome inherited by offspring from a common ancestor are known as runs of homozygosity (ROH). ROH can be used to calculate genomic inbreeding and to identify genomic regions that are potentially under historical selection pressure. The dataset of our study consisted of 254 Azeri (AZ) and 115 Khuzestani (KHZ) river buffalo genotyped for ~ 65,000 SNPs for the following two purposes: 1) to estimate and compare inbreeding calculated using ROH (F_ROH), excess of homozygosity (F_HOM), correlation between uniting gametes (F_UNI), and diagonal elements of the genomic relationship matrix (F_GRM); 2) to identify frequently occurring ROH (i.e. ROH islands) for our selection signature and gene enrichment studies.

Results

In this study, 9102 ROH were identified, with an average number of 21.2 ± 13.1 and 33.2 ± 15.9 segments per animal in AZ and KHZ breeds, respectively. On average in AZ, 4.35% (108.8 ± 120.3 Mb), and in KHZ, 5.96% (149.1 ± 107.7 Mb) of the genome was autozygous. The estimated inbreeding values based on F_HOM, F_UNI and F_GRM were higher in AZ than they were in KHZ, which was in contrast to the F_ROH estimates. We identified 11 ROH islands (four in AZ and seven in KHZ). In the KHZ breed, the genes located in ROH islands were enriched for multiple Gene Ontology (GO) terms (P ≤ 0.05). The genes located in ROH islands were associated with diverse biological functions and traits such as body size and muscle development (BMP2), immune response (CYP27B1), milk production and components (MARS, ADRA1A, and KCTD16), coat colour and pigmentation (PMEL and MYO1A), reproductive traits (INHBC, INHBE, STAT6 and PCNA), and bone development (SUOX).

Conclusion

The calculated F_ROH was in line with expected higher inbreeding in KHZ than in AZ because of the smaller effective population size of KHZ. Thus, we find that F_ROH can be used as a robust estimate of genomic inbreeding. Further, the majority of ROH peaks were overlapped with or in close proximity to the previously reported genomic regions with signatures of selection. This tells us that it is likely that the genes in the ROH islands have been subject to artificial or natural selection.

Biological Network Approach for the Identification of Regulatory Long Non-Coding RNAs Associated With Metabolic Efficiency in Cattle

Background: Genomic regions associated with divergent livestock feed efficiency have been found predominantly outside protein coding sequences. Long non-coding RNAs (lncRNA) can modulate chromatin accessibility, gene expression and act as important metabolic regulators in mammals. By integrating phenotypic, transcriptomic, and metabolomic data with quantitative trait locus data in prioritizing co-expression network analyses, we aimed to identify and functionally characterize lncRNAs with a potential key regulatory role in metabolic efficiency in cattle. Materials and Methods: Crossbred animals (n = 48) of a Charolais x Holstein F₂-population were allocated to groups of high or low metabolic efficiency based on residual feed intake in bulls, energy corrected milk in cows and intramuscular fat content in both genders. Tissue samples from jejunum, liver, skeletal muscle and rumen were subjected to global transcriptomic analysis via stranded total RNA sequencing (RNAseq) and blood plasma samples were used for profiling of 640 metabolites. To identify lncRNAs within the indicated tissues, a project-specific transcriptome annotation was established. Subsequently, novel transcripts were categorized for potential lncRNA status, yielding a total of 7,646 predicted lncRNA transcripts belonging to 3,287 loci. A regulatory impact factor approach highlighted 92, 55, 35, and 73 lncRNAs in jejunum, liver, muscle, and rumen, respectively. Their ensuing high regulatory impact factor scores indicated a potential regulatory key function in a gene set comprising loci displaying differential expression, tissue specificity and loci overlapping with quantitative trait locus regions for residual feed intake or milk production. These were subjected to a partial correlation and information theory analysis with the prioritized gene set. Results and Conclusions: Independent, significant and group-specific correlations (|r| > 0.8) were used to build a network for the high and the low metabolic efficiency group resulting in 1,522 and 1,732 nodes, respectively. Eight lncRNAs displayed a particularly high connectivity (>100 nodes). Metabolites and genes from the partial correlation and information theory networks, which each correlated significantly with the respective lncRNA, were included in an enrichment analysis indicating distinct affected pathways for the eight lncRNAs. LncRNAs associated with metabolic efficiency were classified to be functionally involved in hepatic amino acid metabolism and protein synthesis and in calcium signaling and neuronal nitric oxide synthase signaling in skeletal muscle cells.

In-depth characterization of the pituitary transcriptome in Simmental and Chinese native cattle

The pituitary gland is a key endocrine organ responsible for growth and development. To get closer to understanding important molecular mechanisms at work in the bovine pituitary and identifying the core genes controlling growth, in the present study we have analyzed the transcriptome profiling of the pituitary glands of two cattle breeds (Wenshan and Simmental cattle) with extreme growth differences using high-throughput RNA sequencing. Our study revealed that the overall gene expression levels are quite similar between the two breeds. A total of 18,065 annotated genes were identified, which accounted for 85% of the annotated genes in cattle. The five most abundant hormone genes (GH, PRL, POMC, CGA, and LHB) were relatively stable in both breeds, indicating their pivotal roles in maintaining the basic functions of the pituitary. In addition, 105 genes were identified as differentially expressed between Wenshan and Simmental pituitary glands, including 83 known and 22 novel genes. Functional association analyses of the differentially expressed genes between the breeds revealed 60 enriched gene ontology terms and 3 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Interestingly, the most enriched KEGG pathway, cell adhesion molecules, can modulate pituitary hormone secretion by cell-cell contact. Our findings demonstrated that SYTL2, SLC38A4, and NCAM2 are new candidates for crucial functions in the secretory pathways of the pituitary gland. These results will further understanding of the important molecular mechanisms at work in the bovine pituitary.

Pre- and post-puberty expression of genes and proteins in the uterus of Bos indicus heifers: the luteal phase effect post-puberty

Progesterone signaling and uterine function are crucial in terms of pregnancy establishment. To investigate how the uterine tissue and its secretion changes in relation to puberty, we sampled tissue and uterine fluid from six pre- and six post-pubertal Brahman heifers. Post-pubertal heifers were sampled in the luteal phase. Gene expression of the uterine tissue was investigated with RNA-sequencing, whereas the uterine fluid was used for protein profiling with mass spectrometry. A total of 4034 genes were differentially expressed (DE) at a nominal P-value of 0.05, and 26 genes were significantly DE after Bonferroni correction (P < 3.1 × 10^-6 ). We also identified 79 proteins (out of 230 proteins) that were DE (P < 1 × 10^-5 ) in the uterine fluid. When we compared proteomics and transcriptome results, four DE proteins were identified as being encoded by DE genes: OVGP1, GRP, CAP1 and HBA. Except for CAP1, the other three had lower expression post-puberty. The function of these four genes hypothetically related to preparation of the uterus for a potential pregnancy is discussed in the context of puberty. All DE genes and proteins were also used in pathway and ontology enrichment analyses to investigate overall function. The DE genes were enriched for terms related to ribosomal activity. Transcription factors that were deemed key regulators of DE genes are also reported. Transcription factors ZNF567, ZNF775, RELA, PIAS2, LHX4, SOX2, MEF2C, ZNF354C, HMG20A, TCF7L2, ZNF420, HIC1, GTF3A and two novel genes had the highest regulatory impact factor scores. These data can help to understand how puberty influences uterine function.