Clonal fitness inferred from time-series modelling of single-cell cancer genomes.

Progress in defining genomic fitness landscapes in cancer, especially those defined by copy number alterations (CNAs), has been impeded by lack of time-series single-cell sampling of polyclonal populations and temporal statistical models1-7. Here we generated 42,000 genomes from multi-year time-series single-cell whole-genome sequencing of breast epithelium and primary triple-negative breast cancer (TNBC) patient-derived xenografts (PDXs), revealing the nature of CNA-defined clonal fitness dynamics induced by TP53 mutation and cisplatin chemotherapy. Using a new Wright-Fisher population genetics model8,9 to infer clonal fitness, we found that TP53 mutation alters the fitness landscape, reproducibly distributing fitness over a larger number of clones associated with distinct CNAs. Furthermore, in TNBC PDX models with mutated TP53, inferred fitness coefficients from CNA-based genotypes accurately forecast experimentally enforced clonal competition dynamics. Drug treatment in three long-term serially passaged TNBC PDXs resulted in cisplatin-resistant clones emerging from low-fitness phylogenetic lineages in the untreated setting. Conversely, high-fitness clones from treatment-naive controls were eradicated, signalling an inversion of the fitness landscape. Finally, upon release of drug, selection pressure dynamics were reversed, indicating a fitness cost of treatment resistance. Together, our findings define clonal fitness linked to both CNA and therapeutic resistance in polyclonal tumours.

Radiomics analysis to predict pulmonary nodule malignancy using machine learning approaches.

BACKGROUND: Low-dose CT screening can reduce lung cancer-related mortality. However, most screen-detected pulmonary abnormalities do not develop into cancer and it often remains challenging to identify malignant nodules, particularly among indeterminate nodules. We aimed to develop and assess prediction models based on radiological features to discriminate between benign and malignant pulmonary lesions detected on a baseline screen. METHODS: Using four international lung cancer screening studies, we extracted 2060 radiomic features for each of 16 797 nodules (513 malignant) among 6865 participants. After filtering out low-quality radiomic features, 642 radiomic and 9 epidemiological features remained for model development. We used cross-validation and grid search to assess three machine learning (ML) models (eXtreme Gradient Boosted Trees, random forest, least absolute shrinkage and selection operator (LASSO)) for their ability to accurately predict risk of malignancy for pulmonary nodules. We report model performance based on the area under the curve (AUC) and calibration metrics in the held-out test set. RESULTS: The LASSO model yielded the best predictive performance in cross-validation and was fit in the full training set based on optimised hyperparameters. Our radiomics model had a test-set AUC of 0.93 (95% CI 0.90 to 0.96) and outperformed the established Pan-Canadian Early Detection of Lung Cancer model (AUC 0.87, 95% CI 0.85 to 0.89) for nodule assessment. Our model performed well among both solid (AUC 0.93, 95% CI 0.89 to 0.97) and subsolid nodules (AUC 0.91, 95% CI 0.85 to 0.95). CONCLUSIONS: We developed highly accurate ML models based on radiomic and epidemiological features from four international lung cancer screening studies that may be suitable for assessing indeterminate screen-detected pulmonary nodules for risk of malignancy.

Vibrational spectroscopic profiling of biomolecular interactions between oak powdery mildew and oak leaves.

Oak powdery mildew, caused by the biotrophic fungus Erysiphe alphitoides, is a prevalent disease affecting oak trees, such as English oak (Quercus robur). While mature oak populations are generally less susceptible to this disease, it can endanger young oak seedlings and new leaves on mature trees. Although disruptions of photosynthate and carbohydrate translocation have been observed, accurately detecting and understanding the specific biomolecular interactions between the fungus and the leaves of oak trees is currently lacking. Herein, via hybrid Raman spectroscopy combined with an advanced artificial neural network algorithm, the underpinning biomolecular interactions between biological soft matter, i.e., Quercus robur leaves and Erysiphe alphitoides, are investigated and profiled, generating a spectral library and shedding light on the changes induced by fungal infection and the tree's defence response. The adaxial surfaces of oak leaves are categorised based on either the presence or absence of Erysiphe alphitoides mildew and further distinguishing between covered or not covered infected leaf tissues, yielding three disease classes including healthy controls, non-mildew covered and mildew-covered. By analysing spectral changes between each disease category per tissue type, we identified important biomolecular interactions including disruption of chlorophyll in the non-vein and venule tissues, pathogen-induced degradation of cellulose and pectin and tree-initiated lignification of cell walls in response, amongst others, in lateral vein and mid-vein tissues. Via our developed computational algorithm, the underlying biomolecular differences between classes were identified and allowed accurate and rapid classification of disease with high accuracy of 69.6% for non-vein, 73.5% for venule, 82.1% for lateral vein and 85.6% for mid-vein tissues. Interfacial wetting differences between non-mildew covered and mildew-covered tissue were further analysed on the surfaces of non-vein and venule tissue. The overall results demonstrated the ability of Raman spectroscopy, combined with advanced AI, to act as a powerful and specific tool to probe foliar interactions between forest pathogens and host trees with the simultaneous potential to probe and catalogue molecular interactions between biological soft matter, paving the way for exploring similar relations in broader forest tree-pathogen systems.

Loss of apelin blocks the emergence of sprouting angiogenesis in experimental tumors.

Angiogenesis inhibitor drugs targeting vascular endothelial growth factor (VEGF) signaling to the endothelial cell (EC) are used to treat various cancer types. However, primary or secondary resistance to therapy is common. Clinical and pre-clinical studies suggest that alternative pro-angiogenic factors are upregulated after VEGF pathway inhibition. Therefore, identification of alternative pro-angiogenic pathway(s) is critical for the development of more effective anti-angiogenic therapy. Here we study the role of apelin as a pro-angiogenic G-protein-coupled receptor ligand in tumor growth and angiogenesis. We found that loss of apelin in mice delayed the primary tumor growth of Lewis lung carcinoma 1 and B16F10 melanoma when combined with the VEGF receptor tyrosine kinase inhibitor, sunitinib. Targeting apelin in combination with sunitinib markedly reduced the tumor vessel density, and decreased microvessel remodeling. Apelin loss reduced angiogenic sprouting and tip cell marker gene expression in comparison to the sunitinib-alone-treated mice. Single-cell RNA sequencing of tumor EC demonstrated that the loss of apelin prevented EC tip cell differentiation. Thus, apelin is a potent pro-angiogenic cue that supports initiation of tumor neovascularization. Together, our data suggest that targeting apelin may be useful as adjuvant therapy in combination with VEGF signaling inhibition to inhibit the growth of advanced tumors.

Probabilistic cell-type assignment of single-cell RNA-seq for tumor microenvironment profiling.

Single-cell RNA sequencing has enabled the decomposition of complex tissues into functionally distinct cell types. Often, investigators wish to assign cells to cell types through unsupervised clustering followed by manual annotation or via 'mapping' to existing data. However, manual interpretation scales poorly to large datasets, mapping approaches require purified or pre-annotated data and both are prone to batch effects. To overcome these issues, we present CellAssign, a probabilistic model that leverages prior knowledge of cell-type marker genes to annotate single-cell RNA sequencing data into predefined or de novo cell types. CellAssign automates the process of assigning cells in a highly scalable manner across large datasets while controlling for batch and sample effects. We demonstrate the advantages of CellAssign through extensive simulations and analysis of tumor microenvironment composition in high-grade serous ovarian cancer and follicular lymphoma.

The Landscape of COVID-19 Research in the United States: a Cross-sectional Study of Randomized Trials Registered on ClinicalTrials.Gov.

IMPORTANCE: SARS-CoV-2 has infected over 200 million people worldwide, resulting in more than 4 million deaths. Randomized controlled trials are the single best tool to identify effective treatments against this novel pathogen. OBJECTIVE: To describe the characteristics of randomized controlled trials of treatments for COVID-19 in the United States launched in the first 9 months of the pandemic. Design, Setting, and Participants We conducted a cross-sectional study of all completed or actively enrolling randomized, interventional, clinical trials for the treatment of COVID-19 in the United States registered on www.clinicaltrials.gov as of August 10, 2020. We excluded trials of vaccines and other interventions intended to prevent COVID-19. Main Outcomes and Measures We used descriptive statistics to characterize the clinical trials and the statistical power for the available studies. For the late-phase trials (i.e., phase 3 and 2/3 studies), we compared the geographic distribution of the clinical trials with the geographic distribution of people diagnosed with COVID-19. RESULTS: We identified 200 randomized controlled trials of treatments for people with COVID-19. Across all trials, 87 (43.5%) were single-center, 64 (32.0%) were unblinded, and 80 (40.0%) were sponsored by industry. The most common treatments included monoclonal antibodies (N=46 trials), small molecule immunomodulators (N=28), antiviral medications (N=24 trials), and hydroxychloroquine (N=20 trials). Of the 9 trials completed by August 2020, the median sample size was 450 (IQR 67-1113); of the 191 ongoing trials, the median planned sample size was 150 (IQR 60-400). Of the late-phase trials (N=54), the most common primary outcome was a severity scale (N=23, 42.6%), followed by a composite of mortality and ventilation (N=10, 18.5%), and mortality alone (N=6, 11.1%). Among these late-phase trials, all trials of antivirals, monoclonal antibodies, or chloroquine/hydroxychloroquine had a power of less than 25% to detect a 20% relative risk reduction in mortality. Had the individual trials for a given class of treatments instead formed a single trial, the power to detect that same reduction in mortality would have been greater than 98%. There was large variability in access to trials with the highest number of trials per capita in the Northeast and the lowest in the Midwest. CONCLUSIONS AND RELEVANCE: A large number of randomized trials were launched early in the pandemic to evaluate treatments for COVID-19. However, many trials were underpowered for important clinical endpoints and substantial geographic disparities were observed, highlighting the importance of improving national clinical trial infrastructure.

Modelling hereditary diffuse gastric cancer initiation using transgenic mouse-derived gastric organoids and single-cell sequencing.

Hereditary diffuse gastric cancer (HDGC) is a cancer syndrome caused by germline variants in CDH1, the gene encoding the cell-cell adhesion molecule E-cadherin. Loss of E-cadherin in cancer is associated with cellular dedifferentiation and poor prognosis, but the mechanisms through which CDH1 loss initiates HDGC are not known. Using single-cell RNA sequencing, we explored the transcriptional landscape of a murine organoid model of HDGC to characterize the impact of CDH1 loss in early tumourigenesis. Progenitor populations of stratified squamous and simple columnar epithelium, characteristic of the mouse stomach, showed lineage-specific transcriptional programs. Cdh1 inactivation resulted in shifts along the squamous differentiation trajectory associated with aberrant expression of genes central to gastrointestinal epithelial differentiation. Cytokeratin 7 (CK7), encoded by the differentiation-dependent gene Krt7, was a specific marker for early neoplastic lesions in CDH1 carriers. Our findings suggest that deregulation of developmental transcriptional programs may precede malignancy in HDGC. © 2021 The Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

The structural features that distinguish PD-L2 from PD-L1 emerged in placental mammals.

Programmed cell death protein 1 (PD-1) is an inhibitory receptor on T lymphocytes that is critical for modulating adaptive immunity. As such, it has been successfully exploited for cancer immunotherapy. Programmed death ligand 1 (PD-L1) and PD-L2 are ligands for PD-1; the former is ubiquitously expressed in inflamed tissues, whereas the latter is restricted to antigen-presenting cells. PD-L2 binds to PD-1 with 3-fold stronger affinity compared with PD-L1. To date, this affinity discrepancy has been attributed to a tryptophan (W110PD-L2) that is unique to PD-L2 and has been assumed to fit snuggly into a pocket on the PD-1 surface. Contrary to this model, using surface plasmon resonance to monitor real-time binding of recombinantly-expressed and -purified proteins, we found that W110PD-L2 acts as an "elbow" that helps shorten PD-L2 engagement with PD-1 and therefore lower affinity. Furthermore, we identified a "latch" between the C and D ß-strands of the binding face as the source of the PD-L2 affinity advantage. We show that the 3-fold affinity advantage of PD-L2 is the consequence of these two opposing features, the W110PD-L2 "elbow" and a C-D region "latch." Interestingly, using phylogenetic analysis, we found that these features evolved simultaneously upon the emergence of placental mammals, suggesting that PD-L2-affinity tuning was part of the alterations to the adaptive immune system required for placental gestation.

Quantitative phosphoproteomic analysis reveals involvement of PD-1 in multiple T cell functions.

Programmed cell death protein 1 (PD-1) is a critical inhibitory receptor that limits excessive T cell responses. Cancer cells have evolved to evade these immunoregulatory mechanisms by upregulating PD-1 ligands and preventing T cell-mediated anti-tumor responses. Consequently, therapeutic blockade of PD-1 enhances T cell-mediated anti-tumor immunity, but many patients do not respond and a significant proportion develop inflammatory toxicities. To improve anti-cancer therapy, it is critical to reveal the mechanisms by which PD-1 regulates T cell responses. We performed global quantitative phosphoproteomic interrogation of PD-1 signaling in T cells. By complementing our analysis with functional validation assays, we show that PD-1 targets tyrosine phosphosites that mediate proximal T cell receptor signaling, cytoskeletal organization, and immune synapse formation. PD-1 ligation also led to differential phosphorylation of serine and threonine sites within proteins regulating T cell activation, gene expression, and protein translation. In silico predictions revealed that kinase/substrate relationships engaged downstream of PD-1 ligation. These insights uncover the phosphoproteomic landscape of PD-1-triggered pathways and reveal novel PD-1 substrates that modulate diverse T cell functions and may serve as future therapeutic targets. These data are a useful resource in the design of future PD-1-targeting therapeutic approaches.

Single cell transcriptomes of normal endometrial derived organoids uncover novel cell type markers and cryptic differentiation of primary tumours.

Endometrial carcinoma, the most common gynaecological cancer, develops from endometrial epithelium which is composed of secretory and ciliated cells. Pathologic classification is unreliable and there is a need for prognostic tools. We used single cell sequencing to study organoid model systems derived from normal endometrial endometrium to discover novel markers specific for endometrial ciliated or secretory cells. A marker of secretory cells (MPST) and several markers of ciliated cells (FAM92B, WDR16, and DYDC2) were validated by immunohistochemistry on organoids and tissue sections. We performed single cell sequencing on endometrial and ovarian tumours and found both secretory-like and ciliated-like tumour cells. We found that ciliated cell markers (DYDC2, CTH, FOXJ1, and p73) and the secretory cell marker MPST were expressed in endometrial tumours and positively correlated with disease-specific and overall survival of endometrial cancer patients. These findings suggest that expression of differentiation markers in tumours correlates with less aggressive disease, as would be expected for tumours that retain differentiation capacity, albeit cryptic in the case of ciliated cells. These markers could be used to improve the risk stratification of endometrial cancer patients, thereby improving their management. We further assessed whether consideration of MPST expression could refine the ProMiSE molecular classification system for endometrial tumours. We found that higher expression levels of MPST could be used to refine stratification of three of the four ProMiSE molecular subgroups, and that any level of MPST expression was able to significantly refine risk stratification of the copy number high subgroup which has the worst prognosis. Taken together, this shows that single cell sequencing of putative cells of origin has the potential to uncover novel biomarkers that could be used to guide management of cancers. © 2020 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

A descriptive marker gene approach to single-cell pseudotime inference.

Motivation: Pseudotime estimation from single-cell gene expression data allows the recovery of temporal information from otherwise static profiles of individual cells. Conventional pseudotime inference methods emphasize an unsupervised transcriptome-wide approach and use retrospective analysis to evaluate the behaviour of individual genes. However, the resulting trajectories can only be understood in terms of abstract geometric structures and not in terms of interpretable models of gene behaviour. Results: Here we introduce an orthogonal Bayesian approach termed 'Ouija' that learns pseudotimes from a small set of marker genes that might ordinarily be used to retrospectively confirm the accuracy of unsupervised pseudotime algorithms. Crucially, we model these genes in terms of switch-like or transient behaviour along the trajectory, allowing us to understand why the pseudotimes have been inferred and learn informative parameters about the behaviour of each gene. Since each gene is associated with a switch or peak time the genes are effectively ordered along with the cells, allowing each part of the trajectory to be understood in terms of the behaviour of certain genes. We demonstrate that this small panel of marker genes can recover pseudotimes that are consistent with those obtained using the entire transcriptome. Furthermore, we show that our method can detect differences in the regulation timings between two genes and identify 'metastable' states-discrete cell types along the continuous trajectories-that recapitulate known cell types. Availability and implementation: An open source implementation is available as an R package at http://www.github.com/kieranrcampbell/ouija and as a Python/TensorFlow package at http://www.github.com/kieranrcampbell/ouijaflow. Supplementary information: Supplementary data are available at Bioinformatics online.

Computational modelling in single-cell cancer genomics: methods and future directions.

Single-cell technologies have revolutionized biomedical research by enabling scalable measurement of the genome, transcriptome, proteome, and epigenome of multiple systems at single-cell resolution. Now widely applied to cancer models, these assays offer new insights into tumour heterogeneity, which underlies cancer initiation, progression, and relapse. However, the large quantities of high-dimensional, noisy data produced by single-cell assays can complicate data analysis, obscuring biological signals with technical artifacts. In this review article, we outline the major challenges in analyzing single-cell cancer genomics data and survey the current computational tools available to tackle these. We further outline unsolved problems that we consider major opportunities for future methods development to help interpret the vast quantities of data being generated.

switchde: inference of switch-like differential expression along single-cell trajectories.

Motivation: Pseudotime analyses of single-cell RNA-seq data have become increasingly common. Typically, a latent trajectory corresponding to a biological process of interest-such as differentiation or cell cycle-is discovered. However, relatively little attention has been paid to modelling the differential expression of genes along such trajectories. Results: We present switchde , a statistical framework and accompanying R package for identifying switch-like differential expression of genes along pseudotemporal trajectories. Our method includes fast model fitting that provides interpretable parameter estimates corresponding to how quickly a gene is up or down regulated as well as where in the trajectory such regulation occurs. It also reports a P -value in favour of rejecting a constant-expression model for switch-like differential expression and optionally models the zero-inflation prevalent in single-cell data. Availability and Implementation: The R package switchde is available through the Bioconductor project at https://bioconductor.org/packages/switchde . Contact: kieran.campbell@sjc.ox.ac.uk. Supplementary information: Supplementary data are available at Bioinformatics online.

Scater: pre-processing, quality control, normalization and visualization of single-cell RNA-seq data in R.

Motivation: Single-cell RNA sequencing (scRNA-seq) is increasingly used to study gene expression at the level of individual cells. However, preparing raw sequence data for further analysis is not a straightforward process. Biases, artifacts and other sources of unwanted variation are present in the data, requiring substantial time and effort to be spent on pre-processing, quality control (QC) and normalization. Results: We have developed the R/Bioconductor package scater to facilitate rigorous pre-processing, quality control, normalization and visualization of scRNA-seq data. The package provides a convenient, flexible workflow to process raw sequencing reads into a high-quality expression dataset ready for downstream analysis. scater provides a rich suite of plotting tools for single-cell data and a flexible data structure that is compatible with existing tools and can be used as infrastructure for future software development. Availability and Implementation: The open-source code, along with installation instructions, vignettes and case studies, is available through Bioconductor at http://bioconductor.org/packages/scater . Contact: davis@ebi.ac.uk. Supplementary information: Supplementary data are available at Bioinformatics online.

The use of a quantitative structure-activity relationship (QSAR) model to predict GABA-A receptor binding of newly emerging benzodiazepines.

The illicit market for new psychoactive substances is forever expanding. Benzodiazepines and their derivatives are one of a number of groups of these substances and thus far their number has grown year upon year. For both forensic and clinical purposes it is important to be able to rapidly understand these emerging substances. However as a consequence of the illicit nature of these compounds, there is a deficiency in the pharmacological data available for these 'new' benzodiazepines. In order to further understand the pharmacology of 'new' benzodiazepines we utilised a quantitative structure-activity relationship (QSAR) approach. A set of 69 benzodiazepine-based compounds was analysed to develop a QSAR training set with respect to published binding values to GABAA receptors. The QSAR model returned an R2 value of 0.90. The most influential factors were found to be the positioning of two H-bond acceptors, two aromatic rings and a hydrophobic group. A test set of nine random compounds was then selected for internal validation to determine the predictive ability of the model and gave an R2 value of 0.86 when comparing the binding values with their experimental data. The QSAR model was then used to predict the binding for 22 benzodiazepines that are classed as new psychoactive substances. This model will allow rapid prediction of the binding activity of emerging benzodiazepines in a rapid and economic way, compared with lengthy and expensive in vitro/in vivo analysis. This will enable forensic chemists and toxicologists to better understand both recently developed compounds and prediction of substances likely to emerge in the future.

Order Under Uncertainty: Robust Differential Expression Analysis Using Probabilistic Models for Pseudotime Inference.

Single cell gene expression profiling can be used to quantify transcriptional dynamics in temporal processes, such as cell differentiation, using computational methods to label each cell with a 'pseudotime' where true time series experimentation is too difficult to perform. However, owing to the high variability in gene expression between individual cells, there is an inherent uncertainty in the precise temporal ordering of the cells. Pre-existing methods for pseudotime estimation have predominantly given point estimates precluding a rigorous analysis of the implications of uncertainty. We use probabilistic modelling techniques to quantify pseudotime uncertainty and propagate this into downstream differential expression analysis. We demonstrate that reliance on a point estimate of pseudotime can lead to inflated false discovery rates and that probabilistic approaches provide greater robustness and measures of the temporal resolution that can be obtained from pseudotime inference.

Assessing white matter microstructure of the newborn with multi-shell diffusion MRI and biophysical compartment models.

Brain white matter connections have become a focus of major interest with important maturational processes occurring in newborns. To study the complex microstructural developmental changes in-vivo, it is imperative that non-invasive neuroimaging approaches are developed for this age-group. Multi-b-value diffusion weighted imaging data were acquired in 13 newborns, and the biophysical compartment diffusion models CHARMED-light and NODDI, providing new microstructural parameters such as intra-neurite volume fraction (νin) and neurite orientation dispersion index (ODI), were developed for newborn data. Comparative analysis was performed and twenty ROIs in the white matter were investigated. Diffusion tensor imaging and both biophysical compartment models highlighted the compact and oriented structure of the corpus-callosum with the highest FA and νin values and the smallest ODI values. We could clearly differentiate, using the FA, νin and ODI, the posterior and anterior internal capsule representing similar cellular structure but with different maturation (i.e. partially myelinated and absence of myelin, respectively). Late maturing regions (external capsule and periventricular crossroads of pathways) had lower νin values, but displayed significant differences in ODI. The compartmented models CHARMED-light and NODDI bring new indices corroborating the cellular architectures, with the lowest νin, reflecting the late maturation of areas with thin non-myelinated fibers, and with highest ODI indicating the presence of fiber crossings and fanning. The application of biophysical compartment diffusion models adds new insights to the brain white matter development in vivo.

Improving T2 -weighted imaging at high field through the use of kT -points.

PURPOSE: At high magnetic field strengths (B(0) ≥ 3 T), the shorter radiofrequency wavelength produces an inhomogeneous distribution of the transmit magnetic field. This can lead to variable contrast across the brain which is particularly pronounced in T(2) -weighted imaging that requires multiple radiofrequency pulses. To obtain T(2) -weighted images with uniform contrast throughout the whole brain at 7 T, short (2-3 ms) 3D tailored radiofrequency pulses (kT -points) were integrated into a 3D variable flip angle turbo spin echo sequence. METHODS: The excitation and refocusing "hard" pulses of a variable flip angle turbo spin echo sequence were replaced with kT -point pulses. Spatially resolved extended phase graph simulations and in vivo acquisitions at 7 T, utilizing both single channel and parallel-transmit systems, were used to test different kT -point configurations. RESULTS: Simulations indicated that an extended optimized k-space trajectory ensured a more homogeneous signal throughout images. In vivo experiments showed that high quality T(2) -weighted brain images with uniform signal and contrast were obtained at 7 T by using the proposed methodology. CONCLUSION: This work demonstrates that T(2) -weighted images devoid of artifacts resulting from B(1)(+) inhomogeneity can be obtained at high field through the optimization of extended kT -point pulses.

Dielectric pads and low- B1+ adiabatic pulses: complementary techniques to optimize structural T1 w whole-brain MP2RAGE scans at 7 tesla.

PURPOSE: To evaluate the combination of low-B1 (+) adiabatic pulses and high permittivity (εr ≈ 165) dielectric pads effectiveness to reproducibly improve the inversion efficiency for whole-brain MP2RAGE scans, at ultra-high field. MATERIALS AND METHODS: Two low-B1 (+) adiabatic pulses, HS8 and TR-FOCI, were compared with the conventional HS1 adiabatic pulse in MP2RAGE acquisitions of four subjects at 7 Tesla. The uniform MP2RAGE images were qualitatively assessed for poor inversion artifacts by trained observers. Each subject was rescanned using dielectric pads. Eight further subjects underwent two MP2RAGE scan sessions using dielectric pads and the TR-FOCI adiabatic pulse. RESULTS: The HS8 and TR-FOCI pulses improved inversion coverage in all subjects compared with the HS1 pulse. However, in subjects whose head lengths are large (≥136 mm) relative to the coil's z-coverage, the B1 (+) field over the cerebellum was insufficient to cause inversion. Dielectric pads increase the B1 (+) field, by ∼50%, over the cerebellum, which in conjunction with the TR-FOCI pulse, reproducibly improves the inversion efficiency over the whole brain for subjects with head lengths ≤155 mm. Minor residual inversion artifacts were present in three of eight subjects (head lengths ≥155 mm). CONCLUSION: The complementary techniques of low-B1 (+) adiabatic RF pulses and high permittivity dielectric pads allow whole-brain structural T1 w images to be reliably acquired at ultra-high field. J. Magn. Reson. Imaging 2014;40:804-812. © 2013 Wiley Periodicals, Inc.