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.

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.

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.

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.

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.

Overview of Metabolomic Analysis and the Integration with Multi-Omics for Economic Traits in Cattle.

Metabolomics has been applied to measure the dynamic metabolic responses, to understand the systematic biological networks, to reveal the potential genetic architecture, etc., for human diseases and livestock traits. For example, the current published results include the detected relevant candidate metabolites, identified metabolic pathways, potential systematic networks, etc., for different cattle traits that can be applied for further metabolomic and integrated omics studies. Therefore, summarizing the applications of metabolomics for economic traits is required in cattle. We here provide a comprehensive review about metabolomic analysis and its integration with other omics in five aspects: (1) characterization of the metabolomic profile of cattle; (2) metabolomic applications in cattle; (3) integrated metabolomic analysis with other omics; (4) methods and tools in metabolomic analysis; and (5) further potentialities. The review aims to investigate the existing metabolomic studies by highlighting the results in cattle, integrated with other omics studies, to understand the metabolic mechanisms underlying the economic traits and to provide useful information for further research and practical breeding programs in cattle.

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.

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.

Integrated Analysis of mRNA and MicroRNA Co-expressed Network for the Differentiation of Bovine Skeletal Muscle Cells After Polyphenol Resveratrol Treatment.

Resveratrol (RSV) has been confirmed to benefit human health. Resveratrol supplemented in the feeds of animals improved pork, chicken, and duck meat qualities. In this study, we identified differentially expressed (DE) messenger RNAs (mRNAs) (n = 3,856) and microRNAs (miRNAs) (n = 93) for the weighted gene co-expression network analysis (WGCNA) to investigate the co-expressed DE mRNAs and DE miRNAs in the primary bovine myoblasts after RSV treatment. The mRNA results indicated that RSV treatments had high correlations with turquoise module (0.91, P-value = 0.01) and blue module (0.93, P-value < 0.01), while only the turquoise module (0.96, P-value < 0.01) was highly correlated with the treatment status using miRNA data. After biological enrichment analysis, the 2,579 DE genes in the turquoise module were significantly enriched in the Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. The top two GO terms were actin filament-based process (GO:0030029) and actin cytoskeleton organization (GO:0030036). The top two KEGG pathways were regulation of actin cytoskeleton (bta04810) and tight junction (bta04530). Then, we constructed the DE mRNA co-expression and DE miRNA co-expression networks in the turquoise module and the mRNA-miRNA targeting networks based on their co-expressions in the key module. In summary, the RSV-induced miRNAs participated in the co-expression networks that could affect mRNA expressions to regulate the primary myoblast differentiation. Our study provided a better understanding of the roles of RSV in inducing miRNA and of the characteristics of DE miRNAs in the key co-expressed module in regulation of mRNAs and revealed new candidate regulatory miRNAs and genes for the beef quality traits.

[Post-partum hypoglycaemia in a child to a father with HNF4A diabetes].

Untreated hypoglycaemia in newborns may result in permanent cognitive damage, why early diagnosis and treatment is important. This case report describes a newborn girl, who developed hypoglycaemia, when she was two hours old despite early feeding. The father of the child had maturity-onset diabetes of the young Type 1, which is caused by an autosomal dominant inherited mutation in the HNF4A gene. Due to this, early blood glucose measurements were performed. The child was treated with extra feeding and recovered without any consequences. A later gene test showed, that the child was carrier of the mutation.

Transcriptomic changes in bovine skeletal muscle cells after resveratrol treatment.

Skeletal muscles constitute a high proportion of the cellular mass that is essential for the growth traits in cattle. Resveratrol (RSV) is a natural polyphenol compound involved in pleiotropic biological activities of muscle. Therefore, the aim of our study was to investigate the transcriptome-level effects of RSV on bovine primary myoblast to reveal differentially expressed genes (DEGs). We treated three replicates of primary myoblasts with 20 µM mother solution containing RSV, whereas three other replicates without RSV were used as control group. Then, we conducted genome-wide transcriptome analysis for the two groups. The results of expression analysis identified 3856 DEGs of which 1805 genes were up-regulated and 2051 genes were down-regulated (adjusted P < 0.05). In addition, qRT-PCR analysis of 19 selected DEGs were consistent with the expression levels observed in the transcriptome data. Gene Ontology (GO) and pathway enrichment analysis showed 72 and 66 significant GO terms and KEGG pathways, respectively (adjusted P < 0.05). The most significant GO term was actin cytoskeleton organization (GO:0030036). The top significant KEGG pathway was focal adhesion (bta04510). Predicted protein-protein interactions (PPIs) showed that CDKN1A encoding cyclindependent kinase inhibitor 1A connects several larger protein complexes. In conclusion, our results found a list of DEGs, significant GO terms and pathways, and provided an improved and expanded understanding of the impact of RSV on cattle muscle cells at the transcriptomic level. The study elucidates the potential of using the genes enriched in pathways mediating resveratrol effects as targets in genomic selection for muscle development and growth in beef cattle.

Copy Number Variations and Expression Levels of Guanylate-Binding Protein 6 Gene Associated with Growth Traits of Chinese Cattle.

Association studies have indicated profound effects of copy number variations (CNVs) on various phenotypes in different species. In this study, we identified the CNV distributions and expression levels of guanylate-binding protein 6 (GBP6) associated with the growth traits of Chinese cattle. The results showed that the phenotypic values of body size and weight of Xianan (XN) cattle were higher than those of Nanyang (NY) cattle. The medium CNV types were mostly identified in the XN and NY breeds, but their CNV distributions were significantly different (adjusted p < 0.05). The association analysis revealed that the body weight, cannon circumference and chest circumference of XN cattle had significantly different values in different CNV types (p < 0.05), with CNV gain types (Log22-ΔΔCt > 0.5) displaying superior phenotypic values. We also found that transcription levels varied in different tissues (p < 0.001) and the CNV gain types showed the highest relative gene expression levels in the muscle tissue, consistent with the highest phenotypic values of body weight and cannon circumference among the three CNV types. Consequently, our results suggested that CNV gain types of GBP6 could be used as the candidate markers in the cattle-breeding program for growth traits.

Expression profiles of the MXD3 gene and association of sequence variants with growth traits in Xianan and Qinchuan cattle.

Max dimerization protein 3 (MXD3) belongs to the MYC superfamily of basic helix-loop-helix leucine zipper transcription factors, and MXD3-MAX heterodimers can bind to promoters of target genes to modulate their expression. The aim of this study was to determine the MXD3 mRNA expression levels in various cattle tissues comprising heart, liver, spleen, lung, kidney, Longissimus dorsi muscle and subcutaneous fat in Chinese Qinchuan and Xianan cattle breeds. The RT-qPCR data showed that the MXD3 gene was variably expressed between all tissues and at levels that were significantly different between two breeds (p < .05). We used the polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) method to investigate the possible association between single-nucleotide polymorphisms (SNP) within the MXD3 gene and five different growth traits in cattle. We found two intronic SNPs (g.2694 C>T and g.3801 T>C) and one SNP in 3'untranslated region (3'UTR) (g.6263 G>A) of MXD3 gene. Association analysis revealed strong associations between pairwise and triple SNP combinations and the growth traits. Based on these results, we suggest that MXD3 polymorphisms could be useful as molecular markers in the Chinese beef cattle breeding program.

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.

MicroRNA-guided prioritization of genome-wide association signals reveals the importance of microRNA-target gene networks for complex traits in cattle.

MicroRNAs (miRNA) are key modulators of gene expression and so act as putative fine-tuners of complex phenotypes. Here, we hypothesized that causal variants of complex traits are enriched in miRNAs and miRNA-target networks. First, we conducted a genome-wide association study (GWAS) for seven functional and milk production traits using imputed sequence variants (13~15 million) and >10,000 animals from three dairy cattle breeds, i.e., Holstein (HOL), Nordic red cattle (RDC) and Jersey (JER). Second, we analyzed for enrichments of association signals in miRNAs and their miRNA-target networks. Our results demonstrated that genomic regions harboring miRNA genes were significantly (P < 0.05) enriched with GWAS signals for milk production traits and mastitis, and that enrichments within miRNA-target gene networks were significantly higher than in random gene-sets for the majority of traits. Furthermore, most between-trait and across-breed correlations of enrichments with miRNA-target networks were significantly greater than with random gene-sets, suggesting pleiotropic effects of miRNAs. Intriguingly, genes that were differentially expressed in response to mammary gland infections were significantly enriched in the miRNA-target networks associated with mastitis. All these findings were consistent across three breeds. Collectively, our observations demonstrate the importance of miRNAs and their targets for the expression of complex traits.

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.

Meta-analysis of genome-wide association studies for cattle stature identifies common genes that regulate body size in mammals.

Stature is affected by many polymorphisms of small effect in humans 1 . In contrast, variation in dogs, even within breeds, has been suggested to be largely due to variants in a small number of genes2,3. Here we use data from cattle to compare the genetic architecture of stature to those in humans and dogs. We conducted a meta-analysis for stature using 58,265 cattle from 17 populations with 25.4 million imputed whole-genome sequence variants. Results showed that the genetic architecture of stature in cattle is similar to that in humans, as the lead variants in 163 significantly associated genomic regions (P < 5 × 10-8) explained at most 13.8% of the phenotypic variance. Most of these variants were noncoding, including variants that were also expression quantitative trait loci (eQTLs) and in ChIP-seq peaks. There was significant overlap in loci for stature with humans and dogs, suggesting that a set of common genes regulates body size in mammals.

The first draft reference genome of the American mink (Neovison vison).

The American mink (Neovison vison) is a semiaquatic species of mustelid native to North America. It's an important animal for the fur industry. Many efforts have been made to locate genes influencing fur quality and color, but this search has been impeded by the lack of a reference genome. Here we present the first draft genome of mink. In our study, two mink individuals were sequenced by Illumina sequencing with 797 Gb sequence generated. Assembly yielded 7,175 scaffolds with an N50 of 6.3 Mb and length of 2.4 Gb including gaps. Repeat sequences constitute around 31% of the genome, which is lower than for dog and cat genomes. The alignments of mink, ferret and dog genomes help to illustrate the chromosomes rearrangement. Gene annotation identified 21,053 protein-coding sequences present in mink genome. The reference genome's structure is consistent with the microsatellite-based genetic map. Mapping of well-studied genes known to be involved in coat quality and coat color, and previously located fur quality QTL provide new knowledge about putative candidate genes for fur traits. The draft genome shows great potential to facilitate genomic research towards improved breeding for high fur quality animals and strengthen our understanding on evolution of Carnivora.

Weight-Averaged Anharmonic Vibrational Analysis of Hydration Structures of Polyamide 6.

Structures of polyamide 6 are investigated for different hydration levels using molecular dynamics (MD) simulations and quantum vibrational calculations. The MD simulations have shown that hydration leads to an increase in the diffusion coefficient, accompanied by a growth of water clusters in the polymer. The IR difference spectra upon hydration are calculated using a weight-averaged method incorporating anharmonicity of the potential energy surface. The predicted IR difference spectrum for the amide A band is in quantitative agreement with the experiment [ Iwamoto , R. ; Murase , H. J. Polym. Sci., Part B: Polym. Phys. 2003 , 41 , 1722 - 1729 ]. The proposed method, combined with experimental IR difference spectra, makes it feasible to elucidate the atomistic structure of hydrated polymer materials.