Annotation of novel transcripts putatively relevant for bovine fat metabolism

Genome-wide detection of copy number variations in polled yak using the Illumina BovineHD BeadChip

Background

Copy number variations (CNVs), which are genetic variations present throughout mammalian genomes, are a vital source of phenotypic variation that can lead to the development of unique traits. In this study we used the Illunima BovineHD BeadChip to conduct genome-wide detection of CNVs in 215 polled yaks.

Results

A total of 1066 CNV regions (CNVRs) were detected with a total length of 181.6 Mb, comprising ~ 7.2% of the bovine autosomal genome. The size of these CNVRs ranged from 5.53 kb to 1148.45 kb, with an average size of 170.31 kb. Eight out of nine randomly chosen CNVRs were successfully validated by qPCR. A functional enrichment analysis of the CNVR-associated genes indicated their relationship to a number of molecular adaptations that enable yaks to thrive at high altitudes. One third of the detected CNVRs were mapped to QTLs associated with six classes of economically important traits, indicating that these CNVRs may play an important role in variations of these traits.

Conclusions

Our genome-wide yak CNV map may thus provide valuable insights into both the molecular mechanisms of high altitude adaptation and the potential genomic basis of economically important traits in yak.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5759-1) contains supplementary material, which is available to authorized users.

Genomic structure of nucleotide diversity among Lyon rat models of metabolic syndrome

Background

The metabolic syndrome (MetS), a complex disorder involving hypertension, obesity, dyslipidemia and insulin resistance, is a major risk factor for heart disease, stroke, and diabetes. The Lyon Hypertensive (LH), Lyon Normotensive (LN) and Lyon Low-pressure (LL) rats are inbred strains simultaneously derived from a common outbred Sprague Dawley colony by selection for high, normal, and low blood pressure, respectively. Further studies found that LH is a MetS susceptible strain, while LN is resistant and LL has an intermediate phenotype. Whole genome sequencing determined that, while the strains are phenotypically divergent, they are nearly 98% similar at the nucleotide level. Using the sequence of the three strains, we applied an approach that harnesses the distribution of Observed Strain Differences (OSD), or nucleotide diversity, to distinguish genomic regions of identity-by-descent (IBD) from those with divergent ancestry between the three strains. This information was then used to fine-map QTL identified in a cross between LH and LN rats in order to identify candidate genes causing the phenotypes.

Results

We identified haplotypes that, in total, contain at least 95% of the identifiable polymorphisms between the Lyon strains that are likely of differing ancestral origin. By intersecting the identified haplotype blocks with Quantitative Trait Loci (QTL) previously identified in a cross between LH and LN strains, the candidate QTL regions have been narrowed by 78%. Because the genome sequence has been determined, we were further able to identify putative functional variants in genes that are candidates for causing the QTL.

Conclusions

Whole genome sequence analysis between the LH, LN, and LL strains identified the haplotype structure of these three strains and identified candidate genes with sequence variants predicted to affect gene function. This approach, merged with additional integrative genetics approaches, will likely lead to novel mechanisms underlying complex disease and provide new drug targets and therapies.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-197) contains supplementary material, which is available to authorized users.

Gene expression phenotypes for lipid metabolism and intramuscular fat in skeletal muscle of cattle

Gene expression phenotypes were evaluated for intramuscular fat (IMF) in bovine skeletal muscle as an alternative to traditional estimates of IMF%. Gene expression data from a time course of LM development in high- and low-marbling Bos taurus cattle crosses were compared to identify genes involved in intramuscular adipocyte lipid metabolism with developmentally similar gene expression profiles. Three sets of genes were identified: triacylglyceride (TAG) synthesis and storage, fatty acid (FA) synthesis, and PPARγ-related genes. In an independent analysis in the LM of 48 Bos indicus cattle, TAG and FA gene sets were enriched in the top 100 genes of which expression was most correlated with IMF% (P = 1.2 × 10(-24) and 3.5 × 10(-9), respectively). In general, genes encoding enzymes involved in the synthesis of FA and TAG in the intramuscular adipocytes were present in the top 100 genes. In B. indicus, effects of a steroid hormone growth promotant (HGP), 2 experimental sites [New South Wales (NSW) and Western Australia (WA)], and 3 tenderness genotypes on the expression levels of genes in the TAG gene set and the correlation of gene expression with IMF% were investigated. Although correlation between expression of 12 individual TAG genes and IMF% was observed in HGP-treated animals in both experimental sites (mean r = 0.43), correlation was not observed for untreated animals at the NSW site (mean r = -0.07, P < 3 × 10(-6)). However, TAG genes showed an average 1.6-fold (P < 0.0004) reduction in expression in the LM of HGP-treated cattle relative to untreated cattle, an effect consistent across both experimental sites. Cattle possessing the favored tenderness calpain 1 and 3 and calpastatin alleles exhibited a greater (P = 0.008) reduction in expression in NSW (1.8-fold reduction, P = 0.0002) compared with WA (1.2-fold reduction, P = 0.03). Tenderness genotype had no impact (P > 0.05) on the correlation of TAG genes with IMF%. In general, the interactions among genotype, treatment and location, and TAG gene set gene expression were consistent with the interactions among the same factors and IMF% detected using 255 animals, of which the 48 in this study were a subset. Thus, the TAG gene set constitutes a gene expression phenotype able to predict effects of different genotypes and treatments on IMF% using much smaller groups than current approaches, even in animals with very low IMF%.

Molecular characterization, polymorphism of bovine ZBTB38 gene and association with body measurement traits in native Chinese cattle breeds

Zinc finger and BTB domain containing 38 (ZBTB38), binding to and repressing methylated DNA, is an important candidate gene for selection of body measurement traits through marker-assisted selection (MAS). The expression of ZBTB38 is regulated in human and animal height as well as other stature indexes. Genomic structural analysis shows that bovine ZBTB38 shares much similarity with human ZBTB38. We discovered and evaluated the potential association of the single nucleotide polymorphism (SNP) of the bovine ZBTB38 gene with body measurement traits in 722 individuals. The latest findings demonstrate that the A841G SNP in exon 1 is significantly associated with Body Length (BL), Hip Height (HH) and Heart Girth (HG). Furthermore, the analysis of A841G SNP marker shows that there are significant effects on the BL (P = 0.0389) in 722 individuals, significant effects on the HH (P = 0.0173) and HG (P = 0.0147) in Qinchuan improvement steers (QI) population, as well as significant effects on the WH (P = 0.0094) in Xuelong (XL) population. These results clearly suggest that the ZBTB38 gene is among of target genes for body measurement traits in bovine reproduction and breeding, and thus provide data for establishment of an animal model using cattle to study big animal body type.