Miracle-dz/UniFAR-FTUs huggingface dataset:
Reconstructing genotypes from GWAS Summary Statistics

The emergence of genomewide association studies (GWAS) has led to the creation of large repositories of human genetic variation, creating enormous opportunities for genetic research and worldwide collaboration. Methods that are based on GWAS summary statistics seek to leverage such records, overcoming barriers that often exist in individual-level data access while also oﬀering signiﬁcant computational savings. Here, we propose a novel framework that can reconstruct allelic and genotypic counts/frequencies for each SNP from case-control GWAS summary statistics. Our framework is simple and eﬃcient without the need of any complicated underlying assumptions. Illustrating the great potential of this framework we also propose three summary-statistics-based applications implemented in a new software package (ReACt): GWAS meta-analysis (with and without sample overlap), case-case GWAS, and, for the ﬁrst time, group polygenic risk score (PRS) estimation. We evaluate our methods against the current state-of-the-art on both synthetic data and real genotype data and show high performance in power and error control. Our novel group PRS method based on summary statistics could not be achieved prior to our proposed framework. We demonstrate here the potential applications and advantages of this approach. Our work further highlights the great potential of summary-statistics-based methodologies towards elucidating the genetic background of complex disease and opens up new avenues for research.
Species boundaries among extremely diverse and sexually dimorphic Arrenurus water mites (Acariformes: Hydrachnidia: Arrenuridae)

Abstract Arrenurus (Arrenuridae) is the most species-rich genus of mites with about 1000 named species that inhabit standing, and to a lesser extent, running water habitats around the world. To date, distinguishing species of Arrenurus has been based on male reproductive morphology. Here, we use morphological and molecular approaches to examine species boundaries among 42 named species of Arrenurus, including four named species that have colour variants (red and green A. americanus, and red and blue A. intermedius, A. manubriator and A. apetiolatus), and two unnamed morphospecies. In this study, we examine male genital structures with the use of SEM techniques, and apply mitochondrial (COI barcode region) and nuclear (28S rRNA) gene fragments to test whether male morphology reflects species boundaries in Arrenurus assessed by molecular analyses. Our results reveal that male reproductive morphology parallels species boundaries as judged by molecular data. We discuss the cases of genetically poorly diversified, yet morphologically clearly defined named species. Moreover, we show that based on the species we examined, colour morphs within otherwise morphologically similar specimens represent within-species variation and, in the absence of other diagnostic traits, colour itself can be misleading in distinguishing species. Our outcomes on molecular taxonomy of Arrenurus provide a background for testing hypotheses about speciation rate in water mites.
Relish plays a dynamic role in the niche to modulate Drosophila blood progenitor homeostasis in development and infection

Immune challenges demand the gearing up of basal hematopoiesis to combat infection. Little is known about how during development, this switch is achieved to take care of the insult. Here, we show that the hematopoietic niche of the larval lymph gland of Drosophila senses immune challenge and reacts to it quickly through the nuclear factor-κB (NF-κB), Relish, a component of the immune deficiency (Imd) pathway. During development, Relish is triggered by ecdysone signaling in the hematopoietic niche to maintain the blood progenitors. Loss of Relish causes an alteration in the cytoskeletal architecture of the niche cells in a Jun Kinase dependent manner, resulting in the trapping of Hh implicated in progenitor maintenance. Notably, during infection, downregulation of Relish in the niche tilts the maintenance program towards precocious differentiation, thereby bolstering the cellular arm of the immune response.
Neural mechanisms of distributed value representations and learning strategies

Learning appropriate representations of the reward environment is extremely challenging in the real world where there are many options to learn about and these options have many attributes or features. Despite existence of alternative solutions for this challenge, neural mechanisms underlying emergence and adoption of value representations and learning strategies remain unknown. To address this, we measured learning and choice during a novel multi-dimensional probabilistic learning task in humans and trained recurrent neural networks (RNNs) to capture our experimental observations. We found that participants estimate stimulus-outcome associations by learning and combining estimates of reward probabilities associated with the informative feature followed by those of informative conjunctions. Through analyzing representations, connectivity, and lesioning of the RNNs, we demonstrate this mixed learning strategy relies on a distributed neural code and distinct contributions of inhibitory and excitatory neurons. Together, our results reveal neural mechanisms underlying emergence of complex learning strategies in naturalistic settings.
Enhanced cGAS-STING-dependent interferon signaling associated with mutations in ATAD3A

Mitochondrial DNA (mtDNA) has been suggested to drive immune system activation, but the induction of interferon signaling by mtDNA has not been demonstrated in a Mendelian mitochondrial disease. We initially ascertained two patients, one with a purely neurological phenotype, and one with features suggestive of systemic sclerosis in a syndromic context, and found them both to demonstrate enhanced interferon-stimulated gene (ISG) expression in blood. We determined each to harbor a previously described de novo dominant-negative heterozygous mutation in ATAD3A, encoding ATPase family AAA domain-containing protein 3A (ATAD3A). We identified five further patients with mutations in ATAD3A, and recorded up-regulated ISG expression and interferon alpha protein in four of them. Knockdown of ATAD3A in THP-1 cells resulted in increased interferon signaling, mediated by cyclic GMP-AMP synthase (cGAS) and stimulator of interferon genes (STING). Enhanced interferon signaling was abrogated in THP-1 cells and patient fibroblasts depleted of mtDNA. Thus, mutations in the mitochondrial membrane protein ATAD3A define a novel type I interferonopathy. Summary Dominant-negative mutations in ATAD3A, a ubiquitously expressed mitochondrial protein, cause mitochondrial DNA-dependent up-regulation of type I interferon signaling in the context of neurological disease and autoimmunity, thereby defining a novel type I interferonopathy.
The SPARC DRC: Building a Resource for the Autonomic Nervous System Community

The Data and Resource Center (DRC) of the NIH-funded SPARC program is developing databases, connectivity maps and simulation tools for the mammalian autonomic nervous system. The experimental data and mathematical models supplied to the DRC by the SPARC consortium are curated, annotated and semantically linked via a single knowledgebase. A data portal has been developed that allows discovery of data and models both via semantic search and via an interface that includes Google Map-like 2D flatmaps for displaying connectivity, and 3D anatomical organ scaffolds that provide a common coordinate framework for cross-species comparisons. We discuss examples that illustrate the data pipeline, which includes data upload, curation, segmentation (for image data), registration against the flatmaps and scaffolds, and finally display via the web portal, including the link to freely available online computational facilities that will enable neuromodulation hypotheses to be investigated by the autonomic neuroscience community and device manufacturers.
The m6A landscape of polyadenylated nuclear (PAN) RNA and its related methylome in the context of KSHV replication

Polyadenylated nuclear (PAN) RNA is a non-coding transcript involved in Kaposi’s sarcoma-associated herpesvirus (KSHV) lytic reactivation and regulation of cellular and viral gene expression. We have shown that PAN RNA has a dynamic secondary structure and protein binding profiles that can be influenced by the epitranscriptomic modifications. N6-methyladenosine (m6A) is an abundant signature found in viral and virus-encoded RNAs. Here, we combined an antibody-independent next-generation mapping with direct RNA sequencing to elucidate the m6A landscape of PAN RNA during the KSHV latent and lytic stages of infection. Using a newly developed method, termed Selenium-modified deoxythymidine triphosphate reverse transcription and Ligation Assisted PCR analysis of m6A (SLAP), we gained insight into the fraction of modification at identified sites. Using comprehensive proteomic approaches, we identified writers, erasers, and readers that regulate the m6A status of PAN. We verified the temporal and spatial subcellular availability of the methylome components for PAN modification by performing confocal microscopy analysis. Additionally, the RNA biochemical probing outlined structural alterations invoked by m6A in the context of full-length PAN RNA. This work represents the first comprehensive overview of the dynamic interplay between the cellular epitranscriptomic machinery and a specific viral RNA.
Pkd1 and Wnt5a genetically interact to control lymphatic vascular morphogenesis in mice

Lymphatic vascular development is regulated by well-characterised signalling and transcriptional pathways. These pathways regulate lymphatic endothelial cell (LEC) migration, motility, polarity and and morphogenesis. Canonical and non-canonical WNT signalling pathways are known to control LEC polarity and development of lymphatic vessels and valves. PKD1, encoding Polycystin-1, is the most commonly mutated gene in polycystic kidney disease but has also been shown to be essential in lymphatic vascular morphogenesis. The mechanism by which Pkd1 acts during lymphangiogenesis remains unclear. Here we find that loss of non-canonical WNT signalling components Wnt5a and Ryk phenocopy lymphatic defects seen in Pkd1 knockout mice. To investigate genetic interaction, we generated Pkd1/Wnt5a double knockout mice. Loss of Wnt5a suppressed phenotypes seen in the lymphatic vasculature of Pkd1−/− mice and Pkd1 deletion suppressed phenotypes observed in Wnt5a−/− mice. Thus, we report mutually suppressive roles for Pkd1 and Wnt5a, with developing lymphatic networks restored to a more wild-type state in double mutant mice. This genetic interaction between Pkd1 and the non-canonical WNT signalling pathway ultimately controls LEC polarity and the morphogenesis of developing vessel networks. Our work suggests that Pkd1 acts at least in part by regulating non-canonical WNT signalling during the formation of lymphatic vascular networks.
In silico analysis of cis-elements and identification of transcription factors putatively involved in the regulation of the OAS cluster genes SDI1 and SDI2

A. thaliana sulfur deficiency-induced 1 and sulfur deficiency-induced 2 (SDI1 and SDI2) are involved in partitioning sulfur among metabolite pools during sulfur deficiency and their transcription is strongly induced by this condition. However, little is currently known about the cis- and trans-factors that regulate SDI expression. To identify potential transcription factors and DNA sequence element regulators of SDI expression we performed a comparative in silico analysis of their promoter sequences cataloguing known and potentially new cis-elements. We further screened an arrayed library of Arabidopsis transcription factors (TF) for binding to the SDI1 and SDI2 promoters. In total 14 candidate TF regulators of SDIs were identified with yeast-one-hybrid analyses, of which five bound to both promoters, 4 were specific to SDI1, and 5 were specific SDI2. Direct association between particular cis-elements in these promoter regions and specific TFs was established via electrophoretic mobility shift assays. SLIM1 was shown to bind SURE cis-element(s) in the proximal promoter region of both SDI1 and SDI2. The bZIP core cis-element in the proximal promoter region of SDI2 was shown to be important for bZIP16, bZIP44, and HYH binding. GBF1 was shown to bind the E-box in the proximal promoter region of SDI2. Additionally, we performed a meta-analysis of expression changes of these 14 TF candidates in a variety of conditions that alter SDI expression. These data will allow for more detailed future analysis of the molecular factors required for transcriptional regulation of SDIs under a range of physiological and metabolic conditions, apart from sulfur deficiency.
Stimulation with mycobacterial glycolipids and PPD reveals different innate immune response profiles in active and latent TB

Upon infection with Mycobacterium tuberculosis (Mtb) the host immune response might clear the bacteria, control its growth leading to latent tuberculosis (LTB), or fail to control its growth resulting in active TB (ATB). There is however no clear understanding of the features underlying a more or less effective response. Mtb glycolipids are abundant in the bacterial cell envelope and modulate the immune response to Mtb, but the patterns of response to glycolipids are still underexplored. To identify the CD45+ leukocyte activation landscape induced by Mtb glycolipids in peripheral blood of ATB and LTB, we performed a detailed assessment of the immune response of PBMCs to the Mtb glycolipids lipoarabinomannan (LAM) and its biosynthetic precursor phosphatidyl-inositol mannoside (PIM), and PPD. At 24 h and 5 days of stimulation, cell profiling and secretome analysis was done using mass cytometry and high-multiplex immunoassay. PIM mainly affected antigen-presenting cells to produce both proinflammatory (IL-2, IL-6, IL-17A, TNF-α and GM-CSF), and IL-4 and IL-10 cytokines, but not IFN-γ. LAM triggered a similar, albeit weaker, response. By contrast, PPD induced an increase in IFN-γ-producing cells. Moreover, PPD also led to increased numbers of IL-2, IL-6, IL-10, IL-17A, TNF-α and GrzB-producing cells. Treatment with an anti-TLR2 antibody led to partial inhibition of PIM-induced IL-6 production in myeloid cells, suggesting that PIM induces IL-6 production through TLR2. Expansion of monocyte subsets in response to PIM or LAM was reduced in both ATB and LTB as compared to healthy controls, suggesting a hyporesponsive/tolerance pattern in Mtb-infected individuals.
Micrococcal nuclease sequencing of porcine sperm suggests a nucleosomal involvement on semen quality and early embryo development

Background The mammalian mature spermatozoon has a unique chromatin structure in which the vast majority of histones are replaced by protamines during spermatogenesis and a small fraction of nucleosomes are retained at specific locations of the genome. The chromatin structure of sperm remains unresolved in most livestock species, including the pig. However, its resolution could provide further light into the identification of the genomic regions related to sperm biology and embryo development and it could also help identifying molecular markers for sperm quality and fertility traits. Here, for the first time in swine, we performed Micrococcal Nuclease coupled with high throughput sequencing on pig sperm and characterized the mono-nucleosomal (MN) and sub-nucleosomal (SN) chromatin fractions. Results We identified 25,293 and 4,239 peaks in the mono-nucleosomal and sub-nucleosomal fractions, covering 0.3% and 0.02% of the porcine genome, respectively. A cross-species comparison of nucleosome-associated DNAs in sperm revealed positional conservation of the nucleosome retention between human and pig. Gene ontology analysis of the genes mapping nearby the mono-nucleosomal peaks and identification of putative transcription factor binding motifs within the mono-nucleosomal peaks showed enrichment for sperm function and embryo development related processes. We found motif enrichment for the transcription factor Znf263, which in humans was suggested to be a key regulator of the genes with paternal preferential expression during early embryo development. Moreover, we found enriched co-occupancy between the RNAs present in pig sperm and the RNA related to sperm quality, and the mono-nucleosomal peaks. We also found preferential co-location between GWAS hits for semen quality in swine and the mono-nucleosomal sites identified in this study. Conclusions These results suggest a clear relationship between nucleosome positioning in sperm and sperm and embryo development.
Predicting skeletal stature using ancient DNA

Objectives Ancient DNA provides an opportunity to separate the genetic and environmental bases of complex traits by allowing direct estimation of genetic values in ancient individuals. Here, we test whether genetic scores for height in ancient individuals are predictive of their actual height, as inferred from skeletal remains. We estimate the contributions of genetic and environmental variables to observed phenotypic variation as a first step towards quantifying individual sources of morphological variation. Materials and Methods We collected stature estimates and femur lengths from West Eurasian skeletal remains with published genome-wide ancient DNA data (n=167, dating from 33,000-850 BP). We also recorded genetic sex, genetic ancestry, date and paleoclimate data for each individual, and δ13C and δ15N stable isotope values where available (n=67). Results A polygenic score (PRS) for height predicts 6.8% of the variance in femur length in our data (n=117, SD=0.0068%, p<0.001), controlling for sex, ancestry, and date. This is consistent with the predictive power of height PRS in present-day populations and the low coverage of ancient samples. Comparatively, sex explains about 15% of the variance in femur length in our sample. Environmental effects also likely play a role in variation, independent of genetics, though with considerable uncertainty (longitude: R2=0.0317, SD=0.009, p=0.019). Discussion Polygenic scores explain a small but significant proportion of the variance in height in ancient individuals, though not enough to make useful predictions of individual phenotypes. However, environmental variables also contribute to phenotypic outcomes and understanding their interaction with direct genetic predictions will provide a framework with which to model how plasticity and genetic changes ultimately combine to drive adaptation and evolution.
CRISPRpas: programmable regulation of alternative polyadenylation by dCas9

Well over half of human mRNA genes produce alternative polyadenylation (APA) isoforms that differ in mRNA metabolism due to 3’ UTR size changes or have variable coding potentials when coupled with alternative splicing. Aberrant APA is implicated in a growing number of human diseases. A programmable tool for APA regulation, hence, would be instrumental for understanding how APA events impact biological processes. Here, using a catalytically dead Cas9 (dCas9), we developed a method, named CRISPRpas, to alter cleavage and polyadenylation site (PAS) usage in 3’ UTRs or introns. We present key features that facilitate CRISPRpas, including targeting DNA strand, distance between PAS and targeting sequence, and strength of the PAS. For intronic PAS, we additionally analyze strengths of 5’ splice site and target location in intron. Our analyses implicate a dynamic competition between PAS usage and nascent RNA decay when RNA polymerase II elongation is blocked. We show modulation of APA of multiple endogenous genes including a gene that contains a single nucleotide polymorphism (SNP) that affects APA in the human population. CRISPRpas expands the CRISPR toolkit for perturbation of gene expression.
Bipedalism in Mexican Albian lizard (Squamata) and the locomotion type in other Cretaceous lizards

Representative locomotion types in lizards include terrestrial, arboreal, grass swimmer, sand swimmer and bipedal. Few studies explain the locomotion habit of extinct lizards, and even less asses those of bipedal ones. Here, we use quantitative methods to infer the type of locomotion of two Albian Mexican lizards (Lower Cretaceous) and three Cretaceous lizards from Brazil, North America and Spain, assessing the similarities of the hindlimb-forelimb length ratio amongst extinct and extant species. Additionally, an ancestral character state reconstruction analysis was performed, to evaluate the evolution of lizard locomotion habits. The species Huehuecuetzpalli mixtecus was bipedal while Tijubina pontei was facultative bipedal, Hoyalacerta sanzi, Tepexisaurus tepexii and Polyglyphanodon sternbergi cannot be differentiated amongst terrestrial or arboreal with the approach used in this work. The ancestral character state reconstruction analysis showed a terrestrial ancestral locomotion type, with a basal character state of hindlimbs longer than forelimbs. Equal length between hind and forelimbs appear to be a derivate state that evolved multiple times in lizard evolutionary history.
Simultaneous suppression of ribosome biogenesis and Tor activation by TRIM-NHL proteins promotes terminal differentiation

Proper tissue development and homeostasis depends on the balance between growth and terminal differentiation, but the mechanisms coordinating these processes remain elusive. Accumulating evidence indicates that ribosome biogenesis (RiBi) and protein synthesis, two of the most energy-consuming cellular processes sustaining growth, are tightly regulated and yet can be uncoupled during stem cell fate transitions. Here, using the Drosophila adult female germline stem cell (GSC) and larval neuroblast (NB) systems, we show that Mei-P26 and Brat, two Drosophila TRIM-NHL paralogs of the mammalian TRIM32 protein family, are responsible for uncoupling RiBi and protein synthesis during GSC and NB differentiation, respectively. Mei-P26 and Brat modify the metabolism of differentiating cells by activating the Target of rapamycin (Tor) kinase to promote translation, while concomitantly repressing RiBi. Depletion of Mei-P26 or Brat results in excessive cellular growth and defective terminal differentiation, which can be rescued by ectopic activation of Tor together with suppression of RiBi. Our results indicate that the anabolic reprogramming established by TRIM-NHL activity by uncoupling RiBi and translation activities creates the conditions required for terminal differentiation.
Molecular diagnosis of new isolate of tomato yellow leaf curl virus in Iraq

Tomato yield and quality in Iraq have been threatened by a variable range of infections caused by tomato yellow leaf curl virus. In previous studies, the TYLCV isolates were partially characterized using molecular tests for small fragments not the entire length of the virus. Sample of TYLCV-infected tomato has applied in this study to diagnose complete sequence of TYLCV isolate. Three sets of primers that belong to three well-identified strains in Iraq were used in a PCR technique and interestingly the results were negative. A new Iraqi isolate has been characterized as a first novel Iraqi isolate detected ever using next generation (NGS) and bioinformatics techniques. The NGS platform has produced about 78,232,062 paired reads of the TYLCV-infected tomato var. Oula F1. The complete raw reads of the infected variety have been analyzed using RepeatExplorer pipeline and Map to reference. The full sequence of TYLCV was reconstructed and extracted to consist of 2770 nt and then deposited in Genbank under accession number MT583814. Copy numbers and genome proportion of this sequence have been calculated that were 3523 and 0,086% respectively. Phylogenetic analyses of full nucleotide sequences confirmed close relationship to Iranian isolate (TYLCV-Kahnooj) than other published viruses. Additionally, a and β DNA satellites have not discovered in the TYLCV-infected sample.