Reconstructing Denisovan Anatomy Using DNA Methylation Maps.

Denisovans are an extinct group of humans whose morphology remains unknown. Here, we present a method for reconstructing skeletal morphology using DNA methylation patterns. Our method is based on linking unidirectional methylation changes to loss-of-function phenotypes. We tested performance by reconstructing Neanderthal and chimpanzee skeletal morphologies and obtained >85% precision in identifying divergent traits. We then applied this method to the Denisovan and offer a putative morphological profile. We suggest that Denisovans likely shared with Neanderthals traits such as an elongated face and a wide pelvis. We also identify Denisovan-derived changes, such as an increased dental arch and lateral cranial expansion. Our predictions match the only morphologically informative Denisovan bone to date, as well as the Xuchang skull, which was suggested by some to be a Denisovan. We conclude that DNA methylation can be used to reconstruct anatomical features, including some that do not survive in the fossil record.

Host-Protozoan Interactions Protect from Mucosal Infections through Activation of the Inflammasome.

While conventional pathogenic protists have been extensively studied, there is an underappreciated constitutive protist microbiota that is an integral part of the vertebrate microbiome. The impact of these species on the host and their potential contributions to mucosal immune homeostasis remain poorly studied. Here, we show that the protozoan Tritrichomonas musculis activates the host epithelial inflammasome to induce IL-18 release. Epithelial-derived IL-18 promotes dendritic cell-driven Th1 and Th17 immunity and confers dramatic protection from mucosal bacterial infections. Along with its role as a "protistic" antibiotic, colonization with T. musculis exacerbates the development of T-cell-driven colitis and sporadic colorectal tumors. Our findings demonstrate a novel mutualistic host-protozoan interaction that increases mucosal host defenses at the cost of an increased risk of inflammatory disease.

Identification of constrained sequence elements across 239 primate genomes.

Noncoding DNA is central to our understanding of human gene regulation and complex diseases1,2, and measuring the evolutionary sequence constraint can establish the functional relevance of putative regulatory elements in the human genome3-9. Identifying the genomic elements that have become constrained specifically in primates has been hampered by the faster evolution of noncoding DNA compared to protein-coding DNA10, the relatively short timescales separating primate species11, and the previously limited availability of whole-genome sequences12. Here we construct a whole-genome alignment of 239 species, representing nearly half of all extant species in the primate order. Using this resource, we identified human regulatory elements that are under selective constraint across primates and other mammals at a 5% false discovery rate. We detected 111,318 DNase I hypersensitivity sites and 267,410 transcription factor binding sites that are constrained specifically in primates but not across other placental mammals and validate their cis-regulatory effects on gene expression. These regulatory elements are enriched for human genetic variants that affect gene expression and complex traits and diseases. Our results highlight the important role of recent evolution in regulatory sequence elements differentiating primates, including humans, from other placental mammals.

Current advances in primate genomics: novel approaches for understanding evolution and disease.

Primate genomics holds the key to understanding fundamental aspects of human evolution and disease. However, genetic diversity and functional genomics data sets are currently available for only a few of the more than 500 extant primate species. Concerted efforts are under way to characterize primate genomes, genetic polymorphism and divergence, and functional landscapes across the primate phylogeny. The resulting data sets will enable the connection of genotypes to phenotypes and provide new insight into aspects of the genetics of primate traits, including human diseases. In this Review, we describe the existing genome assemblies as well as genetic variation and functional genomic data sets. We highlight some of the challenges with sample acquisition. Finally, we explore how technological advances in single-cell functional genomics and induced pluripotent stem cell-derived organoids will facilitate our understanding of the molecular foundations of primate biology.

The manifold costs of being a non-native English speaker in science.

The use of English as the common language of science represents a major impediment to maximising the contribution of non-native English speakers to science. Yet few studies have quantified the consequences of language barriers on the career development of researchers who are non-native English speakers. By surveying 908 researchers in environmental sciences, this study estimates and compares the amount of effort required to conduct scientific activities in English between researchers from different countries and, thus, different linguistic and economic backgrounds. Our survey demonstrates that non-native English speakers, especially early in their careers, spend more effort than native English speakers in conducting scientific activities, from reading and writing papers and preparing presentations in English, to disseminating research in multiple languages. Language barriers can also cause them not to attend, or give oral presentations at, international conferences conducted in English. We urge scientific communities to recognise and tackle these disadvantages to release the untapped potential of non-native English speakers in science. This study also proposes potential solutions that can be implemented today by individuals, institutions, journals, funders, and conferences. Please see the Supporting information files (S2-S6 Text) for Alternative Language Abstracts and Figs 5 and 6.

Transcriptome innovations in primates revealed by single-molecule long-read sequencing.

Transcriptomic diversity greatly contributes to the fundamentals of disease, lineage-specific biology, and environmental adaptation. However, much of the actual isoform repertoire contributing to shaping primate evolution remains unknown. Here, we combined deep long- and short-read sequencing complemented with mass spectrometry proteomics in a panel of lymphoblastoid cell lines (LCLs) from human, three other great apes, and rhesus macaque, producing the largest full-length isoform catalog in primates to date. Around half of the captured isoforms are not annotated in their reference genomes, significantly expanding the gene models in primates. Furthermore, our comparative analyses unveil hundreds of transcriptomic innovations and isoform usage changes related to immune function and immunological disorders. The confluence of these evolutionary innovations with signals of positive selection and their limited impact in the proteome points to changes in alternative splicing in genes involved in immune response as an important target of recent regulatory divergence in primates.

Influence of dengue virus serotypes on the abundance of Aedes aegypti insect-specific viruses (ISVs).

A comprehensive understanding of the virome in mosquito vectors is crucial for assessing the potential transmission of viral agents, designing effective vector control strategies, and advancing our knowledge of insect-specific viruses (ISVs). In this study, we utilized Oxford Nanopore Technologies metagenomics to characterize the virome of Aedes aegypti mosquitoes collected in various regions of Colombia, a country hyperendemic for dengue virus (DENV). Analyses were conducted on groups of insects with previous natural DENV infection (DENV-1 and DENV-2 serotypes), as well as mosquito samples that tested negative for virus infection (DENV-negative). Our findings indicate that the Ae. aegypti virome exhibits a similar viral composition at the ISV family and species levels in both DENV-positive and DENV-negative samples across all study sites. However, differences were observed in the relative abundance of viral families such as Phenuiviridae, Partitiviridae, Flaviviridae, Rhabdoviridae, Picornaviridae, Bromoviridae, and Virgaviridae, depending on the serotype of DENV-1 and DENV-2. In addition, ISVs are frequently found in the core virome of Ae. aegypti, such as Phasi Charoen-like phasivirus (PCLV), which was the most prevalent and showed variable abundance in relation to the presence of specific DENV serotypes. Phylogenetic analyses of the L, M, and S segments of the PCLV genome are associated with sequences from different regions of the world but show close clustering with sequences from Brazil and Guadeloupe, indicating a shared evolutionary relationship. The profiling of the Ae. aegypti virome in Colombia presented here improves our understanding of viral diversity within mosquito vectors and provides information that opens the way to possible connections between ISVs and arboviruses. Future studies aimed at deepening our understanding of the mechanisms underlying the interactions between ISVs and DENV serotypes in Ae. aegypti could provide valuable information for the design of effective vector-borne viral disease control and prevention strategies.IMPORTANCEIn this study, we employed a metagenomic approach to characterize the virome of Aedes aegypti mosquitoes, with and without natural DENV infection, in several regions of Colombia. Our findings indicate that the mosquito virome is predominantly composed of insect-specific viruses (ISVs) and that infection with different DENV serotypes (DENV-1 and DENV-2) could lead to alterations in the relative abundance of viral families and species constituting the core virome in Aedes spp. The study also sheds light on the identification of the genome and evolutionary relationships of the Phasi Charoen-like phasivirus in Ae. aegypti in Colombia, a widespread ISV in areas with high DENV incidence.

From contigs towards chromosomes: automatic improvement of long read assemblies (ILRA).

Recent advances in long read technologies not only enable large consortia to aim to sequence all eukaryotes on Earth, but they also allow individual laboratories to sequence their species of interest with relatively low investment. Long read technologies embody the promise of overcoming scaffolding problems associated with repeats and low complexity sequences, but the number of contigs often far exceeds the number of chromosomes and they may contain many insertion and deletion errors around homopolymer tracts. To overcome these issues, we have implemented the ILRA pipeline to correct long read-based assemblies. Contigs are first reordered, renamed, merged, circularized, or filtered if erroneous or contaminated. Illumina short reads are used subsequently to correct homopolymer errors. We successfully tested our approach by improving the genome sequences of Homo sapiens, Trypanosoma brucei, and Leptosphaeria spp., and by generating four novel Plasmodium falciparum assemblies from field samples. We found that correcting homopolymer tracts reduced the number of genes incorrectly annotated as pseudogenes, but an iterative approach seems to be required to correct more sequencing errors. In summary, we describe and benchmark the performance of our new tool, which improved the quality of novel long read assemblies up to 1 Gbp. The pipeline is available at GitHub: https://github.com/ThomasDOtto/ILRA.

New paradigms in the study of the cholinergic system and metabolic diseases: Acetyl-and-butyrylcholinesterase.

Acetylcholinesterase (AChE) and butyrylcholinesterase (BChE) are enzymes that belong to the neuromuscular cholinergic system, their main function is to hydrolyze the neurotransmitter acetylcholine (ACh), through their hydrolysis these enzymes regulate the neuronal and neuromuscular cholinergic system. They have recently attracted considerable attention due to the discovery of new enzymatic and nonenzymatic functions. These discoveries have aroused the interest of numerous scientists, consolidating the relevance of this group of enzymes. Recent investigations have revealed a positive correlation between several risk factors for metabolic syndrome (MetS) and the expression of cholinesterases (ChE's), which underscore the impact of high ChE's activity on the pro-inflammatory state associated with MetS. In addition, the excessive hydrolysis of ACh and other choline esters (succinylcholine, propionylcholine, butyrylcholine, etc.) by both ChE's results in the overproduction of fatty acid precursor metabolites, which facilitate the synthesis of very low-density lipoproteins and triacylglycerols. Participation in these processes may represent the link between ChE's and metabolic disorders. However, further scientific research is required to fully elucidate the involvement of ChE's in metabolic diseases. This review aims to collect recent research studies that contribute to understanding the association between the cholinergic system and metabolic diseases.

Emergence of Dengue Virus Serotype 2 Cosmopolitan Genotype, Colombia.

Using Oxford Nanopore technologies and phylogenetic analyses, we sequenced and identified the cosmopolitan genotype of dengue virus serotype 2 isolated from 2 patients in the city of Villavicencio, Meta department, Colombia. This identification suggests the emergence of this genotype in the country, which warrants further surveillance to identify its epidemic potential.

Wuchereria bancrofti Lymphatic Filariasis, Barrancabermeja, Colombia, 2023.

We describe a recent case of lymphatic filariasis in Colombia caused by Wuchereria bancrofti nematodes. Our study combines clinical-epidemiologic findings with phylogenetic data. Resurgence of lymphatic filariasis may be linked to increasing urbanization trends and migration from previously endemic regions. Fieldwork can be a beneficial tool for screening and containing transmission.

Independent Recruitment of Different Types of Phospholipases A2 to the Venoms of Caenophidian Snakes: The Rise of PLA2-IIE within Pseudoboini (Dipsadidae).

Snake venoms harbor a wide and diverse array of enzymatic and nonenzymatic toxic components, allowing them to exert myriad effects on their prey. However, they appear to trend toward a few optimal compositional scaffolds, dominated by four major toxin classes: SVMPs, SVSPs, 3FTxs, and PLA2s. Nevertheless, the latter appears to be restricted to vipers and elapids, as it has never been reported as a major venom component in rear-fanged species. Here, by investigating the original transcriptomes from 19 species distributed in eight genera from the Pseudoboini tribe (Dipsadidae: Xenodontinae) and screening among seven additional tribes of Dipsadidae and three additional families of advanced snakes, we discovered that a novel type of venom PLA2, resembling a PLA2-IIE, has been recruited to the venom of some species of the Pseudoboini tribe, where it is a major component. Proteomic and functional analyses of these venoms further indicate that these PLA2s play a relevant role in the venoms from this tribe. Moreover, we reconstructed the phylogeny of PLA2s across different snake groups and show that different types of these toxins have been recruited in at least five independent events in caenophidian snakes. Additionally, we present the first compositional profiling of Pseudoboini venoms. Our results demonstrate how relevant phenotypic traits are convergently recruited by different means and from homologous and nonhomologous genes in phylogenetically and ecologically divergent snake groups, possibly optimizing venom composition to overcome diverse adaptative landscapes.

Copy number variation underlies complex phenotypes in domestic dog breeds and other canids.

Extreme phenotypic diversity, a history of artificial selection, and socioeconomic value make domestic dog breeds a compelling subject for genomic research. Copy number variation (CNV) is known to account for a significant part of inter-individual genomic diversity in other systems. However, a comprehensive genome-wide study of structural variation as it relates to breed-specific phenotypes is lacking. We have generated whole genome CNV maps for more than 300 canids. Our data set extends the canine structural variation landscape to more than 100 dog breeds, including novel variants that cannot be assessed using microarray technologies. We have taken advantage of this data set to perform the first CNV-based genome-wide association study (GWAS) in canids. We identify 96 loci that display copy number differences across breeds, which are statistically associated with a previously compiled set of breed-specific morphometrics and disease susceptibilities. Among these, we highlight the discovery of a long-range interaction involving a CNV near MED13L and TBX3, which could influence breed standard height. Integration of the CNVs with chromatin interactions, long noncoding RNA expression, and single nucleotide variation highlights a subset of specific loci and genes with potential functional relevance and the prospect to explain trait variation between dog breeds.

CAAStools: a toolbox to identify and test Convergent Amino Acid Substitutions.

MOTIVATION: Coincidence of Convergent Amino Acid Substitutions (CAAS) with phenotypic convergences allow pinpointing genes and even individual mutations that are likely to be associated with trait variation within their phylogenetic context. Such findings can provide useful insights into the genetic architecture of complex phenotypes. RESULTS: Here we introduce CAAStools, a set of bioinformatics tools to identify and validate CAAS in orthologous protein alignments for predefined groups of species representing the phenotypic values targeted by the user. AVAILABILITY AND IMPLEMENTATION: CAAStools source code is available at http://github.com/linudz/caastools, along with documentation and examples.

A highly resolved nuclear phylogeny uncovers strong phylogenetic conservatism and correlated evolution of fruit color and size in Solanum L.

Mutualisms between plants and fruit-eating animals were key to the radiation of angiosperms. Still, phylogenetic uncertainties limit our understanding of fleshy-fruit evolution, as in the case of Solanum, a genus with remarkable fleshy-fruit diversity, but with unresolved phylogenetic relationships. We used 1786 nuclear genes from 247 species, including 122 newly generated transcriptomes/genomes, to reconstruct the Solanum phylogeny and examine the tempo and mode of the evolution of fruit color and size. Our analysis resolved the backbone phylogeny of Solanum, providing high support for its clades. Our results pushed back the origin of Solanum to 53.1 million years ago (Ma), with most major clades diverging between 35 and 27 Ma. Evolution of Solanum fruit color and size revealed high levels of trait conservatism, where medium-sized berries that remain green when ripe are the likely ancestral form. Our analyses revealed that fruit size and color are evolutionary correlated, where dull-colored fruits are two times larger than black/purple and red fruits. We conclude that the strong phylogenetic conservatism shown in the color and size of Solanum fruits could limit the influences of fruit-eating animals on fleshy-fruit evolution. Our findings highlight the importance of phylogenetic constraints on the diversification of fleshy-fruit functional traits.

An exploratory study of clinical factors associated with IGF-1 and IGFBP-3 in preterm infants.

BACKGROUND: Despite advances in parenteral nutrition, postnatal growth failure in very low birthweight (VLBW) preterm infants is common and associated with chronic health problems. Insulin-like growth factor 1 (IGF-1) is positively associated with improved infant growth, but factors which promote IGF-1 levels in this population have not been clearly identified. The objective of this study was to explore early factors that influence IGF-1 in VLBW preterm infants. METHODS: VLBW infants were enrolled into a prospective, randomized controlled nutrition trial (N = 87). Outcome measures included IGF-1 and IGFBP-3 levels measured at 35 weeks PMA. Linear regression analyses tested the relationships between candidate clinical predictors and levels of IGF-1 and IGFBP-3. RESULTS: Higher protein intake, longer duration of parenteral nutrition, and lower IGFBP-3 levels at 1 week of life were associated with lower IGF-1 levels at 35 weeks PMA. Neither early markers of insulin resistance nor degree of illness were associated with IGF-1 levels at 35 weeks PMA. CONCLUSION: Optimization of early nutrient intake, and attention to route of delivery, may have a lasting influence on IGF-1/IGFBP-3, and in turn, long-term health outcomes. IMPACT: In very low birthweight preterm infants, early protein intake, duration of parenteral nutrition, and insulin-like growth factor binding protein 3 (IGFBP-3) levels at 1 week of life are positively associated with insulin-like growth factor 1 (IGF-1) levels at 35 weeks postmenstrual age. Data from this study highlight the influence of early nutrition on components of the endocrine axis in preterm infants. Strategies aimed at early initiation of enteral nutrition, as well as optimizing composition of parenteral nutrition, may bolster hormones involved in promoting preterm infant growth.

The Impact of Health Inequities on Population-Based Breast Cancer Survival in a Colombian Population, 2008-2015.

OBJECTIVE: To obtain breast cancer survival estimates in Manizales, Colombia, considering socioeconomic level, health insurance regime and residential area, while adjusting for age, histology and stage at diagnosis. METHODS: Analytical cohort study based on breast cancer incident cases recorded by the Population-based Manizales Cancer Registry between 2008-2015. Patients were followed-up for 60 months. Cause-specific survival was calculated using the Kaplan-Meier method for variables of interest, with the Wilcoxon-Breslow-Gehan test for differences. Cox multivariate regression models were fitted. RESULTS: 856 breast cancer cases were included. The 5-year cause-specific survival for the entire cohort was 78.2%. It was higher in women with special/exception health insurance, high socioeconomic level, <50 years old, ductal carcinoma, and stages I and II. Residential area did not impact survival. In Cox models, the subsidized health insurance regime (HR: 4.87 vs contributory) and low socioeconomic level (HR: 2.45 vs high) were predictors of the hazard of death in women with breast cancer, adjusted for age, histology, stage and interactions age-stage and insurance-stage. A positive interaction (synergistic effect modification) between health insurance regime and stage regarding to survival was observed. CONCLUSION: Socioeconomic factors significantly contribute to the inequities in breast cancer survival, independent of the stage at diagnosis. This suggests the need for comprehensive interventions to remove barriers to accessing the health system. This research provides evidence of survival gaps mediated by certain social determinants of health and generates data on the overall performance of the Colombian health system.

Cost-effectiveness of the McGill interactive pediatric oncogenetic guidelines in identifying Li-Fraumeni syndrome in female patients with osteosarcoma.

BACKGROUND: Li-Fraumeni syndrome (LFS) is a penetrant cancer predisposition syndrome (CPS) associated with the development of many tumor types in young people including osteosarcoma and breast cancer (BC). The McGill Interactive Pediatric OncoGenetic Guidelines (MIPOGG) decision-support tool provides a standardized approach to identify patients at risk of CPSs. METHODS: We conducted a cost-utility analysis, from the healthcare payer perspective, to compare MIPOGG-guided, physician-guided, and universal genetic testing strategies to detect LFS in female patients diagnosed at an age of less than 18 years with osteosarcoma. We developed a decision tree and discrete-event simulation model to simulate the clinical and cost outcomes of the three genetic referral strategies on a cohort of female children diagnosed with osteosarcoma, especially focused on BC as subsequent cancer. Outcomes included BC incidence, quality-adjusted life-years (QALYs), healthcare costs, and incremental cost-utility ratios (ICURs). We conducted probabilistic and scenario analyses to assess the uncertainty surrounding model parameters. RESULTS: Compared to the physician-guided testing, the MIPOGG-guided strategy was marginally more expensive by $105 (-$516; $743), but slightly more effective by 0.003 (-0.04; 0.045) QALYs. Compared to MIPOGG, the universal testing strategy was $1333 ($732; $1953) more costly and associated with 0.011 (-0.043; 0.064) additional QALYs. The ICUR for the MIPOGG strategy was $33,947/QALY when compared to the physician strategy; the ICUR for universal testing strategy was $118,631/QALY when compared to the MIPOGG strategy. DISCUSSION: This study provides evidence for clinical and policy decision-making on the cost-effectiveness of genetic referral strategies to identify LFS in the setting of osteosarcoma. MIPOGG-guided strategy was most likely to be cost-effective at a willingness-to-pay threshold value of $50,000/QALY.

Exposure of protected areas in Central America to extreme weather events.

Central America and the Caribbean are regularly battered by megadroughts, heavy rainfall, heat waves, and tropical cyclones. Although 21st-century climate change is expected to increase the frequency, intensity, and duration of these extreme weather events (EWEs), their incidence in regional protected areas (PAs) remains poorly explored. We examined historical and projected EWEs across the region based on 32 metrics that describe distinct dimensions (i.e., intensity, duration, and frequency) of heat waves, cyclones, droughts, and rainfall and compared trends in PAs with trends in unprotected lands. From the early 21st century onward, exposure to EWEs increased across the region, and PAs were predicted to be more exposed to climate extremes than unprotected areas (as shown by autoregressive model coefficients at p < 0.05 significance level). This was particularly true for heat waves, which were projected to have a significantly higher average (tested by Wilcoxon tests at p < 0.01) intensity and duration, and tropical cyclones, which affected PAs more severely in carbon-intensive scenarios. PAs were also predicted to be significantly less exposed to droughts and heavy rainfall than unprotected areas (tested by Wilcoxon tests at p < 0.01). However, droughts that could threaten connectivity between PAs are increasingly common in this region. We estimated that approximately 65% of the study area will experience at least one drought episode that is more intense and longer lasting than previous droughts. Collectively, our results highlight that new conservation strategies adapted to threats associated with EWEs need to be tailored and implemented promptly. Unless urgent action is taken, significant damage may be inflicted on the unique biodiversity of the region.

Ciclones, olas de calor, sequías y lluvias intensas son eventos comunes en Centroamérica y el Caribe, cuya frecuencia, intensidad y duración se espera aumente durante el siglo XXI a causa del cambio climático. Sin embargo, en la actualidad, se desconoce cuál será la incidencia de estos eventos meteorológicos extremos (EME) dentro de las áreas protegidas. En este estudio examinamos la exposición histórica y futura a los extremos climáticos y comparamos el grado de exposición dentro y fuera de las áreas protegidas de toda la región por medio de 32 métricas que describen distintas dimensiones (intensidad, duración y frecuencia) de las olas de calor, los ciclones, las sequías y las precipitaciones. Los resultados indican que a medida que aumente el número de EME, las áreas protegidas estarán más expuestas a los extremos climáticos que las áreas no protegidas. Esto es especialmente cierto en el caso de las olas de calor, que, según las proyecciones, tendrán una intensidad y una duración medias significativamente mayores, y de los ciclones tropicales, que afectarán más gravemente a las zonas protegidas en los escenarios intensivos en carbono. Nuestros resultados también indican que las zonas protegidas estarán significativamente menos expuestas a sequías o lluvias torrenciales que las zonas no protegidas. Sin embargo, las sequías que podrían amenazar la conectividad entre áreas protegidas son cada vez más frecuentes en esta región. Se estima que aproximadamente el 65% del área de estudio experimentará al menos un episodio de sequía más intenso y duradero que las sequías anteriores. En conjunto, nuestros resultados ponen de relieve la necesidad de diseñar y aplicar con prontitud nuevas estrategias de conservación adaptadas a las amenazas asociadas a los EWE. A menos que se tomen medidas urgentes, la biodiversidad única de la región podría sufrir daños considerables.