Massive genomic rearrangement acquired in a single catastrophic event during cancer development.

Cancer is driven by somatically acquired point mutations and chromosomal rearrangements, conventionally thought to accumulate gradually over time. Using next-generation sequencing, we characterize a phenomenon, which we term chromothripsis, whereby tens to hundreds of genomic rearrangements occur in a one-off cellular crisis. Rearrangements involving one or a few chromosomes crisscross back and forth across involved regions, generating frequent oscillations between two copy number states. These genomic hallmarks are highly improbable if rearrangements accumulate over time and instead imply that nearly all occur during a single cellular catastrophe. The stamp of chromothripsis can be seen in at least 2%-3% of all cancers, across many subtypes, and is present in ∼25% of bone cancers. We find that one, or indeed more than one, cancer-causing lesion can emerge out of the genomic crisis. This phenomenon has important implications for the origins of genomic remodeling and temporal emergence of cancer.

Image2Flow: A proof-of-concept hybrid image and graph convolutional neural network for rapid patient-specific pulmonary artery segmentation and CFD flow field calculation from 3D cardiac MRI data.

Computational fluid dynamics (CFD) can be used for non-invasive evaluation of hemodynamics. However, its routine use is limited by labor-intensive manual segmentation, CFD mesh creation, and time-consuming simulation. This study aims to train a deep learning model to both generate patient-specific volume-meshes of the pulmonary artery from 3D cardiac MRI data and directly estimate CFD flow fields. This proof-of-concept study used 135 3D cardiac MRIs from both a public and private dataset. The pulmonary arteries in the MRIs were manually segmented and converted into volume-meshes. CFD simulations were performed on ground truth meshes and interpolated onto point-point correspondent meshes to create the ground truth dataset. The dataset was split 110/10/15 for training, validation, and testing. Image2Flow, a hybrid image and graph convolutional neural network, was trained to transform a pulmonary artery template to patient-specific anatomy and CFD values, taking a specific inlet velocity as an additional input. Image2Flow was evaluated in terms of segmentation, and the accuracy of predicted CFD was assessed using node-wise comparisons. In addition, the ability of Image2Flow to respond to increasing inlet velocities was also evaluated. Image2Flow achieved excellent segmentation accuracy with a median Dice score of 0.91 (IQR: 0.86-0.92). The median node-wise normalized absolute error for pressure and velocity magnitude was 11.75% (IQR: 9.60-15.30%) and 9.90% (IQR: 8.47-11.90), respectively. Image2Flow also showed an expected response to increased inlet velocities with increasing pressure and velocity values. This proof-of-concept study has shown that it is possible to simultaneously perform patient-specific volume-mesh based segmentation and pressure and flow field estimation using Image2Flow. Image2Flow completes segmentation and CFD in ~330ms, which is ~5000 times faster than manual methods, making it more feasible in a clinical environment.

Shape-driven deep neural networks for fast acquisition of aortic 3D pressure and velocity flow fields.

Computational fluid dynamics (CFD) can be used to simulate vascular haemodynamics and analyse potential treatment options. CFD has shown to be beneficial in improving patient outcomes. However, the implementation of CFD for routine clinical use is yet to be realised. Barriers for CFD include high computational resources, specialist experience needed for designing simulation set-ups, and long processing times. The aim of this study was to explore the use of machine learning (ML) to replicate conventional aortic CFD with automatic and fast regression models. Data used to train/test the model consisted of 3,000 CFD simulations performed on synthetically generated 3D aortic shapes. These subjects were generated from a statistical shape model (SSM) built on real patient-specific aortas (N = 67). Inference performed on 200 test shapes resulted in average errors of 6.01% ±3.12 SD and 3.99% ±0.93 SD for pressure and velocity, respectively. Our ML-based models performed CFD in ∼0.075 seconds (4,000x faster than the solver). This proof-of-concept study shows that results from conventional vascular CFD can be reproduced using ML at a much faster rate, in an automatic process, and with reasonable accuracy.

The industrial melanism mutation in British peppered moths is a transposable element.

Discovering the mutational events that fuel adaptation to environmental change remains an important challenge for evolutionary biology. The classroom example of a visible evolutionary response is industrial melanism in the peppered moth (Biston betularia): the replacement, during the Industrial Revolution, of the common pale typica form by a previously unknown black (carbonaria) form, driven by the interaction between bird predation and coal pollution. The carbonaria locus has been coarsely localized to a 200-kilobase region, but the specific identity and nature of the sequence difference controlling the carbonaria-typica polymorphism, and the gene it influences, are unknown. Here we show that the mutation event giving rise to industrial melanism in Britain was the insertion of a large, tandemly repeated, transposable element into the first intron of the gene cortex. Statistical inference based on the distribution of recombined carbonaria haplotypes indicates that this transposition event occurred around 1819, consistent with the historical record. We have begun to dissect the mode of action of the carbonaria transposable element by showing that it increases the abundance of a cortex transcript, the protein product of which plays an important role in cell-cycle regulation, during early wing disc development. Our findings fill a substantial knowledge gap in the iconic example of microevolutionary change, adding a further layer of insight into the mechanism of adaptation in response to natural selection. The discovery that the mutation itself is a transposable element will stimulate further debate about the importance of 'jumping genes' as a source of major phenotypic novelty.

Repeat associated mechanisms of genome evolution and function revealed by the Mus caroli and Mus pahari genomes.

Understanding the mechanisms driving lineage-specific evolution in both primates and rodents has been hindered by the lack of sister clades with a similar phylogenetic structure having high-quality genome assemblies. Here, we have created chromosome-level assemblies of the Mus caroli and Mus pahari genomes. Together with the Mus musculus and Rattus norvegicus genomes, this set of rodent genomes is similar in divergence times to the Hominidae (human-chimpanzee-gorilla-orangutan). By comparing the evolutionary dynamics between the Muridae and Hominidae, we identified punctate events of chromosome reshuffling that shaped the ancestral karyotype of Mus musculus and Mus caroli between 3 and 6 million yr ago, but that are absent in the Hominidae. Hominidae show between four- and sevenfold lower rates of nucleotide change and feature turnover in both neutral and functional sequences, suggesting an underlying coherence to the Muridae acceleration. Our system of matched, high-quality genome assemblies revealed how specific classes of repeats can play lineage-specific roles in related species. Recent LINE activity has remodeled protein-coding loci to a greater extent across the Muridae than the Hominidae, with functional consequences at the species level such as reproductive isolation. Furthermore, we charted a Muridae-specific retrotransposon expansion at unprecedented resolution, revealing how a single nucleotide mutation transformed a specific SINE element into an active CTCF binding site carrier specifically in Mus caroli, which resulted in thousands of novel, species-specific CTCF binding sites. Our results show that the comparison of matched phylogenetic sets of genomes will be an increasingly powerful strategy for understanding mammalian biology.

Reducing Contrast Agent Dose in Cardiovascular MR Angiography with Deep Learning.

BACKGROUND: Contrast-enhanced magnetic resonance angiography (MRA) is used to assess various cardiovascular conditions. However, gadolinium-based contrast agents (GBCAs) carry a risk of dose-related adverse effects. PURPOSE: To develop a deep learning method to reduce GBCA dose by 80%. STUDY TYPE: Retrospective and prospective. POPULATION: A total of 1157 retrospective and 40 prospective congenital heart disease patients for training/validation and testing, respectively. FIELD STRENGTH/SEQUENCE: A 1.5 T, T1-weighted three-dimensional (3D) gradient echo. ASSESSMENT: A neural network was trained to enhance low-dose (LD) 3D MRA using retrospective synthetic data and tested with prospective LD data. Image quality for LD (LD-MRA), enhanced LD (ELD-MRA), and high-dose (HD-MRA) was assessed in terms of signal-to-noise ratio (SNR), contrast-to-noise ratio (CNR), and a quantitative measure of edge sharpness and scored for perceptual sharpness and contrast on a 1-5 scale. Diagnostic confidence was assessed on a 1-3 scale. LD- and ELD-MRA were assessed against HD-MRA for sensitivity/specificity and agreement of vessel diameter measurements (aorta and pulmonary arteries). STATISTICAL TESTS: SNR, CNR, edge sharpness, and vessel diameters were compared between LD-, ELD-, and HD-MRA using one-way repeated measures analysis of variance with post-hoc t-tests. Perceptual quality and diagnostic confidence were compared using Friedman's test with post-hoc Wilcoxon signed-rank tests. Sensitivity/specificity was compared using McNemar's test. Agreement of vessel diameters was assessed using Bland-Altman analysis. RESULTS: SNR, CNR, edge sharpness, perceptual sharpness, and perceptual contrast were lower (P < 0.05) for LD-MRA compared to ELD-MRA and HD-MRA. SNR, CNR, edge sharpness, and perceptual contrast were comparable between ELD and HD-MRA, but perceptual sharpness was significantly lower. Sensitivity/specificity was 0.824/0.921 for LD-MRA and 0.882/0.960 for ELD-MRA. Diagnostic confidence was 2.72, 2.85, and 2.92 for LD, ELD, and HD-MRA, respectively (PLD-ELD , PLD-HD  < 0.05). Vessel diameter measurements were comparable, with biases of 0.238 (LD-MRA) and 0.278 mm (ELD-MRA). DATA CONCLUSION: Deep learning can improve contrast in LD cardiovascular MRA. LEVEL OF EVIDENCE LEVEL: 2 TECHNICAL EFFICACY: Stage 2.

Rapid whole-heart CMR with single volume super-resolution.

BACKGROUND: Three-dimensional, whole heart, balanced steady state free precession (WH-bSSFP) sequences provide delineation of intra-cardiac and vascular anatomy. However, they have long acquisition times. Here, we propose significant speed-ups using a deep-learning single volume super-resolution reconstruction, to recover high-resolution features from rapidly acquired low-resolution WH-bSSFP images. METHODS: A 3D residual U-Net was trained using synthetic data, created from a library of 500 high-resolution WH-bSSFP images by simulating 50% slice resolution and 50% phase resolution. The trained network was validated with 25 synthetic test data sets. Additionally, prospective low-resolution data and high-resolution data were acquired in 40 patients. In the prospective data, vessel diameters, quantitative and qualitative image quality, and diagnostic scoring was compared between the low-resolution, super-resolution and reference high-resolution WH-bSSFP data. RESULTS: The synthetic test data showed a significant increase in image quality of the low-resolution images after super-resolution reconstruction. Prospectively acquired low-resolution data was acquired ~× 3 faster than the prospective high-resolution data (173 s vs 488 s). Super-resolution reconstruction of the low-resolution data took < 1 s per volume. Qualitative image scores showed super-resolved images had better edge sharpness, fewer residual artefacts and less image distortion than low-resolution images, with similar scores to high-resolution data. Quantitative image scores showed super-resolved images had significantly better edge sharpness than low-resolution or high-resolution images, with significantly better signal-to-noise ratio than high-resolution data. Vessel diameters measurements showed over-estimation in the low-resolution measurements, compared to the high-resolution data. No significant differences and no bias was found in the super-resolution measurements in any of the great vessels. However, a small but significant for the underestimation was found in the proximal left coronary artery diameter measurement from super-resolution data. Diagnostic scoring showed that although super-resolution did not improve accuracy of diagnosis, it did improve diagnostic confidence compared to low-resolution imaging. CONCLUSION: This paper demonstrates the potential of using a residual U-Net for super-resolution reconstruction of rapidly acquired low-resolution whole heart bSSFP data within a clinical setting. We were able to train the network using synthetic training data from retrospective high-resolution whole heart data. The resulting network can be applied very quickly, making these techniques particularly appealing within busy clinical workflow. Thus, we believe that this technique may help speed up whole heart CMR in clinical practice.

Quantitative insertion-site sequencing (QIseq) for high throughput phenotyping of transposon mutants.

Genetic screening using random transposon insertions has been a powerful tool for uncovering biology in prokaryotes, where whole-genome saturating screens have been performed in multiple organisms. In eukaryotes, such screens have proven more problematic, in part because of the lack of a sensitive and robust system for identifying transposon insertion sites. We here describe quantitative insertion-site sequencing, or QIseq, which uses custom library preparation and Illumina sequencing technology and is able to identify insertion sites from both the 5' and 3' ends of the transposon, providing an inbuilt level of validation. The approach was developed using piggyBac mutants in the human malaria parasite Plasmodium falciparum but should be applicable to many other eukaryotic genomes. QIseq proved accurate, confirming known sites in >100 mutants, and sensitive, identifying and monitoring sites over a >10,000-fold dynamic range of sequence counts. Applying QIseq to uncloned parasites shortly after transfections revealed multiple insertions in mixed populations and suggests that >4000 independent mutants could be generated from relatively modest scales of transfection, providing a clear pathway to genome-scale screens in P. falciparum QIseq was also used to monitor the growth of pools of previously cloned mutants and reproducibly differentiated between deleterious and neutral mutations in competitive growth. Among the mutants with fitness defects was a mutant with a piggyBac insertion immediately upstream of the kelch protein K13 gene associated with artemisinin resistance, implying mutants in this gene may have competitive fitness costs. QIseq has the potential to enable the scale-up of piggyBac-mediated genetics across multiple eukaryotic systems.

The aorta after coarctation repair - effects of calibre and curvature on arterial haemodynamics.

BACKGROUND: Aortic shape has been proposed as an important determinant of adverse haemodynamics following coarctation repair. However, previous studies have not demonstrated a consistent relationship between shape and vascular load. In this study, 3D aortic shape was evaluated using principal component analysis (PCA), allowing investigation of the relationship between 3D shape and haemodynamics. METHODS: Sixty subjects (38 male, 25.0 ± 7.8 years) with repaired coarctation were recruited. Central aortic haemodynamics including wave intensity analysis were measured noninvasively using a combination of blood pressure and phase contrast cardiovascular magnetic resonance (CMR). 3D curvature and radius data were derived from CMR angiograms. PCA was separately performed on 3D radius and curvature data to assess the role of arch geometry on haemodynamics. Clinical findings were corroborated using 1D vascular models. RESULTS: There were no independent associations between 3D curvature and any hemodynamic parameters. However, the magnitude of the backwards compression wave was related to the 1st (r = - 0.36, p = 0.005), 3rd (r = 0.27, p = 0.036) and 4th (r = - 0.31, p = 0.017) principle components of radius. The 4th principle componentof radius also correlated with central aortic systolic pressure. These aortas had larger aortic roots, more transverse arch hypoplasia and narrower aortic isthmuses. CONCLUSIONS: There are major modes of variation in 3D aortic shape after coarctation repair witha modest association between variation in aortic radius and pathological wave reflections, but not with 3D curvature. Taken together, these data suggest that shape is not the major determinant of vascular load following coarctation repair, and calibre is more important than curvature.

Correction to: The aorta after coarctation repair - effects of calibre and curvature on arterial haemodynamics.

In the original version of this article [1], published on 11 April 2019, there is 1 error in the 'Conclusion' paragraph of the abstract.

Prognostic and functional implications of left atrial late gadolinium enhancement cardiovascular magnetic resonance.

BACKGROUND: Left atrial (LA) late gadolinium enhancement (LGE) on cardiovascular magnetic resonance (CMR) imaging is indicative of fibrosis, and has been correlated with reduced LA function, increased LA volume, and poor procedural outcomes in cohorts with atrial fibrillation (AF). However, the role of LGE as a prognostic biomarker for arrhythmia in cardiac disease has not been examined. METHODS: In this study, we assessed LA LGE using a 3D LGE CMR sequence to examine its relationships with new onset atrial arrhythmia, and LA and left ventricular (LV) mechanical function. RESULTS: LA LGE images were acquired in 111 patients undergoing CMR imaging, including 66 patients with no prior history of an atrial arrhythmia. During the median follow-up of 2.7 years (interquartile range (IQR) 1.8-3.7 years), 15/66 (23%) of patients developed a new atrial arrhythmia. LA LGE ≥10% of LA myocardial volume was significantly associated with an increased rate of new-onset atrial arrhythmia, with a hazard ratio of 3.16 (95% CI 1.14-8.72), p = 0.026. There were significant relationships between LA LGE and both LA ejection fraction (r = - 0.39, p < 0.0005) and echocardiographic LV septal e' (r = - 0.24, p = 0.04) and septal E/e' (r = 0.31, p = 0.007). CONCLUSIONS: Elevated LA LGE is associated with reduced LA function and reduced LV diastolic function. LA LGE is associated with new onset atrial arrhythmia during follow-up.

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.

Copy number variation arising from gene conversion on the human Y chromosome.

We describe the variation in copy number of a ~ 10 kb region overlapping the long intergenic noncoding RNA (lincRNA) gene, TTTY22, within the IR3 inverted repeat on the short arm of the human Y chromosome, leading to individuals with 0-3 copies of this region in the general population. Variation of this CNV is common, with 266 individuals having 0 copies, 943 (including the reference sequence) having 1, 23 having 2 copies, and two having 3 copies, and was validated by breakpoint PCR, fibre-FISH, and 10× Genomics Chromium linked-read sequencing in subsets of 1234 individuals from the 1000 Genomes Project. Mapping the changes in copy number to the phylogeny of these Y chromosomes previously established by the Project identified at least 20 mutational events, and investigation of flanking paralogous sequence variants showed that the mutations involved flanking sequences in 18 of these, and could extend over > 30 kb of DNA. While either gene conversion or double crossover between misaligned sister chromatids could formally explain the 0-2 copy events, gene conversion is the more likely mechanism, and these events include the longest non-allelic gene conversion reported thus far. Chromosomes with three copies of this CNV have arisen just once in our data set via another mechanism: duplication of 420 kb that places the third copy 230 kb proximal to the existing proximal copy. Our results establish gene conversion as a previously under-appreciated mechanism of generating copy number changes in humans and reveal the exceptionally large size of the conversion events that can occur.

G&T-seq: parallel sequencing of single-cell genomes and transcriptomes.

The simultaneous sequencing of a single cell's genome and transcriptome offers a powerful means to dissect genetic variation and its effect on gene expression. Here we describe G&T-seq, a method for separating and sequencing genomic DNA and full-length mRNA from single cells. By applying G&T-seq to over 220 single cells from mice and humans, we discovered cellular properties that could not be inferred from DNA or RNA sequencing alone.

Development and validation of a comprehensive genomic diagnostic tool for myeloid malignancies.

The diagnosis of hematologic malignancies relies on multidisciplinary workflows involving morphology, flow cytometry, cytogenetic, and molecular genetic analyses. Advances in cancer genomics have identified numerous recurrent mutations with clear prognostic and/or therapeutic significance to different cancers. In myeloid malignancies, there is a clinical imperative to test for such mutations in mainstream diagnosis; however, progress toward this has been slow and piecemeal. Here we describe Karyogene, an integrated targeted resequencing/analytical platform that detects nucleotide substitutions, insertions/deletions, chromosomal translocations, copy number abnormalities, and zygosity changes in a single assay. We validate the approach against 62 acute myeloid leukemia, 50 myelodysplastic syndrome, and 40 blood DNA samples from individuals without evidence of clonal blood disorders. We demonstrate robust detection of sequence changes in 49 genes, including difficult-to-detect mutations such as FLT3 internal-tandem and mixed-lineage leukemia (MLL) partial-tandem duplications, and clinically significant chromosomal rearrangements including MLL translocations to known and unknown partners, identifying the novel fusion gene MLL-DIAPH2 in the process. Additionally, we identify most significant chromosomal gains and losses, and several copy neutral loss-of-heterozygosity mutations at a genome-wide level, including previously unreported changes such as homozygosity for DNMT3A R882 mutations. Karyogene represents a dependable genomic diagnosis platform for translational research and for the clinical management of myeloid malignancies, which can be readily adapted for use in other cancers.

A comprehensive characterization of myocardial and vascular phenotype in pediatric chronic kidney disease using cardiovascular magnetic resonance imaging.

BACKGROUND: Children with chronic kidney disease (CKD) have increased cardiovascular mortality. Identifying high-risk children who may benefit from further therapeutic intervention is difficult as cardiovascular abnormalities are subtle. Although transthoracic echocardiography may be used to detect sub-clinical abnormalities, it has well-known problems with reproducibility that limit its ability to accurately detect these changes. Cardiovascular magnetic resonance (CMR) is the reference standard method for assessing blood flow, cardiac structure and function. Furthermore, recent innovations enable the assessment of radial and longitudinal myocardial velocity, such that detection of sub-clinical changes is now possible. Thus, CMR may be ideal for cardiovascular assessment in pediatric CKD. This study aims to comprehensively assess cardiovascular function in pediatric CKD using CMR and determine its relationship with CKD severity. METHODS: A total of 120 children (40 mild, 40 moderate, 20 severe pre-dialysis CKD subjects and 20 healthy controls) underwent CMR with non-invasive blood pressure (BP) measurements. Cardiovascular parameters measured included systemic vascular resistance (SVR), total arterial compliance (TAC), left ventricular (LV) structure, ejection fraction (EF), cardiac timings, radial and longitudinal systolic and diastolic myocardial velocities. Between group comparisons and regression modelling were used to identify abnormalities in CKD and determine the effects of renal severity on myocardial function. RESULTS: The elevation in mean BP in CKD was accompanied by significantly increased afterload (SVR), without evidence of arterial stiffness (TAC) or increased fluid overload. Left ventricular volumes and global function were not abnormal in CKD. However, there was evidence of LV remodelling, prolongation of isovolumic relaxation time and reduced systolic and diastolic myocardial velocities. CONCLUSION: Abnormal cardiovascular function is evident in pre-dialysis pediatric CKD. Novel CMR biomarkers may be useful for the detection of subtle abnormalities in this population. Further studies are needed to determine to prognostic value of these biomarkers.

Analysis of Plasmodium falciparum diversity in natural infections by deep sequencing.

Malaria elimination strategies require surveillance of the parasite population for genetic changes that demand a public health response, such as new forms of drug resistance. Here we describe methods for the large-scale analysis of genetic variation in Plasmodium falciparum by deep sequencing of parasite DNA obtained from the blood of patients with malaria, either directly or after short-term culture. Analysis of 86,158 exonic single nucleotide polymorphisms that passed genotyping quality control in 227 samples from Africa, Asia and Oceania provides genome-wide estimates of allele frequency distribution, population structure and linkage disequilibrium. By comparing the genetic diversity of individual infections with that of the local parasite population, we derive a metric of within-host diversity that is related to the level of inbreeding in the population. An open-access web application has been established for the exploration of regional differences in allele frequency and of highly differentiated loci in the P. falciparum genome.

The TraDIS toolkit: sequencing and analysis for dense transposon mutant libraries.

UNLABELLED: Transposon insertion sequencing is a high-throughput technique for assaying large libraries of otherwise isogenic transposon mutants providing insight into gene essentiality, gene function and genetic interactions. We previously developed the Transposon Directed Insertion Sequencing (TraDIS) protocol for this purpose, which utilizes shearing of genomic DNA followed by specific PCR amplification of transposon-containing fragments and Illumina sequencing. Here we describe an optimized high-yield library preparation and sequencing protocol for TraDIS experiments and a novel software pipeline for analysis of the resulting data. The Bio-Tradis analysis pipeline is implemented as an extensible Perl library which can either be used as is, or as a basis for the development of more advanced analysis tools. This article can serve as a general reference for the application of the TraDIS methodology. AVAILABILITY AND IMPLEMENTATION: The optimized sequencing protocol is included as supplementary information. The Bio-Tradis analysis pipeline is available under a GPL license at https://github.com/sanger-pathogens/Bio-Tradis CONTACT: parkhill@sanger.ac.uk SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

PET-CMR in heart failure - synergistic or redundant imaging?

Imaging in heart failure (HF) provides data for diagnosis, prognosis and disease monitoring. Both MRI and nuclear imaging techniques have been successfully used for this purpose in HF. Positron Emission Tomography-Cardiac Magnetic Resonance (PET-CMR) is an example of a new multimodality diagnostic imaging technique with potential applications in HF. The threshold for adopting a new diagnostic tool to clinical practice must necessarily be high, lest they exacerbate costs without improving care. New modalities must demonstrate clinical superiority, or at least equivalence, combined with another important advantage, such as lower cost or improved patient safety. The purpose of this review is to outline the current status of multimodality PET-CMR with regard to HF applications, and determine whether the clinical utility of this new technology justifies the cost.

Chromosomal rearrangements maintain a polymorphic supergene controlling butterfly mimicry.

Supergenes are tight clusters of loci that facilitate the co-segregation of adaptive variation, providing integrated control of complex adaptive phenotypes. Polymorphic supergenes, in which specific combinations of traits are maintained within a single population, were first described for 'pin' and 'thrum' floral types in Primula and Fagopyrum, but classic examples are also found in insect mimicry and snail morphology. Understanding the evolutionary mechanisms that generate these co-adapted gene sets, as well as the mode of limiting the production of unfit recombinant forms, remains a substantial challenge. Here we show that individual wing-pattern morphs in the polymorphic mimetic butterfly Heliconius numata are associated with different genomic rearrangements at the supergene locus P. These rearrangements tighten the genetic linkage between at least two colour-pattern loci that are known to recombine in closely related species, with complete suppression of recombination being observed in experimental crosses across a 400-kilobase interval containing at least 18 genes. In natural populations, notable patterns of linkage disequilibrium (LD) are observed across the entire P region. The resulting divergent haplotype clades and inversion breakpoints are found in complete association with wing-pattern morphs. Our results indicate that allelic combinations at known wing-patterning loci have become locked together in a polymorphic rearrangement at the P locus, forming a supergene that acts as a simple switch between complex adaptive phenotypes found in sympatry. These findings highlight how genomic rearrangements can have a central role in the coexistence of adaptive phenotypes involving several genes acting in concert, by locally limiting recombination and gene flow.