Cell-type-specific consequences of mosaic structural variants in hematopoietic stem and progenitor cells.

The functional impact and cellular context of mosaic structural variants (mSVs) in normal tissues is understudied. Utilizing Strand-seq, we sequenced 1,133 single-cell genomes from 19 human donors of increasing age, and discovered the heterogeneous mSV landscapes of hematopoietic stem and progenitor cells. While mSVs are continuously acquired throughout life, expanded subclones in our cohort are confined to individuals >60. Cells already harboring mSVs are more likely to acquire additional somatic structural variants, including megabase-scale segmental aneuploidies. Capitalizing on comprehensive single-cell micrococcal nuclease digestion with sequencing reference data, we conducted high-resolution cell-typing for eight hematopoietic stem and progenitor cells. Clonally expanded mSVs disrupt normal cellular function by dysregulating diverse cellular pathways, and enriching for myeloid progenitors. Our findings underscore the contribution of mSVs to the cellular and molecular phenotypes associated with the aging hematopoietic system, and establish a foundation for deciphering the molecular links between mSVs, aging and disease susceptibility in normal tissues.

Prenatal exposure to maternal disadvantage-related inflammatory biomarkers: associations with neonatal white matter microstructure.

Prenatal exposure to heightened maternal inflammation has been associated with adverse neurodevelopmental outcomes, including atypical brain maturation and psychiatric illness. In mothers experiencing socioeconomic disadvantage, immune activation can be a product of the chronic stress inherent to such environmental hardship. While growing preclinical and clinical evidence has shown links between altered neonatal brain development and increased inflammatory states in utero, the potential mechanism by which socioeconomic disadvantage differentially impacts neural-immune crosstalk remains unclear. In the current study, we investigated associations between socioeconomic disadvantage, gestational inflammation, and neonatal white matter microstructure in 320 mother-infant dyads over-sampled for poverty. We analyzed maternal serum levels of four cytokines (IL-6, IL-8, IL-10, TNF-α) over the course of pregnancy in relation to offspring white matter microstructure and socioeconomic disadvantage. Higher average maternal IL-6 was associated with very low socioeconomic status (SES; INR < 200% poverty line) and lower neonatal corticospinal fractional anisotropy (FA) and lower uncinate axial diffusivity (AD). No other cytokine was associated with SES. Higher average maternal IL-10 was associated with lower FA and higher radial diffusivity (RD) in corpus callosum and corticospinal tracts, higher optic radiation RD, lower uncinate AD, and lower FA in inferior fronto-occipital fasciculus and anterior limb of internal capsule tracts. SES moderated the relationship between average maternal TNF-α levels during gestation and neonatal white matter diffusivity. When these interactions were decomposed, the patterns indicated that this association was significant and positive among very low SES neonates, whereby TNF-α was inversely and significantly associated with inferior cingulum AD. By contrast, among the more advantaged neonates (lower-to-higher SES [INR ≥ 200% poverty line]), TNF-α was positively and significantly associated with superior cingulum AD. Taken together, these findings suggest that the relationship between prenatal cytokine exposure and white matter microstructure differs as a function of SES. These patterns are consistent with a scenario where gestational inflammation's effects on white matter development diverge depending on the availability of foundational resources in utero.

Gaps and complex structurally variant loci in phased genome assemblies.

There has been tremendous progress in phased genome assembly production by combining long-read data with parental information or linked-read data. Nevertheless, a typical phased genome assembly generated by trio-hifiasm still generates more than 140 gaps. We perform a detailed analysis of gaps, assembly breaks, and misorientations from 182 haploid assemblies obtained from a diversity panel of 77 unique human samples. Although trio-based approaches using HiFi are the current gold standard, chromosome-wide phasing accuracy is comparable when using Strand-seq instead of parental data. Importantly, the majority of assembly gaps cluster near the largest and most identical repeats (including segmental duplications [35.4%], satellite DNA [22.3%], or regions enriched in GA/AT-rich DNA [27.4%]). Consequently, 1513 protein-coding genes overlap assembly gaps in at least one haplotype, and 231 are recurrently disrupted or missing from five or more haplotypes. Furthermore, we estimate that 6-7 Mbp of DNA are misorientated per haplotype irrespective of whether trio-free or trio-based approaches are used. Of these misorientations, 81% correspond to bona fide large inversion polymorphisms in the human species, most of which are flanked by large segmental duplications. We also identify large-scale alignment discontinuities consistent with 11.9 Mbp of deletions and 161.4 Mbp of insertions per haploid genome. Although 99% of this variation corresponds to satellite DNA, we identify 230 regions of euchromatic DNA with frequent expansions and contractions, nearly half of which overlap with 197 protein-coding genes. Such variable and incompletely assembled regions are important targets for future algorithmic development and pangenome representation.

Inversion polymorphism in a complete human genome assembly.

The telomere-to-telomere (T2T) complete human reference has significantly improved our ability to characterize genome structural variation. To understand its impact on inversion polymorphisms, we remapped data from 41 genomes against the T2T reference genome and compared it to the GRCh38 reference. We find a ~ 21% increase in sensitivity improving mapping of 63 inversions on the T2T reference. We identify 26 misorientations within GRCh38 and show that the T2T reference is three times more likely to represent the correct orientation of the major human allele. Analysis of 10 additional samples reveals novel rare inversions at chromosomes 15q25.2, 16p11.2, 16q22.1-23.1, and 22q11.21.

Age-related differences in resting-state functional connectivity from childhood to adolescence.

The human brain is active at rest, and spontaneous fluctuations in functional MRI BOLD signals reveal an intrinsic functional architecture. During childhood and adolescence, functional networks undergo varying patterns of maturation, and measures of functional connectivity within and between networks differ as a function of age. However, many aspects of these developmental patterns (e.g. trajectory shape and directionality) remain unresolved. In the present study, we characterised age-related differences in within- and between-network resting-state functional connectivity (rsFC) and integration (i.e. participation coefficient, PC) in a large cross-sectional sample of children and adolescents (n = 628) aged 8-21 years from the Lifespan Human Connectome Project in Development. We found evidence for both linear and non-linear differences in cortical, subcortical, and cerebellar rsFC, as well as integration, that varied by age. Additionally, we found that sex moderated the relationship between age and putamen integration where males displayed significant age-related increases in putamen PC compared with females. Taken together, these results provide evidence for complex, non-linear differences in some brain systems during development.

Functional analysis of structural variants in single cells using Strand-seq.

Somatic structural variants (SVs) are widespread in cancer, but their impact on disease evolution is understudied due to a lack of methods to directly characterize their functional consequences. We present a computational method, scNOVA, which uses Strand-seq to perform haplotype-aware integration of SV discovery and molecular phenotyping in single cells by using nucleosome occupancy to infer gene expression as a readout. Application to leukemias and cell lines identifies local effects of copy-balanced rearrangements on gene deregulation, and consequences of SVs on aberrant signaling pathways in subclones. We discovered distinct SV subclones with dysregulated Wnt signaling in a chronic lymphocytic leukemia patient. We further uncovered the consequences of subclonal chromothripsis in T cell acute lymphoblastic leukemia, which revealed c-Myb activation, enrichment of a primitive cell state and informed successful targeting of the subclone in cell culture, using a Notch inhibitor. By directly linking SVs to their functional effects, scNOVA enables systematic single-cell multiomic studies of structural variation in heterogeneous cell populations.

Family income buffers the relationship between childhood adverse experiences and putamen volume.

Adverse experiences and family income in childhood have been associated with altered brain development. While there is a large body of research examining these associations, it has primarily used cross-sectional data sources and studied adverse experiences and family income in isolation. However, it is possible that low family income and adverse experiences represent dissociable and potentially interacting profiles of risk. To address this gap in the literature, we examined brain structure as a function of adverse experiences in childhood and family income in 158 youths with up to five waves of MRI data. Specifically, we assessed the interactive effect of these two risk factors on six regions of interest: hippocampus, putamen, amygdala, nucleus accumbens, caudate, and thalamus. Adverse experiences and family income interacted to predict putamen volume (B = 0.086, p = 0.011) but only in participants with family income one standard deviation below the mean (slope estimate = -0.11, p = 0.03). These results suggest that adverse experiences in childhood result in distinct patterns of brain development across the socioeconomic gradient. Given previous findings implicating the role of the putamen in psychopathology-related behaviors, these results emphasize the importance of considering life events and socioeconomic context when evaluating markers of risk. Future research should include interactive effects of environmental exposures and family income to better characterize risk for psychopathology in diverse samples.

Structural Connectivity and Emotion Recognition Impairment in Children and Adolescents with Chromosome 22q11.2 Deletion Syndrome.

Children with chromosome 22q11.2 deletion syndrome (22q11.2DS) exhibit impaired ability to process and understand emotions in others. We measured structural connectivity in children and adolescents with 22q11.2DS (n = 28) and healthy controls (n = 29). Compared to controls, those with 22q11.2DS had poorer social skills and more difficulty recognizing facial emotions. Children with 22q11.2DS also had higher fractional anisotropic diffusion in right amygdala to fusiform gyrus white matter pathways. Right amygdala to fusiform gyrus fractional anisotropy values partially mediated the relationship between 22q11.2DS and social skills, as well as the relationship between 22q11.2DS and emotion recognition accuracy. These findings provide insight into the neural origins of social skills deficits seen in 22q11.2DS and may serve as a biomarker for risk of future psychiatric problems.

White matter alterations associated with lifetime and current depression in adolescents: Evidence for cingulum disruptions.

INTRODUCTION: Compared to research on adults with depression, relatively little work has examined white matter microstructure differences in depression arising earlier in life. Here we tested hypotheses about disruptions to white matter structure in adolescents with current and past depression, with an a priori focus on the cingulum bundles, uncinate fasciculi, corpus collosum, and superior longitudinal fasciculus. METHODS: One hundred thirty-one children from the Preschool Depression Study were assessed using a Human Connectome Project style diffusion imaging sequence which was processed with HCP pipelines and TRACULA to generate estimates of fractional anisotropy (FA), mean diffusivity (MD), axial diffusivity (AD) and radial diffusivity (RD). RESULTS: We found that reduced FA, reduced AD, and increased RD in the dorsal cingulum bundle were associated with a lifetime diagnosis of major depression and greater cumulative and current depression severity. Reduced FA, reduced AD, and increased RD in the ventral cingulum were associated with greater cumulative depression severity. CONCLUSION: These findings support the emergence of white matter differences detected in adolescence associated with earlier life and concurrent depression. They also highlight the importance of connections of the cingulate to other brain regions in association with depression, potentially relevant to understanding emotion dysregulation and functional connectivity differences in depression.

A high-resolution map of small-scale inversions in the gibbon genome.

Gibbons are the most speciose family of living apes, characterized by a diverse chromosome number and rapid rate of large-scale rearrangements. Here we performed single-cell template strand sequencing (Strand-seq), molecular cytogenetics, and deep in silico analysis of a southern white-cheeked gibbon genome, providing the first comprehensive map of 238 previously hidden small-scale inversions. We determined that more than half are gibbon specific, at least fivefold higher than shown for other primate lineage-specific inversions, with a significantly high number of small heterozygous inversions, suggesting that accelerated evolution of inversions may have played a role in the high sympatric diversity of gibbons. Although the precise mechanisms underlying these inversions are not yet understood, it is clear that segmental duplication-mediated NAHR only accounts for a small fraction of events. Several genomic features, including gene density and repeat (e.g., LINE-1) content, might render these regions more break-prone and susceptible to inversion formation. In the attempt to characterize interspecific variation between southern and northern white-cheeked gibbons, we identify several large assembly errors in the current GGSC Nleu3.0/nomLeu3 reference genome comprising more than 49 megabases of DNA. Finally, we provide a list of 182 candidate genes potentially involved in gibbon diversification and speciation.

Store-operated calcium entry controls innate and adaptive immune cell function in inflammatory bowel disease.

Inflammatory bowel disease (IBD) is characterized by dysregulated intestinal immune responses. Using mass cytometry (CyTOF) to analyze the immune cell composition in the lamina propria (LP) of patients with ulcerative colitis (UC) and Crohn's disease (CD), we observed an enrichment of CD4+ effector T cells producing IL-17A and TNF, CD8+ T cells producing IFNγ, T regulatory (Treg) cells, and innate lymphoid cells (ILC). The function of these immune cells is regulated by store-operated Ca2+ entry (SOCE), which results from the opening of Ca2+ release-activated Ca2+ (CRAC) channels formed by ORAI and STIM proteins. We observed that the pharmacologic inhibition of SOCE attenuated the production of proinflammatory cytokines including IL-2, IL-4, IL-6, IL-17A, TNF, and IFNγ by human colonic T cells and ILCs, reduced the production of IL-6 by B cells and the production of IFNγ by myeloid cells, but had no effect on the viability, differentiation, and function of intestinal epithelial cells. T cell-specific deletion of CRAC channel genes in mice showed that Orai1, Stim1, and Stim2-deficient T cells have quantitatively distinct defects in SOCE, which correlate with gradually more pronounced impairment of cytokine production by Th1 and Th17 cells and the severity of IBD. Moreover, the pharmacologic inhibition of SOCE with a selective CRAC channel inhibitor attenuated IBD severity and colitogenic T cell function in mice. Our data indicate that SOCE inhibition may be a suitable new approach for the treatment of IBD.

Developmental trajectories of cortical thickness by functional brain network: The roles of pubertal timing and socioeconomic status.

The human cerebral cortex undergoes considerable changes during development, with cortical maturation patterns reflecting regional heterogeneity that generally progresses in a posterior-to-anterior fashion. However, the organizing principles that govern cortical development remain unclear. In the current study, we characterized age-related differences in cortical thickness (CT) as a function of sex, pubertal timing, and two dissociable indices of socioeconomic status (i.e., income-to-needs and maternal education) in the context of functional brain network organization, using a cross-sectional sample (n = 789) diverse in race, ethnicity, and socioeconomic status from the Lifespan Human Connectome Project in Development (HCP-D). We found that CT generally followed a linear decline from 5 to 21 years of age, except for three functional networks that displayed nonlinear trajectories. We found no main effect of sex or age by sex interaction for any network. Earlier pubertal timing was associated with reduced mean CT and CT in seven networks. We also found a significant age by maternal education interaction for mean CT across cortex and CT in the dorsal attention network, where higher levels of maternal education were associated with steeper age-related decreases in CT. Taken together, our results suggest that these biological and environmental variations may impact the emerging functional connectome.

The pathology fellowship application crisis: The current state and suggestions for remediation.

Problems within the Pathology fellowship application process in the US have been recognized and reported for years. Recently, members of the Graduate Medical Education Committee (GMEC) of the Association of Pathology Chairs (APC) and collaborators collected survey data from the residents themselves and the fellowship programs, as represented by both the fellowship program directors (members of the Fellowship Directors Ad Hoc Committee, FDAHC) and the program administrators (members of the Graduate Medical Education Administrators Section, GMEAS). These data are presented and discussed, and potential steps to resolve some of the problems around fellowship applications in pathology are presented.

Graded Variation in T1w/T2w Ratio during Adolescence: Measurement, Caveats, and Implications for Development of Cortical Myelin.

Adolescence is characterized by the maturation of cortical microstructure and connectivity supporting complex cognition and behavior. Axonal myelination influences brain connectivity during development by enhancing neural signaling speed and inhibiting plasticity. However, the maturational timing of cortical myelination during human adolescence remains poorly understood. Here, we take advantage of recent advances in high-resolution cortical T1w/T2w mapping methods, including principled correction of B1+ transmit field effects, using data from the Human Connectome Project in Development (HCP-D; N = 628, ages 8-21). We characterize microstructural changes relevant to myelination by estimating age-related differences in T1w/T2w throughout the cerebral neocortex from childhood to early adulthood. We apply Bayesian spline models and clustering analysis to demonstrate graded variation in age-dependent cortical T1w/T2w differences that are correlated with the sensorimotor-association (S-A) axis of cortical organization reported by others. In sensorimotor areas, T1w/T2w ratio measures start at high levels at early ages, increase at a fast pace, and decelerate at later ages (18-21). In intermediate multimodal areas along the S-A axis, T1w/T2w starts at intermediate levels and increases linearly at an intermediate pace. In transmodal/paralimbic association areas, T1w/T2w starts at low levels and increases linearly at the slowest pace. These data provide evidence for graded variation of the T1w/T2w ratio along the S-A axis that may reflect cortical myelination changes during adolescence underlying the development of complex information processing and psychological functioning. We discuss the implications of these results as well as caveats in interpreting magnetic resonance imaging (MRI)-based estimates of myelination.SIGNIFICANCE STATEMENT Myelin is a lipid membrane that is essential to healthy brain function. Myelin wraps axons to increase neural signaling speed, enabling complex neuronal functioning underlying learning and cognition. Here, we characterize the developmental timing of myelination across the cerebral cortex during adolescence using a noninvasive proxy measure, T1w/T2w mapping. Our results provide new evidence demonstrating graded variation across the cortex in the timing of T1w/T2w changes during adolescence, with rapid T1w/T2w increases in lower-order sensory areas and gradual T1w/T2w increases in higher-order association areas. This spatial pattern of microstructural brain development closely parallels the sensorimotor-to-association axis of cortical organization and plasticity during ontogeny.

Recurrent inversion polymorphisms in humans associate with genetic instability and genomic disorders.

Unlike copy number variants (CNVs), inversions remain an underexplored genetic variation class. By integrating multiple genomic technologies, we discover 729 inversions in 41 human genomes. Approximately 85% of inversions <2 kbp form by twin-priming during L1 retrotransposition; 80% of the larger inversions are balanced and affect twice as many nucleotides as CNVs. Balanced inversions show an excess of common variants, and 72% are flanked by segmental duplications (SDs) or retrotransposons. Since flanking repeats promote non-allelic homologous recombination, we developed complementary approaches to identify recurrent inversion formation. We describe 40 recurrent inversions encompassing 0.6% of the genome, showing inversion rates up to 2.7 × 10-4 per locus per generation. Recurrent inversions exhibit a sex-chromosomal bias and co-localize with genomic disorder critical regions. We propose that inversion recurrence results in an elevated number of heterozygous carriers and structural SD diversity, which increases mutability in the population and predisposes specific haplotypes to disease-causing CNVs.

Familial long-read sequencing increases yield of de novo mutations.

Studies of de novo mutation (DNM) have typically excluded some of the most repetitive and complex regions of the genome because these regions cannot be unambiguously mapped with short-read sequencing data. To better understand the genome-wide pattern of DNM, we generated long-read sequence data from an autism parent-child quad with an affected female where no pathogenic variant had been discovered in short-read Illumina sequence data. We deeply sequenced all four individuals by using three sequencing platforms (Illumina, Oxford Nanopore, and Pacific Biosciences) and three complementary technologies (Strand-seq, optical mapping, and 10X Genomics). Using long-read sequencing, we initially discovered and validated 171 DNMs across two children-a 20% increase in the number of de novo single-nucleotide variants (SNVs) and indels when compared to short-read callsets. The number of DNMs further increased by 5% when considering a more complete human reference (T2T-CHM13) because of the recovery of events in regions absent from GRCh38 (e.g., three DNMs in heterochromatic satellites). In total, we validated 195 de novo germline mutations and 23 potential post-zygotic mosaic mutations across both children; the overall true substitution rate based on this integrated callset is at least 1.41 × 10-8 substitutions per nucleotide per generation. We also identified six de novo insertions and deletions in tandem repeats, two of which represent structural variants. We demonstrate that long-read sequencing and assembly, especially when combined with a more complete reference genome, increases the number of DNMs by >25% compared to previous studies, providing a more complete catalog of DNM compared to short-read data alone.

Early Childhood Socioeconomic Status and Cognitive and Adaptive Outcomes at the Transition to Adulthood: The Mediating Role of Gray Matter Development Across Five Scan Waves.

BACKGROUND: Early low socioeconomic status (SES) is associated with poor outcomes in childhood, many of which endure into adulthood. It is critical to determine how early low SES relates to trajectories of brain development and whether these mediate relationships to poor outcomes. We use data from a unique 17-year longitudinal study with five waves of structural brain imaging to prospectively examine relationships between preschool SES and cognitive, social, academic, and psychiatric outcomes in early adulthood. METHODS: Children (n = 216, 50% female, 47.2% non-White) were recruited from a study of early onset depression and followed approximately annually. Family income-to-needs ratios (SES) were assessed when children were ages 3 to 5 years. Volumes of cortical gray and white matter and subcortical gray matter collected across five scan waves were processed using the FreeSurfer Longitudinal pipeline. When youth were ages 16+ years, cognitive function was assessed using the NIH Toolbox, and psychiatric diagnoses, high-risk behaviors, educational function, and social function were assessed using clinician administered and parent/youth report measures. RESULTS: Lower preschool SES related to worse cognitive, high-risk, educational, and social outcomes (|standardized B| = 0.20-0.31, p values < .003). Lower SES was associated with overall lower cortical (standardized B = 0.12, p < .0001) and subcortical gray matter (standardized B = 0.17, p < .0001) volumes, as well as a shallower slope of subcortical gray matter growth over time (standardized B = 0.04, p = .012). Subcortical gray matter mediated the relationship of preschool SES to cognition and high-risk behaviors. CONCLUSIONS: These novel longitudinal data underscore the key role of brain development in understanding the long-lasting relations of early low SES to outcomes in children.

Construction of Whole Genomes from Scaffolds Using Single Cell Strand-Seq Data.

Accurate reference genome sequences provide the foundation for modern molecular biology and genomics as the interpretation of sequence data to study evolution, gene expression, and epigenetics depends heavily on the quality of the genome assembly used for its alignment. Correctly organising sequenced fragments such as contigs and scaffolds in relation to each other is a critical and often challenging step in the construction of robust genome references. We previously identified misoriented regions in the mouse and human reference assemblies using Strand-seq, a single cell sequencing technique that preserves DNA directionality Here we demonstrate the ability of Strand-seq to build and correct full-length chromosomes by identifying which scaffolds belong to the same chromosome and determining their correct order and orientation, without the need for overlapping sequences. We demonstrate that Strand-seq exquisitely maps assembly fragments into large related groups and chromosome-sized clusters without using new assembly data. Using template strand inheritance as a bi-allelic marker, we employ genetic mapping principles to cluster scaffolds that are derived from the same chromosome and order them within the chromosome based solely on directionality of DNA strand inheritance. We prove the utility of our approach by generating improved genome assemblies for several model organisms including the ferret, pig, Xenopus, zebrafish, Tasmanian devil and the Guinea pig.

ASHLEYS: automated quality control for single-cell Strand-seq data.

SUMMARY: Single-cell DNA template strand sequencing (Strand-seq) enables chromosome length haplotype phasing, construction of phased assemblies, mapping sister-chromatid exchange events and structural variant discovery. The initial quality control of potentially thousands of single-cell libraries is still done manually by domain experts. ASHLEYS automates this tedious task, delivers near-expert performance and labels even large datasets in seconds. AVAILABILITY AND IMPLEMENTATION: github.com/friendsofstrandseq/ashleys-qc, MIT license. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

A feasibility trial of parent HPV vaccine reminders and phone-based motivational interviewing.

BACKGROUND: We assessed the feasibility and acceptability of a sequential approach of parent-targeted HPV vaccine reminders and phone-based Motivation Interviewing (MI). METHODS: In 2016, we selected all 11- to 12-year-old boys and girls seen in one clinic whose vaccine records did not include the HPV vaccine (n=286). By gender, we individually randomized parents of adolescents to an interactive text message (74 girls and 45 boys), postcard reminder (46 boys and no girls because of previously demonstrated efficacy), or standard care group (75 girls and 46 boys). Reminders were sent with medical director permission and a HIPAA waiver. Two months after reminders, among the adolescents whose vaccine records still did not include the HPV vaccine, we selected a gender-stratified random sample of 20 parents for phone-based MI. We assessed the percentage of deliverable messages, the percentage of parents' responding to the interactive text message, parent acceptability of receiving a text message, and MI parent responsiveness and interviewer competence (MI Treatment Integrity Coding system). RESULTS: Nearly all messages were deliverable (98% of postcards and 74% of text messages). Six of the 88 parents (7%) receiving text messages scheduled an appointment through our interactive system. The acceptability survey response rate was 37% (38/102). Respondents were favorable toward vaccine reminders for all parents (82%). Among 20 sampled parents, 17 were reached by phone of whom 7 completed MI, 4 had or were getting the HPV vaccine for their child, and 5 expressed disinterest. Across the 7 MI calls, the interviewer was rated 100% MI adherent and scored an average 4.19 rating for Global Spirit. CONCLUSION: Without providing explicit consent to receive vaccine-related messages, parents nonetheless found postcards and interactive text messages acceptable. Centralizing MI to phone calls with trained staff was acceptable to parents and resulted in highly MI-adherent interviews.