Evaluation of Single Nucleotide Polymorphisms (SNPs) Associated with Genetic Resistance to Bovine Paratuberculosis in Marchigiana Beef Cattle, an Italian Native Breed.

Mycobacterium avium ssp. paratuberculosis (MAP) is the causative agent of paratuberculosis (PTB), a widespread chronic enteritis of ruminants. The progression of the infection depends on the containment action of innate and cell-mediated immunity (CMI), and it is related to environmental and genetic factors. In particular, PTB susceptibility seems to be associated with specific genes coding for immune regulators involved in the cell-mediated response during the infection. The aim of this preliminary study was to verify, in Italian beef cattle, an association between MAP infectious status and the presence of single nucleotide polymorphisms (SNPs) in candidate genes. To the best of our knowledge, this is the first investigation conducted on a native beef cattle breed, known as Marchigiana, reared in Central Italy. The present research, based on a longitudinal study, aimed to identify and correlate phenotypic and genetic profiles characteristic of the subjects potentially able to contrast or contain PTB. In a MAP-infected herd, ELISA, IFN-γ tests, qPCR, and cultures were performed at a follow-up, occurring within a period ranging from three to six years, to evaluate the individual state of infection. Animals testing positive for at least one test were considered infected. DNA samples of 112 bovines, with known MAP statuses, were analyzed to verify an association with SNPs in the genes encoding gamma-interferon (BoIFNG), interleukin receptor 10 (IL10RA), interleukin receptor 12 (IL12RB2), and toll-like receptors (TLR1, TLR2, TLR4). Regarding statistical analysis, the differences among target genes and pairs of alleles in the analyzed groups of animals, were evaluated at a significance level of p < 0.05. For IL10RA and for IL12RB2 genes, relevant differences in genotypic frequencies among the considered cattle groups were observed. For all candidate genes studied in this investigation, SNP genotypes already associated with PTB resistance were found more frequently in our population, suggesting potential resistance traits in the Marchigiana breed.

InteraChem: Exploring Excited States in Virtual Reality with Ab Initio Interactive Molecular Dynamics.

InteraChem is an ab initio interactive molecular dynamics (AI-IMD) visualizer that leverages recent advances in virtual reality hardware and software, as well as the graphical processing unit (GPU)-accelerated TeraChem electronic structure package, in order to render quantum chemistry in real time. We introduce the exploration of electronically excited states via AI-IMD using the floating occupation molecular orbital-complete active space configuration interaction method. The optimization tools in InteraChem enable identification of excited state minima as well as minimum energy conical intersections for further characterization of excited state chemistry in small- to medium-sized systems. We demonstrate that finite-temperature Hartree-Fock theory is an efficient method to perform ground state AI-IMD. InteraChem allows users to track electronic properties such as molecular orbitals and bond order in real time, resulting in an interactive visualization tool that aids in the interpretation of excited state chemistry data and makes quantum chemistry more accessible for both research and educational purposes.

Voice-controlled quantum chemistry.

Over the past decade, artificial intelligence has been propelled forward by advances in machine learning algorithms and computational hardware, opening up myriads of new avenues for scientific research. Nevertheless, virtual assistants and voice control have yet to be widely used in the natural sciences. Here, we present ChemVox, an interactive Amazon Alexa skill that uses speech recognition to perform quantum chemistry calculations. This new application interfaces Alexa with cloud computing and returns the results through a capable device. ChemVox paves the way to making computational chemistry routinely accessible to the wider community.

Selective bond formation triggered by short optical pulses: quantum dynamics of a four-center ring closure.

We report bond formation induced by an ultrashort UV pulse. The photochemical process is described by quantum dynamics as coherent electronic and nuclear motions during the ultrashort pulse induced ring closure of norbornadiene to quadricyclane. Norbornadiene consists of two ethylene moieties connected by a rigid (CH2)3 bridge. Upon photoexcitation, two new sigma bonds are formed, resulting in the closure of a four-atom ring. As a medium-sized polyatomic molecule, norbornadiene exhibits a high density of strongly coupled electronic states from about 6 eV above the ground state. We report on inducing the formation of the new bonds using a short femtosecond UV pulse to pump a non-equilibrium electronic density in the open form that evolves towards the closed ring form. As the coherent electronic-nuclear coupled dynamics unfold, the excited states change character through non-adiabatic interactions and become valence states for the two new C-C bonds of quadricyclane. Our three-dimensional fully quantum dynamical grid simulations during the first 200 fs show that short UV pulses of different polarization initiate markedly different initial non-equilibrium electronic densities that follow different dynamical paths to the S0/S1 conical intersection. They lead to different initial relative yields of quadricyclane, thereby opening the way to controlling bond-making with attopulses.

Genome-Wide DNA Methylation and Gene Expression Profiles in Cows Subjected to Different Stress Level as Assessed by Cortisol in Milk.

Dairy cattle health, wellbeing and productivity are deeply affected by stress. Its influence on metabolism and immune response is well known, but the underlying epigenetic mechanisms require further investigation. In this study, we compared DNA methylation and gene expression signatures between two dairy cattle populations falling in the high- and low-variant tails of the distribution of milk cortisol concentration (MC), a neuroendocrine marker of stress in dairy cows. Reduced Representation Bisulfite Sequencing was used to obtain a methylation map from blood samples of these animals. The high and low groups exhibited similar amounts of methylated CpGs, while we found differences among non-CpG sites. Significant methylation changes were detected in 248 genes. We also identified significant fold differences in the expression of 324 genes. KEGG and Gene Ontology (GO) analysis showed that genes of both groups act together in several pathways, such as nervous system activity, immune regulatory functions and glucocorticoid metabolism. These preliminary results suggest that, in livestock, cortisol secretion could act as a trigger for epigenetic regulation and that peripheral changes in methylation can provide an insight into central nervous system functions.

Modern quantum chemistry with [Open]Molcas.

MOLCAS/OpenMolcas is an ab initio electronic structure program providing a large set of computational methods from Hartree-Fock and density functional theory to various implementations of multiconfigurational theory. This article provides a comprehensive overview of the main features of the code, specifically reviewing the use of the code in previously reported chemical applications as well as more recent applications including the calculation of magnetic properties from optimized density matrix renormalization group wave functions.

Gallop Racing Shifts Mature mRNA towards Introns: Does Exercise-Induced Stress Enhance Genome Plasticity?

Physical exercise is universally recognized as stressful. Among the "sport species", the horse is probably the most appropriate model for investigating the genomic response to stress due to the homogeneity of its genetic background. The aim of this work is to dissect the whole transcription modulation in Peripheral Blood Mononuclear Cells (PBMCs) after exercise with a time course framework focusing on unexplored regions related to introns and intergenic portions. PBMCs NGS from five 3 year old Sardinian Anglo-Arab racehorses collected at rest and after a 2000 m race was performed. Apart from differential gene expression ascertainment between the two time points the complexity of transcription for alternative transcripts was identified. Interestingly, we noted a transcription shift from the coding to the non-coding regions. We further investigated the possible causes of this phenomenon focusing on genomic repeats, using a differential expression approach and finding a strong general up-regulation of repetitive elements such as LINE. Since their modulation is also associated with the "exonization", the recruitment of repeats that act with regulatory functions, suggesting that there might be an active regulation of this transcriptional shift. Thanks to an innovative bioinformatic approach, our study could represent a model for the transcriptomic investigation of stress.

Was the Cinta Senese Pig Already a Luxury Food in the Late Middle Ages? Ancient DNA and Archaeozoological Evidence from Central Italy.

The Cinta senese is a pig breed, highly esteemed for its meat and derived products, characterized by a black coat with a typical white "belt" and documented by scant iconography, since the 13th-14th century in Italy. A piece of pottery showing a Cinta pig was found in the Graffignano castle (Northern Latium, Italy) dated 15th-16th centuries, spurring us to investigate the diet of the inhabitants. Ancient DNA analysis was carried out on 21 pig specimens on three nuclear SNPs: (1) g.43597545C>T, on the KIT gene, informative for the identification of the Cinta senese breed; (2) rs81460129, on an intergenic region in chr. 16, which discriminates between domestic pigs and wild boars, and; (3) a SNP on the ZFY/ZFX homologous genes, to determine the sex of the individuals. Our results indicate that the Cinta senese was present in Northern Latium in Late Medieval time, although it was not the only breed, and that pigs, including Cinta, interbred with wild boars, suggesting free-range breeding for all types of pigs. Moreover, the unexpected high proportion of young females may be considered as evidence for the wealth of the family inhabiting the castle.

Web-ARM: A Web-Based Interface for the Automatic Construction of QM/MM Models of Rhodopsins.

This article introduces Web-ARM, a specialized tool, online available, designed to build quantum mechanical/molecular mechanical models of rhodopsins, a widely spread family of light-responsive proteins. Web-ARM allows the rapidly building of models of rhodopsins with a documented quality and the prediction of trends in UV-vis absorption maximum wavelengths, based on their excitation energies computed at the CASPT2//CASSCF/Amber level of theory. Web-ARM builds upon the recently reported, python-based a-ARM protocol [J. Chem. Theory Comput., 2019, 15, 3134-3152] and, as such, necessitates only a crystallographic structure or a comparative model in PDB format and a very basic knowledge of the studied rhodopsin system. The user-friendly web interface uses such input to generate congruous, gas-phase models of rhodopsins and, if requested, their mutants. We present two possible applications of Web-ARM, which showcase how the interface can be employed to assist both research and educational activities in fields at the interface between chemistry and biology. The first application shows how, through Web-ARM, research projects (e.g., rhodopsin and rhodopsin mutant screening) can be carried out in significantly less time with respect to using the required computational photochemistry tools via a command line. The second application documents the use of Web-ARM in a real-life educational/training activity, through a hands-on experience illustrating the concepts of rhodopsin color tuning.

Fluorescence Enhancement of a Microbial Rhodopsin via Electronic Reprogramming.

The engineering of microbial rhodopsins with enhanced fluorescence is of great importance in the expanding field of optogenetics. Here we report the discovery of two mutants (W76S/Y179F and L83Q) of a sensory rhodopsin from the cyanobacterium Anabaena PCC7120 with opposite fluorescence behavior. In fact, while W76S/Y179F displays, with respect to the wild-type protein, a nearly 10-fold increase in red-light emission, the second is not emissive. Thus, the W76S/Y179F, L83Q pair offers an unprecedented opportunity for the investigation of fluorescence enhancement in microbial rhodopsins, which is pursued by combining transient absorption spectroscopy and multiconfigurational quantum chemistry. The results of such an investigation point to an isomerization-blocking electronic effect as the direct cause of instantaneous (subpicosecond) fluorescence enhancement.

New Insights on Water Buffalo Genomic Diversity and Post-Domestication Migration Routes From Medium Density SNP Chip Data.

The domestic water buffalo is native to the Asian continent but through historical migrations and recent importations, nowadays has a worldwide distribution. The two types of water buffalo, i.e., river and swamp, display distinct morphological and behavioral traits, different karyotypes and also have different purposes and geographical distributions. River buffaloes from Pakistan, Iran, Turkey, Egypt, Romania, Bulgaria, Italy, Mozambique, Brazil and Colombia, and swamp buffaloes from China, Thailand, Philippines, Indonesia and Brazil were genotyped with a species-specific medium-density 90K SNP panel. We estimated the levels of molecular diversity and described population structure, which revealed historical relationships between populations and migration events. Three distinct gene pools were identified in pure river as well as in pure swamp buffalo populations. Genomic admixture was seen in the Philippines and in Brazil, resulting from importations of animals for breed improvement. Our results were largely consistent with previous archeological, historical and molecular-based evidence for two independent domestication events for river- and swamp-type buffaloes, which occurred in the Indo-Pakistani region and close to the China/Indochina border, respectively. Based on a geographical analysis of the distribution of diversity, our evidence also indicated that the water buffalo spread out of the domestication centers followed two major divergent migration directions: river buffaloes migrated west from the Indian sub-continent while swamp buffaloes migrated from northern Indochina via an east-south-eastern route. These data suggest that the current distribution of water buffalo diversity has been shaped by the combined effects of multiple migration events occurred at different stages of the post-domestication history of the species.

Genome assembly and transcriptome resource for river buffalo, Bubalus bubalis (2n = 50).

Water buffalo is a globally important species for agriculture and local economies. A de novo assembled, well-annotated reference sequence for the water buffalo is an important prerequisite for studying the biology of this species, and is necessary to manage genetic diversity and to use modern breeding and genomic selection techniques. However, no such genome assembly has been previously reported. There are 2 species of domestic water buffalo, the river (2 n = 50) and the swamp (2 n = 48) buffalo. Here we describe a draft quality reference sequence for the river buffalo created from Illumina GA and Roche 454 short read sequences using the MaSuRCA assembler. The assembled sequence is 2.83 Gb, consisting of 366 983 scaffolds with a scaffold N50 of 1.41 Mb and contig N50 of 21 398 bp. Annotation of the genome was supported by transcriptome data from 30 tissues and identified 21 711 predicted protein coding genes. Searches for complete mammalian BUSCO gene groups found 98.6% of curated single copy orthologs present among predicted genes, which suggests a high level of completeness of the genome. The annotated sequence is available from NCBI at accession GCA_000471725.1.

How Do Methyl Groups Enhance the Triplet Chemiexcitation Yield of Dioxetane?

Chemiluminescence is the emission of light as a result of a nonadiabatic chemical reaction. The present work is concerned with understanding the yield of chemiluminescence, in particular how it dramatically increases upon methylation of 1,2-dioxetane. Both ground-state and nonadiabatic dynamics (including singlet excited states) of the decomposition reaction of various methyl-substituted dioxetanes have been simulated. Methyl-substitution leads to a significant increase in the dissociation time scale. The rotation around the O-C-C-O dihedral angle is slowed; thus, the molecular system stays longer in the "entropic trap" region. A simple kinetic model is proposed to explain how this leads to a higher chemiluminescence yield. These results have important implications for the design of efficient chemiluminescent systems in medical, environmental, and industrial applications.

De novo assembly, functional annotation, and analysis of the giant reed (Arundo donax L.) leaf transcriptome provide tools for the development of a biofuel feedstock.

BACKGROUND: Arundo donax has attracted renewed interest as a potential candidate energy crop for use in biomass-to-liquid fuel conversion processes and biorefineries. This is due to its high productivity, adaptability to marginal land conditions, and suitability for biofuel and biomaterial production. Despite its importance, the genomic resources currently available for supporting the improvement of this species are still limited. RESULTS: We used RNA sequencing (RNA-Seq) to de novo assemble and characterize the A. donax leaf transcriptome. The sequencing generated 1249 million clean reads that were assembled using single-k-mer and multi-k-mer approaches into 62,596 unique sequences (unitranscripts) with an N50 of 1134 bp. TransDecoder and Trinotate software suites were used to obtain putative coding sequences and annotate them by mapping to UniProtKB/Swiss-Prot and UniRef90 databases, searching for known transcripts, proteins, protein domains, and signal peptides. Furthermore, the unitranscripts were annotated by mapping them to the NCBI non-redundant, GO and KEGG pathway databases using Blast2GO. The transcriptome was also characterized by BLAST searches to investigate homologous transcripts of key genes involved in important metabolic pathways, such as lignin, cellulose, purine, and thiamine biosynthesis and carbon fixation. Moreover, a set of homologous transcripts of key genes involved in stomatal development and of genes coding for stress-associated proteins (SAPs) were identified. Additionally, 8364 simple sequence repeat (SSR) markers were identified and surveyed. SSRs appeared more abundant in non-coding regions (63.18%) than in coding regions (36.82%). This SSR dataset represents the first marker catalogue of A. donax. 53 SSRs (PolySSRs) were then predicted to be polymorphic between ecotype-specific assemblies, suggesting genetic variability in the studied ecotypes. CONCLUSIONS: This study provides the first publicly available leaf transcriptome for the A. donax bioenergy crop. The functional annotation and characterization of the transcriptome will be highly useful for providing insight into the molecular mechanisms underlying its extreme adaptability. The identification of homologous transcripts involved in key metabolic pathways offers a platform for directing future efforts in genetic improvement of this species. Finally, the identified SSRs will facilitate the harnessing of untapped genetic diversity. This transcriptome should be of value to ongoing functional genomics and genetic studies in this crop of paramount economic importance.

Optomechanical Control of Quantum Yield in Trans-Cis Ultrafast Photoisomerization of a Retinal Chromophore Model.

The quantum yield of a photochemical reaction is one of the most fundamental quantities in photochemistry, as it measures the efficiency of the transduction of light energy into chemical energy. Nature has evolved photoreceptors in which the reactivity of a chromophore is enhanced by its molecular environment to achieve high quantum yields. The retinal chromophore sterically constrained inside rhodopsin proteins represents an outstanding example of such a control. In a more general framework, mechanical forces acting on a molecular system can strongly modify its reactivity. Herein, we show that the exertion of tensile forces on a simplified retinal chromophore model provokes a substantial and regular increase in the trans-to-cis photoisomerization quantum yield in a counterintuitive way, as these extension forces facilitate the formation of the more compressed cis photoisomer. A rationale for the mechanochemical effect on this photoisomerization mechanism is also proposed.

Toward Automatic Rhodopsin Modeling as a Tool for High-Throughput Computational Photobiology.

We report on a prototype protocol for the automatic and fast construction of congruous sets of QM/MM models of rhodopsin-like photoreceptors and of their mutants. In the present implementation the information required for the construction of each model is essentially a crystallographic structure or a comparative model complemented with information on the protonation state of ionizable side chains and distributions of external counterions. Starting with such information, a model formed by a fixed environment system, a flexible cavity system, and a chromophore system is automatically generated. The results of the predicted vertical excitation energy for 27 different rhodopsins including vertebrate, invertebrate, and microbial pigments indicate that such basic models could be employed for predicting trends in spectral changes and/or correlate the spectral changes with structural variations in large sets of proteins.

Use of canonical discriminant analysis to study signatures of selection in cattle.

BACKGROUND: Cattle include a large number of breeds that are characterized by marked phenotypic differences and thus constitute a valuable model to study genome evolution in response to processes such as selection and domestication. Detection of "signatures of selection" is a useful approach to study the evolutionary pressures experienced throughout history. In the present study, signatures of selection were investigated in five cattle breeds farmed in Italy using a multivariate approach. METHODS: A total of 4094 bulls from five breeds with different production aptitudes (two dairy breeds: Italian Holstein and Italian Brown Swiss; two beef breeds: Piemontese and Marchigiana; and one dual purpose breed: Italian Simmental) were genotyped using the Illumina BovineSNP50 v.1 beadchip. Canonical discriminant analysis was carried out on the matrix of single nucleotide polymorphisms (SNP) genotyping data, separately for each chromosome. Scores for each canonical variable were calculated and then plotted in the canonical space to quantify the distance between breeds. SNPs for which the correlation with the canonical variable was in the 99th percentile for a specific chromosome were considered to be significantly associated with that variable. Results were compared with those obtained using an FST-based approach. RESULTS: Based on the results of the canonical discriminant analysis, a large number of signatures of selection were detected, among which several had strong signals in genomic regions that harbour genes known to have an impact on production and morphological bovine traits, including MSTN, LCT, GHR, SCD, NCAPG, KIT, and ASIP. Moreover, new putative candidate genes were identified, such as GCK, B3GALNT1, MGAT1, GALNTL1, PRNP, and PRND. Similar results were obtained with the FST-based approach. CONCLUSIONS: The use of canonical discriminant analysis on 50 K SNP genotypes allowed the extraction of new variables that maximize the separation between breeds. This approach is quite straightforward, it can compare more than two groups simultaneously, and relative distances between breeds can be visualized. The genes that were highlighted in the canonical discriminant analysis were in concordance with those obtained using the FST index.

Space and Time Evolution of the Electrostatic Potential During the Activation of a Visual Pigment.

Animal and microbial retinal proteins employ the Schiff base of retinal as their chromophore. Here, the possible consequences of the charge translocation associated with the light-induced dynamics of the chromophore of a visual opsin are investigated along a representative semiclassical trajectory. We show that the evolution of the electrostatic potential projected by the chromophore onto the surrounding protein displays intense but topographically localized sudden variations in proximity of the decay region. pKa calculations carried out on selected snapshots used as probes, indicate that the only residue which may be sensitive to the electrostatic potential shift is Glu181. Accordingly, our results suggest that the frail Tyr191/268-Glu181-Wat2-Ser186 hydrogen bond network may be perturbed by the transient variations of the electrostatic potential.

Structural and Dynamic Characterization of the C313Y Mutation in Myostatin Dimeric Protein, Responsible for the "Double Muscle" Phenotype in Piedmontese Cattle.

The knowledge of the molecular effects of the C313Y mutation, responsible for the "double muscle" phenotype in Piedmontese cattle, can help understanding the actual mechanism of phenotype determination and paves the route for a better modulation of the positive effects of this economic important phenotype in the beef industry, while minimizing the negative side effects, now inevitably intersected. The structure and dynamic behavior of the active dimeric form of Myostatin in cattle was analyzed by means of three state-of-the-art Molecular Dynamics simulations, 200-ns long, of wild-type and C313Y mutants. Our results highlight a role for the conserved Arg333 in establishing a network of short and long range interactions between the two monomers in the wild-type protein that is destroyed upon the C313Y mutation even in a single monomer. Furthermore, the native protein shows an asymmetry in residue fluctuation that is absent in the double monomer mutant. Time window analysis on further 200-ns of simulation demonstrates that this is a characteristic behavior of the protein, likely dependent on long range communications between monomers. The same behavior, in fact, has already been observed in other mutated dimers. Finally, the mutation does not produce alterations in the secondary structure elements that compose the characteristic TGF-ß cystine-knot motif.

Molcas 8: New capabilities for multiconfigurational quantum chemical calculations across the periodic table.

In this report, we summarize and describe the recent unique updates and additions to the Molcas quantum chemistry program suite as contained in release version 8. These updates include natural and spin orbitals for studies of magnetic properties, local and linear scaling methods for the Douglas-Kroll-Hess transformation, the generalized active space concept in MCSCF methods, a combination of multiconfigurational wave functions with density functional theory in the MC-PDFT method, additional methods for computation of magnetic properties, methods for diabatization, analytical gradients of state average complete active space SCF in association with density fitting, methods for constrained fragment optimization, large-scale parallel multireference configuration interaction including analytic gradients via the interface to the Columbus package, and approximations of the CASPT2 method to be used for computations of large systems. In addition, the report includes the description of a computational machinery for nonlinear optical spectroscopy through an interface to the QM/MM package Cobramm. Further, a module to run molecular dynamics simulations is added, two surface hopping algorithms are included to enable nonadiabatic calculations, and the DQ method for diabatization is added. Finally, we report on the subject of improvements with respects to alternative file options and parallelization.