Expression of mitochondrial protein genes encoded by nuclear and mitochondrial genomes correlate with energy metabolism in dairy cattle

Nuclear Genome-Encoded Mitochondrial OXPHOS Complex I Genes in Female Buffalo Show Tissue-Specific Differences

Buffalo physiology intricately balances energy, profoundly influencing health, productivity, and reproduction. This study explores nuclear-mitochondrial crosstalk, revealing OXPHOS Complex I gene expression variations in buffalo tissues through high-throughput RNA sequencing. Unveiling tissue-specific disparities, the research elucidates the genomic landscape of crucial energy production genes, with broader implications for veterinary and agricultural progress. Post-slaughter, tissues from post-pubertal female buffaloes underwent meticulous processing and RNA extraction using the TRIzol method. RNA-Seq library preparation and IlluminaHiSeq 2500 sequencing were performed on QC-passed samples. Data underwent stringent filtration, mapping to the Bubalus bubalis genome using HISAT2. DESeq2 facilitated differential expression gene (DEG) analysis focusing on 57 Mitocarta 3-derived genes associated with OXPHOS complex I. Nuclear-encoded mitochondrial protein transcripts of OXPHOS complex 1 exhibited tissue-specific variations, with 51 genes expressing significantly across tissues. DEG analysis emphasized tissue-specific expression patterns, highlighting a balanced OXPHOS complex I subunit expression in the kidney vs. brain. Gene Ontology (GO) enrichment showcased mitochondria-centric terms, revealing distinct proton motive force-driven mitochondrial ATP synthesis regulation. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses emphasized Thermogenesis and OXPHOS pathways, enriching our understanding of tissue-specific energy metabolism. Noteworthy up-regulation of NDUFB10 in the heart and kidney aligned with heightened metabolic activity. Brain-specific up-regulation of NDUFAF6 indicated a focus on mitochondrial function, while variations in NDUFA11 and ACAD9 underscored pivotal roles in the heart and kidney. GO and KEGG analyses highlighted tissue-specific mitochondrial ATP synthesis and NADH dehydrogenase processes, providing molecular insights into organ-specific metabolic demands and regulatory mechanisms. Our study unveils conserved and tissue-specific nuances in nuclear-encoded mitochondrial OXPHOS complex I genes, laying a foundation for understanding diverse energy demands and potential health implications.

Differential Expression of Nuclear-Encoded Mitochondrial Protein Genes of ATP Synthase Across Different Tissues of Female Buffalo

The physiological well-being of buffaloes, encompassing phenotypic traits, reproductive health, and productivity, depends on their energy status. Mitochondria, the architects of energy production, orchestrate a nuanced interplay between nuclear and mitochondrial domains. Oxidative phosphorylation complexes and associated proteins wield significant influence over metabolic functions, energy synthesis, and organelle dynamics, often linked to tissue-specific pathologies. The unexplored role of ATP synthase in buffalo tissues prompted a hypothesis: in-depth exploration of nuclear-derived mitochondrial genes, notably ATP synthase, reveals distinctive tissue-specific diversity. RNA extraction and sequencing of buffalo tissues (kidney, heart, brain, and ovary) enabled precise quantification of nuclear-derived mitochondrial protein gene expression. The analysis unveiled 24 ATP synthase transcript variants, each with unique tissue-specific patterns. Kidney, brain, and heart exhibited elevated gene expression compared to ovaries, with 10, 8, and 19 up-regulated genes, respectively. The kidney showed 3 and 12 down-regulated genes compared to the brain and heart. The heart-brain comparison highlighted ten highly expressed genes in ATP synthase functions. Gene ontology and pathway analyses revealed enriched functions linked to ATP synthesis and oxidative phosphorylation, offering a comprehensive understanding of energy production in buffalo tissues. This analysis enhances understanding of tissue-specific gene expression, emphasizing the influence of energy demands. Revealing intricate links between mitochondrial gene expression and tissue specialization in buffaloes, it provides nuanced insights into tissue-specific expression of nuclear-encoded mitochondrial protein genes, notably ATP synthase, advancing the comprehension of buffalo tissue biology.

Quantifying the effects of the mitochondrial genome on milk production traits in dairy cows: Empirical results and modeling challenges

Substantial advances in livestock traits have been achieved primarily through selection strategies targeting variation in the nuclear genome, with little attention given to mitogenome variation. We analyzed the influence of the mitogenome on milk production traits of Holstein cattle in Croatia based on strategically generated next-generation sequencing data for 109 cows pedigree-linked to 7,115 milk production records (milk, fat, and protein yield) from 3,006 cows (first 5 lactations). Because little is known about the biology of the relationship between mitogenome variation and production traits, our quantitative genetic modeling was complex. Thus, the proportion of total variance explained by mitogenome inheritance was estimated using 5 different models: (1) a cytoplasmic model with maternal lineages (CYTO), (2) a haplotypic model with mitogenome sequences (HAPLO), (3) an amino acid model with unique amino acid sequences (AMINO), (4) an evolutionary model based on a phylogenetic analysis using Bayesian Evolutionary Analysis Sampling Trees phylogenetic analysis (EVOL), and (5) a mitogenome SNP model (SNPmt). The polygenic autosomal and X chromosome additive genetic effects based on pedigree were modeled, together with the effects of herd-year-season interaction, permanent environment, location, and age at first calving. The estimated proportions of phenotypic variance explained by mitogenome in 4 different models (CYTO, HAPLO, AMINO, and SNPmt) were found to be substantial given the size of mitogenome, ranging from 5% to 7% for all 3 milk traits. At the same time, a negligible proportion of the phenotypic variance was explained by mitogenome with the EVOL model. Similarly, in all models, no proportion of phenotypic variance was explained by the X chromosome. Although our results should be confirmed in other dairy cattle populations, including a large number of sequenced mitogenomes and nuclear genomes, the potential of utilizing mitogenome information in animal breeding is promising, especially as the acquisition of complete genome sequences becomes cost-effective.

Mitochondrial sequence variants: testing imputation accuracy and their association with dairy cattle milk traits

BACKGROUND

Mitochondrial genomes differ from the nuclear genome and in humans it is known that mitochondrial variants contribute to genetic disorders. Prior to genomics, some livestock studies assessed the role of the mitochondrial genome but these were limited and inconclusive. Modern genome sequencing provides an opportunity to re-evaluate the potential impact of mitochondrial variation on livestock traits. This study first evaluated the empirical accuracy of mitochondrial sequence imputation and then used real and imputed mitochondrial sequence genotypes to study the role of mitochondrial variants on milk production traits of dairy cattle.

RESULTS

The empirical accuracy of imputation from Single Nucleotide Polymorphism (SNP) panels to mitochondrial sequence genotypes was assessed in 516 test animals of Holstein, Jersey and Red breeds using Beagle software and a sequence reference of 1883 animals. The overall accuracy estimated as the Pearson's correlation squared (R2) between all imputed and real genotypes across all animals was 0.454. The low accuracy was attributed partly to the majority of variants having low minor allele frequency (MAF < 0.005) but also due to variants in the hypervariable D-loop region showing poor imputation accuracy. Beagle software provides an internal estimate of imputation accuracy (DR2), and 10 percent of the total 1927 imputed positions showed DR2 greater than 0.9 (N = 201). There were 151 sites with empirical R2 > 0.9 (of 954 variants segregating in the test animals) and 138 of these overlapped the sites with DR2 > 0.9. This suggests that the DR2 statistic is a reasonable proxy to select sites that are imputed with higher accuracy for downstream analyses. Accordingly, in the second part of the study mitochondrial sequence variants were imputed from real mitochondrial SNP panel genotypes of 9515 Australian Holstein, Jersey and Red dairy cattle. Then, using only sites with DR2 > 0.900 and real genotypes, we undertook a genome-wide association study (GWAS) for milk, fat and protein yields. The GWAS mitochondrial SNP effects were not significant.

CONCLUSION

The accuracy of imputation of mitochondrial genotypes from the SNP panel to sequence was generally low. The Beagle DR2 statistic enabled selection of sites imputed with higher empirical accuracy. We recommend building larger reference populations with mitochondrial sequence to improve the accuracy of imputing less common variants and ensuring that SNP panels include common variants in the D-loop region.

Imprinting at the KBTBD6 locus involves species-specific maternal methylation and monoallelic expression in livestock animals

BACKGROUND

The primary differentially methylated regions (DMRs) which are maternally hypermethylated serve as imprinting control regions (ICRs) that drive monoallelic gene expression, and these ICRs have been investigated due to their implications in mammalian development. Although a subset of genes has been identified as imprinted, in-depth comparative approach needs to be developed for identification of species-specific imprinted genes. Here, we examined DNA methylation status and allelic expression at the KBTBD6 locus across species and tissues and explored potential mechanisms of imprinting.

RESULTS

Using whole-genome bisulfite sequencing and RNA-sequencing on parthenogenetic and normal porcine embryos, we identified a maternally hypermethylated DMR between the embryos at the KBTBD6 promoter CpG island and paternal monoallelic expression of KBTBD6. Also, in analyzed domesticated mammals but not in humans, non-human primates and mice, the KBTBD6 promoter CpG islands were methylated in oocytes and/or allelically methylated in tissues, and monoallelic KBTBD6 expression was observed, indicating livestock-specific imprinting. Further analysis revealed that these CpG islands were embedded within transcripts in porcine and bovine oocytes which coexisted with an active transcription mark and DNA methylation, implying the presence of transcription-dependent imprinting.

CONCLUSIONS

In this study, our comparative approach revealed an imprinted expression of the KBTBD6 gene in domesticated mammals, but not in humans, non-human primates, and mice which implicates species-specific evolution of genomic imprinting.

Allele specific binding of histone modifications and a transcription factor does not predict allele specific expression in correlated ChIP-seq peak-exon pairs

Allele specific expression (ASE) is widespread in many species including cows. Therefore, regulatory regions which control gene expression should show cis-regulatory variation which mirrors this differential expression within the animal. ChIP-seq peaks for histone modifications and transcription factors measure activity at functional regions and the height of some peaks have been shown to correlate across tissues with the expression of particular genes, suggesting these peaks are putative regulatory regions. In this study we identified ASE in the bovine genome in multiple tissues and investigated whether ChIP-seq peaks for four histone modifications and the transcription factor CTCF show allele specific binding (ASB) differences in the same tissues. We then investigate whether peak height and gene expression, which correlates across tissues, also correlates within the animal by investigating whether the direction of ASB in putative regulatory regions, mirrors that of the ASE in the genes they are putatively regulating. We found that ASE and ASB were widespread in the bovine genome but vary in extent between tissues. However, even when the height of a peak was positively correlated across tissues with expression of an exon, ASE of the exon and ASB of the peak were in the same direction only half the time. A likely explanation for this finding is that the correlations between peak height and exon expression do not indicate that the height of the peak causes the extent of exon expression, at least in some cases.

AgAnimalGenomes: browsers for viewing and manually annotating farm animal genomes

Current genome sequencing technologies have made it possible to generate highly contiguous genome assemblies for non-model animal species. Despite advances in genome assembly methods, there is still room for improvement in the delineation of specific gene features in the genomes. Here we present genome visualization and annotation tools to support seven livestock species (bovine, chicken, goat, horse, pig, sheep, and water buffalo), available in a new resource called AgAnimalGenomes. In addition to supporting the manual refinement of gene models, these browsers provide visualization tracks for hundreds of RNAseq experiments, as well as data generated by the Functional Annotation of Animal Genomes (FAANG) Consortium. For species with predicted gene sets from both Ensembl and RefSeq, the browsers provide special tracks showing the thousands of protein-coding genes that disagree across the two gene sources, serving as a valuable resource to alert researchers to gene model issues that may affect data interpretation. We describe the data and search methods available in the new genome browsers and how to use the provided tools to edit and create new gene models.

Tissue heterogeneity of mitochondrial activity, biogenesis and mitochondrial protein gene expression in buffalo

BACKGROUND

Cellular metabolism is most invariant process, occurring in all living organisms, which involves mitochondrial proteins from both nuclear and mitochondrial genomes. The mitochondrial DNA (mtDNA) copy number, protein-coding genes (mtPCGs) expression, and activity vary between various tissues to fulfill specific energy demands across the tissues.

METHODS AND RESULTS

In present study, we investigated the OXPHOS complexes and citrate synthase activity in isolated mitochondria from various tissues of freshly slaughtered buffaloes (n = 3). Further, the evaluation of tissue-specific diversity based on the quantification of mtDNA copy numbers was performed and also comprised an expression study of 13 mtPCGs. We found that the functional activity of individual OXPHOS complex I was significantly higher in the liver compared to muscle and brain. Additionally, OXPHOS complex III and V activities was observed significantly higher levels in liver compared to heart, ovary, and brain. Similarly, CS-specific activity differs between tissues, with the ovary, kidney, and liver having significantly greater. Furthermore, we revealed the mtDNA copy number was strictly tissue-specific, with muscle and brain tissues exhibiting the highest levels. Among 13 PCGs expression analyses, mRNA abundances in all genes were differentially expressed among the different tissue.

CONCLUSIONS

Overall, our results indicate the existence of a tissue-specific variation in mitochondrial activity, bioenergetics, and mtPCGs expression among various types of buffalo tissues. This study serves as a critical first stage in gathering vital comparable data about the physiological function of mitochondria in energy metabolism in distinct tissues, laying the groundwork for future mitochondrial based diagnosis and research.

Mitochondrial Transfer in Cardiovascular Disease: From Mechanisms to Therapeutic Implications

Mitochondrial dysfunction has been proven to play a critical role in the pathogenesis of cardiovascular diseases. The phenomenon of intercellular mitochondrial transfer has been discovered in the cardiovascular system. Studies have shown that cell-to-cell mitochondrial transfer plays an essential role in regulating cardiovascular system development and maintaining normal tissue homeostasis under physiological conditions. In pathological conditions, damaged cells transfer dysfunctional mitochondria toward recipient cells to ask for help and take up exogenous functional mitochondria to alleviate injury. In this review, we summarized the mechanism of mitochondrial transfer in the cardiovascular system and outlined the fate and functional role of donor mitochondria. We also discussed the advantage and challenges of mitochondrial transfer strategies, including cell-based mitochondrial transplantation, extracellular vesicle-based mitochondrial transplantation, and naked mitochondrial transplantation, for the treatment of cardiovascular disorders. We hope this review will provide perspectives on mitochondrial-targeted therapeutics in cardiovascular diseases.

Feather Gene Expression Elucidates the Developmental Basis of Plumage Iridescence in African Starlings

Iridescence is widespread in the living world, occurring in organisms as diverse as bacteria, plants, and animals. Yet, compared to pigment-based forms of coloration, we know surprisingly little about the developmental and molecular bases of the structural colors that give rise to iridescence. Birds display a rich diversity of iridescent structural colors that are produced in feathers by the arrangement of melanin-containing organelles called melanosomes into nanoscale configurations, but how these often unusually shaped melanosomes form, or how they are arranged into highly organized nanostructures, remains largely unknown. Here, we use functional genomics to explore the developmental basis of iridescent plumage using superb starlings (Lamprotornis superbus), which produce both iridescent blue and non-iridescent red feathers. Through morphological and chemical analyses, we confirm that hollow, flattened melanosomes in iridescent feathers are eumelanin-based, whereas melanosomes in non-iridescent feathers are solid and amorphous, suggesting that high pheomelanin content underlies red coloration. Intriguingly, the nanoscale arrangement of melanosomes within the barbules was surprisingly similar between feather types. After creating a new genome assembly, we use transcriptomics to show that non-iridescent feather development is associated with genes related to pigmentation, metabolism, and mitochondrial function, suggesting non-iridescent feathers are more energetically expensive to produce than iridescent feathers. However, iridescent feather development is associated with genes related to structural and cellular organization, suggesting that, while nanostructures themselves may passively assemble, barbules and melanosomes may require active organization to give them their shape. Together, our analyses suggest that iridescent feathers form through a combination of passive self-assembly and active processes.

Putative Causal Variants Are Enriched in Annotated Functional Regions From Six Bovine Tissues

Genetic variants which affect complex traits (causal variants) are thought to be found in functional regions of the genome. Identifying causal variants would be useful for predicting complex trait phenotypes in dairy cows, however, functional regions are poorly annotated in the bovine genome. Functional regions can be identified on a genome-wide scale by assaying for post-translational modifications to histone proteins (histone modifications) and proteins interacting with the genome (e.g., transcription factors) using a method called Chromatin immunoprecipitation followed by sequencing (ChIP-seq). In this study ChIP-seq was performed to find functional regions in the bovine genome by assaying for four histone modifications (H3K4Me1, H3K4Me3, H3K27ac, and H3K27Me3) and one transcription factor (CTCF) in 6 tissues (heart, kidney, liver, lung, mammary and spleen) from 2 to 3 lactating dairy cows. Eighty-six ChIP-seq samples were generated in this study, identifying millions of functional regions in the bovine genome. Combinations of histone modifications and CTCF were found using ChromHMM and annotated by comparing with active and inactive genes across the genome. Functional marks differed between tissues highlighting areas which might be particularly important to tissue-specific regulation. Supporting the cis-regulatory role of functional regions, the read counts in some ChIP peaks correlated with nearby gene expression. The functional regions identified in this study were enriched for putative causal variants as seen in other species. Interestingly, regions which correlated with gene expression were particularly enriched for potential causal variants. This supports the hypothesis that complex traits are regulated by variants that alter gene expression. This study provides one of the largest ChIP-seq annotation resources in cattle including, for the first time, in the mammary gland of lactating cows. By linking regulatory regions to expression QTL and trait QTL we demonstrate a new strategy for identifying causal variants in cattle.

Mitochondrial protein gene expression and the oxidative phosphorylation pathway associated with feed efficiency and energy balance in dairy cattle

Feed efficiency and energy balance are important traits underpinning profitability and environmental sustainability in animal production. They are complex traits, and our understanding of their underlying biology is currently limited. One measure of feed efficiency is residual feed intake (RFI), which is the difference between actual and predicted intake. Variation in RFI among individuals is attributable to the metabolic efficiency of energy utilization. High RFI (H_RFI) animals require more energy per unit of weight gain or milk produced compared with low RFI (L_RFI) animals. Energy balance (EB) is a closely related trait calculated very similarly to RFI. Cellular energy metabolism in mitochondria involves mitochondrial protein (MiP) encoded by both nuclear (NuMiP) and mitochondrial (MtMiP) genomes. We hypothesized that MiP genes are differentially expressed (DE) between H_RFI and L_RFI animal groups and similarly between negative and positive EB groups. Our study aimed to characterize MiP gene expression in white blood cells of H_RFI and L_RFI cows using RNA sequencing to identify genes and biological pathways associated with feed efficiency in dairy cattle. We used the top and bottom 14 cows ranked for RFI and EB out of 109 animals as H_RFI and L_RFI, and positive and negative EB groups, respectively. The gene expression counts across all nuclear and mitochondrial genes for animals in each group were used for differential gene expression analyses, weighted gene correlation network analysis, functional enrichment, and identification of hub genes. Out of 244 DE genes between RFI groups, 38 were MiP genes. The DE genes were enriched for the oxidative phosphorylation (OXPHOS) and ribosome pathways. The DE MiP genes were underexpressed in L_RFI (and negative EB) compared with the H_RFI (and positive EB) groups, suggestive of reduced mitochondrial activity in the L_RFI group. None of the MtMiP genes were among the DE MiP genes between the groups, which suggests a non-rate limiting role of MtMiP genes in feed efficiency and warrants further investigation. The role of MiP, particularly the NuMiP and OXPHOS pathways in RFI, was also supported by our gene correlation network analysis and the hub gene identification. We validated the findings in an independent data set. Overall, our study suggested that differences in feed efficiency in dairy cows may be linked to differences in cellular energy demand. This study broadens our knowledge of the biology of feed efficiency in dairy cattle.