Integrating physics in deep learning algorithms: a force field as a PyTorch module

MOTIVATION

Deep learning algorithms applied to structural biology often struggle to converge to meaningful solutions when limited data is available, since they are required to learn complex physical rules from examples. State-of-the-art force-fields, however, cannot interface with deep learning algorithms due to their implementation.

RESULTS

We present MadraX, a forcefield implemented as a differentiable PyTorch module, able to interact with deep learning algorithms in an end-to-end fashion.

AVAILABILITY AND IMPLEMENTATION

MadraX documentation, together with tutorials and installation guide, is available at madrax.readthedocs.io.

Exercise pulmonary hypertension by the mPAP/CO slope in primary mitral regurgitation

AIMS

Exercise-induced pulmonary hypertension (PH), defined by a mean pulmonary arterial pressure over cardiac output (mPAP/CO) slope >3 mmHg/L/min, has important diagnostic and prognostic implications. The aim of this study is to investigate the value of the mPAP/CO slope in patients with more than moderate primary mitral regurgitation (MR) with preserved ejection fraction and no or discordant symptoms.

METHODS AND RESULTS

A total of 128 consecutive patients were evaluated with exercise echocardiography and cardiopulmonary testing. Clinical outcome was defined as the composite of mitral valve intervention, new-onset atrial fibrillation, cardiovascular hospitalization, and all-cause mortality. The mean age was 63 years, 61% were male, and the mean LVEF was 66 ± 6%. The mPAP/CO slope correlated with peak VO2 (r = -0.52, P < 0.001), while the peak systolic pulmonary artery pressure (sPAP) did not (r = -0.06, P = 0.584). Forty-six per cent (n = 59) had peak exercise sPAP ≥60 mmHg, and 37% (n = 47) had mPAP/CO slope >3 mmHg/L/min. Event-free survival was 55% at 1 year and 46% at 2 years, with reduced survival in patients with mPAP/CO slope >3 mmHg/L/min (hazard ratio, 4.9; 95% confidence interval, 2.9-8.2; P < 0.001). In 53 cases (41%), mPAP/CO slope and peak sPAP were discordant: patients with slope >3 mmHg/L/mmHg and sPAP <60 mmHg (n = 21) had worse outcome vs. peak sPAP ≥60 mmHg and normal slope (n = 32, log-rank P = 0.003). The mPAP/CO slope improved predictive models for outcome, incremental to resting and exercise sPAP, and peak VO2.

CONCLUSION

Exercise PH defined by the mPAP/CO slope >3 mmHg/L/min is associated with decreased exercise capacity and a higher risk of adverse events in significant primary MR and no or discordant symptoms. The slope provides a greater prognostic value than single sPAP measures and peak VO2.

Very high epilepsy prevalence in rural Southern Rwanda: The underestimated burden of epilepsy in sub-Saharan Africa

OBJECTIVES

Up to 85% of people living with epilepsy (PwE) reside in low-and middle-income countries. In sub-Saharan Africa, the lifetime prevalence of epilepsy is 16 per 1000 persons. In Northern rural Rwanda, a 47.7 per 1000 prevalence has been reported. As variations in prevalence across geographical areas have been observed, we studied the prevalence in Southern rural Rwanda using the same robust methodology as applied in the North.

METHODS

We conducted a three-stage, cross-sectional, door-to-door survey in two rural villages in Southern Rwanda from June 2022 to April 2023. First, trained enumerators administered the validated Limoges questionnaire for epilepsy screening. Second, neurologists examined the persons who had screened positively to confirm the epilepsy diagnosis. Third, cases with an inconclusive assessment were separately reexamined by two neurologists to reevaluate the diagnosis.

RESULTS

Enumerators screened 1745 persons (54.4% female, mean age: 24 ± 19.3 years), of whom 304 (17.4%) screened positive. Epilepsy diagnosis was confirmed in 133 (52.6% female, mean age: 30 ± 18.2 years) and active epilepsy in 130 persons. Lifetime epilepsy prevalence was 76.2 per 1000 (95% CI: 64.2-89.7‰). The highest age-specific rate occurred in the 29-49 age group. No gender-specific differences were noted. In 22.6% of the PwE, only non-convulsive seizures occurred. The treatment gap was 92.2%, including a diagnosis gap of 79.4%.

CONCLUSION

We demonstrated a very high epilepsy prevalence in Southern rural Rwanda, with over 20% of cases having only non-convulsive seizures, which are often underdiagnosed in rural Africa. In line with previous Rwandan reports, we reiterate the high burden of the disease in the country. Geographic variation in prevalence throughout Africa may result from differences in risk and aetiological factors. Case-control studies are underway to understand such differences and propose adapted health policies for epilepsy prevention.

Transcriptomic meta-analysis to identify potential antifungal targets in Candida albicans

BACKGROUND

Candida albicans is a fungal pathogen causing human infections. Here we investigated differential gene expression patterns and functional enrichment in C. albicans strains grown under different conditions.

METHODS

A systematic GEO database search identified 239 "Candida albicans" datasets, of which 14 were selected after rigorous criteria application. Retrieval of raw sequencing data from the ENA database was accompanied by essential metadata extraction from dataset descriptions and original articles. Pre-processing via the tailored nf-core pipeline for C. albicans involved alignment, gene/transcript quantification, and diverse quality control measures. Quality assessment via PCA and DESeq2 identified significant genes (FDR < = 0.05, log2-fold change > = 1 or <= -1), while topGO conducted GO term enrichment analysis. Exclusions were made based on data quality and strain relevance, resulting in the selection of seven datasets from the SC5314 strain background for in-depth investigation.

RESULTS

The meta-analysis of seven selected studies unveiled a substantial number of genes exhibiting significant up-regulation (24,689) and down-regulation (18,074). These differentially expressed genes were further categorized into 2,497 significantly up-regulated and 2,573 significantly down-regulated Gene Ontology (GO) IDs. GO term enrichment analysis clustered these terms into distinct groups, providing insights into the functional implications. Three target gene lists were compiled based on previous studies, focusing on central metabolism, ion homeostasis, and pathogenicity. Frequency analysis revealed genes with higher occurrence within the identified GO clusters, suggesting their potential as antifungal targets. Notably, the genes TPS2, TPS1, RIM21, PRA1, SAP4, and SAP6 exhibited higher frequencies within the clusters. Through frequency analysis within the GO clusters, several key genes emerged as potential targets for antifungal therapies. These include RSP5, GLC7, SOD2, SOD5, SOD1, SOD6, SOD4, SOD3, and RIM101 which exhibited higher occurrence within the identified clusters.

CONCLUSION

This comprehensive study significantly advances our understanding of the dynamic nature of gene expression in C. albicans. The identification of genes with enhanced potential as antifungal drug targets underpins their value for future interventions. The highlighted genes, including TPS2, TPS1, RIM21, PRA1, SAP4, SAP6, RSP5, GLC7, SOD2, SOD5, SOD1, SOD6, SOD4, SOD3, and RIM101, hold promise for the development of targeted antifungal therapies.

Nanoscale Chemical Diversity of Coke Deposits on Nanoprinted Metal Catalysts Visualized by Tip-Enhanced Raman Spectroscopy

Coke formation is the prime cause of catalyst deactivation, where undesired carbon wastes block the catalyst surface and hinder further reaction in a broad gamut of industrial chemical processes. Yet, the origins of coke formation and their distribution across the catalyst remain elusive, obstructing the design of coke-resistant catalysts. Here, the first-time application of tip-enhanced Raman spectroscopy (TERS) is demonstrated as a nanoscale chemical probe to localize and identify coke deposits on a post-mortem metal nanocatalyst. Monitoring coke at the nanoscale circumvents bulk averaging and reveals the local nature of coke with unmatched detail. The nature of coke is chemically diverse and ranges from nanocrystalline graphite to disordered and polymeric coke, even on a single nanoscale location of a top-down nanoprinted SiO2 -supported Pt catalyst. Surprisingly, not all Pt is an equal producer of coke, where clear isolated coke "hotspots" are present non-homogeneously on Pt which generate large amounts of disordered coke. After their formation, coke shifts to the support and undergoes long-range transport on the surrounding SiO2 surface, where it becomes more graphitic. The presented results provide novel guidelines to selectively free-up the coked metal surface at more mild rejuvenation conditions, thus securing the long-term catalyst stability.

Achieving ultra-low and -uniform residual magnetic fields in a very large magnetically shielded room for fundamental physics experiments

High-precision searches for an electric dipole moment of the neutron (nEDM) require stable and uniform magnetic field environments. We present the recent achievements of degaussing and equilibrating the magnetically shielded room (MSR) for the n2EDM experiment at the Paul Scherrer Institute. We present the final degaussing configuration that will be used for n2EDM after numerous studies. The optimized procedure results in a residual magnetic field that has been reduced by a factor of two. The ultra-low field is achieved with the full magnetic-field-coil system, and a large vacuum vessel installed, both in the MSR. In the inner volume of ∼ 1.4 m 3 , the field is now more uniform and below 300 pT. In addition, the procedure is faster and dissipates less heat into the magnetic environment, which in turn, reduces its thermal relaxation time from 12 h down to 1.5 h .

Tools to differentiate between Filamin C and Titin truncating variant carriers: value of MRI

Whereas truncating variants of the giant protein Titin (TTNtv) are the main cause of familial dilated cardiomyopathy (DCM), recently Filamin C truncating variants (FLNCtv) were identified as a cause of arrhythmogenic cardiomyopathy (ACM). Our aim was to characterize and compare clinical and MRI features of TTNtv and FLNCtv in the Belgian population. In index patients referred for genetic testing of ACM/DCM, FLNCtv and TTNtv were found in 17 (3.6%) and 33 (12.3%) subjects, respectively. Further family cascade screening yielded 24 and 19 additional truncating variant carriers in FLNC and TTN, respectively. The main phenotype was ACM in FLNCtv carriers whereas TTNtv carriers showed either an ACM or DCM phenotype. Non-sustained Ventricular Tachycardia was frequent in both populations. MRI data, available in 28/40 FLNCtv and 32/52 TTNtv patients, showed lower Left Ventricular (LV) ejection fraction and lower LV strain in TTNtv patients (p < 0.01). Conversely, both the frequency (68% vs 22%) and extent of non-ischemic myocardial late gadolinium enhancement (LGE) was significantly higher in FLNCtv patients (p < 0.01). Hereby, ring-like LGE was found in 16/19 (84%) FLNCtv versus 1/7 (14%) of TTNtv patients (p < 0.01). In conclusion, a large number of FLNCtv and TTNtv patients present with an ACM phenotype but can be separated by cardiac MRI. Whereas FLNCtv patients often have extensive myocardial fibrosis, typically following a ring-like pattern, LV dysfunction without or limited replacement fibrosis is the common TTNtv phenotype.

SNPeffect 5.0: large-scale structural phenotyping of protein coding variants extracted from next-generation sequencing data using AlphaFold models

BACKGROUND

Next-generation sequencing technologies yield large numbers of genetic alterations, of which a subset are missense variants that alter an amino acid in the protein product. These variants can have a potentially destabilizing effect leading to an increased risk of misfolding and aggregation. Multiple software tools exist to predict the effect of single-nucleotide variants on proteins, however, a pipeline integrating these tools while starting from an NGS data output list of variants is lacking.

RESULTS

The previous version SNPeffect 4.0 (De Baets in Nucleic Acids Res 40(D1):D935-D939, 2011) provided an online database containing pre-calculated variant effects and low-throughput custom variant analysis. Here, we built an automated and parallelized pipeline that analyzes the impact of missense variants on the aggregation propensity and structural stability of proteins starting from the Variant Call Format as input. The pipeline incorporates the AlphaFold Protein Structure Database to achieve high coverage for structural stability analyses using the FoldX force field. The effect on aggregation-propensity is analyzed using the established predictors TANGO and WALTZ. The pipeline focuses solely on the human proteome and can be used to analyze proteome stability/damage in a given sample based on sequencing results.

CONCLUSION

We provide a bioinformatics pipeline that allows structural phenotyping from sequencing data using established stability and aggregation predictors including FoldX, TANGO, and WALTZ; and structural proteome coverage provided by the AlphaFold database. The pipeline and installation guide are freely available for academic users on https://github.com/vibbits/snpeffect and requires a computer cluster.

Lifelong endurance exercise and its relation with coronary atherosclerosis

AIMS

The impact of long-term endurance sport participation (on top of a healthy lifestyle) on coronary atherosclerosis and acute cardiac events remains controversial.

METHODS AND RESULTS

The Master@Heart study is a well-balanced prospective observational cohort study. Overall, 191 lifelong master endurance athletes, 191 late-onset athletes (endurance sports initiation after 30 years of age), and 176 healthy non-athletes, all male with a low cardiovascular risk profile, were included. Peak oxygen uptake quantified fitness. The primary endpoint was the prevalence of coronary plaques (calcified, mixed, and non-calcified) on computed tomography coronary angiography. Analyses were corrected for multiple cardiovascular risk factors. The median age was 55 (50-60) years in all groups. Lifelong and late-onset athletes had higher peak oxygen uptake than non-athletes [159 (143-177) vs. 155 (138-169) vs. 122 (108-138) % predicted]. Lifelong endurance sports was associated with having ≥1 coronary plaque [odds ratio (OR) 1.86, 95% confidence interval (CI) 1.17-2.94], ≥ 1 proximal plaque (OR 1.96, 95% CI 1.24-3.11), ≥ 1 calcified plaques (OR 1.58, 95% CI 1.01-2.49), ≥ 1 calcified proximal plaque (OR 2.07, 95% CI 1.28-3.35), ≥ 1 non-calcified plaque (OR 1.95, 95% CI 1.12-3.40), ≥ 1 non-calcified proximal plaque (OR 2.80, 95% CI 1.39-5.65), and ≥1 mixed plaque (OR 1.78, 95% CI 1.06-2.99) as compared to a healthy non-athletic lifestyle.

CONCLUSION

Lifelong endurance sport participation is not associated with a more favourable coronary plaque composition compared to a healthy lifestyle. Lifelong endurance athletes had more coronary plaques, including more non-calcified plaques in proximal segments, than fit and healthy individuals with a similarly low cardiovascular risk profile. Longitudinal research is needed to reconcile these findings with the risk of cardiovascular events at the higher end of the endurance exercise spectrum.

Investigating the development of the autonomic nervous system in infancy through pupillometry

We aim to investigate early developmental trajectories of the autonomic nervous system (ANS) as indexed by the pupillary light reflex (PLR) in infants with (i.e. preterm birth, feeding difficulties, or siblings of children with autism spectrum disorder) and without (controls) increased likelihood for atypical ANS development. We used eye-tracking to capture the PLR in 216 infants in a longitudinal follow-up study spanning 5 to 24 months of age, and linear mixed models to investigate effects of age and group on three PLR parameters: baseline pupil diameter, latency to constriction and relative constriction amplitude. An increase with age was found in baseline pupil diameter (F(3,273.21) = 13.15, p < 0.001, [Formula: see text] = 0.13), latency to constriction (F(3,326.41) = 3.84, p = 0.010, [Formula: see text] = 0.03) and relative constriction amplitude(F(3,282.53) = 3.70, p = 0.012, [Formula: see text] = 0.04). Group differences were found for baseline pupil diameter (F(3,235.91) = 9.40, p < 0.001, [Formula: see text] = 0.11), with larger diameter in preterms and siblings than in controls, and for latency to constriction (F(3,237.10) = 3.48, p = 0.017, [Formula: see text] = 0.04), with preterms having a longer latency than controls. The results align with previous evidence, with development over time that could be explained by ANS maturation. To better understand the cause of the group differences, further research in a larger sample is necessary, combining pupillometry with other measures to further validate its value.

Histopathological correlations and fat replacement imaging patterns in recessive limb-girdle muscular dystrophy type 12

BACKGROUND

Despite the widespread use of proton density fat fraction (PDFF) measurements with magnetic resonance imaging (MRI) to track disease progression in muscle disorders, it is still unclear how these findings relate to histopathological changes in muscle biopsies of patients with limb-girdle muscular dystrophy autosomal recessive type 12 (LGMDR12). Furthermore, although it is known that LGMDR12 leads to a selective muscle involvement distinct from other muscular dystrophies, the spatial distribution of fat replacement within these muscles is unknown.

METHODS

We included 27 adult patients with LGMDR12 and 27 age-matched and sex-matched healthy controls and acquired 6-point Dixon images of the thighs and T1 and short tau inversion recovery (STIR) MR images of the whole body. In 16 patients and 15 controls, we performed three muscle biopsies, one in the semimembranosus, vastus lateralis, and rectus femoris muscles, which are severely, intermediately, and mildly affected in LGMDR12, respectively. We correlated the PDFF to the fat percentage measured on biopsies of the corresponding muscles, as well as to the Rochester histopathology grading scale.

RESULTS

In patients, we demonstrated a strong correlation of PDFF on MRI and muscle biopsy fat percentage for the semimembranosus (r = 0.85, P < 0.001) and vastus lateralis (r = 0.68, P = 0.005). We found similar results for the correlation between PDFF and the Rochester histopathology grading scale. Out of the five patients with inflammatory changes on muscle biopsy, three showed STIR hyperintensities in the corresponding muscle on MRI. By modelling the PDFF on MRI for 18 thigh muscles from origin to insertion, we observed a significantly inhomogeneous proximo-distal distribution of fat replacement in all thigh muscles of patients with LGMDR12 (P < 0.001), and different patterns of fat replacement within each of the muscles.

CONCLUSIONS

We showed a strong correlation of fat fraction on MRI and fat percentage on muscle biopsy for diseased muscles and validated the use of Dixon fat fraction imaging as an outcome measure in LGMDR12. The inhomogeneous fat replacement within thigh muscles on imaging underlines the risk of analysing only samples of muscles instead of the entire muscles, which has important implications for clinical trials.

From DNA Copy Number Gains and Tumor Dependencies to Novel Therapeutic Targets for High-Risk Neuroblastoma

Neuroblastoma is a pediatric tumor arising from the sympatho-adrenal lineage and a worldwide leading cause of childhood cancer-related deaths. About half of high-risk patients die from the disease while survivors suffer from multiple therapy-related side-effects. While neuroblastomas present with a low mutational burden, focal and large segmental DNA copy number aberrations are highly recurrent and associated with poor survival. It can be assumed that the affected chromosomal regions contain critical genes implicated in neuroblastoma biology and behavior. More specifically, evidence has emerged that several of these genes are implicated in tumor dependencies thus potentially providing novel therapeutic entry points. In this review, we briefly review the current status of recurrent DNA copy number aberrations in neuroblastoma and provide an overview of the genes affected by these genomic variants for which a direct role in neuroblastoma has been established. Several of these genes are implicated in networks that positively regulate MYCN expression or stability as well as cell cycle control and apoptosis. Finally, we summarize alternative approaches to identify and prioritize candidate copy-number driven dependency genes for neuroblastoma offering novel therapeutic opportunities.

Cyclin D2 overexpression drives B1a-derived MCL-like lymphoma in mice

Mantle cell lymphoma (MCL) is an aggressive B cell lymphoma with poor long-term overall survival. Currently, MCL research and development of potential cures is hampered by the lack of good in vivo models. MCL is characterized by recurrent translocations of CCND1 or CCND2, resulting in overexpression of the cell cycle regulators cyclin D1 or D2, respectively. Here, we show, for the first time, that hematopoiesis-specific activation of cyclin D2 is sufficient to drive murine MCL-like lymphoma development. Furthermore, we demonstrate that cyclin D2 overexpression can synergize with loss of p53 to form aggressive and transplantable MCL-like lymphomas. Strikingly, cyclin D2-driven lymphomas display transcriptional, immunophenotypic, and functional similarities with B1a B cells. These MCL-like lymphomas have B1a-specific B cell receptors (BCRs), show elevated BCR and NF-κB pathway activation, and display increased MALT1 protease activity. Finally, we provide preclinical evidence that inhibition of MALT1 protease activity, which is essential for the development of early life-derived B1a cells, can be an effective therapeutic strategy to treat MCL.

Quantifying the Likelihood of Structural Models through a Dynamically Enhanced Powder X-Ray Diffraction Protocol

Structurally characterizing new materials is tremendously challenging, especially when single crystal structures are hardly available which is often the case for covalent organic frameworks. Yet, knowledge of the atomic structure is key to establish structure-function relations and enable functional material design. Herein, a new protocol is proposed to unambiguously predict the structure of poorly crystalline materials through a likelihood ordering based on the X-ray diffraction (XRD) pattern. Key of the procedure is the broad set of structures generated from a limited number of building blocks and topologies, which is submitted to operando structural characterization. The dynamic averaging in the latter accounts for the operando conditions and inherent temporal character of experimental measurements, yielding unparalleled agreement with experimental powder XRD patterns. The proposed concept can hence unquestionably identify the structure of experimentally synthesized materials, a crucial step to design next generation functional materials.

Energy metabolism in ALS: an underappreciated opportunity?

Amyotrophic lateral sclerosis (ALS) is a relentlessly progressive and fatal neurodegenerative disorder that primarily affects motor neurons. Despite our increased understanding of the genetic factors contributing to ALS, no effective treatment is available. A growing body of evidence shows disturbances in energy metabolism in ALS. Moreover, the remarkable vulnerability of motor neurons to ATP depletion has become increasingly clear. Here, we review metabolic alterations present in ALS patients and models, discuss the selective vulnerability of motor neurons to energetic stress, and provide an overview of tested and emerging metabolic approaches to treat ALS. We believe that a further understanding of the metabolic biology of ALS can lead to the identification of novel therapeutic targets.

Semi-automated digital measurement as the method of choice for beta cell mass analysis

Pancreas injury by partial duct ligation (PDL) activates beta cell differentiation and proliferation in adult mouse pancreas but remains controversial regarding the anticipated increase in beta cell volume. Several reports unable to show beta cell volume augmentation in PDL pancreas used automated digital image analysis software. We hypothesized that fully automatic beta cell morphometry without manual micrograph artifact remediation introduces bias and therefore might be responsible for reported discrepancies and controversy. However, our present results prove that standard digital image processing with automatic thresholding is sufficiently robust albeit less sensitive and less adequate to demonstrate a significant increase in beta cell volume in PDL versus Sham-operated pancreas. We therefore conclude that other confounding factors such as quality of surgery, selection of samples based on relative abundance of the transcription factor Neurogenin 3 (Ngn3) and tissue processing give rise to inter-laboratory inconsistencies in beta cell volume quantification in PDL pancreas.

Covert deformed wing virus infections have long-term deleterious effects on honeybee foraging and survival

Several studies have suggested that covert stressors can contribute to bee colony declines. Here we provide a novel case study and show using radiofrequency identification tracking technology that covert deformed wing virus (DWV) infections in adult honeybee workers seriously impact long-term foraging and survival under natural foraging conditions. In particular, our experiments show that adult workers injected with low doses of DWV experienced increased mortality rates, that DWV caused workers to start foraging at a premature age, and that the virus reduced the workers' total activity span as foragers. Altogether, these results demonstrate that covert DWV infections have strongly deleterious effects on honeybee foraging and survival. These results are consistent with previous studies that suggested DWV to be an important contributor to the ongoing bee declines in Europe and the USA. Overall, our study underlines the strong impact that covert pathogen infections can have on individual and group-level performance in bees.

Combined Activity of DCL2 and DCL3 Is Crucial in the Defense against Potato Spindle Tuber Viroid

Viroids are self replicating non-coding RNAs capable of infecting a wide range of plant hosts. They do not encode any proteins, thus the mechanism by which they escape plant defenses remains unclear. RNAi silencing is a major defense mechanism against virus infections, with the four DCL proteins being principal components of the pathway. We have used Nicotiana benthamiana as a model to study Potato spindle tuber viroid infection. This viroid is a member of the Pospiviroidae family and replicates in the nucleus via an asymmetric rolling circle mechanism. We have created knock-down plants for all four DCL genes and their combinations. Previously, we showed that DCL4 has a positive effect on PSTVd infectivity since viroid levels drop when DCL4 is suppressed. Here, we show that PSTVd levels remain decreased throughout infection in DCL4 knockdown plants, and that simultaneous knockdown of DCL1, DCL2 or DCL3 together with DCL4 cannot reverse this effect. Through infection of plants suppressed for multiple DCLs we further show that a combined suppression of DCL2 and DCL3 has a major effect in succumbing plant antiviral defense. Based on our results, we further suggest that Pospoviroids may have evolved to be primarily processed by DCL4 as it seems to be a DCL protein with less detrimental effects on viroid infectivity. These findings pave the way to delineate the complexity of the relationship between viroids and plant RNA silencing response.

The zinc finger transcription factor PW1/PEG3 restrains murine beta cell cycling

AIMS/HYPOTHESIS

Pw1 or paternally-expressed gene 3 (Peg3) encodes a zinc finger transcription factor that is widely expressed during mouse embryonic development and later restricted to multiple somatic stem cell lineages in the adult. The aim of the present study was to define Pw1 expression in the embryonic and adult pancreas and investigate its role in the beta cell cycle in Pw1 wild-type and mutant mice.

METHODS

We analysed PW1 expression by immunohistochemistry in pancreas of nonpregant and pregnant mice and following injury by partial duct ligation. Its role in the beta cell cycle was studied in vivo using a novel conditional knockout mouse and in vitro by lentivirus-mediated gene knockdown.

RESULTS

We showed that PW1 is expressed in early pancreatic progenitors at E9.5 but becomes progressively restricted to fully differentiated beta cells as they become established after birth and withdraw from the cell cycle. Notably, PW1 expression declines when beta cells are induced to proliferate and loss of PW1 function activates the beta cell cycle.

CONCLUSIONS/INTERPRETATION

These results indicate that PW1 is a co-regulator of the beta cell cycle and can thus be considered a novel therapeutic target in diabetes.