Hybrid spectral CT system with clinical rapid kVp-switching x-ray tube and dual-layer detector for improved iodine quantification.

Spectral computed tomography (CT) is a powerful diagnostic tool offering quantitative material decomposition results that enhance clinical imaging by providing physiologic and functional insights. Iodine, a widely used contrast agent, improves visualization in various clinical contexts. However, accurately detecting low-concentration iodine presents challenges in spectral CT systems, particularly crucial for conditions like pancreatic cancer assessment. In this study, we present preliminary results from our hybrid spectral CT instrumentation which includes clinical-grade hardware (rapid kVp-switching x-ray tube, dual-layer detector). This combination expands spectral datasets from two to four channels, wherein we hypothesize improved quantification accuracy for low-dose and low-iodine concentration cases. We modulate the system duty cycle to evaluate its impact on quantification noise and bias. We evaluate iodine quantification performance by comparing two hybrid weighting strategies alongside rapid kVp-switching. This evaluation is performed with a polyamide phantom containing seven iodine inserts ranging from 0.5 to 20 mg/mL. In comparison to alternative methodologies, the maximum separation configuration, incorporating data from both the 80 kVp, low photon energy detector layer and the 140 kVp, high photon energy detector layer produces spectral images containing low quantitative noise and bias. This study presents initial evaluations on a hybrid spectral CT system, leveraging clinical hardware to demonstrate the potential for enhanced precision and sensitivity in spectral imaging. This research holds promise for advancing spectral CT imaging performance across diverse clinical scenarios.

Proteomic approaches for protein kinase substrate identification in Apicomplexa.

Apicomplexa is a phylum of protist parasites, notable for causing life-threatening diseases including malaria, toxoplasmosis, cryptosporidiosis, and babesiosis. Apicomplexan pathogenesis is generally a function of lytic replication, dissemination, persistence, host cell modification, and immune subversion. Decades of research have revealed essential roles for apicomplexan protein kinases in establishing infections and promoting pathogenesis. Protein kinases modify their substrates by phosphorylating serine, threonine, tyrosine, or other residues, resulting in rapid functional changes in the target protein. Post-translational modification by phosphorylation can activate or inhibit a substrate, alter its localization, or promote interactions with other proteins or ligands. Deciphering direct kinase substrates is crucial to understand mechanisms of kinase signaling, yet can be challenging due to the transient nature of kinase phosphorylation and potential for downstream indirect phosphorylation events. However, with recent advances in proteomic approaches, our understanding of kinase function in Apicomplexa has improved dramatically. Here, we discuss methods that have been used to identify kinase substrates in apicomplexan parasites, classifying them into three main categories: i) kinase interactome, ii) indirect phosphoproteomics and iii) direct labeling. We briefly discuss each approach, including their advantages and limitations, and highlight representative examples from the Apicomplexa literature. Finally, we conclude each main category by introducing prospective approaches from other fields that would benefit kinase substrate identification in Apicomplexa.

Multi-ancestry genome-wide association study of kidney cancer identifies 63 susceptibility regions.

Here, in a multi-ancestry genome-wide association study meta-analysis of kidney cancer (29,020 cases and 835,670 controls), we identified 63 susceptibility regions (50 novel) containing 108 independent risk loci. In analyses stratified by subtype, 52 regions (78 loci) were associated with clear cell renal cell carcinoma (RCC) and 6 regions (7 loci) with papillary RCC. Notably, we report a variant common in African ancestry individuals ( rs7629500 ) in the 3' untranslated region of VHL, nearly tripling clear cell RCC risk (odds ratio 2.72, 95% confidence interval 2.23-3.30). In cis-expression quantitative trait locus analyses, 48 variants from 34 regions point toward 83 candidate genes. Enrichment of hypoxia-inducible factor-binding sites underscores the importance of hypoxia-related mechanisms in kidney cancer. Our results advance understanding of the genetic architecture of kidney cancer, provide clues for functional investigation and enable generation of a validated polygenic risk score with an estimated area under the curve of 0.65 (0.74 including risk factors) among European ancestry individuals.

The Toxoplasma monocarboxylate transporters are involved in the metabolism within the apicoplast and are linked to parasite survival.

The apicoplast is a four-membrane plastid found in the apicomplexans, which harbors biosynthesis and organelle housekeeping activities in the matrix. However, the mechanism driving the flux of metabolites, in and out, remains unknown. Here, we used TurboID and genome engineering to identify apicoplast transporters in Toxoplasma gondii. Among the many novel transporters, we show that one pair of apicomplexan monocarboxylate transporters (AMTs) appears to have evolved from a putative host cell that engulfed a red alga. Protein depletion showed that AMT1 and AMT2 are critical for parasite growth. Metabolite analyses supported the notion that AMT1 and AMT2 are associated with biosynthesis of isoprenoids and fatty acids. However, stronger phenotypic defects were observed for AMT2, including in the inability to establish T. gondii parasite virulence in mice. This study clarifies, significantly, the mystery of apicoplast transporter composition and reveals the importance of the pair of AMTs in maintaining the apicoplast activity in apicomplexans.

Single-cell profiling of MC1R-inhibited melanocytes.

The human red hair color (RHC) trait is caused by increased pheomelanin (red-yellow) and reduced eumelanin (black-brown) pigment in skin and hair due to diminished melanocortin 1 receptor (MC1R) function. In addition, individuals harboring the RHC trait are predisposed to melanoma development. While MC1R variants have been established as causative of RHC and are a well-defined risk factor for melanoma, it remains unclear mechanistically why decreased MC1R signaling alters pigmentation and increases melanoma susceptibility. Here, we use single-cell RNA sequencing (scRNA-seq) of melanocytes isolated from RHC mouse models to define a MC1R-inhibited Gene Signature (MiGS) comprising a large set of previously unidentified genes which may be implicated in melanogenesis and oncogenic transformation. We show that one of the candidate MiGS genes, TBX3, a well-known anti-senescence transcription factor implicated in melanoma progression, binds both E-box and T-box elements to regulate genes associated with melanogenesis and senescence bypass. Our results provide key insights into further mechanisms by which melanocytes with reduced MC1R signaling may regulate pigmentation and offer new candidates of study toward understanding how individuals with the RHC phenotype are predisposed to melanoma.

3D genomic features across >50 diverse cell types reveal insights into the genomic architecture of childhood obesity.

The prevalence of childhood obesity is increasing worldwide, along with the associated common comorbidities of type 2 diabetes and cardiovascular disease in later life. Motivated by evidence for a strong genetic component, our prior genome-wide association study (GWAS) efforts for childhood obesity revealed 19 independent signals for the trait; however, the mechanism of action of these loci remains to be elucidated. To molecularly characterize these childhood obesity loci we sought to determine the underlying causal variants and the corresponding effector genes within diverse cellular contexts. Integrating childhood obesity GWAS summary statistics with our existing 3D genomic datasets for 57 human cell types, consisting of high-resolution promoter-focused Capture-C/Hi-C, ATAC-seq, and RNA-seq, we applied stratified LD score regression and calculated the proportion of genome-wide SNP heritability attributable to cell type-specific features, revealing pancreatic alpha cell enrichment as the most statistically significant. Subsequent chromatin contact-based fine-mapping was carried out for genome-wide significant childhood obesity loci and their linkage disequilibrium proxies to implicate effector genes, yielded the most abundant number of candidate variants and target genes at the BDNF, ADCY3, TMEM18 and FTO loci in skeletal muscle myotubes and the pancreatic beta-cell line, EndoC-BH1. One novel implicated effector gene, ALKAL2 - an inflammation-responsive gene in nerve nociceptors - was observed at the key TMEM18 locus across multiple immune cell types. Interestingly, this observation was also supported through colocalization analysis using expression quantitative trait loci (eQTL) derived from the Genotype-Tissue Expression (GTEx) dataset, supporting an inflammatory and neurologic component to the pathogenesis of childhood obesity. Our comprehensive appraisal of 3D genomic datasets generated in a myriad of different cell types provides genomic insights into pediatric obesity pathogenesis.

Apicomplexan phosphodiesterases in cyclic nucleotide turnover: conservation, function, and therapeutic potential.

Apicomplexa encompasses a large number of intracellular parasites infecting a wide range of animals. Cyclic nucleotide signaling is crucial for a variety of apicomplexan life stages and cellular processes. The cyclases and kinases that synthesize and respond to cyclic nucleotides (i.e., 3',5'-cyclic guanosine monophosphate and 3',5'-cyclic adenosine monophosphate) are highly conserved and essential throughout the parasite phylum. Growing evidence indicates that phosphodiesterases (PDEs) are also critical for regulating cyclic nucleotide signaling via cyclic nucleotide hydrolysis. Here, we discuss recent advances in apicomplexan PDE biology and opportunities for therapeutic interventions, with special emphasis on the major human apicomplexan parasite genera Plasmodium, Toxoplasma, Cryptosporidium, and Babesia. In particular, we show a highly flexible repertoire of apicomplexan PDEs associated with a wide range of cellular requirements across parasites and lifecycle stages. Despite this phylogenetic diversity, cellular requirements of apicomplexan PDEs for motility, host cell egress, or invasion are conserved. However, the molecular wiring of associated PDEs is extremely malleable suggesting that PDE diversity and redundancy are key for the optimization of cyclic nucleotide turnover to respond to the various environments encountered by each parasite and life stage. Understanding how apicomplexan PDEs are regulated and integrating multiple signaling systems into a unified response represent an untapped avenue for future exploration.

Loss of MC1R signaling implicates TBX3 in pheomelanogenesis and melanoma predisposition.

The human Red Hair Color (RHC) trait is caused by increased pheomelanin (red-yellow) and reduced eumelanin (black-brown) pigment in skin and hair due to diminished melanocortin 1 receptor (MC1R) function. In addition, individuals harboring the RHC trait are predisposed to melanoma development. While MC1R variants have been established as causative of RHC and are a well-defined risk factor for melanoma, it remains unclear mechanistically why decreased MC1R signaling alters pigmentation and increases melanoma susceptibility. Here, we use single-cell RNA-sequencing (scRNA-seq) of melanocytes isolated from RHC mouse models to reveal a Pheomelanin Gene Signature (PGS) comprising genes implicated in melanogenesis and oncogenic transformation. We show that TBX3, a well-known anti-senescence transcription factor implicated in melanoma progression, is part of the PGS and binds both E-box and T-box elements to regulate genes associated with melanogenesis and senescence bypass. Our results provide key insights into mechanisms by which MC1R signaling regulates pigmentation and how individuals with the RHC phenotype are predisposed to melanoma.

Oncogenic CDK13 mutations impede nuclear RNA surveillance.

RNA surveillance pathways detect and degrade defective transcripts to ensure RNA fidelity. We found that disrupted nuclear RNA surveillance is oncogenic. Cyclin-dependent kinase 13 (CDK13) is mutated in melanoma, and patient-mutated CDK13 accelerates zebrafish melanoma. CDK13 mutation causes aberrant RNA stabilization. CDK13 is required for ZC3H14 phosphorylation, which is necessary and sufficient to promote nuclear RNA degradation. Mutant CDK13 fails to activate nuclear RNA surveillance, causing aberrant protein-coding transcripts to be stabilized and translated. Forced aberrant RNA expression accelerates melanoma in zebrafish. We found recurrent mutations in genes encoding nuclear RNA surveillance components in many malignancies, establishing nuclear RNA surveillance as a tumor-suppressive pathway. Activating nuclear RNA surveillance is crucial to avoid accumulation of aberrant RNAs and their ensuing consequences in development and disease.

Massively parallel reporter assays and variant scoring identified functional variants and target genes for melanoma loci and highlighted cell-type specificity.

The most recent genome-wide association study (GWAS) of cutaneous melanoma identified 54 risk-associated loci, but functional variants and their target genes for most have not been established. Here, we performed massively parallel reporter assays (MPRAs) by using malignant melanoma and normal melanocyte cells and further integrated multi-layer annotation to systematically prioritize functional variants and susceptibility genes from these GWAS loci. Of 1,992 risk-associated variants tested in MPRAs, we identified 285 from 42 loci (78% of the known loci) displaying significant allelic transcriptional activities in either cell type (FDR < 1%). We further characterized MPRA-significant variants by motif prediction, epigenomic annotation, and statistical/functional fine-mapping to create integrative variant scores, which prioritized one to six plausible candidate variants per locus for the 42 loci and nominated a single variant for 43% of these loci. Overlaying the MPRA-significant variants with genome-wide significant expression or methylation quantitative trait loci (eQTLs or meQTLs, respectively) from melanocytes or melanomas identified candidate susceptibility genes for 60% of variants (172 of 285 variants). CRISPRi of top-scoring variants validated their cis-regulatory effect on the eQTL target genes, MAFF (22q13.1) and GPRC5A (12p13.1). Finally, we identified 36 melanoma-specific and 45 melanocyte-specific MPRA-significant variants, a subset of which are linked to cell-type-specific target genes. Analyses of transcription factor availability in MPRA datasets and variant-transcription-factor interaction in eQTL datasets highlighted the roles of transcription factors in cell-type-specific variant functionality. In conclusion, MPRAs along with variant scoring effectively prioritized plausible candidates for most melanoma GWAS loci and highlighted cellular contexts where the susceptibility variants are functional.

ezQTL: A Web Platform for Interactive Visualization and Colocalization of QTLs and GWAS Loci.

Genome-wide association studies (GWAS) have identified thousands of genomic loci associated with complex diseases and traits, including cancer. The vast majority of common trait-associated variants identified via GWAS fall in non-coding regions of the genome, posing a challenge in elucidating the causal variants, genes, and mechanisms involved. Expression quantitative trait locus (eQTL) and other molecular QTL studies have been valuable resources in identifying candidate causal genes from GWAS loci through statistical colocalization methods. While QTL colocalization is becoming a standard analysis in post-GWAS investigation, an easy web tool for users to perform formal colocalization analyses with either user-provided or public GWAS and eQTL datasets has been lacking. Here, we present ezQTL, a web-based bioinformatic application to interactively visualize and analyze genetic association data such as GWAS loci and molecular QTLs under different linkage disequilibrium (LD) patterns (1000 Genomes Project, UK Biobank, or user-provided data). This application allows users to perform data quality control for variants matched between different datasets, LD visualization, and two-trait colocalization analyses using two state-of-the-art methodologies (eCAVIAR and HyPrColoc), including batch processing. ezQTL is a free and publicly available cross-platform web tool, which can be accessed online at https://analysistools.cancer.gov/ezqtl.

Investigating the genetic architecture of eye colour in a Canadian cohort.

Eye color is highly variable in populations with European ancestry, ranging from low to high quantities of melanin in the iris. Polymorphisms in the HERC2/OCA2 locus have the largest effect on eye color in these populations, although other genomic regions also influence eye color. We performed genome-wide association studies of eye color in a Canadian cohort of European ancestry (N = 5,641) and investigated candidate causal variants. We uncovered several candidate causal signals in the HERC2/OCA2 region, whereas other loci likely harbor a single causal signal. We observed colocalization of eye color signals with the expression or methylation profiles of cultured primary melanocytes. Genetic correlations of eye and hair color suggest high genome-wide pleiotropy, but locus-level differences in the genetic architecture of both traits. Overall, we provide a better picture of the polymorphisms underpinning eye color variation, which may be a consequence of specific molecular processes in the iris melanocytes.

Rare germline deleterious variants increase susceptibility for lung cancer.

Although multiple common susceptibility loci for lung cancer (LC) have been identified by genome-wide association studies, they can explain only a small portion of heritability. The etiological contribution of rare deleterious variants (RDVs) to LC risk is not fully characterized and may account for part of the missing heritability. Here, we sequenced the whole exomes of 2777 participants from the Environment and Genetics in Lung cancer Etiology study, a homogenous population including 1461 LC cases and 1316 controls. In single-variant analyses, we identified a new RDV, rs77187983 [EHBP1, odds ratio (OR) = 3.13, 95% confidence interval (CI) = 1.34-7.30, P = 0.008] and replicated two previously reported RDVs, rs11571833 (BRCA2, OR = 2.18; 95% CI = 1.25-3.81, P = 0.006) and rs752672077 (MPZL2, OR = 3.70, 95% CI = 1.04-13.15, P = 0.044). In gene-based analyses, we confirmed BRCA2 (P = 0.007) and ATM (P = 0.014) associations with LC risk and identified TRIB3 (P = 0.009), involved in maintaining genome stability and DNA repair, as a new candidate susceptibility gene. Furthermore, cases were enriched with RDVs in homologous recombination repair [carrier frequency (CF) = 22.9% versus 19.5%, P = 0.017] and Fanconi anemia (CF = 12.5% versus 10.2%, P = 0.036) pathways. Our results were not significant after multiple testing corrections but were enriched in cases versus controls from large scale public biobank resources, including The Cancer Genome Atlas, FinnGen and UK Biobank. Our study identifies novel candidate genes and highlights the importance of RDVs in DNA repair-related genes for LC susceptibility. These findings improve our understanding of LC heritability and may contribute to the development of risk stratification and prevention strategies.

Defining novel causal SNPs and linked phenotypes at melanoma-associated loci.

A number of genomic regions have been associated with melanoma risk through genome-wide association studies; however, the causal variants underlying the majority of these associations remain unknown. Here, we sequenced either the full locus or the functional regions including exons of 19 melanoma-associated loci in 1959 British melanoma cases and 737 controls. Variant filtering followed by Fisher's exact test analyses identified 66 variants associated with melanoma risk. Sequential conditional logistic regression identified the distinct haplotypes on which variants reside, and massively parallel reporter assays provided biological insights into how these variants influence gene function. We performed further analyses to link variants to melanoma risk phenotypes and assessed their association with melanoma-specific survival. Our analyses replicate previously known associations in the melanocortin 1 receptor (MC1R) and tyrosinase (TYR) loci, while identifying novel potentially causal variants at the MTAP/CDKN2A and CASP8 loci. These results improve our understanding of the architecture of melanoma risk and outcome.

Functional Analysis of the Expanded Phosphodiesterase Gene Family in Toxoplasma gondii Tachyzoites.

Toxoplasma motility is both activated and suppressed by 3',5'-cyclic nucleotide signaling. Cyclic GMP (cGMP) signaling through Toxoplasma gondii protein kinase G (TgPKG) activates motility, whereas cyclic AMP (cAMP) signaling through TgPKAc1 inhibits motility. Despite their importance, it remains unclear how cGMP and cAMP levels are maintained in Toxoplasma. Phosphodiesterases (PDEs) are known to inactivate cyclic nucleotides and are highly expanded in the Toxoplasma genome. Here, we analyzed the expression and function of the 18-member TgPDE family in tachyzoites, the virulent life stage of Toxoplasma. We detected the expression of 11 of 18 TgPDEs, confirming prior expression studies. A knockdown screen of the TgPDE family revealed four TgPDEs that contribute to lytic Toxoplasma growth (TgPDE1, TgPDE2, TgPDE5, and TgPDE9). Depletion of TgPDE1 or TgPDE2 caused severe growth defects, prompting further investigation. While TgPDE1 was important for extracellular motility, TgPDE2 was important for host cell invasion, parasite replication, host cell egress, and extracellular motility. TgPDE1 displayed a plasma membrane/cytomembranous distribution, whereas TgPDE2 displayed an endoplasmic reticulum/cytomembranous distribution. Biochemical analysis of TgPDE1 and TgPDE2 purified from Toxoplasma lysates revealed that TgPDE1 hydrolyzes both cGMP and cAMP, whereas TgPDE2 was cAMP specific. Interactome studies of TgPDE1 and TgPDE2 indicated that they do not physically interact with each other or other TgPDEs but may be regulated by kinases and proteases. Our studies have identified TgPDE1 and TgPDE2 as central regulators of tachyzoite cyclic nucleotide levels and enable future studies aimed at determining how these enzymes are regulated and cooperate to control Toxoplasma motility and growth. IMPORTANCE Apicomplexan parasites require motility to actively infect host cells and cause disease. Cyclic nucleotide signaling governs apicomplexan motility, but it is unclear how cyclic nucleotide levels are maintained in these parasites. In search of novel regulators of cyclic nucleotides in the model apicomplexan Toxoplasma, we identified and characterized two catalytically active phosphodiesterases, TgPDE1 and TgPDE2, that are important for Toxoplasma's virulent tachyzoite life cycle. Enzymes that generate, sense, or degrade cyclic nucleotides make attractive targets for therapies aimed at paralyzing and killing apicomplexan parasites.

Integrated Analysis of Coexpression and Exome Sequencing to Prioritize Susceptibility Genes for Familial Cutaneous Melanoma.

The application of whole-exome sequencing has led to the identification of high- and moderate-risk variants that contribute to cutaneous melanoma susceptibility. However, confirming disease-causing variants remains challenging. We applied a gene coexpression network analysis to prioritize the candidate genes identified from whole-exome sequencing of 34 melanoma-prone families, with at least three affected members sequenced per family (N = 119 cases). A coexpression network was constructed from genotype-tissue expression project, skin melanoma from the cancer genome atlas, and primary melanocyte cultures. We performed module-specific enrichment and focused on modules associated with pigmentation processes because they are the best-studied and most well-known risk factors for melanoma susceptibility. We found that pigmentation-associated modules across the four expression datasets examined were enriched for well-known melanoma susceptibility genes plus genes associated with pigmentation. We also used network properties to prioritize genes within pigmentation modules as candidate susceptibility genes. Integrating information from coexpression network analysis and variant prioritization, we identified 36 genes (such as DCT, TPCN2, TRPM1, ATP10A, and EPHA5) as potential melanoma risk genes in the families. Our approach also allowed us to link families with private gene mutations on the basis of gene coexpression patterns and thereby may provide an innovative perspective in gene identification in high-risk families.

Functional annotation and investigation of the 10q24.33 melanoma risk locus identifies a common variant that influences transcriptional regulation of OBFC1.

The 10q24.33 locus is known to be associated with susceptibility to cutaneous malignant melanoma (CMM), but the mechanisms underlying this association have been not extensively investigated. We carried out an integrative genomic analysis of 10q24.33 using epigenomic annotations and in vitro reporter gene assays to identify regulatory variants. We found two putative functional single nucleotide polymorphisms (SNPs) in an enhancer and in the promoter of OBFC1, respectively, in neural crest and CMM cells, one, rs2995264, altering enhancer activity. The minor allele G of rs2995264 correlated with lower OBFC1 expression in 470 CMM tumors and was confirmed to increase the CMM risk in a cohort of 484 CMM cases and 1801 controls of Italian origin. Hi-C and chromosome conformation capture (3C) experiments showed the interaction between the enhancer-SNP region and the promoter of OBFC1 and an isogenic model characterized by CRISPR-Cas9 deletion of the enhancer-SNP region confirmed the potential regulatory effect of rs2995264 on OBFC1 transcription. Moreover, the presence of G-rs2995264 risk allele reduced the binding affinity of the transcription factor MEOX2. Biologic investigations showed significant cell viability upon depletion of OBFC1, specifically in CMM cells that were homozygous for the protective allele. Clinically, high levels of OBFC1 expression associated with histologically favorable CMM tumors. Finally, preliminary results suggested the potential effect of decreased OBFC1 expression on telomerase activity in tumorigenic conditions. Our results support the hypothesis that reduced expression of OBFC1 gene through functional heritable DNA variation can contribute to malignant transformation of normal melanocytes.

Novel MAPK/AKT-impairing germline NRAS variant identified in a melanoma-prone family.

While several high-penetrance melanoma risk genes are known, variation in these genes fail to explain melanoma susceptibility in a large proportion of high-risk families. As part of a melanoma family sequencing study, including 435 families from Mediterranean populations we identified a novel NRAS variant (c.170A > C, p.D57A) in an Italian melanoma-prone family. This variant is absent in exomes in gnomAD, ESP, UKBiobank, and the 1000 Genomes Project, as well as in 11,273 Mediterranean individuals and 109 melanoma-prone families from the US and Australia. This variant occurs in the GTP-binding pocket of NRAS. Differently from other RAS activating alterations, NRAS D57A expression is unable to activate MAPK-pathway both constitutively and after stimulation but enhances EGF-induced PI3K-pathway signaling in serum starved conditions in vitro. Consistent with in vitro data demonstrating that NRAS D57A does not enrich GTP binding, molecular modeling suggests that the D57A substitution would be expected to impair Mg2 + binding and decrease nucleotide-binding and GTPase activity of NRAS. While we cannot firmly establish NRAS c.170A > C (p.D57A) as a melanoma susceptibility variant, further investigation of NRAS as a familial melanoma gene is warranted.

Toxoplasma bradyzoites exhibit physiological plasticity of calcium and energy stores controlling motility and egress.

Toxoplasma gondii has evolved different developmental stages for disseminating during acute infection (i.e., tachyzoites) and establishing chronic infection (i.e., bradyzoites). Calcium ion (Ca2+) signaling tightly regulates the lytic cycle of tachyzoites by controlling microneme secretion and motility to drive egress and cell invasion. However, the roles of Ca2+ signaling pathways in bradyzoites remain largely unexplored. Here, we show that Ca2+ responses are highly restricted in bradyzoites and that they fail to egress in response to agonists. Development of dual-reporter parasites revealed dampened Ca2+ responses and minimal microneme secretion by bradyzoites induced in vitro or harvested from infected mice and tested ex vivo. Ratiometric Ca2+ imaging demonstrated lower Ca2+ basal levels, reduced magnitude, and slower Ca2+ kinetics in bradyzoites compared with tachyzoites stimulated with agonists. Diminished responses in bradyzoites were associated with downregulation of Ca2+-ATPases involved in intracellular Ca2+ storage in the endoplasmic reticulum (ER) and acidocalcisomes. Once liberated from cysts by trypsin digestion, bradyzoites incubated in glucose plus Ca2+ rapidly restored their intracellular Ca2+ and ATP stores, leading to enhanced gliding. Collectively, our findings indicate that intracellular bradyzoites exhibit dampened Ca2+ signaling and lower energy levels that restrict egress, and yet upon release they rapidly respond to changes in the environment to regain motility.

A large Canadian cohort provides insights into the genetic architecture of human hair colour.

Hair colour is a polygenic phenotype that results from differences in the amount and ratio of melanins located in the hair bulb. Genome-wide association studies (GWAS) have identified many loci involved in the pigmentation pathway affecting hair colour. However, most of the associated loci overlap non-protein coding regions and many of the molecular mechanisms underlying pigmentation variation are still not understood. Here, we conduct GWAS meta-analyses of hair colour in a Canadian cohort of 12,741 individuals of European ancestry. By performing fine-mapping analyses we identify candidate causal variants in pigmentation loci associated with blonde, red and brown hair colour. Additionally, we observe colocalization of several GWAS hits with expression and methylation quantitative trait loci (QTLs) of cultured melanocytes. Finally, transcriptome-wide association studies (TWAS) further nominate the expression of EDNRB and CDK10 as significantly associated with hair colour. Our results provide insights on the mechanisms regulating pigmentation biology in humans.