Mitochondrial haplotypes influence metabolic traits across bovine inter- and intra-species cybrids

The relationship between mitochondrial DNA haplotype and its copy number on body weight and morphological traits of Wuliangshan black-bone chickens

Mitochondria play a pivotal role as carriers of genetic information through their circular DNA molecules. The rapid evolution of the D-loop region in mitochondria makes it an ideal molecular marker for exploring genetic differentiation among individuals within species and populations with close kinship. However, the influence of mtDNA D-loop region haplotypes and mtDNA copy numbers on phenotypic traits, particularly production traits in chickens, remains poorly understood. In this comprehensive study, we conducted D-loop region amplification and sequencing in the blood mitochondria of 232 female Wuliangshan black-bone chickens. Our investigation identified a total of 38 haplotypes, with a focus on 10 haplotypes that included more than five individuals. We meticulously analyzed the correlations between these haplotypes and a range of traits, encompassing body weight, tibial length, tibial circumference, body oblique length, chest width, and chest depth. The results unveiled significant disparities in specific tested traits across different haplotypes, indicating a tangible association between mtDNA haplotypes and traits in chickens. These findings underscore the potential impact of mitochondrial DNA variations on energy metabolism, ultimately leading to divergent chicken phenotypes. Furthermore, our examination revealed positive correlations between mtDNA copy numbers and tested traits for select haplotypes, while other haplotypes exhibited non-uniform relationships between traits and mtDNA copy numbers. In addition, phylogenetic analysis disclosed the involvement of two subspecies of red jungle chicken in the origin of Wuliangshan black-bone chickens. Consequently, our research contributes novel insights into mitochondrial genomic selection, augments comprehension of the roles played by haplotypes and mtDNA copy numbers in chicken population genetics and phylogenetic analysis, and furnishes fundamental data crucial for the preservation and provenance determination of black-bone chickens.

Maternal Genetic Diversity Analysis of Guanling Cattle by Mitochondrial Genome Sequencing

Owing to the unique geography and the isolated environment, Guanling cattle, which is one of five local cattle breeds in Guizhou, China, has developed unique characteristics. The number of pure Guanling cattle decreased markedly because of the hybridization with foreign breeds. In the present study, the maternal genetic diversity of 58 Guanling bulls was assessed by whole mitochondrial genome sequencing. Genetic polymorphisms and phylogenetic analyses classified Guanling cattle into two main lineages, where 43.10% of Guanling cattle were closely related to the foreign breeds and 56.90% displayed distinct features in mitochondrial genomic diversity. PCA analysis further separated Guanling cattle into four populations, one of which was clustered with the foreign breeds. The result of the structure plot and genetic polymorphisms revealed high genetic diversity within two populations that have a long genetic distance from the foreign breeds. Overall, our findings suggest that the whole mitochondrial genome sequencing analysis is a useful and reliable tool to study maternal genetic diversity and to identify the pure population of Guanling cattle. The results will be beneficial to the breeding management of Guanling cattle.

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.

Mitochondrial DNA D-loop SNPs unveil molecular signatures of milk production variation in Murrah buffalo

BACKGROUND

The Murrah buffalo, pivotal in Asian agriculture, faces challenges in maximizing milk production despite significant breeding efforts. Recognizing its economic importance, this study investigates mtDNA D-loop variations in Murrah buffalo as potential indicators of milk production variability, addressing challenges in maximizing yield despite significant breeding efforts.

METHODS AND RESULTS

Analyzing mtDNA D-loop sequences from 50 buffaloes, we categorized them into Low (Group 1), Medium (Group 2), and High ECM (Group 3) groups based on milk yields, fat and protein percentage of a 30-day period data. Somatic cell mtDNA D-loop analysis revealed distinct genetic variations, with significant differences among ECM groups. Group 2 showed higher SNP prevalence, group 3 had more insertions/deletions, and Group 1 exhibited the highest transition frequency. Notably, a consistent "C" deletion at the 714th position occurred in Groups 1 and 3, prevalent in 68% of Group 2. A G-A variation at the 93rd position was specific to the medium ECM group. Negative Tajima D values indicated unique variations in each group, with Group 1 having the highest number, and a specific SNP linked to Group 2 was identified. These SNPs in the D-loop region could impact mtDNA replication, influencing mitochondrial content among animals. Our results provide valuable insights into the role of mtDNA D-loop polymorphisms in milk production traits in Murrah buffalo.

CONCLUSIONS

Our research highlights the potential for valuable markers of cellular energy efficiency in Murrah buffalo. Exploring diverse cytoplasmic backgrounds opens avenues for mtDNA-based selection strategies, enhancing milk production and optimizing genetic traits for the dairy industry.

Primary Cell Culture as a Model System for Evolutionary Molecular Physiology

Primary cell culture is a powerful model system to address fundamental questions about organismal physiology at the cellular level, especially for species that are difficult, or impossible, to study under natural or semi-natural conditions. Due to their ease of use, primary fibroblast cultures are the dominant model system, but studies using both somatic and germ cells are also common. Using these models, genome evolution and phylogenetic relationships, the molecular and biochemical basis of differential longevities among species, and the physiological consequences of life history evolution have been studied in depth. With the advent of new technologies such as gene editing and the generation of induced pluripotent stem cells (iPSC), the field of molecular evolutionary physiology will continue to expand using both descriptive and experimental approaches.

Association analysis of mitochondrial genome polymorphisms with backfat thickness in pigs

Mitochondrial DNA (mtDNA) variations and associated effects on economic traits have been widely reported in farm animals, as these genetic polymorphisms can affect the efficiency of energy production and cell metabolism. In studies related to metabolism, the deposition of fat was highly correlated with mitochondria. However, the effect of mtDNA polymorphisms on porcine backfat thickness (BFT) remained unclear. In this study, 243 pigs were collected to analyse the relationship between BFT and mtDNA polymorphisms. There were considerable differences in BFT, ranging from 5 mm to 18 mm. MtDNA D-loop sequencing discovered 48 polymorphic sites. Association analysis revealed that 30 variations were associated with BFT (P < 0.05). The polymorphism m.794A > G showed the maximum difference in BFT between A and G carriers, which differed at ∼2.5 mm (P < 0.001). The 48 polymorphic sites generated 22 haplotypes (H1-H22), which clustered into 4 haplogroups (HG1-HG4). HG1 had a lower BFT value than other three haplogroups (P < 0.01), whereas H4 in HG1 exhibited the lowest BFT of all haplotypes analyzed (P < 0.01). The results of this study highlight an association between mtDNA polymorphisms and BFT, and suggest the potential application of mtDNA in pig molecular breeding practices.

Mitochondrial oxygen respiration is associated with loneliness in a sample of community-dwelling adults

Loneliness, a distressing perception of insufficient social support, is associated with physical illness and premature mortality that may be explained by reduced mitochondrial efficiency. In this human study (n = 50), loneliness was associated with less efficient mitochondrial functioning, indicated by the bioenergetic health index (BHI; r = -.39, p = .009), coupling efficiency (r = -.34, p = .021), and phosphorylating respiration (r = .39, p = .009). These findings remained significant when controlling for age, sex, and body mass index. The association between loneliness and mitochondrial functioning is important, given that both have been associated with age-related diseases.

Genome-wide local ancestry and evidence for mitonuclear coadaptation in African hybrid cattle populations (Bos taurus/indicus)

The phenotypic diversity of African cattle reflects adaptation to a wide range of agroecological conditions, human-mediated selection preferences, and complex patterns of admixture between the humpless Bos taurus (taurine) and humped Bos indicus (zebu) subspecies, which diverged 150-500 thousand years ago. Despite extensive admixture, all African cattle possess taurine mitochondrial haplotypes, even populations with significant zebu biparental and male uniparental nuclear ancestry. This has been interpreted as the result of human-mediated dispersal ultimately stemming from zebu bulls imported from South Asia during the last three millennia. Here, we assess whether ancestry at mitochondrially targeted nuclear genes in African admixed cattle is impacted by mitonuclear functional interactions. Using high-density SNP data, we find evidence for mitonuclear coevolution across hybrid African cattle populations with a significant increase of taurine ancestry at mitochondrially targeted nuclear genes. Our results, therefore, support the hypothesis of incompatibility between the taurine mitochondrial genome and the zebu nuclear genome.

Recovery of mitogenomes from whole genome sequences to infer maternal diversity in 1883 modern taurine and indicine cattle

Maternal diversity based on a sub-region of mitochondrial genome or variants were commonly used to understand past demographic events in livestock. Additionally, there is growing evidence of direct association of mitochondrial genetic variants with a range of phenotypes. Therefore, this study used complete bovine mitogenomes from a large sequence database to explore the full spectrum of maternal diversity. Mitogenome diversity was evaluated among 1883 animals representing 156 globally important cattle breeds. Overall, the mitogenomes were diverse: presenting 11 major haplogroups, expanding to 1309 unique haplotypes, with nucleotide diversity 0.011 and haplotype diversity 0.999. A small proportion of African taurine (3.5%) and indicine (1.3%) haplogroups were found among the European taurine breeds and composites. The haplogrouping was largely consistent with the population structure derived from alternate clustering methods (e.g. PCA and hierarchical clustering). Further, we present evidence confirming a new indicine subgroup (I1a, 64 animals) mainly consisting of breeds originating from China and characterised by two private mutations within the I1 haplogroup. The total genetic variation was attributed mainly to within-breed variance (96.9%). The accuracy of the imputation of missing genotypes was high (99.8%) except for the relatively rare heteroplasmic genotypes, suggesting the potential for trait association studies within a breed.

Mitonuclear incompatibility as a hidden driver behind the genome ancestry of African admixed cattle

Background

Africa is an important watershed in the genetic history of domestic cattle, as two lineages of modern cattle, Bos taurus and B. indicus, form distinct admixed cattle populations. Despite the predominant B. indicus nuclear ancestry of African admixed cattle, B. indicus mitochondria have not been found on the continent. This discrepancy between the mitochondrial and nuclear genomes has been previously hypothesized to be driven by male-biased introgression of Asian B. indicus into ancestral African B. taurus. Given that this hypothesis mandates extreme demographic assumptions relying on random genetic drift, we propose a novel hypothesis of selection induced by mitonuclear incompatibility and assess these hypotheses with regard to the current genomic status of African admixed cattle.

Results

By analyzing 494 mitochondrial and 235 nuclear genome sequences, we first confirmed the genotype discrepancy between mitochondrial and nuclear genome in African admixed cattle: the absence of B. indicus mitochondria and the predominant B. indicus autosomal ancestry. We applied approximate Bayesian computation (ABC) to assess the posterior probabilities of two selection hypotheses given this observation. The results of ABC indicated that the model assuming both male-biased B. indicus introgression and selection induced by mitonuclear incompatibility explains the current genomic discrepancy most accurately. Subsequently, we identified selection signatures at autosomal loci interacting with mitochondria that are responsible for integrity of the cellular respiration system. By contrast with B. indicus-enriched genome ancestry of African admixed cattle, local ancestries at these selection signatures were enriched with B. taurus alleles, concurring with the key expectation of selection induced by mitonuclear incompatibility.

Conclusions

Our findings support the current genome status of African admixed cattle as a potential outcome of male-biased B. indicus introgression, where mitonuclear incompatibility exerted selection pressure against B. indicus mitochondria. This study provides a novel perspective on African cattle demography and supports the role of mitonuclear incompatibility in the hybridization of mammalian species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01206-x.

Cybrid Model Supports Mitochondrial Genetic Effect on Pig Litter Size

In pigs, mitochondrial DNA (mtDNA) polymorphism and the correlation to reproductive performance across breeds and individuals have been largely reported, however, experimental proof has never been provided. In this study, we analyzed 807 sows for correlation of total number born (TNB) and mitotype, which presented the maximum of 1.73 piglets for mtDNA contribution. Cybrid models representing different mitotypes were generated for identification of the mtDNA effect. Results indicated significant differences on cellular and molecular characteristics among cybrids, including energy metabolic traits, mtDNA copy numbers and transcriptions, mRNA and protein expressions on mitochondrial biogenesis genes and reproduction-related genes. Referring to mitotypes, the cybrids with prolific mitotypes presented significantly higher oxygen consumption rate (OCR) productions, mtDNA transcriptions and copy numbers than those with common mitotypes, while both mRNA and protein expressions of PPARA, TFAM, ER1, ER2, and ESRRG in prolific cybrids were significantly higher than those with common mitotypes. Cybrid models reflected the mtDNA effect on pig litter size, suggesting the potential application of mtDNA polymorphism in pig selection and breeding practices.

Genome-Wide Identification of Discriminative Genetic Variations in Beef and Dairy Cattle via an Information-Theoretic Approach

Analyzing the associations between genotypic changes and phenotypic traits on a genome-wide scale can contribute to understanding the functional roles of distinct genetic variations during breed development. We performed a whole-genome analysis of Angus and Jersey cattle breeds using conditional mutual information, which is an information-theoretic method estimating the conditional independency among multiple factor variables. The proposed conditional mutual information-based approach allows breed-discriminative genetic variations to be explicitly identified from tens of millions of SNP (single nucleotide polymorphism) positions on a genome-wide scale while minimizing the usage of prior knowledge. Using this data-driven approach, we identified biologically relevant functional genes, including breed-specific variants for cattle traits such as beef and dairy production. The identified lipid-related genes were shown to be significantly associated with lipid and triglyceride metabolism, fat cell differentiation, and muscle development. In addition, we confirmed that milk-related genes are involved in mammary gland development, lactation, and mastitis-associated processes. Our results provide the distinct properties of Angus and Jersey cattle at a genome-wide level. Moreover, this study offers important insights into discovering unrevealed genetic variants for breed-specific traits and the identification of genetic signatures of diverse cattle breeds with respect to target breed-specific properties.

Mitochondrial DNA Haplotypes Influence Energy Metabolism across Chicken Transmitochondrial Cybrids

The association between mitochondrial DNA haplotype and productive performances has been widely reported in chicken breeds. However, there has not been physiological evidence of this seen previously. In this study, chicken transmitochondrial cells were generated using the nucleus of the DF-1 cell line and mitochondria of primary cell lines derived from two native chicken breeds, Tibetan chicken and Shouguang chicken. Generally, Tibetan chicken primary cells showed a stronger metabolic capacity than Shouguang chicken primary cells. However, the Tibetan chicken cybrids had a dramatic drop in relative mtDNA copies and oxygen consumption. Higher rates of oxygen consumption (OCR) and expression levels of mitochondrial biogenesis and fusion genes were observed in Shouguang chicken cybrids, potentially reflecting that the mitochondrial DNA haplotype of Shouguang chicken had better coordination with the DF-1 nucleus than others. Meanwhile, mitonuclear incompatibility occurred in Tibetan chicken cybrids. The results demonstrate functional differences among mitochondrial DNA haplotypes and may shed light on the interaction between the mitochondria and nucleus in Gallus gallus domesticus.

Inter-individual variation in adaptive capacity at onset of lactation: Linking metabolic phenotype with mitochondrial DNA haplotype in Holstein dairy cows

Modern-day dairy cows express great variation in metabolic capacity to adapt to the onset of lactation. Although breeding programs increased the breeding value for longevity and robustness in the sires, a respective phenotype in female offspring has not been improving as predicted. Fundamental energy generating pathways such as mitochondrial fatty acid oxidation might have a crucial role for robustness and metabolic efficiency in dairy cows. Therefore, mitochondrial enzyme expression was examined in liver samples of one set of animals before and after calving. Furthermore, the mitochondrial DNA sequence was determined for each individual of a second set of animals using liver samples. Results from the first trial indicate that the expression and extent of phosphorylation of acetyl-CoA carboxylase (ACC) is the major key step for modulating fatty acid flux into the mitochondria at early onset of lactation in Holstein dairy cows. In the second trial, mitochondrial DNA sequencing and identification of mutation patterns yielded three major haplotypes. Haplotype H2 was closely associated with liver fat content, plasma glycerol and acyl-carnitine concentrations. The mitochondrial DNA haplotype, which is a feature of the maternal lines, might be related to the inter-individual variation in metabolic capacity of Holstein dairy cows.

Mitochondrial Recombination Reveals Mito-Mito Epistasis in Yeast

Genetic variation in mitochondrial DNA (mtDNA) provides adaptive potential although the underlying genetic architecture of fitness components within mtDNAs is not known. To dissect functional variation within mtDNAs, we first identified naturally occurring mtDNAs that conferred high or low fitness in Saccharomyces cerevisiae by comparing growth in strains containing identical nuclear genotypes but different mtDNAs. During respiratory growth under temperature and oxidative stress conditions, mitotype effects were largely independent of nuclear genotypes even in the presence of mito-nuclear interactions. Recombinant mtDNAs were generated to determine fitness components within high- and low-fitness mtDNAs. Based on phenotypic distributions of isogenic strains containing recombinant mtDNAs, we found that multiple loci contributed to mitotype fitness differences. These mitochondrial loci interacted in epistatic, nonadditive ways in certain environmental conditions. Mito-mito epistasis (i.e., nonadditive interactions between mitochondrial loci) influenced fitness in progeny from four different crosses, suggesting that mito-mito epistasis is a widespread phenomenon in yeast and other systems with recombining mtDNAs. Furthermore, we found that interruption of coadapted mito-mito interactions produced recombinant mtDNAs with lower fitness. Our results demonstrate that mito-mito epistasis results in functional variation through mitochondrial recombination in fungi, providing modes for adaptive evolution and the generation of mito-mito incompatibilities.