Genome-wide association study identifies functional genomic variants associated with young stock survival in Nordic Red Dairy Cattle.

Identifying quantitative trait loci (QTL) associated with calf survival is essential for both reducing economic loss in cattle industry and understanding the genetic basis of the trait. To identify mutations and genes underlying young stock survival (YSS), we performed GWAS using de-regressed estimated breeding values of a YSS index and its component traits defined by sex and age in 3,077 Nordic Red Dairy Cattle (RDC) bulls and 2 stillbirth traits (first lactation and later lactations) in 5,141 RDC bulls. Two associated QTL regions on Bos taurus autosome (BTA) 4 and 6 were identified for the YSS index. The results of 4 YSS component traits indicate that same QTL regions were associated with bull and heifer calf mortality, but the effects were different over the growing period and suggested an additional QTL on BTA23. The GWAS on stillbirth identified 3 additional QTL regions on BTA5, 14, and 24 compared with YSS and its component traits. The conditional test of BTA6 showed at least 2 closely located QTL segregating for YSS component traits and stillbirth. We found 2 independent QTL for stillbirth on BTA23. The post-GWAS revealed LCORL, PPM1K, SSP1, MED28, and LAP3 are putative causal genes on BTA6, and a frame shift variant within LCORL, BTA6:37401770 (rs384548488) could be the putative causal variant. On BTA4, the GRB10 gene is the putative causal gene and BTA4:5296018 is the putative causal variant. In addition, NDUFA9 and FGF23 on BTA5, LYN on BTA14, and KCNK5 on BTA23 are putative causal genes for QTL for stillbirth. The gene analysis also proposed several candidate genes. Our findings shed new light on the candidate genes affecting calf survival, and the knowledge could be utilized to reduce calf mortality and thereby enhance welfare of dairy cattle.

Girl Power in Glaucoma: The Role of Estrogen in Primary Open Angle Glaucoma.

Estrogen is essential in maintaining various physiological features in women, and a decline in estrogen levels are known to give rise to numerous unfortunate symptoms associated with menopause. To alleviate these symptoms hormone replacement therapy with estrogen is often used, and has been shown to be fruitful in improving quality of life in women suffering from postmenopausal discomforts. An often forgotten condition associated with menopause is the optic nerve disorder, glaucoma. Thus, estrogen may also have an impact in maintaining the retinal ganglion cells (RGCs), which make up the optic nerve, thereby preventing glaucomatous neurodegeneration. This review aims to provide an overview of possible associations of estrogen and the glaucoma subtype, primary open-angle glaucoma (POAG), by evaluating the current literature through a PubMed-based literature search. Multiple in vitro and in vivo studies of RGC protection, as well as clinical and epidemiological data concerning the well-defined retinal neurodegenerative disorder POAG have been reviewed. Over all, deficiencies in retinal estrogen may potentially instigate RGC loss, visual disability, and eventual blindness. Estrogen replacement therapy may therefore be a beneficial future treatment. However, more studies are needed to confirm the relevance of estrogen in glaucoma prevention.

Structures of liquid and aqueous water isotopologues at ambient temperature from ab initio path integral simulations.

The heavy hydrogen isotopes D and T are found in trace amounts in water, and when their concentration increases they can play an intricate role in modulating the physical properties of the liquid. We present an analysis of the microscopic structures of ambient light water (H2O(l)), heavy water (D2O(l)), T2O(l), HDO(aq) and HTO(aq) studied by ab initio path integral molecular dynamics (PIMD). Unlike previous ab initio PIMD investigations of H2O(l) and D2O(l) [Chen et al., Phys. Rev. Lett., 2003, 91, 215503] [Machida et al., J. Chem. Phys., 2017, 148, 102324] we find that D2O(l) is more structured than H2O(l), as is predicted by the experiment. The agreement between the experiment and our simulation for H2O(l) and D2O(l) allows us to accurately predict the intra- and intermolecular structures of T2O(l) HDO(aq) and HTO(aq). T2O(l) is found to have a similar intermolecular structure to that of D2O(l), while the intramolecular structure is more compact, giving rise to a smaller dipole moment than those of H2O(l) and D2O(l). For the mixed isotope species, HDO(aq) and HTO(aq), we find smaller dipole moments and fewer hydrogen bonds when compared with the pure species H2O and D2O. We can attribute this effect to the relative compactness of the mixed isotope species, which results in a lower dipole moment than that of the pure species.

Ab initio study of nuclear quantum effects on sub- and supercritical water.

The structures of water in the ambient, subcritical, and supercritical conditions at various densities were studied systematically by ab initio path integral molecular dynamics simulations. It was found that the nuclear quantum effects (NQEs) have a significant impact on the structure of hydrogen bonds in close contact, not only in the ambient condition but also in the sub- and supercritical conditions. The NQEs on the structure beyond the hydrogen bond contact are important in ambient water, but not much for water in the sub- and supercritical conditions. The NQEs are furthermore important for determining the number of hydrogen bonds in the ambient conditions, and this role is, however, diminished in the sub- and supercritical conditions. The NQEs do, nevertheless, show their importance in determining the intramolecular structure of water and the close contact structures of the hydrogen bonds, even at sub- and supercritical conditions. Using the RPBE-D3 functional, the computed radial distribution functions for ambient water are in excellent agreement with experimental data, upgrading our previous results using the BLYP-D2 functional [Machida et al., J. Chem. Phys. 148, 102324 (2018)]. The computed radial distribution functions for water in the sub- and supercritical conditions were carefully compared with experiment. In particular, we found that the first peak in hydrogen pair distribution functions matches only when the NQEs are taken into account.

Nuclear quantum effects on autoionization of water isotopologs studied by ab initio path integral molecular dynamics.

In this study, we investigate the nuclear quantum effects (NQEs) on the acidity constant (pKA) of liquid water isotopologs under the ambient condition by path integral molecular dynamics (PIMD) simulations. We compared simulations using a fully explicit solvent model with a classical polarizable force field, density functional tight binding, and ab initio density functional theory, which correspond to empirical, semiempirical, and ab initio PIMD simulations, respectively. The centroid variable with respect to the proton coordination number of a water molecule was restrained to compute the gradient of the free energy, which measures the reversible work of the proton abstraction for the quantum mechanical system. The free energy curve obtained by thermodynamic integration was used to compute the pKA value based on probabilistic determination. This technique not only reproduces the pKA value of liquid D2O experimentally measured (14.86) but also allows for a theoretical prediction of the pKA values of liquid T2O and aqueous HDO and HTO, which are unknown due to their scarcity. It is also shown that the NQEs on the free energy curve can result in a downshift of 4.5 ± 0.9 pKA units in the case of liquid water, which indicates that the NQEs plays an indispensable role in the absolute determination of pKA. The results of this study can help inform further extensions into the calculation of the acidity constants of isotope substituted species with high accuracy.

High-Amylose Maize, Potato, and Butyrylated Starch Modulate Large Intestinal Fermentation, Microbial Composition, and Oncogenic miRNA Expression in Rats Fed A High-Protein Meat Diet.

High red meat intake is associated with the risk of colorectal cancer (CRC), whereas dietary fibers, such as resistant starch (RS) seemed to protect against CRC. The aim of this study was to determine whether high-amylose potato starch (HAPS), high-amylose maize starch (HAMS), and butyrylated high-amylose maize starch (HAMSB)-produced by an organocatalytic route-could oppose the negative effects of a high-protein meat diet (HPM), in terms of fermentation pattern, cecal microbial composition, and colonic biomarkers of CRC. Rats were fed a HPM diet or an HPM diet where 10% of the maize starch was substituted with either HAPS, HAMS, or HAMSB, for 4 weeks. Feces, cecum digesta, and colonic tissue were obtained for biochemical, microbial, gene expression (oncogenic microRNA), and immuno-histochemical (O6-methyl-2-deoxyguanosine (O6MeG) adduct) analysis. The HAMS and HAMSB diets shifted the fecal fermentation pattern from protein towards carbohydrate metabolism. The HAMSB diet also substantially increased fecal butyrate concentration and the pool, compared with the other diets. All three RS treatments altered the cecal microbial composition in a diet specific manner. HAPS and HAMSB showed CRC preventive effects, based on the reduced colonic oncogenic miR17-92 cluster miRNA expression, but there was no significant diet-induced differences in the colonic O6MeG adduct levels. Overall, HAMSB consumption showed the most potential for limiting the negative effects of a high-meat diet.

Employing general fit-bases for construction of potential energy surfaces with an adaptive density-guided approach.

We present an approach to treat sets of general fit-basis functions in a single uniform framework, where the functional form is supplied on input, i.e., the use of different functions does not require new code to be written. The fit-basis functions can be used to carry out linear fits to the grid of single points, which are generated with an adaptive density-guided approach (ADGA). A non-linear conjugate gradient method is used to optimize non-linear parameters if such are present in the fit-basis functions. This means that a set of fit-basis functions with the same inherent shape as the potential cuts can be requested and no other choices with regards to the fit-basis functions need to be taken. The general fit-basis framework is explored in relation to anharmonic potentials for model systems, diatomic molecules, water, and imidazole. The behaviour and performance of Morse and double-well fit-basis functions are compared to that of polynomial fit-basis functions for unsymmetrical single-minimum and symmetrical double-well potentials. Furthermore, calculations for water and imidazole were carried out using both normal coordinates and hybrid optimized and localized coordinates (HOLCs). Our results suggest that choosing a suitable set of fit-basis functions can improve the stability of the fitting routine and the overall efficiency of potential construction by lowering the number of single point calculations required for the ADGA. It is possible to reduce the number of terms in the potential by choosing the Morse and double-well fit-basis functions. These effects are substantial for normal coordinates but become even more pronounced if HOLCs are used.

Infrared Spectra of Protonated Water Clusters, H+(H2O)4, in Eigen and Zundel Forms Studied by Vibrational Quasi-Degenerate Perturbation Theory.

The infrared spectrum of H+(H2O)4 recently observed in a wide spectral range has shown a series of bands in a range of 1700-2500 cm-1, which can not be understood by the standard harmonic normal mode analysis. Here, we theoretically investigate the origin of these bands with a focus on (1) the possibility of coexistence of multiple isomers in the Eigen [H3O+(H2O)3] and Zundel [H5O2+(H2O)2] forms and (2) the effect of anharmonic coupling that gives rise to nonzero intensities for overtones and combination bands. Anharmonic vibrational calculations are carried out for the Eigen and Zundel clusters by the second-order vibrational quasi-degenerate perturbation theory (VQDPT2) based on optimized coordinates. The anharmonic potential energy surface and the dipole moment surfaces are generated by a multiresolution approach combining one-dimensional (1D) grid potential functions derived from CCSD(T)-F12, 2D and 3D grid potential functions derived from B3LYP for important coupling terms, and a quartic force field derived from B3LYP for less important terms. The spectrum calculated for the Eigen cluster is in excellent agreement with the experiment, assigning the bands in the range of 1700-2500 cm-1 to overtones and combination bands of a H3O+ moiety in line with recent reports [ J. Phys. Chem. A 2015 , 119 , 9425 ; Science 2016 , 354 , 1131 ]. On the other hand, characteristic OH stretching bands of the Zundel cluster is found to be absent in the experimental spectrum. We therefore conclude that the experimental spectrum originates solely from the Eigen cluster. Nonetheless, the present calculation for the Eigen cluster poorly reproduces a band observed at 1765 cm-1. A possible nature of this band is discussed.

Partitioning of genomic variance reveals biological pathways associated with udder health and milk production traits in dairy cattle.

BACKGROUND: We have used a linear mixed model (LMM) approach to examine the joint contribution of genetic markers associated with a biological pathway. However, with these markers being scattered throughout the genome, we are faced with the challenge of modelling the contribution from several, sometimes even all, chromosomes at once. Due to linkage disequilibrium (LD), all markers may be assumed to account for some genomic variance; but the question is whether random sets of markers account for the same genomic variance as markers associated with a biological pathway? RESULTS: We applied the LMM approach to identify biological pathways associated with udder health and milk production traits in dairy cattle. A random gene sampling procedure was applied to assess the biological pathways in a dataset that has an inherently complex genetic correlation pattern due to the population structure of dairy cattle, and to linkage disequilibrium within the bovine genome and within the genes associated to the biological pathway. CONCLUSIONS: Several biological pathways that were significantly associated with health and production traits were identified in dairy cattle; i.e. the markers linked to these pathways explained more of the genomic variance and provided a better model fit than 95 % of the randomly sampled gene groups. Our results show that immune related pathways are associated with production traits, and that pathways that include a causal marker for production traits are identified with our procedure. We are confident that the LMM approach provides a general framework to exploit and integrate prior biological information and could potentially lead to improved understanding of the genetic architecture of complex traits and diseases.

A random set scoring model for prioritization of disease candidate genes using protein complexes and data-mining of GeneRIF, OMIM and PubMed records.

BACKGROUND: Prioritizing genetic variants is a challenge because disease susceptibility loci are often located in genes of unknown function or the relationship with the corresponding phenotype is unclear. A global data-mining exercise on the biomedical literature can establish the phenotypic profile of genes with respect to their connection to disease phenotypes. The importance of protein-protein interaction networks in the genetic heterogeneity of common diseases or complex traits is becoming increasingly recognized. Thus, the development of a network-based approach combined with phenotypic profiling would be useful for disease gene prioritization. RESULTS: We developed a random-set scoring model and implemented it to quantify phenotype relevance in a network-based disease gene-prioritization approach. We validated our approach based on different gene phenotypic profiles, which were generated from PubMed abstracts, OMIM, and GeneRIF records. We also investigated the validity of several vocabulary filters and different likelihood thresholds for predicted protein-protein interactions in terms of their effect on the network-based gene-prioritization approach, which relies on text-mining of the phenotype data. Our method demonstrated good precision and sensitivity compared with those of two alternative complex-based prioritization approaches. We then conducted a global ranking of all human genes according to their relevance to a range of human diseases. The resulting accurate ranking of known causal genes supported the reliability of our approach. Moreover, these data suggest many promising novel candidate genes for human disorders that have a complex mode of inheritance. CONCLUSION: We have implemented and validated a network-based approach to prioritize genes for human diseases based on their phenotypic profile. We have devised a powerful and transparent tool to identify and rank candidate genes. Our global gene prioritization provides a unique resource for the biological interpretation of data from genome-wide association studies, and will help in the understanding of how the associated genetic variants influence disease or quantitative phenotypes.

Structured RNAs and synteny regions in the pig genome.

BACKGROUND: Annotating mammalian genomes for noncoding RNAs (ncRNAs) is nontrivial since far from all ncRNAs are known and the computational models are resource demanding. Currently, the human genome holds the best mammalian ncRNA annotation, a result of numerous efforts by several groups. However, a more direct strategy is desired for the increasing number of sequenced mammalian genomes of which some, such as the pig, are relevant as disease models and production animals. RESULTS: We present a comprehensive annotation of structured RNAs in the pig genome. Combining sequence and structure similarity search as well as class specific methods, we obtained a conservative set with a total of 3,391 structured RNA loci of which 1,011 and 2,314, respectively, hold strong sequence and structure similarity to structured RNAs in existing databases. The RNA loci cover 139 cis-regulatory element loci, 58 lncRNA loci, 11 conflicts of annotation, and 3,183 ncRNA genes. The ncRNA genes comprise 359 miRNAs, 8 ribozymes, 185 rRNAs, 638 snoRNAs, 1,030 snRNAs, 810 tRNAs and 153 ncRNA genes not belonging to the here fore mentioned classes. When running the pipeline on a local shuffled version of the genome, we obtained no matches at the highest confidence level. Additional analysis of RNA-seq data from a pooled library from 10 different pig tissues added another 165 miRNA loci, yielding an overall annotation of 3,556 structured RNA loci. This annotation represents our best effort at making an automated annotation. To further enhance the reliability, 571 of the 3,556 structured RNAs were manually curated by methods depending on the RNA class while 1,581 were declared as pseudogenes. We further created a multiple alignment of pig against 20 representative vertebrates, from which RNAz predicted 83,859 de novo RNA loci with conserved RNA structures. 528 of the RNAz predictions overlapped with the homology based annotation or novel miRNAs. We further present a substantial synteny analysis which includes 1,004 lineage specific de novo RNA loci and 4 ncRNA loci in the known annotation specific for Laurasiatheria (pig, cow, dolphin, horse, cat, dog, hedgehog). CONCLUSIONS: We have obtained one of the most comprehensive annotations for structured ncRNAs of a mammalian genome, which is likely to play central roles in both health modelling and production. The core annotation is available in Ensembl 70 and the complete annotation is available at http://rth.dk/resources/rnannotator/susscr102/version1.02.

Pairwise comparisons of ten porcine tissues identify differential transcriptional regulation at the gene, isoform, promoter and transcription start site level.

The transcriptome is the absolute set of transcripts in a tissue or cell at the time of sampling. In this study RNA-Seq is employed to enable the differential analysis of the transcriptome profile for ten porcine tissues in order to evaluate differences between the tissues at the gene and isoform expression level, together with an analysis of variation in transcription start sites, promoter usage, and splicing. Totally, 223 million RNA fragments were sequenced leading to the identification of 59,930 transcribed gene locations and 290,936 transcript variants using Cufflinks with similarity to approximately 13,899 annotated human genes. Pairwise analysis of tissues for differential expression at the gene level showed that the smallest differences were between tissues originating from the porcine brain. Interestingly, the relative level of differential expression at the isoform level did generally not vary between tissue contrasts. Furthermore, analysis of differential promoter usage between tissues, revealed a proportionally higher variation between cerebellum (CBE) versus frontal cortex and cerebellum versus hypothalamus (HYP) than in the remaining comparisons. In addition, the comparison of differential transcription start sites showed that the number of these sites is generally increased in comparisons including hypothalamus in contrast to other pairwise assessments. A comprehensive analysis of one of the tissue contrasts, i.e. cerebellum versus heart for differential variation at the gene, isoform, and transcription start site (TSS), and promoter level showed that several of the genes differed at all four levels. Interestingly, these genes were mainly annotated to the "electron transport chain" and neuronal differentiation, emphasizing that "tissue important" genes are regulated at several levels. Furthermore, our analysis shows that the "across tissue approach" has a promising potential when screening for possible explanations for variations, such as those observed at the gene expression levels.

Tensor decomposition and vibrational coupled cluster theory.

The use of tensor decomposition in the calculation of anharmonic vibrational wave functions is discussed. The correlation amplitudes of vibrational coupled cluster (VCC) and vibrational configuration interaction (VCI) theories are considered as tensors and decomposed. A pilot code is implemented allowing a numerical study of the performance of the canonical decomposition/parallel factors (CP) for three and higher mode couplings in computations on water, formaldehyde, and 1,2,5-thiadiazole. The results show that there is a significant perspective in applying tensor decomposition in the context of anharmonic vibrational wave functions, with the CP tensor decomposition providing compression of data and a computational convenient representation. The calculations also illustrate how the multiplicative separability of the VCC ansatz with respect to noninteracting degrees of freedom goes well together with a tensor decomposition approach. Tensor decomposition opens for adjusting the computational effort spent on a particular mode-coupling according to the significance of that particular coupling, which is guaranteed to decrease to zero in the case of VCC in the limit of noninteracting subsystems.

MicroRNA sequence variation can impact interactions with target mRNA in cattle.

BACKGROUND: Small non-coding microRNAs (miRNAs) are important modulators at post-transcriptional levels. Single-nucleotide polymorphisms (SNPs) located in miRNA genes can alter the secondary structure of pre-miRNA to either impair or promote the miRNA maturation processes. Furthermore, SNPs located in the miRNA seed regions can stabilize or disturb miRNA-target interactions, thereby, quantitatively influence the expression of target genes. Therefore, the main objective of this study was to detect SNPs in bovine miRNAs using the whole-genome re-sequencing datasets of 1632 cattle of five breeds from the 1000 bull genomes project. RESULTS: In total, our study identified 1109, 334, and 130 SNPs in the miRNA precursor, mature, and seed regions, respectively. The heterozygosity values were generally less than 0.3, and the minor allele frequencies (MAFs) were mainly less than 0.1. Most SNPs were in Hardy-Weinberg equilibrium (HWE) (HWE-P > 0.05). Furthermore, we found that the majority of SNPs (MAF > 0.1 and HWE-P > 0.05) in the miRNA seed regions altered the repertoire of target genes, which in turn were enriched in different KEGG pathways or GO terms. Thus target prediction for bta-miR-2888 revealed loss of 309 target genes and gain of 691 target genes. The 691 gained target genes were significantly enriched in 60 KEGG pathways and 21 GO terms. CONCLUSION: In summary, our study identified candidate SNPs in miRNA precursor, mature, and seed regions that are likely to affect target RNA interactions, thereby potentially influencing cattle phenotypic traits.

Gene prioritization for livestock diseases by data integration.

Identifying causal genes that underlie complex traits such as susceptibility to disease is a primary aim of genetic and biomedical studies. Genetic mapping of quantitative trait loci (QTL) and gene expression profiling based on high-throughput technologies are common first approaches toward identifying associations between genes and traits; however, it is often difficult to assess whether the biological function of a putative candidate gene is consistent with a particular phenotype. Here, we have implemented a network-based disease gene prioritization approach for ranking genes associated with quantitative traits and diseases in livestock species. The approach uses ortholog mapping and integrates information on disease or trait phenotypes, gene-associated phenotypes, and protein-protein interactions. It was used for ranking all known genes present in the cattle genome for their potential roles in bovine mastitis. Gene-associated phenome profile and transcriptome profile in response to Escherichia coli infection in the mammary gland were integrated to make a global inference of bovine genes involved in mastitis. The top ranked genes were highly enriched for pathways and biological processes underlying inflammation and immune responses, which supports the validity of our approach for identifying genes that are relevant to animal health and disease. These gene-associated phenotypes were used for a local prioritization of candidate genes located in a QTL affecting the susceptibility to mastitis. Our study provides a general framework for prioritizing genes associated with various complex traits in different species. To our knowledge this is the first time that gene expression, ortholog mapping, protein interactions, and biomedical text data have been integrated systematically for ranking candidate genes in any livestock species.

Gene expression profiling of porcine skeletal muscle in the early recovery phase following acute physical activity.

Acute physical activity elicits changes in gene expression in skeletal muscles to promote metabolic changes and to repair exercise-induced muscle injuries. In the present time-course study, pigs were submitted to an acute bout of treadmill running until near exhaustion to determine the impact of unaccustomed exercise on global transcriptional profiles in porcine skeletal muscles. Using a combined microarray and candidate gene approach, we identified a suite of genes that are differentially expressed in muscles during postexercise recovery. Several members of the heat shock protein family and proteins associated with proteolytic events, such as the muscle-specific E3 ubiquitin ligase atrogin-1, were significantly upregulated, suggesting that protein breakdown, prevention of protein aggregation and stabilization of unfolded proteins are important processes for restoration of cellular homeostasis. We also detected an upregulation of genes that are associated with muscle cell proliferation and differentiation, including MUSTN1, ASB5 and CSRP3, possibly reflecting activation, differentiation and fusion of satellite cells to facilitate repair of muscle damage. In addition, exercise increased expression of the orphan nuclear hormone receptor NR4A3, which regulates metabolic functions associated with lipid, carbohydrate and energy homeostasis. Finally, we observed an unanticipated induction of the long non-coding RNA transcript NEAT1, which has been implicated in RNA processing and nuclear retention of adenosine-to-inosine edited mRNAs in the ribonucleoprotein bodies called paraspeckles. These findings expand the complexity of pathways affected by acute contractile activity of skeletal muscle, contributing to a better understanding of the molecular processes that occur in muscle tissue in the recovery phase.

Vibrational absorption spectra from vibrational coupled cluster damped linear response functions calculated using an asymmetric Lanczos algorithm.

We report the theory and implementation of vibrational coupled cluster (VCC) damped response functions. From the imaginary part of the damped VCC response function the absorption as function of frequency can be obtained, requiring formally the solution of the now complex VCC response equations. The absorption spectrum can in this formulation be seen as a matrix function of the characteristic VCC Jacobian response matrix. The asymmetric matrix version of the Lanczos method is used to generate a tridiagonal representation of the VCC response Jacobian. Solving the complex response equations in the relevant Lanczos space provides a method for calculating the VCC damped response functions and thereby subsequently the absorption spectra. The convergence behaviour of the algorithm is discussed theoretically and tested for different levels of completeness of the VCC expansion. Comparison is made with results from the recently reported [P. Seidler, M. B. Hansen, W. Györffy, D. Toffoli, and O. Christiansen, J. Chem. Phys. 132, 164105 (2010)] vibrational configuration interaction damped response function calculated using a symmetric Lanczos algorithm. Calculations of IR spectra of oxazole, cyclopropene, and uracil illustrate the usefulness of the new VCC based method.

Global assessment of genomic variation in cattle by genome resequencing and high-throughput genotyping.

BACKGROUND: Integration of genomic variation with phenotypic information is an effective approach for uncovering genotype-phenotype associations. This requires an accurate identification of the different types of variation in individual genomes. RESULTS: We report the integration of the whole genome sequence of a single Holstein Friesian bull with data from single nucleotide polymorphism (SNP) and comparative genomic hybridization (CGH) array technologies to determine a comprehensive spectrum of genomic variation. The performance of resequencing SNP detection was assessed by combining SNPs that were identified to be either in identity by descent (IBD) or in copy number variation (CNV) with results from SNP array genotyping. Coding insertions and deletions (indels) were found to be enriched for size in multiples of 3 and were located near the N- and C-termini of proteins. For larger indels, a combination of split-read and read-pair approaches proved to be complementary in finding different signatures. CNVs were identified on the basis of the depth of sequenced reads, and by using SNP and CGH arrays. CONCLUSIONS: Our results provide high resolution mapping of diverse classes of genomic variation in an individual bovine genome and demonstrate that structural variation surpasses sequence variation as the main component of genomic variability. Better accuracy of SNP detection was achieved with little loss of sensitivity when algorithms that implemented mapping quality were used. IBD regions were found to be instrumental for calculating resequencing SNP accuracy, while SNP detection within CNVs tended to be less reliable. CNV discovery was affected dramatically by platform resolution and coverage biases. The combined data for this study showed that at a moderate level of sequencing coverage, an ensemble of platforms and tools can be applied together to maximize the accurate detection of sequence and structural variants.

MicroRNA bta-miR-365-3p inhibits proliferation but promotes differentiation of primary bovine myoblasts by targeting the activin A receptor type I.

BACKGROUND: MicroRNAs act as post-transcriptional regulators that repress translation or degrade mRNA transcripts. Each microRNA has many mRNA targets and each mRNA may be targeted by several microRNAs. Skeletal muscles express a plethora of microRNA genes that regulate muscle development and function by controlling the expression of protein-coding target genes. To expand our understanding of the role of microRNA, specifically bta-miR-365-3p, in muscle biology, we investigated its functions in regulating primary bovine myoblast proliferation and differentiation. RESULTS: Firstly, we found that bta-miR-365-3p was predominantly expressed in skeletal muscle and heart tissue in Chinese Qinchuan beef cattle. Quantitative PCR and western blotting results showed that overexpression of bta-miR-365-3p significantly reduced the expression levels of cyclin D1 (CCND1), cyclin dependent kinase 2 (CDK2) and proliferating cell nuclear antigen (PCNA) but stimulated the expression levels of muscle differentiation markers, i.e., MYOD1, MYOG at both mRNA and protein level. Moreover, downregulation of bta-miR-365-3p increased the expression of CCND1, CDK2 and PCNA but decreased the expression of MYOD1 and MYOG at both mRNA and protein levels. Furthermore, flow cytometry, EdU proliferation assays and immunostaining results showed that increased levels of bta-miR-365-3p suppressed cell proliferation but promoted myotube formation, whereas decreased levels of bta-miR-365-3p resulted in the opposite consequences. Finally, we identified that activin A receptor type I (ACVR1) could be a direct target of bta-miR-365-3p. It was demonstrated that bta-miR-365-3p can bind to the 3'UTR of ACVR1 gene to regulate its expression based on dual luciferase gene reporter assays. Consistently, knock-down of ACVR1 was associated with decreased expressions of CDK2, CCND1 and PCNA but increased expression of MYOG and MYOD1 both at mRNA and protein level. CONCLUSION: Collectively, these data suggested that bta-miR-365-3p represses proliferation but promotes differentiation of bovine myoblasts through several biological mechanisms involving downregulation of ACVR1.