A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3

Modulators of MAPK pathway activity during filamentous growth in Saccharomyces cerevisiae

Mitogen-activated protein kinase (MAPK) pathways control the response to intrinsic and extrinsic stimuli. In the budding yeast Saccharomyces cerevisiae, cells undergo filamentous growth, which is regulated by the fMAPK pathway. To better understand the regulation of the fMAPK pathway, a genetic screen was performed to identify spontaneous mutants with elevated activity of an fMAPK-pathway dependent growth reporter (ste4 FUS1-HIS3). In total, 159 mutants were isolated and analyzed by secondary screens for invasive growth by the plate-washing assay, and filament formation by microscopy. Thirty-two mutants were selected for whole-genome sequencing, which identified new alleles in genes encoding known regulators of the fMAPK pathway. These included gain-of-function alleles in STE11, which encodes the MAPKKK, as well as loss-of-function alleles in KSS1, which encodes the MAP kinase, and RGA1, which encodes a GTPase activating protein (GAP) for CDC42. New alleles in previously identified pathway modulators were also uncovered in ALY1, AIM44, RCK2, IRA2, REG1 and in genes that regulate protein folding (KAR2), glycosylation (MNN4), and turnover (BLM10). C-terminal truncations in the transcription factor Ste12p were also uncovered that resulted in elevated reporter activity, presumably identifying an inhibitory domain in the C-terminus of the protein. We also show that a wide variety of filamentous growth phenotypes result from mutations in different regulators of the response. The alleles identified here expand the connections surrounding MAPK pathway regulation and reveal new features of proteins that function in the signaling cascade.

Profiling the genome and proteome of metabolic dysfunction-associated steatotic liver disease identifies potential therapeutic targets

BACKGROUND & AIMS

Metabolic dysfunction-associated steatotic liver disease (MASLD) affects over 25% of the population and currently has no effective treatments. Plasma proteins with causal evidence may represent promising drug targets. We aimed to identify plasma proteins in the causal pathway of MASLD and explore their interaction with obesity.

METHODS

We analysed 2,941 plasma proteins in 43,978 European participants from UK Biobank. We performed genome-wide association study (GWAS) for all MASLD-associated proteins and created the largest MASLD GWAS (109,885 cases/1,014,923 controls). We performed Mendelian Randomization (MR) and integrated proteins and their encoding genes in MASLD ranges to identify candidate causal proteins. We then validated them through independent replication, exome sequencing, liver imaging, bulk and single-cell gene expression, liver biopsies, pathway, and phenome-wide data. We explored the role of obesity by MR and multivariable MR across proteins, body mass index, and MASLD.

RESULTS

We found 929 proteins associated with MASLD, reported five novel genetic loci associated with MASLD, and identified 17 candidate MASLD protein targets. We identified four novel targets for MASLD (CD33, GRHPR, HMOX2, and SCG3), provided protein evidence supporting roles of AHCY, FCGR2B, ORM1, and RBKS in MASLD, and validated nine previously known targets. We found that CD33, FCGR2B, ORM1, RBKS, and SCG3 mediated the association of obesity and MASLD, and HMOX2, ORM1, and RBKS had effect on MASLD independent of obesity.

CONCLUSIONS

This study identified new protein targets in the causal pathway of MASLD, providing new insights into the multi-omics architecture and pathophysiology of MASLD. These findings advise further therapeutic interventions for MASLD.

Transcriptional immune suppression and upregulation of double stranded DNA damage and repair repertoires in ecDNA-containing tumors

Extrachromosomal DNA is a common cause of oncogene amplification in cancer. The non-chromosomal inheritance of ecDNA enables tumors to rapidly evolve, contributing to treatment resistance and poor outcome for patients. The transcriptional context in which ecDNAs arise and progress, including chromosomally-driven transcription, is incompletely understood. We examined gene expression patterns of 870 tumors of varied histological types, to identify transcriptional correlates of ecDNA. Here we show that ecDNA containing tumors impact four major biological processes. Specifically, ecDNA containing tumors upregulate DNA damage and repair, cell cycle control, and mitotic processes, but downregulate global immune regulation pathways. Taken together, these results suggest profound alterations in gene regulation in ecDNA containing tumors, shedding light on molecular processes that give rise to their development and progression.

A familial missense variant in the Alzheimer's Disease gene SORL1 impairs its maturation and endosomal sorting

The SORL1 gene has recently emerged as a strong Alzheimer's Disease (AD) risk gene. Over 500 different variants have been identified in the gene and the contribution of individual variants to AD development and progression is still largely unknown. Here, we describe a family consisting of 2 parents and 5 offspring. Both parents were affected with dementia and one had confirmed AD pathology with an age of onset >75 years. All offspring were affected with AD with ages at onset ranging from 53yrs-74yrs. DNA was available from the parent with confirmed AD and 5 offspring. We identified a coding variant, p.(Arg953Cys), in SORL1 in 5 of 6 individuals affected by AD. Notably, variant carriers had severe AD pathology, and the SORL1 variant segregated with TDP-43 pathology (LATE-NC). We further characterized this variant and show that this Arginine substitution occurs at a critical position in the YWTD-domain of the SORL1 translation product, SORL1. Functional studies further show that the p.R953C variant leads to retention of the SORL1 protein in the endoplasmic reticulum which leads to decreased maturation and shedding of the receptor and prevents its normal endosomal trafficking. Together, our analysis suggests that p.R953C is a pathogenic variant of SORL1 and sheds light on mechanisms of how missense SORL1 variants may lead to AD.

Distinguishing recrudescence from reinfection in lymphatic filariasis: a genomics-based approach for monitoring worm burden

Background

The Global Program to Eliminate Lymphatic Filariasis is the largest public health program based on mass drug administration (MDA). Despite decades of MDA, ongoing transmission in some countries remains a challenge. To optimize interventions, it is essential to differentiate between recrudescence (poor drug response and persistent infection) and new infections (ongoing transmission). Since adult filariae are inaccessible in humans, an approach that relies on genotyping the offspring microfilariae (mf) is required.

Methods

We utilized Brugia malayi adults and mf obtained from gerbils with a known pedigree to develop and validate our whole-genome amplification and kinship analysis approach. We then sequenced the genomes of Wuchereria bancrofti mf from infected humans from Côte d'Ivoire (CDI), characterized the population genetic diversity, and made inferences about the adult breeders. We developed a whole-exome capture panel for W. bancrofti to enrich parasite nuclear DNA from lower-quality samples contaminated with host DNA.

Results

We established a robust analysis pipeline using B. malayi adult and mf. We estimated the pre-treatment genetic diversity in W. bancrofti from 269 mf collected from 18 individuals, and further analyzed 1-year post-treatment samples of 74 mf from 4 individuals. By reconstructing and temporally tracking sibling relationships across pre- and post-treatment samples, we differentiated between new and established maternal families, suggesting reinfection in one subject and recrudescence in three subjects. Estimated reproductively active adult females ranged between 3 and 9 in the studied subjects. Hemizygosity of the male X-chromosome allowed for direct inference of haplotypes, facilitating robust maternal parentage inference, even when the genetic diversity was low. Population structure analysis revealed genetically distinct parasites among our CDI samples. Sequence composition and variant analysis of whole-exome libraries showed that the hybridization capture approach can effectively enrich parasite nuclear DNA and identify protein-coding variants with ∼95% genotype concordance rate.

Conclusions

We have generated resources to facilitate development of field-deployable genotyping tools that can estimate worm burdens and monitor parasite populations. These tools are essential for the success of lymphatic filariasis MDA programs. With further expansion of the databases to include geographically diverse samples, we will be able to spatially track parasite movement associated with host/vector migration.

Hypoxia-inducible factor induces cysteine dioxygenase and promotes cysteine homeostasis in Caenorhabditis elegans

Dedicated genetic pathways regulate cysteine homeostasis. For example, high levels of cysteine activate cysteine dioxygenase, a key enzyme in cysteine catabolism in most animal and many fungal species. The mechanism by which cysteine dioxygenase is regulated is largely unknown. In an unbiased genetic screen for mutations that activate cysteine dioxygenase (cdo-1) in the nematode C. elegans, we isolated loss-of-function mutations in rhy-1 and egl-9, which encode proteins that negatively regulate the stability or activity of the oxygen-sensing hypoxia inducible transcription factor (hif-1). EGL-9 and HIF-1 are core members of the conserved eukaryotic hypoxia response. However, we demonstrate that the mechanism of HIF-1-mediated induction of cdo-1 is largely independent of EGL-9 prolyl hydroxylase activity and the von Hippel-Lindau E3 ubiquitin ligase, the classical hypoxia signaling pathway components. We demonstrate that C. elegans cdo-1 is transcriptionally activated by high levels of cysteine and hif-1. hif-1-dependent activation of cdo-1 occurs downstream of an H2S-sensing pathway that includes rhy-1, cysl-1, and egl-9. cdo-1 transcription is primarily activated in the hypodermis where it is also sufficient to drive sulfur amino acid metabolism. Thus, the regulation of cdo-1 by hif-1 reveals a negative feedback loop that maintains cysteine homeostasis. High levels of cysteine stimulate the production of an H2S signal. H2S then acts through the rhy-1/cysl-1/egl-9 signaling pathway to increase HIF-1-mediated transcription of cdo-1, promoting degradation of cysteine via CDO-1.

Females translate male mRNA transferred during mating

Although RNA is found in the seminal fluid of diverse organisms, it is unknown whether this RNA is functional within females. Here, we develop an experimental proteomic method called VESPA (Variant Enabled SILAC Proteomic Analysis) to test the hypothesis that Drosophila male seminal fluid RNA is translated by females. We find strong evidence for 67 male-derived, female-translated proteins (mdFTPs) in female lower reproductive tracts at six hours postmating, many with predicted functions relevant to reproduction. Gene knockout experiments indicate that genes coding for mdFTPs play diverse roles in postmating interactions, with effects on fertilization efficiency, and the formation and persistence of the insemination reaction mass, a trait hypothesized to be involved in sexual conflict. These findings advance our understanding of reproduction by revealing a novel mechanism of postmating molecular interactions between the sexes that strengthens and extends male influences on reproductive outcomes in previously unrecognized ways. Given the diverse species known to carry RNA in seminal fluid, this discovery has broad significance for understanding molecular mechanisms of cooperation and conflict during reproduction.

A mutation in Themis contributes to peanut-induced oral anaphylaxis in CC027 mice

Background

The development of peanut allergy is due to a combination of genetic and environmental factors, although specific genes have proven difficult to identify. Previously, we reported that peanut-sensitized CC027/GeniUnc (CC027) mice develop anaphylaxis upon oral challenge to peanut, unlike C3H/HeJ (C3H) mice.

Objective

To determine the genetic basis of orally-induced anaphylaxis to peanut in CC027 mice.

Methods

A genetic mapping population between CC027 and C3H mice was designed to identify the genetic factors that drive oral anaphylaxis. A total of 356 CC027xC3H backcrossed mice were generated, sensitized to peanut, then challenged to peanut by oral gavage. Anaphylaxis and peanut-specific IgE were quantified for all mice. T-cell phenotyping was conducted on CC027 and five additional CC strains.

Results

Anaphylaxis to peanut was absent in 77% of backcrossed mice, with 19% showing moderate anaphylaxis, and 4% having severe anaphylaxis. A total of eight genetic loci were associated with variation in response to peanut challenge, six associated with anaphylaxis (temperature decrease) and two associated with peanut-specific IgE levels. There were two major loci that impacted multiple aspects of the severity of acute anaphylaxis, at which the CC027 allele was associated with worse outcome. At one of these loci, CC027 has a private genetic variant in the Themis (thymocyte-expressed molecule involved in selection) gene. Consistent with Themis' described functions, we found that CC027 have more immature T cells with fewer CD8+, CD4+, and CD4+CD25+CD127- regulatory T cells.

Conclusion

Our results demonstrate a key role for Themis in the orally-reactive CC027 mouse model of peanut allergy.

Cross-species transmission and PB2 mammalian adaptations of highly pathogenic avian influenza A/H5N1 viruses in Chile

H5N1 highly pathogenic avian influenza viruses (HPAIV) emerged in wild birds in Chile in December 2022 and spilled over into poultry, marine mammals, and one human. Between December 9, 2022 - March 14, 2023, a coordinated government/academic response detected HPAIV by real-time RT-PCR in 8.5% (412/4735) of samples from 23 avian and 3 mammal orders. Whole-genome sequences obtained from 77 birds and 8 marine mammals revealed that all Chilean H5N1 viruses belong to lineage 2.3.4.4b and cluster monophyletically with viruses from Peru, indicating a single introduction from North America into Peru/Chile. Mammalian adaptations were identified in the PB2 segment: D701N in two sea lions, one human, and one shorebird, and Q591K in the human and one sea lion. Minor variant analysis revealed that D701N was present in 52.9 - 70.9% of sequence reads, indicating the presence of both genotypes within hosts. Further surveillance of spillover events is warranted to assess the emergence and potential onward transmission of mammalian adapted H5N1 HPAIV in South America.

Characterizing common and rare variations in non-traditional glycemic biomarkers using multivariate approaches on multi-ancestry ARIC study

Glycated hemoglobin, fasting glucose, glycated albumin, and fructosamine are biomarkers that reflect different aspects of the glycemic process. Genetic studies of these glycemic biomarkers can shed light on unknown aspects of type 2 diabetes genetics and biology. While there exists several GWAS of glycated hemoglobin and fasting glucose, very few GWAS have focused on glycated albumin or fructosamine. We performed a multi-phenotype GWAS of glycated albumin and fructosamine from 7,395 White and 2,016 Black participants in the Atherosclerosis Risk in Communities (ARIC) study on the common variants from genotyped/imputed data. We found 2 genome-wide significant loci, one mapping to known type 2 diabetes gene (ARAP1/STARD10, p = 2.8 × 10-8) and another mapping to a novel gene (UGT1A, p = 1.4 × 10-8) using multi-omics gene mapping strategies in diabetes-relevant tissues. We identified additional loci that were ancestry-specific (e.g., PRKCA from African ancestry individuals, p = 1.7 × 10-8) and sex-specific (TEX29 locus in males only, p = 3.0 × 10-8). Further, we implemented multi-phenotype gene-burden tests on whole-exome sequence data from 6,590 White and 2,309 Black ARIC participants. Eleven genes across different rare variant aggregation strategies were exome-wide significant only in multi-ancestry analysis. Four out of 11 genes had notable enrichment of rare predicted loss of function variants in African ancestry participants despite smaller sample size. Overall, 8 out of 15 loci/genes were implicated to influence these biomarkers via glycemic pathways. This study illustrates improved locus discovery and potential effector gene discovery by leveraging joint patterns of related biomarkers across entire allele frequency spectrum in multi-ancestry analyses. Most of the loci/genes we identified have not been previously implicated in studies of type 2 diabetes, and future investigation of the loci/genes potentially acting through glycemic pathways may help us better understand risk of developing type 2 diabetes.

Identifying Aging and Alzheimer Disease-Associated Somatic Variations in Excitatory Neurons From the Human Frontal Cortex

Background and Objectives

With age, somatic mutations accumulated in human brain cells can lead to various neurologic disorders and brain tumors. Because the incidence rate of Alzheimer disease (AD) increases exponentially with age, investigating the association between AD and the accumulation of somatic mutation can help understand the etiology of AD.

Methods

We designed a somatic mutation detection workflow by contrasting genotypes derived from whole-genome sequencing (WGS) data with genotypes derived from scRNA-seq data and applied this workflow to 76 participants from the Religious Order Study and the Rush Memory and Aging Project (ROSMAP) cohort. We focused only on excitatory neurons, the dominant cell type in the scRNA-seq data.

Results

We identified 196 sites that harbored at least 1 individual with an excitatory neuron-specific somatic mutation (ENSM), and these 196 sites were mapped to 127 genes. The single base substitution (SBS) pattern of the putative ENSMs was best explained by signature SBS5 from the Catalogue of Somatic Mutations in Cancer (COSMIC) mutational signatures, a clock-like pattern correlating with the age of the individual. The count of ENSMs per individual also showed an increasing trend with age. Among the mutated sites, we found 2 sites tend to have more mutations in older individuals (16:6899517 [RBFOX1], p = 0.04; 4:21788463 [KCNIP4], p < 0.05). In addition, 2 sites were found to have a higher odds ratio to detect a somatic mutation in AD samples (6:73374221 [KCNQ5], p = 0.01 and 13:36667102 [DCLK1], p = 0.02). Thirty-two genes that harbor somatic mutations unique to AD and the KCNQ5 and DCLK1 genes were used for gene ontology (GO)-term enrichment analysis. We found the AD-specific ENSMs enriched in the GO-term "vocalization behavior" and "intraspecies interaction between organisms." Of interest we observed both age-specific and AD-specific ENSMs enriched in the K⁺ channel-associated genes.

Discussion

Our results show that combining scRNA-seq and WGS data can successfully detect putative somatic mutations. The putative somatic mutations detected from ROSMAP data set have provided new insights into the association of AD and aging with brain somatic mutagenesis.

Towards robust clinical genome interpretation: developing a consistent terminology to characterize disease-gene relationships - allelic requirement, inheritance modes and disease mechanisms

PURPOSE

The terminology used for gene-disease curation and variant annotation to describe inheritance, allelic requirement, and both sequence and functional consequences of a variant is currently not standardized. There is considerable discrepancy in the literature and across clinical variant reporting in the derivation and application of terms. Here we standardize the terminology for the characterization of disease-gene relationships to facilitate harmonized global curation, and to support variant classification within the ACMG/AMP framework.

METHODS

Terminology for inheritance, allelic requirement, and both structural and functional consequences of a variant used by Gene Curation Coalition (GenCC) members and partner organizations was collated and reviewed. Harmonized terminology with definitions and use examples was created, reviewed, and validated.

RESULTS

We present a standardized terminology to describe gene-disease relationships, and to support variant annotation. We demonstrate application of the terminology for classification of variation in the ACMG SF 2.0 genes recommended for reporting of secondary findings. Consensus terms were agreed and formalized in both sequence ontology (SO) and human phenotype ontology (HPO) ontologies. GenCC member groups intend to use or map to these terms in their respective resources.

CONCLUSION

The terminology standardization presented here will improve harmonization, facilitate the pooling of curation datasets across international curation efforts and, in turn, improve consistency in variant classification and genetic test interpretation.

The impact of damaging epilepsy and cardiac genetic variant burden in sudden death in the young

Background

Sudden unexpected death in children is a tragic event. Understanding the genetics of sudden death in the young (SDY) enables family counseling and cascade screening. The objective of this study was to characterize genetic variation in an SDY cohort using whole genome sequencing.

Methods

The SDY Case Registry is a National Institutes of Health/Centers for Disease Control surveillance effort to discern the prevalence, causes, and risk factors for SDY. The SDY Case Registry prospectively collected clinical data and DNA biospecimens from SDY cases <20 years of age. SDY cases were collected from medical examiner and coroner offices spanning 13 US jurisdictions from 2015-2019. The cohort included 211 children (mean age 1 year; range 0-20 years), determined to have died suddenly and unexpectedly and in whom DNA biospecimens and next-of-kin consent were ascertained. A control cohort consisted of 211 randomly sampled, sex-and ancestry-matched individuals from the 1000 Genomes Project. Genetic variation was evaluated in epilepsy, cardiomyopathy and arrhythmia genes in the SDY and control cohorts. American College of Medical Genetics/Genomics guidelines were used to classify variants as pathogenic or likely pathogenic. Additionally, genetic variation predicted to be damaging was identified using a Bayesian-based artificial intelligence (AI) tool.

Results

The SDY cohort was 42% European, 30% African, 17% Hispanic, and 11% with mixed ancestries, and 39% female. Six percent of the cohort was found to harbor a pathogenic or likely pathogenic genetic variant in an epilepsy, cardiomyopathy or arrhythmia gene. The genomes of SDY cases, but not controls, were enriched for rare, damaging variants in epilepsy, cardiomyopathy and arrhythmia-related genes. A greater number of rare epilepsy genetic variants correlated with younger age at death.

Conclusions

While damaging cardiomyopathy and arrhythmia genes are recognized contributors to SDY, we also observed an enrichment in epilepsy-related genes in the SDY cohort, and a correlation between rare epilepsy variation and younger age at death. These findings emphasize the importance of considering epilepsy genes when evaluating SDY.

Integrated metabolic and genetic analysis reveals distinct features of primary differentiated thyroid cancer and its metastatic potential in humans

Differentiated thyroid cancer (DTC) affects thousands of lives worldwide every year. Typically, DTC is a treatable disease with a good prognosis. Yet, some patients are subjected to partial or total thyroidectomy and radioiodine therapy to prevent local disease recurrence and metastasis. Unfortunately, thyroidectomy and/or radioiodine therapy often worsen(s) the quality of life and might be unnecessary in indolent DTC cases. This clinical setting highlights the unmet need for a precise molecular diagnosis of DTC, which should dictate appropriate therapy. Here we propose a differential multi-omics model approach to distinguish normal gland from thyroid tumor and to indicate potential metastatic diseases in papillary thyroid cancer (PTC), a sub-class of DTC. Based on PTC patient samples, our data suggest that elevated nuclear and mitochondrial DNA mutational burden, intratumor heterogeneity, shortened telomere length, and altered metabolic profile reflect the potential for metastatic disease. Specifically, normal and tumor thyroid tissues from these patients had a distinct yet well-defined metabolic profile with high levels of anabolic metabolites and/or other metabolites associated with the energy maintenance of tumor cells. Altogether, this work indicates that a differential and integrated multi-omics approach might improve DTC management, perhaps preventing unnecessary thyroid gland removal and/or radioiodine therapy. Well-designed, prospective translational clinical trials will ultimately show the value of this targeted molecular approach.

TRANSLATIONAL RELEVANCE

In this article, we propose a new integrated metabolic, genomic, and cytopathologic methods to diagnose Differentiated Thyroid Cancer when the conventional methods failed. Moreover, we suggest metabolic and genomic markers to help predict high-risk Papillary Thyroid Cancer. Both might be important tools to avoid unnecessary surgery and/or radioiodine therapy that can worsen the quality of life of the patients more than living with an indolent Thyroid nodule.

Germline pathogenic variants in 786 neuroblastoma patients

Importance

Neuroblastoma accounts for 12% of childhood cancer deaths. The genetic contribution of rare pathogenic germline variation in patients without a family history remains unclear.

Objective

To define the prevalence, spectrum, and clinical significance of pathogenic germline variation in cancer predisposition genes (CPGs) in neuroblastoma patients.

Design Setting and Participants

Germline DNA sequencing was performed on the peripheral blood from 786 neuroblastoma patients unselected for family history. Rare variants mapping to CPGs were evaluated for pathogenicity and the percentage of cases harboring pathogenic (P) or likely pathogenic (LP) variants was quantified. The frequency of CPG P-LP variants in neuroblastoma cases was compared to two distinct cancer-free control cohorts to assess enrichment. Matched tumor DNA sequencing was evaluated for "second hits" at CPGs and germline DNA array data from 5,585 neuroblastoma cases and 23,505 cancer-free control children was analyzed to identify rare germline copy number variants (CNVs) affecting genes with an excess burden of P-LP variants in neuroblastoma. Neuroblastoma patients with germline P-LP variants were compared to those without P-LP variants to test for association with clinical characteristics, tumor features, and patient survival.

Main Outcomes and Measures

Rare variant prevalence, pathogenicity, enrichment, and association with clinical characteristics, tumor features, and patient survival.

Results

We observed 116 P-LP variants in CPGs involving 13.9% (109/786) of patients, representing a significant excess burden of P-LP variants compared to controls (9.1%; P = 5.14 × 10-5, Odds Ratio: 1.60, 95% confidence interval: 1.27-2.00). BARD1 harbored the most significant burden of P-LP variants compared to controls (1.0% vs. 0.03%; P = 8.18 × 10-7; Odds Ratio: 32.30, 95% confidence interval: 6.44-310.35). Rare germline CNVs disrupting BARD1 were also identified in neuroblastoma patients (0.05%) but absent in controls (P = 7.08 × 10-3; Odds Ratio: 29.47, 95% confidence interval: 1.52 - 570.70). Overall, P-LP variants in DNA repair genes in this study were enriched in cases compared to controls (8.1% vs. 5.7%; P = 0.01; Odds Ratio: 1.45, 95% confidence interval: 1.08-1.92). Neuroblastoma patients harboring a germline P-LP variant had a worse overall survival when compared to patients without P-LP variants (P = 8.6 × 10-3), and this remained significant in a multivariate Cox proportional-hazards model (P = 0.01).

Conclusions and Relevance

Neuroblastoma patients harboring germline P-LP variants in CPGs have worse overall survival and BARD1 is an important predisposition gene affected by both common and rare pathogenic variation. Germline sequencing should be performed for all neuroblastoma patients at diagnosis to inform genetic counseling and support future longitudinal and mechanistic studies. Patients with a germline P-LP variant should be closely monitored, regardless of risk group assignment.

Cell-type-specific co-expression inference from single cell RNA-sequencing data

The inference of gene co-expressions from microarray and RNA-sequencing data has led to rich insights on biological processes and disease mechanisms. However, the bulk samples analyzed in most studies are a mixture of different cell types. As a result, the inferred co-expressions are confounded by varying cell type compositions across samples and only offer an aggregated view of gene regulations that may be distinct across different cell types. The advancement of single cell RNA-sequencing (scRNA-seq) technology has enabled the direct inference of co-expressions in specific cell types, facilitating our understanding of cell-type-specific biological functions. However, the high sequencing depth variations and measurement errors in scRNA-seq data present significant challenges in inferring cell-type-specific gene co-expressions, and these issues have not been adequately addressed in the existing methods. We propose a statistical approach, CS-CORE, for estimating and testing cell-type-specific co-expressions, built on a general expression-measurement model that explicitly accounts for sequencing depth variations and measurement errors in the observed single cell data. Systematic evaluations show that most existing methods suffer from inflated false positives and biased co-expression estimates and clustering analysis, whereas CS-CORE has appropriate false positive control, unbiased co-expression estimates, good statistical power and satisfactory performance in downstream co-expression analysis. When applied to analyze scRNA-seq data from postmortem brain samples from Alzheimer’s disease patients and controls and blood samples from COVID-19 patients and controls, CS-CORE identified cell-type-specific co-expressions and differential co-expressions that were more reproducible and/or more enriched for relevant biological pathways than those inferred from other methods.

Somatic Variants in SVIL in Cerebral Aneurysms

Background and Objectives

While somatic mutations have been well-studied in cancer, their roles in other complex traits are much less understood. Our goal is to identify somatic variants that may contribute to the formation of saccular cerebral aneurysms.

Methods

We performed whole-exome sequencing on aneurysm tissues and paired peripheral blood. RNA sequencing and the CRISPR/Cas9 system were then used to perform functional validation of our results.

Results

Somatic variants involved in supervillin (SVIL) or its regulation were found in 17% of aneurysm tissues. In the presence of a mutation in theSVILgene, the expression level of SVIL was downregulated in the aneurysm tissue compared with normal control vessels. Downstream signaling pathways that were induced by knockdown ofSVILvia the CRISPR/Cas9 system in vascular smooth muscle cells (vSMCs) were determined by evaluating changes in gene expression and protein kinase phosphorylation. We found thatSVILregulated the phenotypic modulation of vSMCs to the synthetic phenotype via Krüppel-like factor 4 and platelet-derived growth factor and affected cell migration of vSMCs via the RhoA/ROCK pathway.

Discussion

We propose that somatic variants form a novel mechanism for the development of cerebral aneurysms. Specifically, somatic variants inSVILresult in the phenotypic modulation of vSMCs, which increases the susceptibility to aneurysm formation. This finding suggests a new avenue for the therapeutic intervention and prevention of cerebral aneurysms.

Natural Antibodies Alert the Adaptive Immune System of the Presence of Transformed Cells in Early Tumorigenesis

Recent studies have revealed a critical role for natural Abs (NAbs) in antitumor immune responses. However, the role of NAbs in cancer immunosurveillance remains unexplored, mainly because of the lack of in vivo models that mimic the early recognition and elimination of transforming cells. In this article, we propose a role for NAbs in alerting the immune system against precancerous neoantigen-expressing cells immediately after they escape intrinsic tumor suppression mechanisms. We identify four distinct reproducible, trackable, MHC-matched neoantigen-expressing cell models that do not form tumors as the end point. This amplified readout in the critical window prior to tumor formation allows investigation of new mediators of cancer immunosurveillance. We found that neoantigen-expressing cells adoptively transferred in NAb-deficient mice persisted, whereas they were eliminated in wild-type mice, indicating that the circulating NAb repertoire alerts the immune system to the presence of transformed cells. Moreover, immunity is mounted against immunogenic and nonimmunogenic neoantigens contained in the NAb-tagged cells, regardless of whether the NAb directly recognizes the neoantigens. Beyond these neoantigen-expressing model systems, we observed a significantly greater tumor burden in chemically and virally induced tumor models in NAb-deficient mice compared with wild-type mice. Restoration of the NAb repertoire in NAb-deficient mice elicited the recognition and elimination of neoantigen-expressing cells and cancer. These data show that NAbs are required and sufficient for elimination of transformed cells early in tumorigenesis. These models can now be used to investigate how NAbs stimulate immunity via recognition receptors to eliminate precancerous cells.

Genomic epidemiology and emergence of SARS-CoV-2 variants of concern in the United Arab Emirates

Since the declaration of SARS-CoV-2 outbreak as a pandemic, the United Arab Emirates (UAE) public health authorities have adopted strict measures to reduce transmission as early as March 2020. As a result of these measures, flight suspension, nationwide RT-PCR and surveillance of viral sequences were extensively implemented. This study aims to characterize the epidemiology, transmission pattern, and emergence of variants of concerns (VOCs) and variants of interests (VOIs) of SARS-CoV-2 in the UAE, followed by the investigation of mutations associated with hospitalized cases. A total of 1274 samples were collected and sequenced from all seven emirates between the period of 25 April 2020 to 15 February 2021. Phylogenetic analysis demonstrated multiple introductions of SARS-CoV-2 into the UAE in the early pandemic, followed by a local spread of root clades (A, B, B.1 and B.1.1). As the international flight resumed, the frequencies of VOCs surged indicating the January peak of positive cases. We observed that the hospitalized cases were significantly associated with the presence of B.1.1.7 (p < 0.001), B.1.351 (p < 0.001) and A.23.1 (p = 0.009). Deceased cases are more likely to occur in the presence of B.1.351 (p < 0.001) and A.23.1 (p = 0.022). Logistic and ridge regression showed that 51 mutations are significantly associated with hospitalized cases with the highest proportion originated from S and ORF1a genes (31% and 29% respectively). Our study provides an epidemiological insight of the emergence of VOCs and VOIs following the borders reopening and worldwide travels. It provides reassurance that hospitalization is markedly more associated with the presence of VOCs. This study can contribute to understand the global transmission of SARS-CoV-2 variants.

Conservation genetics as a management tool: The five best-supported paradigms to assist the management of threatened species

About 50 y ago, Crow and Kimura [An Introduction to Population Genetics Theory (1970)] and Ohta and Kimura [Genet. Res. 22, 201-204 (1973)] laid the foundations of conservation genetics by predicting the relationship between population size and genetic marker diversity. This work sparked an enormous research effort investigating the importance of population dynamics, in particular small population size, for population mean performance, population viability, and evolutionary potential. In light of a recent perspective [J. C. Teixeira, C. D. Huber, Proc. Natl. Acad. Sci. U.S.A. 118, 10 (2021)] that challenges some fundamental assumptions in conservation genetics, it is timely to summarize what the field has achieved, what robust patterns have emerged, and worthwhile future research directions. We consider theory and methodological breakthroughs that have helped management, and we outline some fundamental and applied challenges for conservation genetics.

Genomic basis of fishing-associated selection varies with population density

Fisheries-associated selection is recognized as one of the strongest potential human drivers of contemporary evolution in natural populations. The results of this study show that while simulated commercial fishing techniques consistently remove fish with traits associated with growth, metabolism, and social behavior, the specific genes under fishing selection differ depending on the density of the targeted population. This finding suggests that different fish populations of varying sizes will respond differently to fishing selection at the genetic level. Furthermore, as a population is fished over time, the genes under selection may change as the population diminishes. This could have repercussions on population resilience. This study highlights the importance of selection but also environmental and density effects on harvested fish populations.

Fisheries induce one of the strongest anthropogenic selective pressures on natural populations, but the genetic effects of fishing remain unclear. Crucially, we lack knowledge of how capture-associated selection and its interaction with reductions in population density caused by fishing can potentially shift which genes are under selection. Using experimental fish reared at two densities and repeatedly harvested by simulated trawling, we show consistent phenotypic selection on growth, metabolism, and social behavior regardless of density. However, the specific genes under selection—mainly related to brain function and neurogenesis—varied with the population density. This interaction between direct fishing selection and density could fundamentally alter the genomic responses to harvest. The evolutionary consequences of fishing are therefore likely context dependent, possibly varying as exploited populations decline. These results highlight the need to consider environmental factors when predicting effects of human-induced selection and evolution.

Parallel genomic responses to historical climate change and high elevation in East Asian songbirds

Parallel evolution can be expected among closely related taxa exposed to similar selective pressures. However, parallelism is typically stronger at the phenotypic level, while genetic solutions to achieve these phenotypic similarities may differ. For polygenic traits, the availability of standing genetic variation (i.e., heterozygosity) may influence such genetic nonparallelism. Here, we examine the extent to which high-elevation adaptation is parallel-and whether the level of parallelism is affected by heterozygosity-by analyzing genomes of 19 Paridae species distributed across East Asia with a dramatic east-west elevation gradient. We find that western highlands endemic parids have consistently lower levels of heterozygosity-likely the result of late-Pleistocene demographic contraction-than do parids found exclusively in eastern lowlands, which remained unglaciated during the late Pleistocene. Three widespread species (east to west) have high levels of heterozygosity similar to that observed in eastern species, although their western populations are less variable than eastern ones. Comparing genomic responses to extreme environments of the Qinghai-Tibet Plateau, we find that the most differentiated genomic regions between each high-elevation taxon and its low-elevation relative are significantly enriched for genes potentially related to the oxygen transport cascade and/or thermogenesis. Despite no parallelism at particular genes, high similarity in gene function is found among comparisons. Furthermore, parallelism is not higher in more heterozygous widespread parids than in highland endemics. Thus, in East Asian parids, parallel functional response to extreme elevation appears to rely on different genes, with differences in heterozygosity having no effect on the degree of genetic parallelism.

Dendritic cell paucity in mismatch repair-proficient colorectal cancer liver metastases limits immune checkpoint blockade efficacy

Liver metastasis is a major cause of mortality for patients with colorectal cancer (CRC). Mismatch repair-proficient (pMMR) CRCs make up about 95% of metastatic CRCs, and are unresponsive to immune checkpoint blockade (ICB) therapy. Here we show that mouse models of orthotopic pMMR CRC liver metastasis accurately recapitulate the inefficacy of ICB therapy in patients, whereas the same pMMR CRC tumors are sensitive to ICB therapy when grown subcutaneously. To reveal local, nonmalignant components that determine CRC sensitivity to treatment, we compared the microenvironments of pMMR CRC cells grown as liver metastases and subcutaneous tumors. We found a paucity of both activated T cells and dendritic cells in ICB-treated orthotopic liver metastases, when compared with their subcutaneous tumor counterparts. Furthermore, treatment with Feline McDonough sarcoma (FMS)-like tyrosine kinase 3 ligand (Flt3L) plus ICB therapy increased dendritic cell infiltration into pMMR CRC liver metastases and improved mouse survival. Lastly, we show that human CRC liver metastases and microsatellite stable (MSS) primary CRC have a similar paucity of T cells and dendritic cells. These studies indicate that orthotopic tumor models, but not subcutaneous models, should be used to guide human clinical trials. Our findings also posit dendritic cells as antitumor components that can increase the efficacy of immunotherapies against pMMR CRC.

SARS-CoV-2 evolution in animals suggests mechanisms for rapid variant selection

SARS-CoV-2 spillback from humans into domestic and wild animals has been well documented, and an accumulating number of studies illustrate that human-to-animal transmission is widespread in cats, mink, deer, and other species. Experimental inoculations of cats, mink, and ferrets have perpetuated transmission cycles. We sequenced full genomes of Vero cell-expanded SARS-CoV-2 inoculum and viruses recovered from cats (n = 6), dogs (n = 3), hamsters (n = 3), and a ferret (n = 1) following experimental exposure. Five nonsynonymous changes relative to the USA-WA1/2020 prototype strain were near fixation in the stock used for inoculation but had reverted to wild-type sequences at these sites in dogs, cats, and hamsters within 1- to 3-d postexposure. A total of 14 emergent variants (six in nonstructural genes, six in spike, and one each in orf8 and nucleocapsid) were detected in viruses recovered from animals. This included substitutions in spike residues H69, N501, and D614, which also vary in human lineages of concern. Even though a live virus was not cultured from dogs, substitutions in replicase genes were detected in amplified sequences. The rapid selection of SARS-CoV-2 variants in vitro and in vivo reveals residues with functional significance during host switching. These observations also illustrate the potential for spillback from animal hosts to accelerate the evolution of new viral lineages, findings of particular concern for dogs and cats living in households with COVID-19 patients. More generally, this glimpse into viral host switching reveals the unrealized rapidity and plasticity of viral evolution in experimental animal model systems.

Dysfunctional Homozygous VRK1-D263G Variant Impairs the Assembly of Cajal Bodies and DNA Damage Response in Hereditary Spastic Paraplegia

Background and Objectives

To conduct a genetic and molecular functional study of a family with members affected of hereditary spastic paraplegia (HSP) of unknown origin and carrying a novel pathogenic vaccinia-related kinase 1 (VRK1) variant.

Methods

Whole-exome sequencing was performed in 2 patients, and their parents diagnosed with HSP. The novel VRK1 variant was detected by whole-exome sequencing, molecularly modeled and biochemically characterized in kinase assays. Functionally, we studied the role of this VRK1 variant in DNA damage response and its effect on the assembly of Cajal bodies (CBs).

Results

We have identified a very rare homozygous variant VRK1-D263G with a neurologic phenotype associated with HSP and moderate intellectual disability. The molecular modeling of this VRK1 variant protein predicted an alteration in the folding of a loop that interferes with the access to the kinase catalytic site. The VRK1-D263G variant is kinase inactive and does not phosphorylate histones H2AX and H3, transcription factors activating transcription factor 2 and p53, coilin needed for assembly of CBs, and p53 binding protein 1, a DNA repair protein. Functionally, this VRK1 variant protein impairs CB formation and the DNA damage response.

Discussion

This report expands the neurologic spectrum of neuromotor syndromes associated with a new and rare VRK1 variant, representing a novel pathogenic participant in complicated HSP and demonstrates that CBs and the DNA damage response are impaired in these patients.

Transmission event of SARS-CoV-2 Delta variant reveals multiple vaccine breakthrough infections

Importance

Vaccine breakthrough by an emergent SARS-CoV-2 variant poses a great risk to global public health.

Objective

To determine the SARS-CoV-2 variant responsible for 6 cases of vaccine breakthrough.

Design

Nasopharyngeal swabs from suspected vaccine breakthrough cases were tested for SARS-CoV-2 by qPCR for Wuhan-Hu1 and Alpha variant. Positive samples were then sequenced by Swift Normalase Amplicon Panels to determine the causal variant.

Setting

Transmission event occurred at events surrounding a wedding outside of Houston, TX. Two patients from India, likely transmitted the Delta variant to other guests.

Participants

Following a positive SARS-CoV-2 qPCR test at a third-party site, six fully vaccinated patients were investigated. Three males and three females ranged from 53 to 69 years old. One patient suffered from diabetes while three others were classified as overweight. No significant other comorbidities were identified. None of the patients had a history of failed vaccination.