Complex trait associations in rare diseases and impacts on Mendelian variant interpretation.

Emerging evidence implicates common genetic variation - aggregated into polygenic scores (PGS) - in the onset and phenotypic presentation of rare diseases. Here, we comprehensively map individual polygenic liability for 1102 open-source PGS in a cohort of 3059 probands enrolled in the Genomic Answers for Kids (GA4K) rare disease study, revealing widespread associations between rare disease phenotypes and PGSs for common complex diseases and traits, blood protein levels, and brain and other organ morphological measurements. Using this resource, we demonstrate increased polygenic liability in probands with an inherited candidate disease variant (VUS) compared to unaffected carrier parents. Further, we show an enrichment for large-effect rare variants in putative core PGS genes for associated complex traits. Overall, our study supports and expands on previous findings of complex trait associations in rare diseases, implicates polygenic liability as a potential mechanism underlying variable penetrance of candidate causal variants, and provides a framework for identifying novel candidate rare disease genes.

Mapping structural variants to rare disease genes using long-read whole genome sequencing and trait-relevant polygenic scores.

Recent studies have revealed the pervasive landscape of rare structural variants (rSVs) present in human genomes. rSVs can have extreme effects on the expression of proximal genes and, in a rare disease context, have been implicated in patient cases where no diagnostic single nucleotide variant (SNV) was found. Approaches for integrating rSVs to date have focused on targeted approaches in known Mendelian rare disease genes. This approach is intractable for rare diseases with many causal loci or patients with complex, multi-phenotype syndromes. We hypothesized that integrating trait-relevant polygenic scores (PGS) would provide a substantial reduction in the number of candidate disease genes in which to assess rSV effects. We further implemented a method for ranking PGS genes to define a set of core/key genes where a rSV has the potential to exert relatively larger effects on disease risk. Among a subset of patients enrolled in the Genomic Answers for Kids (GA4K) rare disease program (N=497), we used PacBio HiFi long-read whole genome sequencing (lrWGS) to identify rSVs intersecting genes in trait-relevant PGSs. Illustrating our approach in Autism (N=54 cases), we identified 22, 019 deletions, 2,041 duplications, 87,826 insertions, and 214 inversions overlapping putative core/key PGS genes. Additionally, by integrating genomic constraint annotations from gnomAD, we observed that rare duplications overlapping putative core/key PGS genes were frequently in higher constraint regions compared to controls (P = 1×10-03). This difference was not observed in the lowest-ranked gene set (P = 0.15). Overall, our study provides a framework for the annotation of long-read rSVs from lrWGS data and prioritization of disease-linked genomic regions for downstream functional validation of rSV impacts. To enable reuse by other researchers, we have made SV allele frequencies and gene associations freely available.

Complex trait associations in rare diseases and impacts on Mendelian variant interpretation.

Emerging evidence implicates common genetic variation - aggregated into polygenic scores (PGS) - impacting the onset and phenotypic presentation of rare diseases. In this study, we quantified individual polygenic liability for 1,151 previously published PGS in a cohort of 2,374 probands enrolled in the Genomic Answers for Kids (GA4K) rare disease study, revealing widespread associations between rare disease phenotypes and PGSs for common complex diseases and traits, blood protein levels, and brain and other organ morphological measurements. We observed increased polygenic burden in probands with variants of unknown significance (VUS) compared to unaffected carrier parents. We further observed an enrichment in overlap between diagnostic and candidate rare disease genes and large-effect PGS genes. Overall, our study supports and expands on previous findings of complex trait associations in rare disease phenotypes and provides a framework for identifying novel candidate rare disease genes and in understanding variable penetrance of candidate Mendelian disease variants.

Pangenome graphs improve the analysis of structural variants in rare genetic diseases.

Rare DNA alterations that cause heritable diseases are only partially resolvable by clinical next-generation sequencing due to the difficulty of detecting structural variation (SV) in all genomic contexts. Long-read, high fidelity genome sequencing (HiFi-GS) detects SVs with increased sensitivity and enables assembling personal and graph genomes. We leverage standard reference genomes, public assemblies (n = 94) and a large collection of HiFi-GS data from a rare disease program (Genomic Answers for Kids, GA4K, n = 574 assemblies) to build a graph genome representing a unified SV callset in GA4K, identify common variation and prioritize SVs that are more likely to cause genetic disease (MAF < 0.01). Using graphs, we obtain a higher level of reproducibility than the standard reference approach. We observe over 200,000 SV alleles unique to GA4K, including nearly 1000 rare variants that impact coding sequence. With improved specificity for rare SVs, we isolate 30 candidate SVs in phenotypically prioritized genes, including known disease SVs. We isolate a novel diagnostic SV in KMT2E, demonstrating use of personal assemblies coupled with pangenome graphs for rare disease genomics. The community may interrogate our pangenome with additional assemblies to discover new SVs within the allele frequency spectrum relevant to genetic diseases.

Intestinal lesions and parasites associated with senescence and prespawn mortality in Chinook Salmon (Oncorhynchus tshawytscha).

Prespawn mortality (PSM) presents a major problem for the recovery of spring Chinook Salmon (Oncorhynchus tshawytscha) populations. In the Willamette River, Oregon, PSM exceeds 90% in some years but factors explaining it are not well understood. We examined intestinal tissue samples using histological slides from over 783 spring Chinook Salmon collected between 2009 and 2021, which included tissues from PSM fish, artificially spawned captive broodstock (BS) and normal river run fish, comprised of trapped (Live) and naturally post-spawned river (RPS) fish collected from the river. We observed degeneration of the intestinal epithelium and loss of villous structure, with concurrent severe enteritis. A natural progression of decline in epithelial integrity (EI) through the summer and fall until spawning and subsequent death was also observed. Live fish exhibited high EI scores (mean = 68%), BS exhibited variable EI scores (35%) and RPS exhibited severe loss of EI (14%). PSM fish exhibited prominent loss of intestinal epithelium with EI scores (13%), very similar to RPS fish, despite having been collected earlier in the year. Hence, we argue that low EI scores are strongly linked with PSM. Ceratonova shasta and Enterocytozoon schreckii were common in all groups, but neither were linked to either PSM or a decline in EI.

Aeromonas salmonicida, causative agent of salmonid furunculosis, isolated from the freshwater parasitic copepod, Salmincola californiensis.

Here, we provide evidence that the freshwater parasitic copepod, Salmincola californiensis, acts as a vector for Aeromonas salmonicida. While investigating the effects of S. californiensis on Chinoook salmon (Oncorhynchus tshawytscha), we tangentially observed that fish infected with the copepod developed furunculosis, caused by A. salmonicida. This occurred despite being reared in pathogen-free well water in a research facility with no prior history of spontaneous infection. We further investigated the possibility of S. californiensis to serve as a vector for the bacterium via detection of fluorescently labelled A. salmonicida inside the egg sacs from copepods in which the fish hosts were experimentally infected with GFP-A449 A. salmonicida. We then evaluated copepod egg sacs that were collected from adult Chinook salmon from a freshwater hatchery with A. salmonicida infections confirmed by either culture or PCR. The bacterium was cultured on tryptic soy agar plates from 75% of the egg sacs, and 61% were positive by PCR. These three separate experiments indicate an alternative tactic of transmission in addition to direct transmission of A. salmonicida in captivity. The copepod may play an important role in transmission of the bacterium when fish are more dispersed, such as in the wild.

Extravillous trophoblast cell lineage development is associated with active remodeling of the chromatin landscape.

The extravillous trophoblast cell lineage is a key feature of placentation and successful pregnancy. Knowledge of transcriptional regulation driving extravillous trophoblast cell development is limited. Here, we map the transcriptome and epigenome landscape as well as chromatin interactions of human trophoblast stem cells and their transition into extravillous trophoblast cells. We show that integrating chromatin accessibility, long-range chromatin interactions, transcriptomic, and transcription factor binding motif enrichment enables identification of transcription factors and regulatory mechanisms critical for extravillous trophoblast cell development. We elucidate functional roles for TFAP2C, SNAI1, and EPAS1 in the regulation of extravillous trophoblast cell development. EPAS1 is identified as an upstream regulator of key extravillous trophoblast cell transcription factors, including ASCL2 and SNAI1 and together with its target genes, is linked to pregnancy loss and birth weight. Collectively, we reveal activation of a dynamic regulatory network and provide a framework for understanding extravillous trophoblast cell specification in trophoblast cell lineage development and human placentation.

Direct haplotype-resolved 5-base HiFi sequencing for genome-wide profiling of hypermethylation outliers in a rare disease cohort.

Long-read HiFi genome sequencing allows for accurate detection and direct phasing of single nucleotide variants, indels, and structural variants. Recent algorithmic development enables simultaneous detection of CpG methylation for analysis of regulatory element activity directly in HiFi reads. We present a comprehensive haplotype resolved 5-base HiFi genome sequencing dataset from a rare disease cohort of 276 samples in 152 families to identify rare (~0.5%) hypermethylation events. We find that 80% of these events are allele-specific and predicted to cause loss of regulatory element activity. We demonstrate heritability of extreme hypermethylation including rare cis variants associated with short (~200 bp) and large hypermethylation events (>1 kb), respectively. We identify repeat expansions in proximal promoters predicting allelic gene silencing via hypermethylation and demonstrate allelic transcriptional events downstream. On average 30-40 rare hypermethylation tiles overlap rare disease genes per patient, providing indications for variation prioritization including a previously undiagnosed pathogenic allele in DIP2B causing global developmental delay. We propose that use of HiFi genome sequencing in unsolved rare disease cases will allow detection of unconventional diseases alleles due to loss of regulatory element activity.

IRF7 and UNC93B1 variants in an infant with recurrent herpes simplex virus infection.

Neonatal herpes simplex virus (HSV) infection is a devastating disease with substantial morbidity and mortality. The genetic basis of susceptibility to HSV in neonates remains undefined. We evaluated a male infant with neonatal skin/eye/mouth (SEM) HSV-1 disease, who had complete recovery after acyclovir but developed HSV-1 encephalitis at 1 year of age. An immune workup showed an anergic PBMC cytokine response to TLR3 stimulation but no other TLRs. Exome sequencing identified rare missense variants in IFN-regulatory factor 7 (IRF7) and UNC-93 homolog B1 (UNC93B1). PBMC single-cell RNA-Seq done during childhood revealed decreased expression of several innate immune genes and a repressed TLR3 pathway signature at baseline in several immune cell populations, including CD14 monocytes. Functional studies in fibroblasts and human leukemia monocytic THP1 cells showed that both variants individually suppressed TLR3-driven IRF3 transcriptional activity and the type I IFN response in vitro. Furthermore, fibroblasts expressing the IRF7 and UNC93B1 variants had higher intracellular viral titers with blunting of the type I IFN response upon HSV-1 challenge. This study reports an infant with recurrent HSV-1 disease complicated by encephalitis associated with deleterious variants in the IRF7 and UNC93B1 genes. Our results suggest that TLR3 pathway mutations may predispose neonates to recurrent, severe HSV.

Impact of crowding on the diversity of expanding populations.

Crowding effects critically impact the self-organization of densely packed cellular assemblies, such as biofilms, solid tumors, and developing tissues. When cells grow and divide, they push each other apart, remodeling the structure and extent of the population's range. Recent work has shown that crowding has a strong impact on the strength of natural selection. However, the impact of crowding on neutral processes, which controls the fate of new variants as long as they are rare, remains unclear. Here, we quantify the genetic diversity of expanding microbial colonies and uncover signatures of crowding in the site frequency spectrum. By combining Luria-Delbrück fluctuation tests, lineage tracing in a novel microfluidic incubator, cell-based simulations, and theoretical modeling, we find that the majority of mutations arise behind the expanding frontier, giving rise to clones that are mechanically "pushed out" of the growing region by the proliferating cells in front. These excluded-volume interactions result in a clone-size distribution that solely depends on where the mutation first arose relative to the front and is characterized by a simple power law for low-frequency clones. Our model predicts that the distribution depends on a single parameter-the characteristic growth layer thickness-and hence allows estimation of the mutation rate in a variety of crowded cellular populations. Combined with previous studies on high-frequency mutations, our finding provides a unified picture of the genetic diversity in expanding populations over the whole frequency range and suggests a practical method to assess growth dynamics by sequencing populations across spatial scales.

Gut microbiome composition associates with corticosteroid treatment, morbidity, and senescence in Chinook salmon (Oncorhynchus tshawytscha).

Pacific salmon experience prolonged elevation in corticosteroid hormones during important life history events including migration, reproduction, and senescence. These periods of elevated corticosteroids correspond with changes to immunity and energy metabolism; therefore, fish may be particularly vulnerable to mortality at these times. Recent studies found that stress-induced cortisol release associated with microbial community shifts in salmonids, raising the question of how longer-term corticosteroid dynamics that accompany life history transitions affect salmonid microbiomes. In this work, we experimentally evaluated the relationships between gut microbiome composition, chronically elevated corticosteroids, and mortality in juvenile Chinook salmon (Oncorhynchus tshawytscha). We found that treatment with slow-release implants of the corticosteroids cortisol or dexamethasone resulted in changes to the gut microbiome. Morbidity was also associated with microbiome composition, suggesting that the gut microbiome reflects individual differences in susceptibility to opportunistic pathogens. Additionally, we analyzed a small number of samples from adult fish at various stages of senescence. Results from these samples suggest that microbiome composition associated with gut integrity, and that the microbial communities of corticosteroid treated juveniles shift in composition toward those of senescent adults. Overall, findings from this work indicate that the gut microbiome correlates with mortality risk during periods of chronic corticosteroid elevation.

Earlier onset of North Atlantic hurricane season with warming oceans.

Numerous Atlantic basin tropical cyclones have recently developed prior to the official start of hurricane season, including several pre-season landfalls in the continental United States. Pre-season and early-season tropical cyclones disproportionately affect populated landmasses, often producing outsized precipitation impacts. Here we show a significant trend towards earlier onset of tropical cyclone activity in the North Atlantic basin, with threshold dates of the first three percentiles of accumulated cyclone energy shifting earlier at a rate exceeding five days decade-1 since 1979, even correcting for biases in climatology due to increased detection of short-lived storms. Initial threshold dates of continental United States named storm landfalls have trended earlier by two days decade-1 since 1900. The trend towards additional pre-season and early-season activity is linked to spring thermodynamic conditions becoming more conducive for tropical cyclone formation. Genesis potential index value increases in the western Atlantic basin are primarily driven by warming ocean temperatures.

Bacterial Filamentation Drives Colony Chirality.

Chirality is ubiquitous in nature, with consequences at the cellular and tissue scales. As Escherichia coli colonies expand radially, an orthogonal component of growth creates a pinwheel-like pattern that can be revealed by fluorescent markers. To elucidate the mechanistic basis of this colony chirality, we investigated its link to left-handed, single-cell twisting during E. coli elongation. While chemical and genetic manipulation of cell width altered single-cell twisting handedness, colonies ceased to be chiral rather than switching handedness, and anaerobic growth altered colony chirality without affecting single-cell twisting. Chiral angle increased with increasing temperature even when growth rate decreased. Unifying these findings, we discovered that colony chirality was associated with the propensity for cell filamentation. Inhibition of cell division accentuated chirality under aerobic growth and generated chirality under anaerobic growth. Thus, regulation of cell division is intrinsically coupled to colony chirality, providing a mechanism for tuning macroscale spatial patterning. IMPORTANCE Chiral objects, such as amino acids, are distinguishable from their mirror image. For living systems, the fundamental mechanisms relating cellular handedness to chirality at the multicellular scale remain largely mysterious. Here, we use chemical, genetic, and environmental perturbations of Escherichia coli to investigate whether pinwheel patterns in bacterial colonies are directly linked to single-cell growth behaviors. We discover that chirality can be abolished without affecting single-cell twisting; instead, the degree of chirality was linked to the proportion of highly elongated cells at the colony edge. Inhibiting cell division boosted the degree of chirality during aerobic growth and even introduced chirality to otherwise achiral colonies during anaerobic growth. These findings reveal a fascinating connection between cell division and macroscopic colony patterning.

An assessment of terminology for intraspecific diversity in fishes, with a focus on "ecotypes" and "life histories".

Understanding and preserving intraspecific diversity (ISD) is important for species conservation. However, ISD units do not have taxonomic standards and are not universally recognized. The terminology used to describe ISD is varied and often used ambiguously. We compared definitions of terms used to describe ISD with use in recent studies of three fish taxa: sticklebacks (Gasterosteidae), Pacific salmon and trout (Oncorhynchus spp., "PST"), and lampreys (Petromyzontiformes). Life history describes the phenotypic responses of organisms to environments and includes biological parameters that affect population growth or decline. Life-history pathway(s) are the result of different organismal routes of development that can result in different life histories. These terms can be used to describe recognizable life-history traits. Life history is generally used in organismal- and ecology-based journals. The terms paired species/species pairs have been used to describe two different phenotypes, whereas in some species and situations a continuum of phenotypes may be expressed. Our review revealed overlapping definitions for race and subspecies, and subspecies and ecotypes. Ecotypes are genotypic adaptations to particular environments, and this term is often used in genetic- and evolution-based journals. "Satellite species" is used for situations in which a parasitic lamprey yields two or more derived, nonparasitic lamprey species. Designatable Units, Evolutionary Significant Units (ESUs), and Distinct Population Segments (DPS) are used by some governments to classify ISD of vertebrate species within distinct and evolutionary significant criteria. In situations where the genetic or life-history components of ISD are not well understood, a conservative approach would be to call them phenotypes.

Laboratory infection rates and associated mortality of juvenile Chinook Salmon (Oncorhynchus tshawytscha) from parasitic copepod (Salmincola californiensis).

Pacific salmon (Oncorhynchus spp.) rearing in lakes and reservoirs above dams have been known to become heavily infected with an ectoparasitic copepod (Salmincola californiensis). Little is known about the factors that affect the parasite infection prevalence and intensity. However, previous research suggests that the parasite may negatively affect the fitness and survival of the host fish. The effect of water temperature, confinement and the density of the free-swimming infectious stage of S. californiensis, the copepodid, on infection prevalence and intensity was evaluated by experimentally exposing juvenile Chinook Salmon (O. tshawytscha). Infection rates observed in wild populations were achieved under warm water (15-16°C) and high copepodid density (150-300/L) treatment conditions. Infection prevalence and intensity were also significantly higher in larger fish. During the infection experiment, 4.5% of infected fish died within 54 days with mortality significantly related to copepod infection intensity. The potential for autoinfection was compared to cross-infection by cohabitation of infected fish with naïve fish. Previously infected fish had significantly greater infection intensity compared with naïve fish, indicating that infected fish can be reinfected and that they may be more susceptible than naïve fish.

Biomechanical Feedback Strengthens Jammed Cellular Packings.

Growth in confined spaces can drive cellular populations through a jamming transition from a fluidlike state to a solidlike state. Experiments have found that jammed budding yeast populations can build up extreme compressive pressures (over 1 MPa), which in turn feed back onto cellular physiology by slowing or even stalling cell growth. Using numerical simulations, we investigate how this feedback impacts the mechanical properties of model jammed cell populations. We find that feedback directs growth toward poorly coordinated regions, resulting in an excess number of cell-cell contacts that rigidify cell packings. Cell packings possess anomalously large shear and bulk moduli that depend sensitively on the strength of feedback. These results demonstrate that mechanical feedback on the single-cell level is a simple mechanism by which living systems may tune their population-level mechanical properties.

Rearing environment affects spatial learning in juvenile Chinook salmon Oncorhynchus tshawytscha.

We tested the prediction that a complex physical rearing environment would enhance short-term spatial memory as assessed by learning ability in a spatial navigation task in juvenile Chinook salmon Oncorhynchus tshawytscha. We reared fish in two low-density treatments, where fish were either in bare fiberglass tanks (bare) or in tanks with physical structure (complex). We also tested conventionally reared high-density hatchery fish to compare with these other experimental treatments. Our reason for including this third hatchery treatment is that the two low-density treatments, aside from the manipulation of structure, followed a rearing programme that is designed to produce fish with more wild-like characteristics. We tested individually marked fish for seven consecutive days and recorded movement and time to exit a testing maze. Stimulus conspecific fish outside the exit of the maze provided positive reinforcement for test fish. Fish from the bare treatment were less likely to exit the start box compared with fish in the complex and hatchery treatments. However, fish in the hatchery treatment were significantly more likely to exit the maze on their own compared with both the bare and complex treatments. Hatchery fish effectively learned the task as shown by a decrease in the number of mistakes over time, but the number of mistakes was significantly greater on the first day of trials. Increasing habitat complexity with structure may not necessarily promote spatial learning ability, but differences between hatchery and experimental treatments in rearing density and motivation to be near conspecifics likely led to observed behavioural differences.

PRECIPITATION AND FATAL MOTOR VEHICLE CRASHES: Continental Analysis with High-Resolution Radar Data.

Precipitation, even at light intensity, contributes a significant risk of fatal motor vehicle crashes across the United States, at nearly all times of day, and in all seasons.

Collective motion conceals fitness differences in crowded cellular populations.

Many cellular populations are tightly packed, such as microbial colonies and biofilms, or tissues and tumours in multicellular organisms. The movement of one cell in these crowded assemblages requires motion of others, so that cell displacements are correlated over many cell diameters. Whenever movement is important for survival or growth, these correlated rearrangements could couple the evolutionary fate of different lineages. However, little is known about the interplay between mechanical forces and evolution in dense cellular populations. Here, by tracking slower-growing clones at the expanding edge of yeast colonies, we show that the collective motion of cells prevents costly mutations from being weeded out rapidly. Joint pushing by neighbouring cells generates correlated movements that suppress the differential displacements required for selection to act. This mechanical screening of fitness differences allows slower-growing mutants to leave more descendants than expected under non-mechanical models, thereby increasing their chance for evolutionary rescue. Our work suggests that, in crowded populations, cells cooperate with surrounding neighbours through inevitable mechanical interactions. This effect has to be considered when predicting evolutionary outcomes, such as the emergence of drug resistance or cancer evolution.

Egg size and growth in steelhead Oncorhynchus mykiss.

The effects of egg size on early development and growth of steelhead Oncorhynchus mykiss were recorded for more than 200 days following hatching. Fish from smaller eggs hatched sooner and at a smaller size than fish from larger eggs, but fish from smaller eggs showed consistently higher growth rates than fish from larger eggs. Since many life-history attributes appear to be determined by size or growth rate at age during the first year, egg size could be a significant predictor of important changes in the life history of individuals.