Miracle-dz/UniFAR-FTUs huggingface dataset:
microRNA-1 represses signaling pathway components to impact embryonic structures derived from all three germ layers

microRNAs are evolutionarily conserved non-coding RNAs that direct post-transcriptional regulation of target transcripts. We use the sea urchin embryo to achieve a comprehensive understanding of miR-1’s function in a developing embryo. Results indicate that miR-1 regulates gut contractions, specification, and positioning of serotonergic neurons, as well as mesodermally-derived muscles, pigment cells, and skeletogenic cells. Gain-of-function of miR-1 generally leads to more severe developmental defects than its loss-of-function. We identified that miR-1 directly suppresses Ets1/2, Tbr, and VegfR7 of the skeletogenic gene regulatory network, and Notch, Nodal, and Wnt1 signaling components. We found that miR-1’s direct suppression of Nodal may indirectly regulate FoxQ2 to impact serotonergic neurons. Excess miR-1 may lead to decreased Nodal and Notch that result in decreased circumpharnygeal muscle fibers and the number of pigment cells. The striking ectopic skeletal branching induced by miR-1 mimic injections may be explained by miR-1’s direct suppression of Nodal that leads to expression changes of Vegf3, and Fgfa that mediate skeletogenesis. This work demonstrates that miR-1 plays a diverse regulatory role that impacts tissues derived from all germ layers. Summary statement This study identifies wide-ranging regulatory roles and regulatory mechanisms of miR-1 that impact structures derived from all three germ layers during embryonic development.
BayesRB: a markov chain Monte Carlo-based polygenic genetic risk score algorithm for dichotomous traits

Identifying high-risk individuals with diseases through reliable prediction models guides screening and preventive treatment. Most complex diseases have a genetic basis influenced by multiple genes and so disease risk can be estimated using polygenic risk score (PRS) algorithms. Many PRS algorithms have been developed so far. Among them, BayesR shows good characteristics of unbiasedness, accuracy, sparseness, and robustness. It detects the associated SNPs, estimates the SNP effects, and makes prediction of disease risks based on all SNPs simultaneously. However, this method assumes that the phenotypes follow a Gaussian distribution, which cannot be met in case-control studies. Here, we made an extension of the BayesR method, called BayesRB, by adding auxiliary variables to the BayesR model. We explored the characteristics, efficacy, and accuracy of BayesRB when estimating SNP effects and predicting disease risks compared with three traditional algorithms under different conditions using both simulated data and real data from the Welcome Trust Case Control Consortium (WTCCC). For SNP effect estimation, BayesRB shows unbiasedness and sparseness for big and small effect SNPs, respectively. For disease risk prediction, BayesRB had the best performance among the methods. This study provides a theoretical basis for complex disease risk prediction and disease prevention
Studies for searching factors relating rumination time around the heat season in Holstein milking cows

To identify factors affecting rumination time (RT), factorial analyses in a dairy farm were performed. In univariate analyses, differences in distribution were observed between low (< 45 kg/d) and high (≥ 45 kg/d) milk yield (MY) before the heat season, early and peak lactation periods and mid and late at the start of the heat season, and 1-2 and ≥ 3 parity at the start of the heat season. Multiple regression analysis confirmed that RT was affected by low MY, early to peak and 1-2 parity cows, low MY, mid to late and 1-2 parity cows, and low MY, early to peak and ≥ 3 parity cows (coefficients = −37.039, −25.353, −44.805 respectively, P < 0.05). When the cows were classified by MY before the heat season, the RT and MY of the high-MY group remained higher than the low-MY group until before the severe heat (THI ≥ 84 continued). However, when cows were classified by RT before the heat season, there was no difference in MY between the low-RT (< 485 min/d) and the high-RT group (> 485 min/d). In conclusions, MY is a factor affecting RT up to moderate heat but RT is not a sufficient condition affecting MY.
In vivo chromatic and spatial tuning of foveolar retinal ganglion cells in Macaca fascicularis

The primate fovea is specialized for high acuity chromatic vision, with the highest density of cone photoreceptors and a disproportionately large representation in visual cortex. The unique visual properties conferred by the fovea are conveyed to the brain by retinal ganglion cells, the somas of which lie at the margin of the foveal pit. Microelectrode recordings of these centermost retinal ganglion cells have been challenging due to the fragility of the fovea in the excised retina. Here we overcome this challenge by combining high resolution fluorescence adaptive optics ophthalmoscopy with calcium imaging to optically record functional responses of foveal retinal ganglion cells in the living eye. We use this approach to study the chromatic responses and spatial transfer functions of retinal ganglion cells using spatially uniform fields modulated in different directions in color space and monochromatic drifting gratings. We repeated these recordings in over 350 cells across three Macaca fascicularis primates over a time frame of weeks to months. We find that the majority of the L vs. M cone opponent cells serving the most central foveolar cones have spatial transfer functions that peak at high spatial frequencies (25-40 c/deg), reflecting strong surround inhibition that sacrifices sensitivity at low frequencies but preserves the transmission of fine detail in the retinal image. In addition, we fit to the drifting grating data a detailed model of how ganglion cell responses draw on the cone mosaic to derive receptive field properties of L vs. M cone opponent cells at the very center of the foveola. The fits are consistent with the hypothesis that midget ganglion cells are specialized to preserve information at the resolution of the cone mosaic. By characterizing the functional properties of retinal ganglion cells in vivo through adaptive optics, we characterize the response characteristics of these cells in situ.
Shaping the genome via lengthwise compaction, phase separation, and lamina adhesion

The link between genomic structure and biological function is yet to be consolidated, it is, however, clear that physical manipulation of the genome, driven by the activity of a variety of proteins, is a crucial step. To understand the consequences of the physical forces underlying genome organization, we build a coarse-grained polymer model of the genome, featuring three fundamentally distinct classes of interactions: lengthwise compaction, i.e., compaction of chromosomes along its contour, self-adhesion among epigenetically similar genomic segments, and adhesion of chromosome segments to the nuclear envelope or lamina. We postulate that these three types of interactions sufficiently represent the concerted action of the different proteins organizing the genome architecture and show that an interplay among these interactions can recapitulate the architectural variants observed across the tree of life. The model elucidates how an interplay of forces arising from the three classes of genomic interactions can drive drastic, yet predictable, changes in the global genome architecture, and makes testable predictions. We posit that precise control over these interactions in vivo is key to the regulation of genome architecture.
Insights on the regulation of photosynthesis in pea leaves exposed to oscillating light

Plants grow in nature often in fluctuating irradiance. In the laboratory, the dynamics of photosynthesis is usually explored by instantaneously exposing dark-adapted plants to constant light and detecting the dark-to-light transient, which is only a poor approximation of natural phenomena. Aiming at a better approximation, we exposed pea leaves to oscillating light and measured, during oscillations, changes in function of photosystem I and II and of the proton-motive force at the thylakoid membrane. The dynamics depends on the oscillation period, leaving information about the regulatory networks. As demonstrated for selected period of the oscillation of 60 s, the regulations try to keep reactions centres of photosystems I and II open. A possible evaluation of obtained data is presented and involvement of particular processes in regulation of photosynthesis is discussed. The forced oscillations provide information-rich fingerprint of complex regulatory networks. Further progress in understanding the networks is expected from experiments involving chemical interventions, plant mutants, and by using mathematical modelling and the system identification and system control tools, as already applied in other parts of science. Highlight Measurement of photosynthetic signals during illumination of plants by light, whose intensity oscillates as sinus function provides information about regulation of photosynthesis in fluctuating light.
Regulation strategies for two-output biomolecular networks

Feedback control theory facilitates the development of self-regulating systems with desired performance which are predictable and insensitive to disturbances. Feedback regulatory topologies are found in many natural systems and have been of key importance in the design of reliable synthetic bio-devices operating in complex biological environments. Here, we study control schemes for biomolecular processes with two outputs of interest, expanding previous traditional concepts describing one-output systems. Regulation of such processes may unlock new design possibilities but can be challenging due to coupling interactions; also potential disturbances applied on one of the outputs may affect both. We therefore propose architectures for robustly manipulating the ratio and linear combinations of the outputs as well as each of the outputs independently. To demonstrate their characteristics, we apply these architectures to a simple process of two mutually activated biomolecular species. We also highlight the potential for experimental implementation by exploring synthetic realizations both in vivo and in vitro. This work presents an important step forward in building bio-devices capable of sophisticated functions.
Development of Dry Powder Formulations of AS01B containing vaccines using Thin-Film Freeze-Drying

AS01B is a liposomal formulation of two immunostimulants namely 3-O-desacyl-4’-monophosphoryl lipid A (MPL) and QS-21. The liposomal formulation of AS01B reduces the endotoxicity of MPL and the lytic activity of QS-21; however, it renders the adjuvant sensitive to accidental slow freezing. The liposomal formulation also represents a major challenge towards the formulation of dry powders of vaccines containing AS01B. In the present study, we tested the feasibility of applying thin-film freeze-drying (TFFD) to engineer dry powders of the AS01B liposomal adjuvant alone or vaccines containing AS01B as an adjuvant. Initially, we showed that after the AS01B liposomal adjuvant was subjected to TFFD using sucrose as a stabilizer at 4% w/v, the particle size distribution of AS01B liposomes reconstituted from the dry powder was identical to the liquid adjuvant before drying. We then showed using ovalbumin (OVA) as a model antigen adjuvanted with AS01B (AS01B/OVA) that subjecting the AS01B/OVA vaccine to TFFD and subsequent reconstitution did not negatively affect the AS01B liposome integrity, nor the immunogenicity of the vaccine. Importantly, the thin-film freeze-dried vaccine was not sensitive to repeated freezing-and-thawing. Finally, the feasibility of using TFFD to prepare dry powders of AS01B-adjuvanted vaccines was further confirmed using AS01B-adjuvanted Fluzone Quadrivalent and Shingrix, which contains AS01B. It is concluded that the TFFD technology can enable the formulation of AS01B-adjuvanted vaccines as freezing-insensitive dry powders in single-vial presentation.
Experimental evidence that animal societies vary in size due to sex differences in cooperation

Cooperative breeding societies allow the costs of reproduction to be shared. However, as groups become larger, such benefits often decrease and competition increases. This is predicted to select for an optimal group size, yet variation in groups is a ubiquitous feature of cooperative breeding animals. Here we experimentally established groups (ngroups=96) of cooperative breeding ostriches, Struthio camelus, with different numbers of males and females and manipulated the potential for cooperation over incubation. There was a clear optimal group size for males. Their reproductive success was maximized in groups with four or more females and no other males, irrespective of cooperation over incubation. Conversely, female reproductive success was strongly dependent on the benefits of cooperating over incubation, being maximized in groups with either many males or many females. In intermediate sized groups, both male and female reproductive success was reduced by sexual conflict over the timing of mating and incubation. Our experiments show that sex differences in the opposing forces of cooperation and competition can explain why variation in cooperative groups is widespread.
Genetic Variation and Disease Severity of Respiratory Syncytial Viruses

Respiratory Syncytial Virus (RSV) disease in newborns ranges from mild symptoms to severe disease requiring hospitalization. RSV is classified into two subtypes (RSVA and RSVB) based on antigenic and genetic differences. The role these genomic variations play in disease severity remains unknown. Genome sequences were obtained using next-generation RNA sequencing on archived frozen nasal swabs of young children (< 8 months-old) infected with RSV in Rochester, NY between 1977-1998. Samples were chosen from both children hospitalized with severe RSV disease (inpatient) and those presenting with mild symptoms (outpatient) during their first cold-season. Both A and B subtypes demonstrated significant differences in the phylogeny and primary-protein structure during this time period. We found a significant association between RSV phylogeny over this time period and disease severity. For both subtypes, the G-protein demonstrated the greatest amino acid substitutions, although the number of amino acid substitutions was higher in the RSVA subtype. We found a significant association between G-protein variation and disease severity for RSVA, but not RSVB. For both subtypes, variation in the M2-2 protein was significantly associated with disease severity. These results suggest that the genetic variability of RSV proteins may contribute to disease severity in humans. Importance Each cold-season Respiratory Syncytial Virus (RSV) infects thousands of children in the US. Some will display mild cold symptoms while others develop severe disease, sometimes resulting in lifelong lung problems or fatality. RSV initiates infection and replicates in the nasopharynx. Substitutions in the RSV genome can be found in clinically isolated nasal-swab samples of RSV infected children. Whether these genome variations contribute to severe disease is unknown. Here we found a statistically significant association between RSV phylogeny and disease severity. Furthermore, we found specific RSV proteins (G and M2-2) whose amino acid variation was statistically associated with severe disease, although which protein was associated depended on subtype. Taken together, our results suggest that RSV genotype contributed to disease severity over this time period.
Dark kinase annotation, mining, and visualization using the Protein Kinase Ontology

The Protein Kinase Ontology (ProKinO) is an integrated knowledge graph that conceptualizes the complex relationships connecting protein kinase sequence, structure, function, and disease in a human and machine-readable format. Here we extend the scope of ProKinO as a discovery tool by including new classes and relationships capturing information on kinase ligand binding sites, expression patterns, and functional features, and demonstrate its application in uncovering new knowledge regarding understudied members of the protein kinase family. Specifically, through graph mining and aggregate SPARQL queries, we identify the p21-activated protein kinase 5 (PAK5) as one of the most frequently mutated dark kinase in human cancers with abnormal expression in multiple cancers, including an unappreciated role in acute myeloid leukemia. We identify recurrent oncogenic mutations in the PAK5 activation loop predicted to alter substrate binding and phosphorylation and identify common ligand/drug binding residues in PAK family kinases, highlighting the potential application of ProKinO in drug discovery. The updated ontology browser and a web component, ProtVista, which allows interactive mining of kinase sequence annotations in 3D structures and Alphafold models, provide a valuable resource for the signaling community. The updated ProKinO database is accessible at http://prokino.uga. edu/browser/.
A novel quantitative trait locus implicates Msh3 in the propensity for genome-wide short tandem repeat expansions in mice

Short tandem repeats (STRs) are a class of rapidly mutating genetic elements characterized by repeated units of 1 or more nucleotides. We leveraged whole genome sequencing data for 152 recombinant inbred (RI) strains from the BXD family derived from C57BL/6J and DBA/2J mice to study the effects of genetic background on genome-wide patterns of new mutations at STRs. We defined quantitative phenotypes describing the numbers and types of germline STR mutations in each strain and identified a locus on chromosome 13 associated with the propensity of STRs to expand. Several dozen genes lie in the QTL region, including Msh3, a known modifier of STR stability at pathogenic repeat expansions in mice and humans. Detailed analysis of the locus revealed a cluster of variants near the 5’ end of Msh3, including multiple protein-coding variants within the DNA mismatch recognition domain of MSH3, and a retrotransposon insertion overlapping an annotated exon. Additionally, gene expression analysis demonstrates co-localization of this QTL with expression QTLs for multiple nearby genes, including Msh3. Our results suggest a novel role for Msh3 in regulating genome-wide patterns of germline STR mutations and demonstrate that inherited genetic variation can contribute to variability in accumulation of new mutations across individuals.
Assessment of dynamics of global DNA methylation during the cryopreservation process of Pinus pinaster embryogenic lines

Long-term in vitro maintenance of embryogenic lines of Pinus species has been associated with lower maturation capacity, because of this, cryopreservation protocols for the embryogenic lines are needed to maintain valuable genotypes. Since cryopreservation may induce epigenetic variations, we evaluate changes in DNA methylation levels through the course of the cryopreservation of maritime pine embryogenic lines, as compared to those lines maintained by repeated subcultures. Six maritime pine embryogenic lines were cryopreserved following a protocol that includes pre-treatments inducing osmotic stress in liquid media. The percentage of methylated cytosines (%5-mC) in total DNA was determined by using a colorimetric assay. Reactive oxygen species (ROS) accumulation in cell lines was also determined by quantifying dihydroethidium (DHE) fluorescence under a confocal laser-scanning microscope. In the first experiment, we found that global DNA methylation was significantly reduced during the cryopreservation protocol. Subsequently, we evaluated the methylation status of both cryopreserved and no cryopreserved lines (maintained by subcultures) and found differences among embryogenic lines but overall, cryopreservation did not alter %5-mC of the recovered lines while periodical subcultures increased methylation rates. In addition, maltose pretreatment did not increase significantly ROS production in embryogenic lines. Our results demonstrate that the genetic stability during cryopreservation highly depends on the embryogenic line studied, but the protocol allows maintaining methylation DNA rates in most of the recovered lines.
Reduced gene templates for supervised analysis of scale-limited CRISPR-Cas9 fitness screens

Pooled genome-wide CRISPR-Cas9 screens are furthering our mechanistic understanding of human biology and have allowed us to identify new oncology therapeutic targets. Scale-limited CRISPR-Cas9 screens – typically employing guide RNA libraries targeting subsets of functionally related genes, individual biological pathways, or portions of the druggable genome – constitute an optimal setting for investigating narrow hypotheses and they are easier to execute on complex models, such as organoids and in vivo models. Different supervised methods are used for the computational analysis of genome-wide CRISPR-Cas9 screens; most are not well suited for scale-limited screens as they require large sets of positive/negative control genes (gene templates) to be included among the screened ones. We have developed a computational framework identifying optimal subsets of known essential and nonessential genes (at different subsampling percentages) that can be used as templates for supervised analyses of scale-limited CRISPR-Cas9 screens, while having a reduced impact on the size of the employed library. Highlights Scale-limited CRISPR-Cas9 screens are experimentally easier than genome-wide screens Reference gene templates are used for supervised analyses of genome-wide screens Reduced templates allow supervised analyses of scale-limited CRISPR-Cas9 screens We present optimal reduced templates and a computational method to assemble them Graphical abstract