Genome sequencing of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is increasingly important to monitor the transmission and adaptive evolution of the virus. The accessibility of high-throughput methods and polymerase chain reaction (PCR) has facilitated a growing ecosystem of protocols. Two differing protocols are tiling multiplex PCR and bait capture enrichment. Each method has advantages and disadvantages but a direct comparison with different viral RNA concentrations has not been performed to assess the performance of these approaches. Here we compare Liverpool amplification, ARTIC amplification, and bait capture using clinical diagnostics samples. All libraries were sequenced using an Illumina MiniSeq with data analyzed using a standardized bioinformatics workflow (SARS-CoV-2 Illumina GeNome Assembly Line; SIGNAL). One sample showed poor SARS-CoV-2 genome coverage and consensus, reflective of low viral RNA concentration. In contrast, the second sample had a higher viral RNA concentration, which yielded good genome coverage and consensus. ARTIC amplification showed the highest depth of coverage results for both samples, suggesting this protocol is effective for low concentrations. Liverpool amplification provided a more even read coverage of the SARS-CoV-2 genome, but at a lower depth of coverage. Bait capture enrichment of SARS-CoV-2 cDNA provided results on par with amplification. While only two clinical samples were examined in this comparative analysis, both the Liverpool and ARTIC amplification methods showed differing efficacy for high and low concentration samples. In addition, amplification-free bait capture enriched sequencing of cDNA is a viable method for generating a SARS-CoV-2 genome sequence and for identification of amplification artifacts.

Sister Marija Krucifiksa Kozulić (1852–1922) was a Croatian nun who is in consideration for beatification by the Vatican, which is facilitated by the identification of her 20th-century remains. Sister Marija was buried in a tomb in Rijeka, Croatia, along with other nuns including her biological sister, Tereza Kozulić (1861–1933). When the remains were exhumed in 2011, they were found in a deteriorated state and commingled with several other sets of remains. Thus, mitochondrial genome sequencing of the long bones was performed to sort the remains by mitochondrial haplotype. Two similar but unique haplotypes belonging to haplogroup H1bu were identified, and samples from these bones were subjected to autosomal short tandem repeat (STR) and single nucleotide polymorphism (SNP) sequencing. Although only partial profiles were obtained, the data were sufficient for kinship analysis with the profile of a paternal niece of Sister Marija (Fides Kozulić). The data indicate that it is 574,195-fold more likely that the two sets of skeletal remains represent 2nd-degree relatives of Fides than sisters who are unrelated to Fides. Although it is impossible to discern which set of remains belongs to Marija and which belongs to Tereza, forensic genomics methods have enabled identification of the sisters.

Orthohantaviruses are re-emerging rodent-borne pathogens distributed all over the world. Here, we report the isolation of a Puumala orthohantavirus (PUUV) strain from bank voles caught in a highly endemic region around the city Osnabrück, north-west Germany. Coding and non-coding sequences of all three segments (S, M, and L) were determined from original lung tissue, after isolation and after additional passaging in VeroE6 cells and a bank vole-derived kidney cell line. Different single amino acid substitutions were observed in the RNA-dependent RNA polymerase (RdRP) of the two stable PUUV isolates. The PUUV strain from VeroE6 cells showed a lower titer when propagated on bank vole cells compared to VeroE6 cells. Additionally, glycoprotein precursor (GPC)-derived virus-like particles of a German PUUV sequence allowed the generation of monoclonal antibodies that allowed the reliable detection of the isolated PUUV strain in the immunofluorescence assay. In conclusion, this is the first isolation of a PUUV strain from Central Europe and the generation of glycoprotein-specific monoclonal antibodies for this PUUV isolate. The obtained virus isolate and GPC-specific antibodies are instrumental tools for future reservoir host studies.

Summary

Background

Direct whole-genome sequencing of Mycobacterium tuberculosis from clinical specimens will be a major breakthrough in tuberculosis diagnosis and control. To date, direct whole-genome sequencing has never been used in genomic epidemiology, and its accuracy in transmission inference remains unknown. We investigated the technical challenges imposed by direct whole-genome sequencing, and used it to infer transmission clusters and predict drug resistance.

Methods

Using an optimised workflow, we did direct whole-genome sequencing for 37 clinical specimens from 23 tuberculosis patients. Nine sputum samples from nine patients who were infected with different non-tuberculous mycobacteria and culture-negative for tuberculosis were used as controls in the qPCR assays and pre-sequencing runs. Additionally, 780 clinical isolates in the region of Comunidad Valenciana (Spain) were whole-genome sequenced between Jan 1, 2014, and Dec 31, 2016. We analysed the genomic variants to build a tuberculosis transmission network for the region, including the clinical specimens, and to predict drug susceptibility profiles.

Findings

After sequencing 37 clinical specimens, 28 specimens (22 [85%] of 26 smear-positive and six [55%] of 11 smear-negative) met the quality criteria for downstream analysis. All 28 clinical specimens clustered with their matching culture isolates, with a median distance of 0 single nucleotide polymorphisms. Of the 28 clinical specimens, 16 (57%) were accurately assigned to ten transmission clusters in the region, and 12 (43%) were unique cases. Transmission inferences and drug-susceptibility predictions from direct whole-genome sequencing data were concordant with sequences from corresponding cultures and phenotypic drug-susceptibility testing. Complete genomic analysis, within a week of specimen receipt, cost €217 per sample (excluding personnel costs).

Interpretation

Direct whole-genome sequencing could be used to accurately delineate transmission clusters of tuberculosis and conduct culture-independent surveillance. Compared with conventional approaches, direct whole-genome sequencing allows researchers to do real-time genomic epidemiology and drug resistance surveillance in settings where culture and drug susceptibility testing are not available.

Funding

European Research Council; Ministerio de Ciencia, Innovación y Universidades (Spanish Government).

The phylogenetic relationships for Ipomoea species are incongruent in previous studies. Comparative karyotype analysis can provide valuable information for phylogenetic relationships among species. A reliable and efficient system for chromosome identification is the foundation for karyotype analysis. However, for the Ipomoea species, individual chromosomes were not identified in all previous reported karyotypes due to their small size, high number and similar morphology. Fluorescence in situ hybridization (FISH) using oligonucleotides (oligos) as probes is a new strategy for chromosome identification and karyotype analysis. Here, we developed the first set of oligo-based probes based on the reference genome of Ipomoea nil, a model species in the genus Ipomoea. In all, we developed four oligo-FISH probes. By a combined use of four oligo probes, sequential FISH analysis were conducted on the same metaphase cells with each round including two probes, which permitted simultaneous identification of all I. nil chromosomes and anchoring 15 pseudochromosomes to individual cytological chromosomes. Moreover, 45S and 5S rDNA were mapped to specific chromosomes. A karyotype based on individually identified chromosomes was established, which was the first FISH-based molecular cytogenetic karyotype of I. nil. Our study has created the basis for studying chromosome variation and evolution in genus Ipomoea by comparative karyotype analysis using this set of oligo probes.

The cave lion is an extinct felid that was widespread across the Holarctic throughout the Late Pleistocene. Its closest extant relative is the lion (Panthera leo), but the timing of the divergence between these two taxa, as well as their taxonomic ranking are contentious. In this study we analyse 31 mitochondrial genome sequences from cave lion individuals that, through a combination of 14C and genetic tip dating, are estimated to be from dates extending well into the mid-Pleistocene. We identified two deeply diverged and well-supported reciprocally monophyletic mitogenome clades in the cave lion, and an additional third distinct lineage represented by a single individual. One of these clades was restricted to Beringia while the other was prevalent across western Eurasia. These observed clade distributions are in line with previous observations that Beringian and European cave lions were morphologically distinct. The divergence dates for these lineages are estimated to be far older than those between extant lions subspecies. By combining our radiocarbon tip-dates with a split time prior that takes into account the most up-to-date fossil stem calibrations, we estimated the mitochondrial DNA divergence between cave lions and lions to be 1.85 Million ya (95% 0.52– 2.91 Mya). Taken together, these results support previous hypotheses that cave lions existed as at least two subspecies during the Pleistocene, and that lions and cave lions were distinct species.

Abstract Relationships of the Neotropical finches in the genera Euphonia and Chlorophonia (Fringillidae: Euphoniinae) have been clarified by recent molecular studies, but species-level relationships within this group have not been thoroughly addressed. In this study, we sampled specimens representing every recognized species of these genera, in addition to 2 outgroup taxa, and used target enrichment to sequence thousands of ultraconserved element (UCE) loci, as well as mitochondrial DNA reconstructed from off-target reads, from each specimen to infer these relationships. We constructed both concatenation and coalescent-based estimates of phylogeny from this dataset using matrices of varying levels of completeness, and we generated a time-scaled ultrametric tree using a recently published fossil-based external calibration. We found uniformly strong support for a monophyletic subfamily Euphoniinae and genus Chlorophonia, but a paraphyletic Euphonia across UCEs and mitochondrial genomes. Otherwise, our inferred relationships were largely concordant with previous studies. Our time-tree indicated a stem divergence time of 13.8 million years ago for this lineage, followed by a relatively young crown age of only 7.1 myr. Reconstructions of biogeographic history based on this tree suggest a South American origin for crown Euphoniinae, possibly resulting from a transoceanic dispersal event from the Eastern Hemisphere, followed by 2 dispersal events into the Caribbean and as many as 6 invasions of North America coinciding with recent estimates of the age at which the Isthmus of Panama had completely formed. We recommend splitting Euphonia and resurrecting the genus Cyanophonia for the 3 blue-hooded species more closely related to Chlorophonia. Based on our results, we suspect that there is undescribed species-level diversity in at least one, possibly many, widespread and phenotypically diverse species.

Vaccination has transformed public health, most notably including the eradication of smallpox. Despite its profound historical importance, little is known of the origins and diversity of the viruses used in smallpox vaccination. Prior to the twentieth century, the method, source and origin of smallpox vaccinations remained unstandardised and opaque. We reconstruct and analyse viral vaccine genomes associated with smallpox vaccination from historical artefacts. Significantly, we recover viral molecules through non-destructive sampling of historical materials lacking signs of biological residues. We use the authenticated ancient genomes to reveal the evolutionary relationships of smallpox vaccination viruses within the poxviruses as a whole.

Polycomb group (PcG) proteins silence gene expression by chemically and physically modifying chromatin. A subset of PcG target loci are compacted and cluster in the nucleus; a conformation that is thought to contribute to gene silencing. However, how these interactions influence gross nuclear organization and their relationship with transcription remains poorly understood. Here we examine the role of Polycomb-repressive complex 1 (PRC1) in shaping 3D genome organization in mouse embryonic stem cells (mESCs). Using a combination of imaging and Hi-C analyses, we show that PRC1-mediated long-range interactions are independent of CTCF and can bridge sites at a megabase scale. Impairment of PRC1 enzymatic activity does not directly disrupt these interactions. We demonstrate that PcG targets coalesce in vivo, and that developmentally induced expression of one of the target loci disrupts this spatial arrangement. Finally, we show that transcriptional activation and the loss of PRC1-mediated interactions are separable events. These findings provide important insights into the function of PRC1, while highlighting the complexity of this regulatory system.