Mitochondrial DNA T7719G in tRNA-Lys gene affects litter size in Small-tailed Han sheep

Mitogenome-wide association study on body measurement traits of Wenshang Barred chickens

Mitochondria are best known for synthesizing ATP through the tricarboxylic acid cycle and oxidative phosphorylation. The cytoplasmic mitochondrial DNA (mtDNA) is important for maintaining the function. This study was designed to reveal the effect of mtDNA on chicken body measurement traits (BMTs). A population of 605 Wenshang Barred chickens were recorded BMTs, including body slope length, keel length, chest width, etc. The single-nucleotide polymorphisms (SNPs) of their mitogenomes were detected by PCR amplification and DNA sequencing. Totally 69 mutations in mitogenome were discovered, including 18 in noncoding region and 51 in coding region. By multi-sequence alignment and haplotype construction, the chickens were clustered into eight haplotypes and further three haplogroups. The association between BMTs and mtDNA SNPs, haplotypes and haplogroups were analyzed in the linear model by ASReml, respectively. Among them, the SNP mt11086 T/C in ND3 was found to significantly affect chest dept (p < .05) and was highly conservative by phylogenetic conservation analyses, which reflected the genetic effect on body size and growth of chickens. No significant association between the mitochondrial haplotypes or haplogroups and BMTs was found. The polymorphic site reflecting body size could be put into chicken breeding programs as the genetic marker.

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.

Comparative Hypothalamic Transcriptome Analysis Reveals Crucial mRNAs, lncRNAs, and circRNAs Affecting Litter Size in Goats

Litter size is an important indicator to measure the reproductive performance of goats, which is affected by the reproductive function of animals. The hypothalamus, as the regulatory center of the endocrine system, plays an important role in the reproduction of female animals. Here, we performed high-throughput RNA sequencing using hypothalamic tissue from high-fecundity and low-fecundity Leizhou goats to explore critical functional genes associated with litter size. Differentially expressed mRNA, lncRNA, and circRNAs were screened using DESeq and were enriched, and then analyzed by Gene Ontology and Kyoto Encyclopedia of Gene and Genome. Results showed that some of these differentially expressed mRNAs could be enriched in reproductive processes, jak-STAT, prolactin signaling pathway, and other signaling pathways related to reproduction, such as SOCS3. Furthermore, the central proteins POSTN, MFAP5, and DCN from protein-protein interaction may regulate animal reproductive activity by affecting cell proliferation and apoptosis. lncRNA MSTRG.33887.2 as well as circRNAs chicirc_098002, chicirc_072583, and chicirc_053531 may be able to influence animal reproduction by participating in folate metabolism and energy metabolism homeostasis through their respective target genes. Our results expand the molecular mechanism of hypothalamic regulation on animal reproduction.

The ND1 and CYTB genes polymorphisms associated with in vitro early embryo development of Sanjabi sheep

This study aimed to investigate the association between polymorphisms of ND1 and CYTB genes and in vitro early embryo development of Sanjabi sheep. Blood and ovarian samples were collected from a local slaughterhouse. The cumulus-oocyte complexes with a diameter greater than 3 mm were aspirated from follicles, and in vitro maturation (IVM) and in vitro culture (IVC) rates of them were recorded. A respective 1200 bp and 980 bp fragments of ND1 and CYTB genes were genotyped using a modified single strand conformation polymorphism (SSCP) method. The results of this study revealed that four different patterns, named as A, B, C, and D were observed for both ND1 and CYTB genes. The ND1 gene polymorphisms had significant effects on the IVM and IVC rate (p < 0.05). The pattern C of the ND1 gene significantly increased the IVM rate compared to the patterns A, B and D. For the IVC, the highest and lowest means were related to the C and B patterns, respectively. The CYTB gene polymorphisms also had significant effects on IVC (p < 0.01), but the IVM did not affected (p = 0.07). Here, the pattern D had the highest and the pattern C had the lowest means for both IVM and IVC rates.

Breed-specific reference sequence optimized mapping accuracy of NGS analyses for pigs

Background

Reference sequences play a vital role in next-generation sequencing (NGS), impacting mapping quality during genome analyses. However, reference genomes usually do not represent the full range of genetic diversity of a species as a result of geographical divergence and independent demographic events of different populations. For the mitochondrial genome (mitogenome), which occurs in high copy numbers in cells and is strictly maternally inherited, an optimal reference sequence has the potential to make mitogenome alignment both more accurate and more efficient. In this study, we used three different types of reference sequences for mitogenome mapping, i.e., the commonly used reference sequence (CU-ref), the breed-specific reference sequence (BS-ref) and the sample-specific reference sequence (SS-ref), respectively, and compared the accuracy of mitogenome alignment and SNP calling among them, for the purpose of proposing the optimal reference sequence for mitochondrial DNA (mtDNA) analyses of specific populations

Results

Four pigs, representing three different breeds, were high-throughput sequenced, subsequently mapping reads to the reference sequences mentioned above, resulting in a largest mapping ratio and a deepest coverage without increased running time when aligning reads to a BS-ref. Next, single nucleotide polymorphism (SNP) calling was carried out by 18 detection strategies with the three tools SAMtools, VarScan and GATK with different parameters, using the bam results mapping to BS-ref. The results showed that all eighteen strategies achieved the same high specificity and sensitivity, which suggested a high accuracy of mitogenome alignment by the BS-ref because of a low requirement for SNP calling tools and parameter choices.

Conclusions

This study showed that different reference sequences representing different genetic relationships to sample reads influenced mitogenome alignment, with the breed-specific reference sequences being optimal for mitogenome analyses, which provides a refined processing perspective for NGS data.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08030-1.

The relationship between mitochondrial ND5 gene polymorphisms and in vitro embryo production in Sanjabi sheep

The aim of this study was to investigate mitochondrial ND5 gene polymorphisms and their relationship with in vitro maturation (IVM) and in vitro culture (IVC) of Sanjabi sheep. Blood and ovarian samples of adult ewes were obtained from a local slaughterhouse. For each ovarian sample, cumulus-oocyte complexes larger than 3 mm in diameter were aspirated from follicles, and their IVM and IVC rates were recorded. A 666-bp fragment of the ND5 gene was amplified using the polymerase chain reaction. The samples were genotyped using a modified single-stranded conformation polymorphism (SSCP) method, and an association study was conducted with IVM and IVC rates. Six different SSCP patterns, designated A, B, C, D, E and F with respective frequencies of 8, 47, 4, 4, 32 and 5%, respectively, were observed. According to the results of association analysis, there was no significant association between the ND1 gene polymorphisms and the IVM and IVC rates (P > 0.05).

Phylogenetic differentiation between Awassi and Hamdani sheep using the mitochondrial 12S rRNA

This study was conducted to assess the role of the mitochondrial 12S rRNA variations in the phylogenetic discrimination between two Iraqi breeds of sheep that differ in geographical distribution. A total of 122 animals (68 Awassi and 54 Hamdani) were included in the study. Direct sequencing of amplicons followed by the construction of a median-joining network and several trees were performed to identify the possible phylogenetic differences between both involved breeds. Genetic diversity, relative frequencies, and analysis of molecular variance (AMOVA) were performed to assess the genetic correlation between both populations. The median-joining network and minimized tree values showed that all observed haplotypes were separated into two groups according to their breed. Comprehensive phylogenetic data revealed only one Asian ancestor for all observed haplotypes. As indicated by AMOVA, the observed diversity was mostly due to between-population variation (1.24836%), while within-population variation (0.91221%) accounted for much less. The currently investigated rRNA amplicons exhibited different mitochondrial manifestations between Awassi and Hamdani breeds. Due to the ability of these 12S rRNA amplicons to mimic the geographical diversity for the currently investigated breeds, it is highly recommended to be used as potent mitochondrial genetic markers among broader ovine sequences.

Association analysis of mitochondrial DNA polymorphisms with oocyte number in pigs

In pigs, correlations between mitochondrial (mt) DNA polymorphisms and economic traits have been widely reported across and within swine breeds. In fecundity studies, the number of oocytes within ovaries was highly correlated with litter size. However, the effect of mitochondrial polymorphisms on porcine oocyte number remained unclear. In this study, 181 porcine ovaries were collected to analyse the relationship between oocyte number and mtDNA polymorphisms. There were considerable differences in oocyte numbers among different ovaries from commercial pig breeds, ranging from 2.7×105 to 1.3×106. Mitochondrial D-loop sequencing discovered 53 polymorphic sites. Association analysis revealed that 13 variations were associated with the number of oocytes (P<0.05). A C323T polymorphism showed the largest value between the C and T carriers, which differed at 105 oocytes (P<0.05). The 53 polymorphic sites generated 45 haplotypes, which clustered into two haplogroups, A and B. Haplogroup A had a higher number of oocytes than Haplogroup B (P<0.05), whereas Haplotype H6 in Haplogroup A had the highest number of oocytes (~7.5×105) of all haplotypes studied (P<0.05). The results of this study highlight a correlation between mtDNA polymorphisms and oocyte number, and suggest the potential application of mtDNA polymorphism analyses in pig selection and breeding practices.

Identifying Pig Mitochondrial TSS: Structure and Functional Features

The transcription start sites (TSSs) of porcine mitochondrial genome were firstly identified in this study, including heavy-strand promoter 1 and 2 (HSP1 and HSP2) harbored at nt 903 and nt 1369 in H strand, respectively, and light-strand promoter (LSP) located at nt 166 in L strand. HSP1 structure and expression features were investigated by analyzing mtDNA copy number, expression of 11 nucleoplasmic genes, mtDNA methylation levels, and gene expression levels of methyl-modifying enzymes, DNMT1 and TETs. The mtDNA copy number presented large differences among 15 organs/tissues, and the largest disparity, nearly 17 times, was found between pancreas (~1890 relative copy numbers) and spleen (~110 relative copy numbers, P < .01). The expression levels of HSP1 strand in these organs/tissues presented similar trends with mtDNA copy number (P < .05), and all of 11 nucleoplasmic genes (POLG, POLRMT, TERT, TFAM, TFB1M, TFB2M, NRF-1, PPARα, ESRRA, SP1 and TUFM) detected in this study displayed significantly higher expression values in pancreas than those in spleen (P < .05). Besides, bisulfite sequencing showed that all cytosine residues in the detected region (D-loop) existed methylation with different levels, and the methylation level in spleen was significantly higher than that in pancreas (P < .05). Unlike nuclear DNA, the tested region contained four types of methylation mode (CA, CC, CT, and CG). In addition, the expression of TET1 in pancreas was significantly higher than that in spleen (P < .05). Collectively, our findings indicated that mtDNA TSSs had correlation to mtDNA copy number, expression of nucleoplasmic gene, and mtDNA methylation level.

Analysis of the impact of cytoplasmic and mitochondrial inheritance on litter size and carcass in rabbits

<p>The effects of mitogenome variation on economically important traits have been reported in a number of domestic animal species. In this study, the first of its kind on rabbits, we have performed the estimation of the contribution of cytoplasmic and D-loop mitochondrial DNA (mtDNA) sequence effects on the litter size and carcass traits in three Pannon rabbit breeds (Pannon Ka, Pannon Large and Pannon White). The observed effects of both estimates, coming from cytoplasmic or D-loop mtDNA variation, were negligible. The most likely explanation for the results obtained is the lack of mitogenome polymorphism in all three populations, as suggested from the analysis performed on the D-loop mtDNA sequence, here assigned to the two most frequent rabbit haplotypes. The extent of potential benefits of the introduction, or alteration, of mitogenome variation in rabbit breeding remains an open question for future research.</p>

Polymorphism of mitochondrial tRNA genes associated with the number of pigs born alive

Background

Mutations in mitochondrial tRNA genes have been widely reported association with human reproductions. It is also important to explore the effect on the number of piglets born alive (NBA). Here, 1017 sows were used to investigate the association between polymorphisms in mitochondrial tRNA genes and NBA.

Results

In total, 16 mutations were found in mitochondrial tRNA genes, of which 13 mutations were significantly associated with NBA (P < 0.05). The reproductions of mutant carriers were significantly greater than that of wild carriers by 0.989 piglets born alive/sow farrowing. To test whether the mutations altered the structure of mitochondrial tRNAs, the secondary and tertiary structures were predicted. In result, C2255T changed the secondary structure of tRNA-Val by elongating the T stem and shrinking the T loop, and C2255T and G2259A in the tRNA-Val gene, C6217T and T6219C in the tRNA-Ala gene, and T15283C in the tRNA-Glu gene altered the tertiary structure of their tRNAs, respectively by changing the folding form of the T arm, and C16487T in the tRNA-Thr gene changed the tertiary structure of mitochondrial tRNA-Thr by influencing the folding form of the acceptor arm.

Conclusions

Results highlight the effect of mitochondrial tRNA genes on the number of piglets born alive, and suggest that polymorphic sites of the tRNA genes be genetic markers for selection of pig reproduction.

Electronic supplementary material

The online version of this article (10.1186/s40104-018-0299-0) contains supplementary material, which is available to authorized users.

Single Nucleotide Polymorphisms in the HIRA Gene Affect Litter Size in Small Tail Han Sheep

Simple Summary

Litter size is one of the most important reproductive traits in sheep. Two single nucleotide polymorphisms (SNPs), g.71874104G>A and g.71833755T>C, in the Histone Cell Cycle Regulator (HIRA) gene, were identified by whole-genome sequencing (WGS) and may be correlated with litter size in sheep. The two SNPs were genotyped and expression patterns of HIRA was determined in sheep breeds with different fecundity and in groups of Small Tail Han sheep producing large or small litters. Association analysis indicated that both SNPs were significantly correlated with litter size in Small Tail Han sheep. Furthermore, high levels of HIRA expression may have a negative effect on litter size in Small Tail Han sheep.

Abstract

Maintenance of appropriate levels of fecundity is critical for efficient sheep production. Opportunities to increase sheep litter size include identifying single gene mutations with major effects on ovulation rate and litter size. Whole-genome sequencing (WGS) data of 89 Chinese domestic sheep from nine different geographical locations and ten Australian sheep were analyzed to detect new polymorphisms affecting litter size. Comparative genomic analysis of sheep with contrasting litter size detected a novel set of candidate genes. Two SNPs, g.71874104G>A and g.71833755T>C, were genotyped in 760 Small Tail Han sheep and analyzed for association with litter size. The two SNPs were significantly associated with litter size, being in strong linkage disequilibrium in the region 71.80–71.87 Mb. This haplotype block contains one gene that may affect litter size, Histone Cell Cycle Regulator (HIRA). HIRA mRNA levels in sheep with different lambing ability were significantly higher in ovaries of Small Tail Han sheep (high fecundity) than in Sunite sheep (low fecundity). Moreover, the expression levels of HIRA in eight tissues of uniparous Small Tail Han sheep were significantly higher than in multiparous Small Tail Han sheep (p < 0.05). HIRA SNPs significantly affect litter size in sheep and are useful as genetic markers for litter size.