Genomic Epidemiology of Severe Acute Respiratory Syndrome Coronavirus 2, Colombia.

Coronavirus disease (COVID-19) in Colombia was first diagnosed in a traveler arriving from Italy on February 26, 2020. However, limited data are available on the origins and number of introductions of COVID-19 into the country. We sequenced the causative agent of COVID-19, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), from 43 clinical samples we collected, along with another 79 genome sequences available from Colombia. We investigated the emergence and importation routes for SARS-CoV-2 into Colombia by using epidemiologic, historical air travel, and phylogenetic observations. Our study provides evidence of multiple introductions, mostly from Europe, and documents >12 lineages. Phylogenetic findings validate the lineage diversity, support multiple importation events, and demonstrate the evolutionary relationship of epidemiologically linked transmission chains. Our results reconstruct the early evolutionary history of SARS-CoV-2 in Colombia and highlight the advantages of genome sequencing to complement COVID-19 outbreak investigations.

Substitutions in Spike and Nucleocapsid proteins of SARS-CoV-2 circulating in South America.

SARS-CoV-2 is a new member of the genus Betacoronavirus, responsible for the COVID-19 pandemic. The virus crossed the species barrier and established in the human population taking advantage of the spike protein high affinity for the ACE receptor to infect the lower respiratory tract. The Nucleocapsid (N) and Spike (S) are highly immunogenic structural proteins and most commercial COVID-19 diagnostic assays target these proteins. In an unpredictable epidemic, it is essential to know about their genetic variability. The objective of this study was to describe the substitution frequency of the S and N proteins of SARS-CoV-2 in South America. A total of 504 amino acid and nucleotide sequences of the S and N proteins of SARS-CoV-2 from seven South American countries (Argentina, Brazil, Chile, Ecuador, Peru, Uruguay, and Colombia), reported as of June 3, and corresponding to samples collected between March and April 2020, were compared through substitution matrices using the Muscle algorithm. Forty-three sequences from 13 Colombian departments were obtained in this study using the Oxford Nanopore and Illumina MiSeq technologies, following the amplicon-based ARTIC network protocol. The substitutions D614G in S and R203K/G204R in N were the most frequent in South America, observed in 83% and 34% of the sequences respectively. Strikingly, genomes with the conserved position D614 were almost completely replaced by genomes with the G614 substitution between March to April 2020. A similar replacement pattern was observed with R203K/G204R although more marked in Chile, Argentina and Brazil, suggesting similar introduction history and/or control strategies of SARS-CoV-2 in these countries. It is necessary to continue with the genomic surveillance of S and N proteins during the SARS-CoV-2 pandemic as this information can be useful for developing vaccines, therapeutics and diagnostic tests.

Molecular analysis of several in-house rRT-PCR protocols for SARS-CoV-2 detection in the context of genetic variability of the virus in Colombia.

The COVID-19 pandemic caused by SARS-CoV-2 is a public health problem unprecedented in the recent history of humanity. Different in-house real-time RT-PCR (rRT-PCR) methods for SARS-CoV-2 diagnosis and the appearance of genomes with mutations in primer regions have been reported. Hence, whole-genome data from locally-circulating SARS-CoV-2 strains contribute to the knowledge of its global variability and the development and fine tuning of diagnostic protocols. To describe the genetic variability of Colombian SARS-CoV-2 genomes in hybridization regions of oligonucleotides of the main in-house methods for SARS-CoV-2 detection, RNA samples with confirmed SARS-CoV-2 molecular diagnosis were processed through next-generation sequencing. Primers/probes sequences from 13 target regions for SARS-CoV-2 detection suggested by 7 institutions and consolidated by WHO during the early stage of the pandemic were aligned with Muscle tool to assess the genetic variability potentially affecting their performance. Finally, the corresponding codon positions at the 3' end of each primer, the open reading frame inspection was identified for each gene/protein product. Complete SARS-CoV-2 genomes were obtained from 30 COVID-19 cases, representative of the current epidemiology in the country. Mismatches between at least one Colombian sequence and five oligonucleotides targeting the RdRP and N genes were observed. The 3' end of 4 primers aligned to the third codon position, showed high risk of nucleotide substitution and potential mismatches at this critical position. Genetic variability was detected in Colombian SARS-CoV-2 sequences in some of the primer/probe regions for in-house rRT-PCR diagnostic tests available at WHO COVID-19 technical guidelines; its impact on the performance and rates of false-negative results should be experimentally evaluated. The genomic surveillance of SARS-CoV-2 is highly recommended for the early identification of mutations in critical regions and to issue recommendations on specific diagnostic tests to ensure the coverage of locally-circulating genetic variants.

Rapid mitochondrial genome sequencing based on Oxford Nanopore Sequencing and a proxy for vertebrate species identification.

Molecular information is crucial for species identification when facing challenging morphology-based specimen identifications. The use of DNA barcodes partially solves this problem, but in some cases when PCR is not an option (i.e., primers are not available, problems in reaction standardization), amplification-free approaches could be an optimal alternative. Recent advances in DNA sequencing, like the MinION device from Oxford Nanopore Technologies (ONT), allow to obtain genomic data with low laboratory and technical requirements, and at a relatively low cost. In this study, we explore ONT sequencing for molecular species identification from a total DNA sample obtained from a neotropical rodent and we also test the technology for complete mitochondrial genome reconstruction via genome skimming. We were able to obtain "de novo" the complete mitogenome of a specimen from the genus Melanomys (Cricetidae: Sigmodontinae) with average depth coverage of 78X using ONT-only data and by combining multiple assembly routines. Our pipeline for an automated species identification was able to identify the sample using unassembled sequence data (raw) in a reasonable computing time, which was substantially reduced when a priori information related to the organism identity was known. Our findings suggest ONT sequencing as a suitable candidate to solve species identification problems in metazoan nonmodel organisms and generate complete mtDNA datasets.

Bacillus subtilis EA-CB0575 genome reveals clues for plant growth promotion and potential for sustainable agriculture.

Bacillus subtilis is a remarkably diverse bacterial species that displays many ecological functions. Given its genomic diversity, the strain Bacillus subtilis EA-CB0575, isolated from the rhizosphere of a banana plant, was sequenced and assembled to determine the genomic potential associated with its plant growth promotion potential. The genome was sequenced by Illumina technology and assembled using Velvet 1.2.10, resulting in a whole genome of 4.09 Mb with 4332 genes. Genes involved in the production of indoles, siderophores, lipopeptides, volatile compounds, phytase, bacilibactin, and nitrogenase were predicted by gene annotation or by metabolic pathway prediction by RAST. These potential traits were determined using in vitro biochemical tests, finding that B. subtilis EA-CB0575 produces two families of lipopeptides (surfactin and fengycin), solubilizes phosphate, fixes nitrogen, and produces indole and siderophores compounds. Finally, strain EA-CB0575 increased 34.60% the total dry weight (TDW) of tomato plants with respect to non-inoculated plants at greenhouse level. These results suggest that the identification of strain-specific genes and predicted metabolic pathways might explain the strain potential to promote plant growth by several mechanisms of action, accelerating the development of plant biostimulants for sustainable agricultural.

Complete mitogenome of the biocontroller fungus Purpureocillium sp. (Ascomycota, Ophiocordycipitaceae, Hypocreales).

The strain Purpureocillium sp. UdeA0106 is an antagonist of nematodes, fungi, and garden symphylans from crops with high economic importance in Colombia (Salazar 2013; Salazar et al. 2014; Cardona et al. 2014; Gallego et al. 2014) and is being studied to be proposed as new species. It was included on the 1000 fungal genomes project to elucidate its phylogenetic relationships with other fungi. Purpureocillium's mitogenome has 23,495 bp of circular size. It contains 15 protein-coding genes without duplications (PCGs), corresponding to the 60% of its total length, 23 transfer genes (7.6% tRNA), two of them duplicated (trnR and trnM), and two ribosomal genes (17.6% rRNA) and a GC content of 28.44%. A phylogenetic tree was proposed using their 14 PCGs mitochondrial genes and was compared with other fungi of the Subphylum Pezizomycotina. Phylogenetics relationships showed UdeA0106 to be close to P. chlamydosporia and M. anisopliae forming a cluster with other fungal biocontrol agents and separated the strain of plant pathogenic fungi.

Mitochondrial genome characterization of Tecia solanivora (Lepidoptera: Gelechiidae) and its phylogenetic relationship with other lepidopteran insects.

The complete mitogenome of the potato tuber moth Tecia solanivora (Lepidoptera: Gelechiidae) was sequenced, annotated, characterized and compared with 140 species of the order Lepidoptera. The circular genome is 15,251 bp, containing 37 genes (13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes and an A+T-rich region). The gene arrangement was identical to other lepidopteran mitogenomes but different from the ancestral arrangement found in most insects for the tRNA-Met gene (A+T-region, tRNA-I, tRNA-Q, tRNA-M). The mitogenome of T. solanivora is highly A+T-biased (78.2%) and exhibits negative AT- and GC-skews. All PCGs are initiated by canonical ATN start codons, except for Cytochrome Oxidase subunit 1 (COI), which is initiated by CGA. Most PCGs have a complete typical stop codon (TAA). Only NAD1 has a TAG stop codon and the COII and NAD5 genes have an incomplete stop codon consisting of just a T. The A+T-rich region is 332 bp long and contains common features found in lepidopteran mitogenomes, including the 'ATAGA' motif, a 17 bp poly (T) stretch and a (AT)8 element preceded by the 'ATTTA' motif. Other tandem repeats like (TAA)4 and (TAT)7 were found, as well as (T)6 and (A)10 mononucleotide repeat elements. Finally, this mitogenome has 20 intergenic spacer regions. The phylogenetic relationship of T. solanivora with 28 other lepidopteran families (12 superfamilies) showed that taxonomic classification by morphological features coincides with the inferred phylogeny. Thus, the Gelechiidae family represents a monophyletic group, suggesting that T. solanivora and Pectinophora gossypiella have a recent common ancestor.