om as , F e r na n d o G on zález Candelas, SeqCOVID-SPAIN consortium, Tanja Stadler & Richard A. NeherThis is a PDF file of a peer-reviewed paper that has been accepted for publication. Although unedited, the content has been subjected to preliminary formatting. Nature is providing this early version of the typeset paper as a service to our authors and readers. The text and figures will undergo copyediting and a proof review before the paper is published in its final form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers apply.

The coronavirus disease 2019 (COVID-19) pandemic has affected the world radically since 2020. Spain was one of the European countries with the highest incidence during the first wave. As a part of a consortium to monitor and study the evolution of the epidemic, we sequenced 2,170 samples, diagnosed mostly before lockdown measures. Here, we identified at least 500 introductions from multiple international sources and documented the early rise of two dominant Spanish epidemic clades (SECs), probably amplified by superspreading events. Both SECs were related closely to the initial Asian variants of SARS-CoV-2 and spread widely across Spain. We inferred a substantial reduction in the effective reproductive number of both SECs due to public-health interventions ( R e < 1), also reflected in the replacement of SECs by a new variant over the summer of 2020. In summary, we reveal a notable difference in the initial genetic makeup of SARS-CoV-2 in Spain compared with other European countries and show evidence to support the effectiveness of lockdown measures in controlling virus spread, even for the most successful genetic variants.

Mycobacterium tuberculosis (Mtb) infection results in a spectrum of clinical and histopathologic manifestations. It has been proposed that the environmental and immune pressures associated with different contexts of infection have different consequences for the associated bacterial populations, affecting drug susceptibility and the emergence of resistance. However, there is little concrete evidence for this model. We prospectively collected sputum samples from 18 newly diagnosed and treatment-naïve patients with tuberculosis and sequenced 795 colony-derived Mtb isolates. Mutant accumulation rates varied considerably between different bacilli isolated from the same individual, and where high rates of mutation were observed, the mutational spectrum was consistent with reactive oxygen species–induced mutagenesis. Elevated bacterial mutation rates were identified in isolates from HIV-negative but not HIV-positive individuals, suggesting that they were immune-driven. These results support the model that mutagenesis of Mtb in vivo is modulated by the host environment, which could drive the emergence of variants associated with drug resistance in a host-dependent manner.

SummaryIt is expected that the next generation of biotech crops displaying enhanced quality traits with benefits to both farmers and consumers will have a better acceptance than first generation biotech crops and will improve public perception of genetic engineering. This will only be true if they are proven to be as safe as traditionally bred crops. In contrast with the first generation of biotech crops where only a single trait is modified, the next generation of biotech crops will add a new level of complexity inherent to the mechanisms underlying their output traits. In this study, a comprehensive evaluation of the comparative safety approach on a qualityimproved biotech crop with metabolic modifications is presented. Three genetically engineered potato lines with silenced polyphenol oxidase (Ppo) transcripts and reduced tuber browning were characterized at both physiological and molecular levels and showed to be equivalent to wild-type (WT) plants when yield-associated traits and photosynthesis were evaluated. Analysis of the primary metabolism revealed several unintended metabolic modifications in the engineered tubers, providing evidence for potential compositional inequivalence between transgenic lines and WT controls. The silencing construct sequence was in silico analysed for potential allergenic cross-reactivity, and no similarities to known allergenic proteins were identified. Moreover, in vivo intake safety evaluation showed no adverse effects in physiological parameters. Taken together, these results provide the first evidence supporting that the safety of next generation biotech crops can be properly assessed following the current evaluation criterion, even if the transgenic and WT crops are not substantially equivalent.

Tuberculosis (TB) surveillance is scarce in most African countries, even though it is the continent with the greatest disease incidence according to the World Health Organization. Liberia is within the 30 countries with the highest TB burden, probably as a consequence of the long civil war and the recent Ebola outbreak, both crippling the health system and depreciating the TB prevention and control programmes. Due to difficulties working in the country, there is a lack of resistance surveys and bacillus characterization. Here, we use genome sequencing of Mycobacterium tuberculosis clinical isolates to fill this gap. Our results highlight that the bacillus population structure is dominated by lineage 4 strains that harbour an outstanding genetic diversity, higher than in the rest of Africa as a whole. Coalescent analyses demonstrate that strains currently circulating in Liberia were introduced several times beginning in the early year 600 CE until very recently coinciding with migratory movements associated with the civil war and Ebola epidemics. A higher multidrug-resistant (MDR)-TB frequency (23.5 %) than current estimates was obtained together with non-catalogued drug-resistance mutations. Additionally, 39 % of strains were in genomic clusters revealing that ongoing transmission is a major contribution to the TB burden in the country. Our report emphasizes the importance of TB surveillance and control in African countries where bacillus diversity, MDR-TB prevalence and transmission are coalescing to jeopardize TB control programmes.

Significance Previous attempts to identify the action of natural selection in the Mycobacterium tuberculosis complex (MTBC) were limited by sample size and averaging across time and lineages. We investigate changes in selective pressures across time for every single gene of the MTBC. We developed a methodology to analyze temporal signals of selection in a large dataset (∼5,000 complete genomes) and showed that 1) almost half of the genes seem to have been under positive selection at some point in time; 2) experimentally confirmed epitopes tend to accumulate more mutations in deeper branches than in external branches; and 3) temporal signals identify genes that were conserved in the past but under positive selection in the present, suggesting ongoing adaptation to the host.

BackgroundModern biology uses experimental systems that involve the exploration of phenotypic variation as a result of the recombination of several genomes. Such systems are useful to investigate the functional evolution of metabolic networks. One such approach is the analysis of transcript and metabolite profiles. These kinds of studies generate a large amount of data, which require dedicated computational tools for their analysis.ResultsThis paper presents a novel software named *omeSOM (transcript/metabol-ome Self Organizing Map) that implements a neural model for biological data clustering and visualization. It allows the discovery of relationships between changes in transcripts and metabolites of crop plants harboring introgressed exotic alleles and furthermore, its use can be extended to other type of omics data. The software is focused on the easy identification of groups including different molecular entities, independently of the number of clusters formed. The *omeSOM software provides easy-to-visualize interfaces for the identification of coordinated variations in the co-expressed genes and co-accumulated metabolites. Additionally, this information is linked to the most widely used gene annotation and metabolic pathway databases.Conclusions*omeSOM is a software designed to give support to the data mining task of metabolic and transcriptional datasets derived from different databases. It provides a user-friendly interface and offers several visualization features, easy to understand by non-expert users. Therefore, *omeSOM provides support for data mining tasks and it is applicable to basic research as well as applied breeding programs. The software and a sample dataset are available free of charge at http://sourcesinc.sourceforge.net/omesom/.

The systematic identity of Senecio madagascariensis is ratified against the opinion that it is conspecific with Senecio inaequidens. Both species are native to South Africa and have been merged in the 'Senecio inaequidens complex', a group of entities difficult to distinguish from each other. Senecio madagascariensis is widespread in South America and Australia, where it is an invasive weed. Mitotic and meiotic studies were conducted on Argentinian material; chromosome counts solved the chromosome number controversy, validating 2n = 20. The karyotype was symmetrical, composed of ten pairs of metacentric chromosomes varying from 1.62 to 2.38 mm in length. The most frequent number of satellited chromosomes was three, but their position was difficult to assign. Meiosis was regular, with a configuration of ten predominantly open bivalents. Univalents and quadrivalents were rarely observed. High frequencies of secondary associations of bivalents, chromosome asynchrony and bivalent grouping were documented, reinforcing the hypothesis that x = 5 is the basic chromosome number. Pollen stainability ranged from 94 to 99%. The relevance of chromosomal studies in the circumscription of S. madagascariensis is discussed. Hybridization and polyploidy, as principal evolutionary forces in this genus, explain the systematic difficulties.