Transposable element landscapes illuminate past evolutionary events in the endangered fern Vandenboschia speciosa

Vandenboschia speciosa is an endangered tetraploid fern species with a large genome (10.5 Gb). Its geographical distribution is characterized by disjoined tertiary flora refuges, with relict populations that survived past climate crises. Here we analyze the transposable elements (TEs) and found that they comprise about 76% of the V. speciosa genome, thus being the most abundant kind of DNA sequences in this gigantic genome. V. speciosa genome is composed of 51% and 5.6% of Class I and Class II elements, respectively. LTR retrotransposons were the most abundant TEs in this species (at least 42% of the genome), followed by non-LTR retrotransposons that constituted at least 8.7% of the genome of this species. We introduce an additional analysis to identify the nature of non-annotated elements (19% of the genome). A BLAST search of the non-annotated contigs against the V. speciosa TE database allowed determining the identity of almost half of them, which were most likely diverged sequence variants of the annotated TEs. In general, TE composition in V. speciosa resembles TE composition in seed plants. In addition, repeat landscapes revealed three episodes of amplification for all TEs, most likely due to demographic changes associated to past climate crises

Eight Million Years of Satellite DNA Evolution in Grasshoppers of the Genus Schistocerca Illuminate the Ins and Outs of the Library Hypothesis

This study was supported by the Fundação de Amparo à Pesquisa do Estado de São Paulo-FAPESP (process number 2014/11763-8) and Coordenadoria de Aperfeiçoamento de Pessoal de Nível Superior-CAPES. The authors are grateful to anonymous reviewers for valuable comments. O.M.P.-G. and F.J.R.-R. acknowledge the scholarship obtained from the Lawski Foundation (Sweden), H.S. to the U.S. National Science Foundation (Grant No. IOS-1253493) and the United State Department of Agriculture (Hatch Grant TEX0-1-6584), D.A.M. was supported by Consejo Nacional de Investigaciones Científicas y Técnicas-CONICET from Argentina. D.C.C.-d.-M. is recipient of a research productivity fellowship from the Conselho Nacional de Desenvolvimento Científico e Tecnológico-CNPq (process number 304758/2014-0). Computing ran on resources provided by the Swedish National Infrastructure for Computing (SNIC) through Uppsala Multidisciplinary Center for Advanced Computational Science (UPPMAX). We also thank Alexander Suh (Uppsala University, Uppsala/Sweden) for helpful comments on an earlier version of this manuscript and to Cristiane Mileo (UNESP-Rio Claro/SP) for helping with some figures.Satellite DNA (satDNA) is an abundant class of tandemly repeated noncoding sequences, showing high rate of change in sequence, abundance, and physical location. However, the mechanisms promoting these changes are still controversial. The library model was put forward to explain the conservation of some satDNAs for long periods, predicting that related species share a common collection of satDNAs, which mostly experience quantitative changes. Here, we tested the library model by analyzing three satDNAs in ten species of Schistocerca grasshoppers. This group represents a valuable material because it diversified during the last 7.9 Myr across the American continent from the African desert locust (Schistocerca gregaria), and this thus illuminates the direction of evolutionary changes. By combining bioinformatic and cytogenetic, we tested whether these three satDNA families found in S. gregaria are also present in nine American species, and whether differential gains and/or losses have occurred in the lineages. We found that the three satDNAs are present in all species but display remarkable interspecies differences in their abundance and sequences while being highly consistent with genus phylogeny. The number of chromosomal loci where satDNA is present was also consistent with phylogeny for two satDNA families but not for the other. Our results suggest eminently chance events for satDNA evolution. Several evolutionary trends clearly imply either massive amplifications or contractions, thus closely fitting the library model prediction that changes are mostly quantitative. Finally, we found that satDNA amplifications or contractions may influence the evolution of monomer consensus sequences and by chance playing a major role in driftlike dynamics.Fundação de Amparo à Pesquisa do Estado de São Paulo-FAPESP 2014/11763-8Coordenadoria de Aperfeiçoamento de Pessoal de Nível Superior-CAPESLawski Foundation (Sweden)U.S. National Science Foundation IOS-1253493United State Department of Agriculture TEX0-1-6584CONICET from ArgentinaCNPq 304758/2014-

Satellitome comparison of two oedipodine grasshoppers highlights the contingent nature of satellite DNA evolution

Background: The full catalog of satellite DNA (satDNA) within a same genome constitutes the satellitome. The Library Hypothesis predicts that satDNA in relative species reflects that in their common ancestor, but the evolutionary mechanisms and pathways of satDNA evolution have never been analyzed for full satellitomes. We compare here the satellitomes of two Oedipodine grasshoppers (Locusta migratoria and Oedaleus decorus) which shared their most recent common ancestor about 22.8 Ma ago. Results: We found that about one third of their satDNA families (near 60 in every species) showed sequence homology and were grouped into 12 orthologous superfamilies. The turnover rate of consensus sequences was extremely variable among the 20 orthologous family pairs analyzed in both species. The satDNAs shared by both species showed poor association with sequence signatures and motives frequently argued as functional, except for short inverted repeats allowing short dyad symmetries and non-B DNA conformations. Orthologous satDNAs frequently showed different FISH patterns at both intra- and interspecific levels. We defined indices of homogenization and degeneration and quantified the level of incomplete library sorting between species. Conclusions: Our analyses revealed that satDNA degenerates through point mutation and homogenizes through partial turnovers caused by massive tandem duplications (the so-called satDNA amplification). Remarkably, satDNA amplification increases homogenization, at intragenomic level, and diversification between species, thus constituting the basis for concerted evolution. We suggest a model of satDNA evolution by means of recursive cycles of amplification and degeneration, leading to mostly contingent evolutionary pathways where concerted evolution emerges promptly after lineages split

Integrated flow cytometry and sequencing to reconstruct evolutionary patterns from dysplasia to acute myeloid leukemia

Clonal evolution in acute myeloid leukemia (AML) originates long before diagnosis and is a dynamic process that may affect survival. However, it remains uninvestigated during routine diagnostic workups. We hypothesized that the mutational status of bone marrow dysplastic cells and leukemic blasts, analyzed at the onset of AML using integrated multidimensional flow cytometry (MFC) immunophenotyping and fluorescence-activated cell sorting (FACS) with next-generation sequencing (NGS), could reconstruct leukemogenesis. Dysplastic cells were detected by MFC in 285 of 348 (82%) newly diagnosed patients with AML. Presence of dysplasia according to MFC and World Health Organization criteria had no prognostic value in older adults. NGS of dysplastic cells and blasts isolated at diagnosis identified 3 evolutionary patterns: stable (n = 12 of 21), branching (n = 4 of 21), and clonal evolution (n = 5 of 21). In patients achieving complete response (CR), integrated MFC and FACS with NGS showed persistent measurable residual disease (MRD) in phenotypically normal cell types, as well as the acquisition of genetic traits associated with treatment resistance. Furthermore, whole-exome sequencing of dysplastic and leukemic cells at diagnosis and of MRD uncovered different clonal involvement in dysplastic myelo-erythropoiesis, leukemic transformation, and chemoresistance. Altogether, we showed that it is possible to reconstruct leukemogenesis in ∼80% of patients with newly diagnosed AML, using techniques other than single-cell multiomics.This work was supported by grants from the Área de Oncología del Instituto de Salud Carlos III, Centro de Investigación Biomédica en Red (CIBER-ONC) (CB16/12/00369, CB16/12/00233, CB16/12/00489, and CB16/12/00284), Instituto de Salud Carlos III/Subdirección General de Investigación Sanitaria (FIS numbers PI16/01661, PI16/00517, and PI19/01518), and the Plan de Investigación de la Universidad de Navarra (PIUNA 2014-18). This work was supported internationally by the Cancer Research UK, FCAECC, and AIRC under the Accelerator Award Program (EDITOR)

The Hoopoe's Uropygial Gland Hosts a Bacterial Community Influenced by the Living Conditions of the Bird

Molecular methods have revealed that symbiotic systems involving bacteria are mostly based on whole bacterial communities. Bacterial diversity in hoopoe uropygial gland secretion is known to be mainly composed of certain strains of enterococci, but this conclusion is based solely on culture-dependent techniques. This study, by using culture-independent techniques (based on the 16S rDNA and the ribosomal intergenic spacer region) shows that the bacterial community in the uropygial gland secretion is more complex than previously thought and its composition is affected by the living conditions of the bird. Besides the known enterococci, the uropygial gland hosts other facultative anaerobic species and several obligated anaerobic species (mostly clostridia). The bacterial assemblage of this community was largely invariable among study individuals, although differences were detected between captive and wild female hoopoes, with some strains showing significantly higher prevalence in wild birds. These results alter previous views on the hoopoe-bacteria symbiosis and open a new window to further explore this system, delving into the possible sources of symbiotic bacteria (e.g. nest environments, digestive tract, winter quarters) or the possible functions of different bacterial groups in different contexts of parasitism or predation of their hoopoe host.This work was supported by the Ministerio de Ciencia y Tecnología (projects CGL2005-06975/BOSFEDER; CGL2007-61251/BOSFEDER), the Ministerio de Ciencia e Innovación (projects CGL2009-14006/BOSFEDER; CGL2010-19233-C03-01/BOSFEDER; CGL2010-19233-C03-03/BOSFEDER), the Ministerio de Economía y Competitividad (projects CGL2013-48193-C3-1-P/BOSFEDER; CGL2013-48193-C3-2-P/BOSFEDER), and the Junta de Andalucía (RNM 345, P09-RNM-4557). SMRR received a grant from the Ministerio de Ciencia e Innovación (FPI program, BES-2011-047677)

Preferential Occupancy of R2 Retroelements on the B Chromosomes of the Grasshopper Eyprepocnemis plorans

R2 non-LTR retrotransposons exclusively insert into the 28S rRNA genes of their host, and are expressed by co-transcription with the rDNA unit. The grasshopper Eyprepocnemis plorans contains transcribed rDNA clusters on most of its A chromosomes, as well as non-transcribed rDNA clusters on the parasitic B chromosomes found in many populations. Here the structure of the E. plorans R2 element, its abundance relative to the number of rDNA units and its retrotransposition activity were determined. Animals screened from five populations contained on average over 12,000 rDNA units on their A chromosomes, but surprisingly only about 100 R2 elements. Monitoring the patterns of R2 insertions in individuals from these populations revealed only low levels of retrotransposition. The low rates of R2 insertion observed in E. plorans differ from the high levels of R2 insertion previously observed in insect species that have many fewer rDNA units. It is proposed that high levels of R2 are strongly selected against in E. plorans, because the rDNA transcription machinery in this species is unable to differentiate between R2-inserted and uninserted units. The B chromosomes of E. plorans contain an additional 7,000 to 15,000 rDNA units, but in contrast to the A chromosomes, from 150 to over 1,500 R2 elements. The higher concentration of R2 in the inactive B chromosomes rDNA clusters suggests these chromosomes can act as a sink for R2 insertions thus further reducing the level of insertions on the A chromosomes. These studies suggest an interesting evolutionary relationship between the parasitic B chromosomes and R2 elements.This study was supported by grants from the Spanish Ministerio de Ciencia y Tecnología (CGL2009-11917) and Plan Andaluz de Investigacion (CVI-6649), and was partially performed by FEDER funds and a grant from the National Institutes of Health (GM42790)

Delimiting the Origin of a B Chromosome by FISH Mapping, Chromosome Painting and DNA Sequence Analysis in Astyanax paranae (Teleostei, Characiformes)

Supernumerary (B) chromosomes have been shown to contain a wide variety of repetitive sequences. For this reason, fluorescent in situ hybridisation (FISH) is a useful tool for ascertaining the origin of these genomic elements, especially when combined with painting from microdissected B chromosomes. In order to investigate the origin of B chromosomes in the fish species Astyanax paranae, these two approaches were used along with PCR amplification of specific DNA sequences obtained from the B chromosomes and its comparison with those residing in the A chromosomes. Remarkably, chromosome painting with the one-arm metacentric B chromosome probe showed hybridization signals on entire B chromosome, while FISH mapping revealed the presence of H1 histone and 18S rDNA genes symmetrically placed in both arms of the B chromosome. These results support the hypothesis that the B chromosome of A. paranae is an isochromosome. Additionally, the chromosome pairs Nos. 2 or 23 are considered the possible B chromosome ancestors since both contain syntenic H1 and 18S rRNA sequences. The analysis of DNA sequence fragments of the histone and rRNA genes obtained from the microdissected B chromosomes showed high similarity with those obtained from 0B individuals, which supports the intraspecific origin of B chromosomes in A. paranae. Finally, the population hereby analysed showed a female-biased B chromosome presence suggesting that B chromosomes in this species could influence sex determinism.This research was funded by grants from the State of São Paulo Research Foundation (FAPESP) to DMZAS (2011/16825-3) and CO (2010/17009-2), grants from National Council for Research and Development (CNPq) to FF and by Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Treatment with tocilizumab or corticosteroids for COVID-19 patients with hyperinflammatory state: a multicentre cohort study (SAM-COVID-19)

Objectives: The objective of this study was to estimate the association between tocilizumab or corticosteroids and the risk of intubation or death in patients with coronavirus disease 19 (COVID-19) with a hyperinflammatory state according to clinical and laboratory parameters. Methods: A cohort study was performed in 60 Spanish hospitals including 778 patients with COVID-19 and clinical and laboratory data indicative of a hyperinflammatory state. Treatment was mainly with tocilizumab, an intermediate-high dose of corticosteroids (IHDC), a pulse dose of corticosteroids (PDC), combination therapy, or no treatment. Primary outcome was intubation or death; follow-up was 21 days. Propensity score-adjusted estimations using Cox regression (logistic regression if needed) were calculated. Propensity scores were used as confounders, matching variables and for the inverse probability of treatment weights (IPTWs). Results: In all, 88, 117, 78 and 151 patients treated with tocilizumab, IHDC, PDC, and combination therapy, respectively, were compared with 344 untreated patients. The primary endpoint occurred in 10 (11.4%), 27 (23.1%), 12 (15.4%), 40 (25.6%) and 69 (21.1%), respectively. The IPTW-based hazard ratios (odds ratio for combination therapy) for the primary endpoint were 0.32 (95%CI 0.22-0.47; p < 0.001) for tocilizumab, 0.82 (0.71-1.30; p 0.82) for IHDC, 0.61 (0.43-0.86; p 0.006) for PDC, and 1.17 (0.86-1.58; p 0.30) for combination therapy. Other applications of the propensity score provided similar results, but were not significant for PDC. Tocilizumab was also associated with lower hazard of death alone in IPTW analysis (0.07; 0.02-0.17; p < 0.001). Conclusions: Tocilizumab might be useful in COVID-19 patients with a hyperinflammatory state and should be prioritized for randomized trials in this situatio