The Bioinformatic Applications of Hi-C and Linked Reads.

Long-range sequencing grants insight into additional genetic information beyond that which can be accessed by both short reads and modern long-read technology. Several new sequencing technologies are available for long-range datasets such as "Hi-C" and "Linked Reads" with high-throughput and high-resolution genome analysis, and are rapidly advancing the field of genome assembly, genome scaffolding, and more comprehensive variant identification. In this article, we focused on five major long-range sequencing technologies: high-throughput chromosome conformation capture (Hi-C), 10x Genomics Linked Reads, haplotagging, transposase enzyme linked long-read sequencing (TELL-seq), and single tube long fragment read (stLFR). We detailed the mechanisms and data products of the five platforms, introduced several of the most important applications, evaluated the quality of sequencing data from different platforms, and discussed the currently available bioinformatics tools. We hope this work will benefit the selection of appropriate long-range technology for specific biological studies.

Identifying the best PCR enzyme for library amplification in NGS.

Background. PCR amplification is a necessary step in many next-generation sequencing (NGS) library preparation methods [1, 2]. Whilst many PCR enzymes are developed to amplify single targets efficiently, accurately and with specificity, few are developed to meet the challenges imposed by NGS PCR, namely unbiased amplification of a wide range of different sizes and GC content. As a result PCR amplification during NGS library prep often results in bias toward GC neutral and smaller fragments. As NGS has matured, optimized NGS library prep kits and polymerase formulations have emerged and in this study we have tested a wide selection of available enzymes for both short-read Illumina library preparation and long fragment amplification ahead of long-read sequencing.We tested over 20 different hi-fidelity PCR enzymes/NGS amplification mixes on a range of Illumina library templates of varying GC content and composition, and find that both yield and genome coverage uniformity characteristics of the commercially available enzymes varied dramatically. Three enzymes Quantabio RepliQa Hifi Toughmix, Watchmaker Library Amplification Hot Start Master Mix (2X) 'Equinox' and Takara Ex Premier were found to give a consistent performance, over all genomes, that mirrored closely that observed for PCR-free datasets. We also test a range of enzymes for long-read sequencing by amplifying size fractionated S. cerevisiae DNA of average size 21.6 and 13.4 kb, respectively.The enzymes of choice for short-read (Illumina) library fragment amplification are Quantabio RepliQa Hifi Toughmix, Watchmaker Library Amplification Hot Start Master Mix (2X) 'Equinox' and Takara Ex Premier, with RepliQa also being the best performing enzyme from the enzymes tested for long fragment amplification prior to long-read sequencing.

Haplotype-specific assembly of shattered chromosomes in esophageal adenocarcinomas.

The epigenetic landscape of cancer is regulated by many factors, but primarily it derives from the underlying genome sequence. Chromothripsis is a catastrophic localized genome shattering event that drives, and often initiates, cancer evolution. We characterized five esophageal adenocarcinoma organoids with chromothripsis using long-read sequencing and transcriptome and epigenome profiling. Complex structural variation and subclonal variants meant that haplotype-aware de novo methods were required to generate contiguous cancer genome assemblies. Chromosomes were assembled separately and scaffolded using haplotype-resolved Hi-C reads, producing accurate assemblies even with up to 900 structural rearrangements. There were widespread differences between the chromothriptic and wild-type copies of chromosomes in topologically associated domains, chromatin accessibility, histone modifications, and gene expression. Differential epigenome peaks were most enriched within 10 kb of chromothriptic structural variants. Alterations in transcriptome and higher-order chromosome organization frequently occurred near differential epigenetic marks. Overall, chromothripsis reshapes gene regulation, causing coordinated changes in epigenetic landscape, transcription, and chromosome conformation.

Somatic evolution of marine transmissible leukemias in the common cockle, Cerastoderma edule.

Transmissible cancers are malignant cell lineages that spread clonally between individuals. Several such cancers, termed bivalve transmissible neoplasia (BTN), induce leukemia-like disease in marine bivalves. This is the case of BTN lineages affecting the common cockle, Cerastoderma edule, which inhabits the Atlantic coasts of Europe and northwest Africa. To investigate the evolution of cockle BTN, we collected 6,854 cockles, diagnosed 390 BTN tumors, generated a reference genome and assessed genomic variation across 61 tumors. Our analyses confirmed the existence of two BTN lineages with hemocytic origins. Mitochondrial variation revealed mitochondrial capture and host co-infection events. Mutational analyses identified lineage-specific signatures, one of which likely reflects DNA alkylation. Cytogenetic and copy number analyses uncovered pervasive genomic instability, with whole-genome duplication, oncogene amplification and alkylation-repair suppression as likely drivers. Satellite DNA distributions suggested ancient clonal origins. Our study illuminates long-term cancer evolution under the sea and reveals tolerance of extreme instability in neoplastic genomes.

The evolution of two transmissible cancers in Tasmanian devils.

Tasmanian devils have spawned two transmissible cancer lineages, named devil facial tumor 1 (DFT1) and devil facial tumor 2 (DFT2). We investigated the genetic diversity and evolution of these clones by analyzing 78 DFT1 and 41 DFT2 genomes relative to a newly assembled, chromosome-level reference. Time-resolved phylogenetic trees reveal that DFT1 first emerged in 1986 (1982 to 1989) and DFT2 in 2011 (2009 to 2012). Subclone analysis documents transmission of heterogeneous cell populations. DFT2 has faster mutation rates than DFT1 across all variant classes, including substitutions, indels, rearrangements, transposable element insertions, and copy number alterations, and we identify a hypermutated DFT1 lineage with defective DNA mismatch repair. Several loci show plausible evidence of positive selection in DFT1 or DFT2, including loss of chromosome Y and inactivation of MGA, but none are common to both cancers. This study reveals the parallel long-term evolution of two transmissible cancers inhabiting a common niche in Tasmanian devils.

The genome sequence of the Norway rat, Rattus norvegicus Berkenhout 1769.

We present a genome assembly from an individual male Rattus norvegicus (the Norway rat; Chordata; Mammalia; Rodentia; Muridae). The genome sequence is 2.44 gigabases in span. The majority of the assembly is scaffolded into 20 chromosomal pseudomolecules, with both X and Y sex chromosomes assembled. This genome assembly, mRatBN7.2, represents the new reference genome for R. norvegicus and has been adopted by the Genome Reference Consortium.

The complete genome sequence of Eimeria tenella (Tyzzer 1929), a common gut parasite of chickens.

We present a genome assembly from a clonal population of Eimeria tenella Houghton parasites (Apicomplexa; Conoidasida; Eucoccidiorida; Eimeriidae). The genome sequence is 53.25 megabases in span. The entire assembly is scaffolded into 15 chromosomal pseudomolecules, with complete mitochondrion and apicoplast organellar genomes also present.

The genome sequence of the European golden eagle, Aquila chrysaetos chrysaetos Linnaeus 1758.

We present a genome assembly from an individual female Aquila chrysaetos chrysaetos (the European golden eagle; Chordata; Aves; Accipitridae). The genome sequence is 1.23 gigabases in span. The majority of the assembly is scaffolded into 28 chromosomal pseudomolecules, including the W and Z sex chromosomes.

Patterns of within-host genetic diversity in SARS-CoV-2.

Monitoring the spread of SARS-CoV-2 and reconstructing transmission chains has become a major public health focus for many governments around the world. The modest mutation rate and rapid transmission of SARS-CoV-2 prevents the reconstruction of transmission chains from consensus genome sequences, but within-host genetic diversity could theoretically help identify close contacts. Here we describe the patterns of within-host diversity in 1181 SARS-CoV-2 samples sequenced to high depth in duplicate. 95.1% of samples show within-host mutations at detectable allele frequencies. Analyses of the mutational spectra revealed strong strand asymmetries suggestive of damage or RNA editing of the plus strand, rather than replication errors, dominating the accumulation of mutations during the SARS-CoV-2 pandemic. Within- and between-host diversity show strong purifying selection, particularly against nonsense mutations. Recurrent within-host mutations, many of which coincide with known phylogenetic homoplasies, display a spectrum and patterns of purifying selection more suggestive of mutational hotspots than recombination or convergent evolution. While allele frequencies suggest that most samples result from infection by a single lineage, we identify multiple putative examples of co-infection. Integrating these results into an epidemiological inference framework, we find that while sharing of within-host variants between samples could help the reconstruction of transmission chains, mutational hotspots and rare cases of superinfection can confound these analyses.

The COVID-19 pandemic has had major health impacts across the globe. The scientific community has focused much attention on finding ways to monitor how the virus responsible for the pandemic, SARS-CoV-2, spreads. One option is to perform genetic tests, known as sequencing, on SARS-CoV-2 samples to determine the genetic code of the virus and to find any differences or mutations in the genes between the viral samples. Viruses mutate within their hosts and can develop into variants that are able to more easily transmit between hosts. Genetic sequencing can reveal how genetically similar two SARS-CoV-2 samples are. But tracking how SARS-CoV-2 moves from one person to the next through sequencing can be tricky. Even a sample of SARS-CoV-2 viruses from the same individual can display differences in their genetic material or within-host variants. Could genetic testing of within-host variants shed light on factors driving SARS-CoV-2 to evolve in humans? To get to the bottom of this, Tonkin-Hill, Martincorena et al. probed the genetics of SARS-CoV-2 within-host variants using 1,181 samples. The analyses revealed that 95.1% of samples contained within-host variants. A number of variants occurred frequently in many samples, which were consistent with mutational hotspots in the SARS-CoV-2 genome. In addition, within-host variants displayed mutation patterns that were similar to patterns found between infected individuals. The shared within-host variants between samples can help to reconstruct transmission chains. However, the observed mutational hotspots and the detection of multiple strains within an individual can make this challenging. These findings could be used to help predict how SARS-CoV-2 evolves in response to interventions such as vaccines. They also suggest that caution is needed when using information on within-host variants to determine transmission between individuals.

Nanopore Sequencing Indicates That Tandem Amplification of Chromosome 20q11.21 in Human Pluripotent Stem Cells Is Driven by Break-Induced Replication.

Copy number variants (CNVs) are genomic rearrangements implicated in numerous congenital and acquired diseases, including cancer. The appearance of culture-acquired CNVs in human pluripotent stem cells (PSCs) has prompted concerns for their use in regenerative medicine. A particular problem in PSC is the frequent occurrence of CNVs in the q11.21 region of chromosome 20. However, the exact mechanism of origin of this amplicon remains elusive due to the difficulty in delineating its sequence and breakpoints. Here, we have addressed this problem using long-read Nanopore sequencing of two examples of this CNV, present as duplication and as triplication. In both cases, the CNVs were arranged in a head-to-tail orientation, with microhomology sequences flanking or overlapping the proximal and distal breakpoints. These breakpoint signatures point to a mechanism of microhomology-mediated break-induced replication in CNV formation, with surrounding Alu sequences likely contributing to the instability of this genomic region.

The genome sequence of the European water vole, Arvicola amphibius Linnaeus 1758.

We present a genome assembly from an individual male Arvicola amphibius (the European water vole; Chordata; Mammalia; Rodentia; Cricetidae). The genome sequence is 2.30 gigabases in span. The majority of the assembly is scaffolded into 18 chromosomal pseudomolecules, including the X sex chromosome. Gene annotation of this assembly on Ensembl has identified 21,394 protein coding genes.

A preoperative estimate of central venous pressure is associated with early Fontan failure.

OBJECTIVE: Early Fontan failure is a serious complication after total cavopulmonary connection, characterized by high central venous pressure, low cardiac output, and resistance to medical therapy. This study aimed to estimate postoperative central venous pressure in patients with total cavopulmonary connection using data routinely collected during preoperative assessment. We sought to determine if this metric correlated with measured postoperative central venous pressure and if it was associated with early Fontan failure. METHODS: In this retrospective study, central venous pressure in total cavopulmonary connection was estimated in 131 patients undergoing pre-total cavopulmonary connection assessment by cardiac magnetic resonance imaging and central venous pressure measurement under general anesthesia. Postoperative central venous pressure during the first 24 hours in the intensive care unit was collected from electronic patient records in a subset of patients. Early Fontan failure was defined as death, transplantation, total cavopulmonary connection takedown, or emergency fenestration within the first 30 days. RESULTS: Estimated central venous pressure in total cavopulmonary connection correlated significantly with central venous pressure during the first 24 hours in the intensive care unit (r = 0.26, P = .03), particularly in patients without a fenestration (r = 0.45, P = .01). Central venous pressure in total cavopulmonary connection was significantly associated with early Fontan failure (odds ratio, 1.1; 95% confidence interval, 1.01-1.21; P = .03). A threshold of central venous pressure in total cavopulmonary connection 33 mm Hg or greater was found to have the highest specificity (90%) and sensitivity (58%) for identifying early Fontan failure (area under receiver operating curve = 0.73; odds ratio, 12.4; 95% confidence interval, 2.5-62.3; P = .002). This association was stronger in patients with single superior vena cava. CONCLUSIONS: Estimated central venous pressure in total cavopulmonary connection is an easily calculated metric combining preoperative pressure and flow data. Higher central venous pressure in total cavopulmonary connection is associated with an increased risk of early Fontan failure and is correlated with directly measured post-total cavopulmonary connection pressure. Identification of patients at risk of early Fontan failure has the potential to guide risk-mitigation strategies.

The genome sequence of the European turtle dove, Streptopelia turtur Linnaeus 1758.

We present a genome assembly from an individual female Streptopelia turtur (the European turtle dove; Chordata; Aves; Columbidae). The genome sequence is 1.18 gigabases in span. The majority of the assembly is scaffolded into 35 chromosomal pseudomolecules, with the W and Z sex chromosomes assembled.

The genome sequence of the eastern grey squirrel, Sciurus carolinensis Gmelin, 1788.

We present a genome assembly from an individual male Sciurus carolinensis (the eastern grey squirrel; Vertebrata; Mammalia; Eutheria; Rodentia; Sciuridae). The genome sequence is 2.82 gigabases in span. The majority of the assembly (92.3%) is scaffolded into 21 chromosomal-level scaffolds, with both X and Y sex chromosomes assembled.

Genomic and transcriptomic variation defines the chromosome-scale assembly of Haemonchus contortus, a model gastrointestinal worm.

Haemonchus contortus is a globally distributed and economically important gastrointestinal pathogen of small ruminants and has become a key nematode model for studying anthelmintic resistance and other parasite-specific traits among a wider group of parasites including major human pathogens. Here, we report using PacBio long-read and OpGen and 10X Genomics long-molecule methods to generate a highly contiguous 283.4 Mbp chromosome-scale genome assembly including a resolved sex chromosome for the MHco3(ISE).N1 isolate. We show a remarkable pattern of conservation of chromosome content with Caenorhabditis elegans, but almost no conservation of gene order. Short and long-read transcriptome sequencing allowed us to define coordinated transcriptional regulation throughout the parasite's life cycle and refine our understanding of cis- and trans-splicing. Finally, we provide a comprehensive picture of chromosome-wide genetic diversity both within a single isolate and globally. These data provide a high-quality comparison for understanding the evolution and genomics of Caenorhabditis and other nematodes and extend the experimental tractability of this model parasitic nematode in understanding helminth biology, drug discovery and vaccine development, as well as important adaptive traits such as drug resistance.

The mutational signature profile of known and suspected human carcinogens in mice.

Epidemiological studies have identified many environmental agents that appear to significantly increase cancer risk in human populations. By analyzing tumor genomes from mice chronically exposed to 1 of 20 known or suspected human carcinogens, we reveal that most agents do not generate distinct mutational signatures or increase mutation burden, with most mutations, including driver mutations, resulting from tissue-specific endogenous processes. We identify signatures resulting from exposure to cobalt and vinylidene chloride and link distinct human signatures (SBS19 and SBS42) with 1,2,3-trichloropropane, a haloalkane and pollutant of drinking water, and find these and other signatures in human tumor genomes. We define the cross-species genomic landscape of tumors induced by an important compendium of agents with relevance to human health.

Serial Assessment of Coronary Flow Reserve by Rubidium-82 Positron Emission Tomography Predicts Mortality in Heart Transplant Recipients.

OBJECTIVES: This study aimed to evaluate the long-term prognostic value of serial assessment of coronary flow reserve (CFR) by rubidium Rb 82 (82Rb) positron emission tomography (PET) in heart transplantation (HT) patients. BACKGROUND: Cardiac allograft vasculopathy is a major determinant of late mortality in HT recipients. The long-term prognostic value of serial CFR quantification by PET imaging in HT patients is unknown. METHODS: A total of 89 patients with history of HT (71% men, 7.0 ± 5.7 years post-HT, age 57 ± 11 years) scheduled for dynamic rest and stress (dipyridamole) 82Rb PET between March 1, 2008 and July 31, 2009 (PET-1) were prospectively enrolled in a single-center study. PET myocardial perfusion studies were reprocessed using U.S. Food and Drug Administration-approved software (Corridor 4DM, version 2017) for calculation of CFR. Follow-up PET (PET-2) imaging was performed in 69 patients at 1.9 ± 0.3 years following PET-1. Patients were categorized based on CFR values considering CFR ≤1.5 as low and CFR >1.5 as high CFR. RESULTS: Forty deaths occurred during the median follow-up time of 8.6 years. Low CFR at PET-1 was associated with a 2.77-fold increase in all-cause mortality (95% confidence interval [CI]: 1.34 to 5.74; p = 0.004). CFR decreased over time in patients with follow-up imaging (PET-1: 2.11 ± 0.74 vs. PET-2: 1.81 ± 0.61; p = 0.003). Twenty-five patients were reclassified based on PET-1 and PET-2 (high to low CFR: n = 18, low to high CFR: n = 7). Overall survival was similar in patients reclassified from high to low as patients with low to low CFR, whereas patients reclassified from low to high had similar survival as patients with high to high CFR. In multivariate Cox regression of patients with PET-2, higher baseline CFR (hazard ratio [HR] for a 0.73 unit (one SD) increase: 0.36, 95% CI: 0.16 to 0.82) and reduction in CFR from PET-1 to PET-2 (HR for a 0.79 unit (one SD) decrease: 1.50 to 7.84) were independent predictors of all-cause mortality. CONCLUSIONS: Serial assessment of CFR by 82Rb PET independently predicts long-term mortality in HT patients.

Best Practices for Illumina Library Preparation.

In this unit, we describe a set of protocols and recommendations for Illumina library preparation. We review best practices in template quantitation methods; template fragmentation methodologies; solid-phase reverse-immobilization cleanup, including buffer exchange and size selection; end repair, A-tailing, and adapter ligation; indexing strategies; considerations regarding whether to use polymerase chain reaction; final library quantification methodologies; and normalization and pooling strategies. These workflows are applicable to both high-throughput and low-throughput Illumina library preparation and should help reduce bias, increase cost effectiveness, and produce high library yields. This is an extensive update of the previous version of this unit. © 2019 by John Wiley & Sons, Inc.

The cardiovascular phenotype of childhood hypertension: a cardiac magnetic resonance study.

BACKGROUND: The cardiovascular phenotype is poorly characterized in treated pediatric hypertension. Cardiovascular magnetic resonance imaging (MRI) can be used to better characterize both cardiac and vascular phenotype in children with hypertension. OBJECTIVE: To use MRI to determine the cardiac and vascular phenotypes of different forms of treated hypertension and compare the results with those of healthy children. MATERIALS AND METHODS: Sixty children (15 with chronic renal disease with hypertension, 15 with renovascular hypertension, 15 with essential hypertension and 15 healthy subjects) underwent MRI with noninvasive blood pressure measurements. Cardiovascular parameters measured include systemic vascular resistance, total arterial compliance, left ventricular mass and volumetric data, ejection fraction and myocardial velocity. Between-group comparisons were used to investigate differences in the hypertension types. RESULTS: Renal hypertension was associated with elevated vascular resistance (P≤0.007) and normal arterial compliance. Conversely, children with essential hypertension had normal resistance but increased compliance (P=0.001). Renovascular hypertension was associated with both increased resistance and compliance (P≤0.03). There was no difference in ventricular volumes, mass or cardiac output between groups. Children with renal hypertension also had lower systolic and diastolic myocardial velocities. CONCLUSION: Cardiovascular MRI may identify distinct vascular and cardiac phenotypes in different forms of treated childhood hypertension. Future studies are needed to investigate how this may inform further optimisation of blood pressure treatment in different types of hypertension.