Integrating single cell expression quantitative trait loci summary statistics to understand complex trait risk genes.

Transcriptome-wide association study (TWAS) is a popular approach to dissect the functional consequence of disease associated non-coding variants. Most existing TWAS use bulk tissues and may not have the resolution to reveal cell-type specific target genes. Single-cell expression quantitative trait loci (sc-eQTL) datasets are emerging. The largest bulk- and sc-eQTL datasets are most conveniently available as summary statistics, but have not been broadly utilized in TWAS. Here, we present a new method EXPRESSO (EXpression PREdiction with Summary Statistics Only), to analyze sc-eQTL summary statistics, which also integrates 3D genomic data and epigenomic annotation to prioritize causal variants. EXPRESSO substantially improves existing methods. We apply EXPRESSO to analyze multi-ancestry GWAS datasets for 14 autoimmune diseases. EXPRESSO uniquely identifies 958 novel gene x trait associations, which is 26% more than the second-best method. Among them, 492 are unique to cell type level analysis and missed by TWAS using whole blood. We also develop a cell type aware drug repurposing pipeline, which leverages EXPRESSO results to identify drug compounds that can reverse disease gene expressions in relevant cell types. Our results point to multiple drugs with therapeutic potentials, including metformin for type 1 diabetes, and vitamin K for ulcerative colitis.

Genomic Insights for Personalized Care: Motivating At-Risk Individuals Toward Evidence-Based Health Practices.

Lung cancer and tobacco use pose significant global health challenges and require a comprehensive translational roadmap for improved prevention strategies. We propose the GREAT care paradigm ( G enomic Informed Care for Motivating High R isk Individuals E ligible for Evidence-b a sed Prevention), which employs polygenic risk scores (PRSs) to stratify disease risk and personalize interventions, such as lung cancer screening and tobacco treatment. We developed PRSs using large-scale multi-ancestry genome-wide association studies and adjusted for genetic ancestry for standardized risk stratification across diverse populations. We applied our PRSs to over 340,000 individuals of diverse ethnic background and found significant odds ratios for lung cancer and difficulty quitting smoking. These findings enable the evaluation of PRS-based interventions in ongoing trials aimed at motivating health behavior changes in high-risk patients. This pioneering approach enhances primary care with genomic insights, promising improved outcomes in cancer prevention and tobacco treatment, and is currently under assessment in clinical trials.

Genetic imputation of kidney transcriptome, proteome and multi-omics illuminates new blood pressure and hypertension targets.

Genetic mechanisms of blood pressure (BP) regulation remain poorly defined. Using kidney-specific epigenomic annotations and 3D genome information we generated and validated gene expression prediction models for the purpose of transcriptome-wide association studies in 700 human kidneys. We identified 889 kidney genes associated with BP of which 399 were prioritised as contributors to BP regulation. Imputation of kidney proteome and microRNAome uncovered 97 renal proteins and 11 miRNAs associated with BP. Integration with plasma proteomics and metabolomics illuminated circulating levels of myo-inositol, 4-guanidinobutanoate and angiotensinogen as downstream effectors of several kidney BP genes (SLC5A11, AGMAT, AGT, respectively). We showed that genetically determined reduction in renal expression may mimic the effects of rare loss-of-function variants on kidney mRNA/protein and lead to an increase in BP (e.g., ENPEP). We demonstrated a strong correlation (r = 0.81) in expression of protein-coding genes between cells harvested from urine and the kidney highlighting a diagnostic potential of urinary cell transcriptomics. We uncovered adenylyl cyclase activators as a repurposing opportunity for hypertension and illustrated examples of BP-elevating effects of anticancer drugs (e.g. tubulin polymerisation inhibitors). Collectively, our studies provide new biological insights into genetic regulation of BP with potential to drive clinical translation in hypertension.

Research Progress of Caspase in Endometriosis.

Endometriosis, a common chronic gynecological disease, refers to the presence and proliferation of endometrial tissue in locations other than the uterine cavity. Approximately 6 to 10% of the population of women of childbearing age are known to have endometriosis; the most common clinical signs are pelvic pain and infertility. Although endometriosis is a benign disease, it exhibits some typical features of malignant tumors, such as proliferation, invasion, metastasis, and recurrence. Endometriosis is considered a chronic, inflammatory, and estrogen-dependent disease, and multiple factors contribute to its occurrence and development. In recent years, increasing attention has been given to the role of apoptosis in the pathogenesis of this disease. Some researchers believe that spontaneous apoptosis of the endometrium is critical in maintaining its normal structure and function, and abnormal apoptosis can promote the occurrence and development of endometriosis. Inflammation is another likely process in the pathogenesis of endometriosis. Inflammation mediates the adhesion, proliferation, differentiation, and invasion of ectopic lesions of endometriosis, primarily by regulating the function of immune cells and increasing the level of proinflammatory cytokines in body fluids. The ultimate initiators of apoptosis and inflammatory cell death (pyroptosis) are the caspase family proteases. In this article, we review the progress in recent years in caspase function as well as the possible role of these enzymes in the pathogenesis of endometriosis, indicating potential treatment strategies.

Association of surfactant protein A2 with acute respiratory failure in children.

BACKGROUND: Interactions among single nucleotide polymorphisms (SNPs) of surfactant protein (SP) are associated with acute respiratory failure (ARF) and its short-term outcome, pulmonary dysfunction at discharge (PDAD) in children. However, genetic association studies using individual SNPs have not been conducted before. We hypothesize that SP genetic variants are associated with pediatric ARF and its short-term complications by themselves. METHODS: We used available genotype and clinical data in the Floros biobank consisting of 248 children aged ≤24 months with ARF; 86 developed PDAD. A logistic regression analysis was performed for each of the 14 selected SNPs, SP-A1 and SP-A2 genotypes. A p-value less than the Bonferroni correction threshold was considered significant. A likelihood ratio test was done to compare two models (one with demographic data and another with genetic variants). RESULTS: Before Bonferroni correction, female sex is associated with a decreased risk of ARF. Black race and the rs721917 of the SFTPD are associated with increased risk of ARF. After Bonferroni correction, the 1A0 1A1 genotype of SFTPA2 was associated with decreased risk of ARF. The likelihood ratio test showed that the model of the genotype information with demographic data was a better fit to predict ARF risk. None of the SP SNPs and SP-A1, SP-A2 genotypes were associated with PDAD. CONCLUSION: Our results indicate that SNPs and genotypes of SPs involved in innate immunity and host defense play an important role in ARF and, in the future, may be used as biomarkers.

The pivotal role of the X-chromosome in the genetic architecture of the human brain.

Genes on the X-chromosome are extensively expressed in the human brain, resulting in substantial influences on brain development, intellectual disability, and other brain-related disorders. To comprehensively investigate the X-chromosome's impact on the cerebral cortex, white matter tract microstructures, and intrinsic and extrinsic brain functions, we examined 2,822 complex brain imaging traits obtained from n=34,000 subjects in the UK Biobank. We unveiled potential autosome-X-chromosome interaction, while proposing an atlas of dosage compensation (DC) for each set of traits. We observed a pronounced X-chromosome impact on the corticospinal tract and the functional amplitude and connectivity of visual networks. In association studies, we identified 50 genome-wide significant trait-locus pairs enriched in Xq28, 22 of which replicated in independent datasets (n=4,900). Notably, 13 newly identified pairs were in the X-chromosome's non-pseudo-autosomal regions (NPR). The volume of the right ventral diencephalon shared genetic architecture with schizophrenia and educational attainment in a locus indexed by rs2361468 (located ~3kb upstream of PJA1, a conserved and ubiquitously expressed gene implicated in multiple psychiatric disorders). No significant associations were identified in the pseudo-autosomal regions (PAR) or the Y-chromosome. Finally, we explored sex-specific associations on the X-chromosome and compared differing genetic effects between sexes. We found much more associations can be identified in males (33 versus 9) given a similar sample size. In conclusion, our research provides invaluable insights into the X-chromosome's role in the human brain, contributing to the observed sex differences in brain structure and function.

Fluorine spillover for ceria- vs silica-supported palladium nanoparticles: A MD study using machine learning potentials.

Supported metallic nanoparticles play a central role in catalysis. However, predictive modeling is particularly challenging due to the structural and dynamic complexity of the nanoparticle and its interface with the support, given that the sizes of interest are often well beyond those accessible via traditional ab initio methods. With recent advances in machine learning, it is now feasible to perform MD simulations with potentials retaining near-density-functional theory (DFT) accuracy, which can elucidate the growth and relaxation of supported metal nanoparticles, as well as reactions on those catalysts, at temperatures and time scales approaching those relevant to experiments. Furthermore, the surfaces of the support materials can also be modeled realistically through simulated annealing to include effects such as defects and amorphous structures. We study the adsorption of fluorine atoms on ceria and silica supported palladium nanoparticles using machine learning potential trained by DFT data using the DeePMD framework. We show defects on ceria and Pd/ceria interfaces are crucial for the initial adsorption of fluorine, while the interplay between Pd and ceria and the reverse oxygen migration from ceria to Pd control spillover of fluorine from Pd to ceria at later stages. In contrast, silica supports do not induce fluorine spillover from Pd particles.

Native American genetic ancestry and pigmentation allele contributions to skin color in a Caribbean population.

Our interest in the genetic basis of skin color variation between populations led us to seek a Native American population with genetically African admixture but low frequency of European light skin alleles. Analysis of 458 genomes from individuals residing in the Kalinago Territory of the Commonwealth of Dominica showed approximately 55% Native American, 32% African, and 12% European genetic ancestry, the highest Native American genetic ancestry among Caribbean populations to date. Skin pigmentation ranged from 20 to 80 melanin units, averaging 46. Three albino individuals were determined to be homozygous for a causative multi-nucleotide polymorphism OCA2NW273KV contained within a haplotype of African origin; its allele frequency was 0.03 and single allele effect size was -8 melanin units. Derived allele frequencies of SLC24A5A111T and SLC45A2L374F were 0.14 and 0.06, with single allele effect sizes of -6 and -4, respectively. Native American genetic ancestry by itself reduced pigmentation by more than 20 melanin units (range 24-29). The responsible hypopigmenting genetic variants remain to be identified, since none of the published polymorphisms predicted in prior literature to affect skin color in Native Americans caused detectable hypopigmentation in the Kalinago.

The variation in skin colour of modern humans is a product of thousands of years of natural selection. All human ancestry can be traced back to African populations, which were dark-skinned to protect them from the intense UV rays of the sun. Over time, humans spread to other parts of the world, and people in the northern latitudes with lower UV developed lighter skin through natural selection. This was likely driven by a need for vitamin D, which requires UV rays for production. Separate genetic mechanisms were involved in the evolution of lighter skin in each of the two main branches of human migration: the European branch (which includes peoples on the Indian subcontinent and Europe) and the East Asian branch (which includes East Asia and the Americas). A variant of the gene SLC24A5 is the primary contributor to lighter skin colour in the European branch, but a corresponding variant driving light skin colour evolution in the East Asian branch remains to be identified. One obstacle to finding such variants is the high prevalence of European ancestry in most people groups, which makes it difficult to separate the influence of European genes from those of other populations. To overcome this issue, Ang et al. studied a population that had a high proportion of Native American and African ancestors, but a relatively small proportion of European ancestors, the Kalinago people. The Kalinago live on the island of Dominica, one of the last Caribbean islands to be colonised by Europeans. Ang et al. were able to collect hundreds of skin pigmentation measurements and DNA samples of the Kalinago, to trace the effect of Native American ancestry on skin colour. Genetic analysis confirmed their oral history records of primarily Native American (55%) ­ one of the highest of any Caribbean population studied to date ­ compared with African (32%) and European (12%) ancestries. Native American ancestry had the highest effect on pigmentation and reduced it by more than 20 melanin units, while the European mutations in the genes SLC24A5 and SLC45A2 and an African gene variant for albinism only contributed 5, 4 and 8 melanin units, respectively. However, none of the so far published gene candidates responsible for skin lightening in Native Americans caused a detectable effect. Therefore, the gene responsible for lighter skin in Native Americans/East Asians has yet to be identified. The work of Ang et al. represents an important step in deciphering the genetic basis of lighter skin colour in Native Americans or East Asians. A better understanding of the genetics of skin pigmentation may help to identify why, for example, East Asians are less susceptible to melanoma than Europeans, despite both having a lighter skin colour. It may also further acceptance of how variations in human skin tones are the result of human migration, random genetic variation, and natural selection for pigmentation in different solar environments.

Future trends in incidence and long-term survival of metastatic cancer in the United States.

BACKGROUND: Previous studies have demonstrated epidemiological trends in individual metastatic cancer subtypes; however, research forecasting long-term incidence trends and projected survivorship of metastatic cancers is lacking. We assess the burden of metastatic cancer to 2040 by (1) characterizing past, current, and forecasted incidence trends, and (2) estimating odds of long-term (5-year) survivorship. METHODS: This retrospective, serial cross-sectional, population-based study used registry data from the Surveillance, Epidemiology, and End Results (SEER 9) database. Average annual percentage change (AAPC) was calculated to describe cancer incidence trends from 1988 to 2018. Autoregressive integrating moving average (ARIMA) models were used to forecast the distribution of primary metastatic cancer and metastatic cancer to specific sites from 2019 to 2040 and JoinPoint models were fitted to estimate mean projected annual percentage change (APC). RESULTS: The average annual percent change (AAPC) in incidence of metastatic cancer decreased by 0.80 per 100,000 individuals (1988-2018) and we forecast an APC decrease by 0.70 per 100,000 individuals (2018-2040). Analyses predict a decrease in metastases to liver (APC = -3.40, 95% CI [-3.50, -3.30]), lung (APC (2019-2030) = -1.90, 95% CI [-2.90, -1.00]); (2030-2040) = -3.70, 95% CI [-4.60, -2.80]), bone (APC = -4.00, 95% CI [-4.30, -3.70]), and brain (APC = -2.30, 95% CI [-2.60, -2.00]). By 2040, patients with metastatic cancer are predicted to have 46.7% greater odds of long-term survivorship, driven by increasing plurality of patients with more indolent forms of metastatic disease. CONCLUSIONS: By 2040, the distribution of metastatic cancer patients is predicted to shift in predominance from invariably fatal to indolent cancers subtypes. Continued research on metastatic cancers is important to guide health policy and clinical intervention efforts, and direct allocations of healthcare resources.

Cancer that has spread beyond the area where it originated and into different organs is called metastatic cancer. This study analyzed trends in metastatic cancer incidence, the proportion of those with metastatic cancer surviving 5 years after diagnosis and the locations in the body each cancer had spread to. The incidence of metastatic cancer decreased between 1988 and 2018 and is expected to continue to decrease until 2040. Some of the most common locations cancer spreads to is the lung, liver, brain, and bone. Metastatic cancer incidence to these areas is predicted to decrease. Also, the likelihood of surviving for more than 5 years after diagnosis with metastatic cancer is predicted to increase by 2040. This research should facilitate optimal planning of future healthcare resources and policy.

Generic two-phase coexistence in a type-2 Schloegl model for autocatalysis on a square lattice: Analysis via heterogeneous master equations.

Schloegl's second model (also known as the quadratic contact process) on a square lattice involves spontaneous annihilation of particles at lattice sites at rate p, and their autocatalytic creation at unoccupied sites with n≥2 occupied neighbors at rate k_{n}. Kinetic Monte Carlo (KMC) simulation reveals that these models exhibit a nonequilibrium discontinuous phase transition with generic two-phase coexistence: the p value for equistability of coexisting populated and vacuum states, p_{eq}(S), depends on the orientation or slope, S, of a planar interface separating those phases. The vacuum state displaces the populated state for p>p_{eq}(S), and the opposite applies for p

Nucleation-mediated reshaping of facetted metallic nanocrystals: Breakdown of the classical free energy picture.

Shape stability is key to avoiding degradation of performance for metallic nanocrystals synthesized with facetted non-equilibrium shapes to optimize properties for catalysis, plasmonics, and so on. Reshaping of facetted nanocrystals is controlled by the surface diffusion-mediated nucleation and growth of new outer layers of atoms. Kinetic Monte Carlo (KMC) simulation of a realistic stochastic atomistic-level model is applied to precisely track the reshaping of Pd octahedra and nanocubes. Unexpectedly, separate constrained equilibrium Monte Carlo analysis of the free energy profile during reshaping reveals a fundamental failure of the classical nucleation theory (CNT) prediction for the reshaping barrier and rate. Why? Nucleation barriers can be relatively low for these processes, so the system is not locally equilibrated before crossing the barrier, as assumed in CNT. This claim is supported by an analysis of a first-passage problem for reshaping within a master equation framework for the model that reasonably captures the behavior in KMC simulations.

Multi-ancestry and multi-trait genome-wide association meta-analyses inform clinical risk prediction for systemic lupus erythematosus.

Systemic lupus erythematosus is a heritable autoimmune disease that predominantly affects young women. To improve our understanding of genetic etiology, we conduct multi-ancestry and multi-trait meta-analysis of genome-wide association studies, encompassing 12 systemic lupus erythematosus cohorts from 3 different ancestries and 10 genetically correlated autoimmune diseases, and identify 16 novel loci. We also perform transcriptome-wide association studies, computational drug repurposing analysis, and cell type enrichment analysis. We discover putative drug classes, including a histone deacetylase inhibitor that could be repurposed to treat lupus. We also identify multiple cell types enriched with putative target genes, such as non-classical monocytes and B cells, which may be targeted for future therapeutics. Using this newly assembled result, we further construct polygenic risk score models and demonstrate that integrating polygenic risk score with clinical lab biomarkers improves the diagnostic accuracy of systemic lupus erythematosus using the Vanderbilt BioVU and Michigan Genomics Initiative biobanks.

1H chemical shift anisotropy: a high sensitivity solid-state NMR dynamics probe for surface studies?

Dynamics play significant roles in chemistry and biochemistry-molecular motions impact both large- and small-scale chemical reactions in addition to biochemical processes. In many systems, including heterogeneous catalysts, the characterization of dynamics remains a challenge. The most common approaches involve the solid-state NMR measurement of anisotropic interactions, in particular 2H quadrupolar coupling and 1H-X dipolar coupling, which generally require isotope enrichment. Due to the high sensitivity of 1H NMR, 1H chemical shift anisotropy (CSA) is a particularly enticing, and underexplored, dynamics probe. We carried out 1H CSA and 1H-13C dipolar coupling measurements in a series of model supported complexes to understand how 1H CSA can be leveraged to gain dynamic information for heterogeneous catalysts. Mathematical descriptions are given for the dynamic averaging of the CSA tensor, and its dependence on orientation and asymmetry. The variability of the orientation of the tensor in the molecular frame, in addition to its magnitude and asymmetry, negatively impacts attempts to extract quantitative dynamic information. Nevertheless, 1H CSA measurements can reveal useful qualitative insights into the motions of a particularly dilute site, such as from a surface species.

Spatial arrangement of dynamic surface species from solid-state NMR and machine learning-accelerated MD simulations.

The surface arrangement of motional organic functionalities is explored by experimental dipolar coupling measurements and the prediction of motionally-averaged coupling constant from molecular dynamics simulations. The use of machine learning potentials was key to reaching the timescale required. The distance between dynamic surface species are important in cooperative heterogeneous catalysis.

Analysis of KIR gene variants in The Cancer Genome Atlas and UK Biobank using KIRCLE.

BACKGROUND: Natural killer (NK) cells represent a critical component of the innate immune system's response against cancer and viral infections, among other diseases. To distinguish healthy host cells from infected or tumor cells, killer immunoglobulin receptors (KIR) on NK cells bind and recognize Human Leukocyte Antigen (HLA) complexes on their target cells. However, NK cells exhibit great diversity in their mechanism of activation, and the outcomes of their activation are not yet understood fully. Just like the HLAs they bind, KIR receptors exhibit high allelic diversity in the human population. Here we provide a method to identify KIR allele variants from whole exome sequencing data and uncover novel associations between these variants and various molecular and clinical correlates. RESULTS: In order to better understand KIRs, we have developed KIRCLE, a novel method for genotyping individual KIR genes from whole exome sequencing data, and used it to analyze approximately sixty-thousand patient samples in The Cancer Genome Atlas (TCGA) and UK Biobank. We were able to assess population frequencies for different KIR alleles and demonstrate that, similar to HLA alleles, individuals' KIR alleles correlate strongly with their ethnicities. In addition, we observed associations between different KIR alleles and HLA alleles, including HLA-B*53 with KIR3DL2*013 (Fisher's exact FDR = 7.64e-51). Finally, we showcased statistically significant associations between KIR alleles and various clinical correlates, including peptic ulcer disease (Fisher's exact FDR = 0.0429) and age of onset of atopy (Mann-Whitney U FDR = 0.0751). CONCLUSIONS: We show that KIRCLE is able to infer KIR variants accurately and consistently, and we demonstrate its utility using data from approximately sixty-thousand individuals from TCGA and UK Biobank to discover novel molecular and clinical correlations with KIR germline variants. Peptic ulcer disease and atopy are just two diseases in which NK cells may play a role beyond their "classical" realm of anti-tumor and anti-viral responses. This tool may be used both as a benchmark for future KIR-variant-inference algorithms, and to better understand the immunogenomics of and disease processes involving KIRs.

MB-SupCon: Microbiome-based Predictive Models via Supervised Contrastive Learning.

Human microbiome consists of trillions of microorganisms. Microbiota can modulate the host physiology through molecule and metabolite interactions. Integrating microbiome and metabolomics data have the potential to predict different diseases more accurately. Yet, most datasets only measure microbiome data but without paired metabolome data. Here, we propose a novel integrative modeling framework, Microbiome-based Supervised Contrastive Learning Framework (MB-SupCon). MB-SupCon integrates microbiome and metabolome data to generate microbiome embeddings, which can be used to improve the prediction accuracy in datasets that only measure microbiome data. As a proof of concept, we applied MB-SupCon on 720 samples with paired 16S microbiome data and metabolomics data from patients with type 2 diabetes. MB-SupCon outperformed existing prediction methods and achieved high average prediction accuracies for insulin resistance status (84.62%), sex (78.98%), and race (80.04%). Moreover, the microbiome embeddings form separable clusters for different covariate groups in the lower-dimensional space, which enhances data visualization. We also applied MB-SupCon on a large inflammatory bowel disease study and observed similar advantages. Thus, MB-SupCon could be broadly applicable to improve microbiome prediction models in multi-omics disease studies.

Hydrophilic But Not Hydrophobic Surfactant Protein Genetic Variants Are Associated With Severe Acute Respiratory Syncytial Virus Infection in Children.

Respiratory syncytial virus (RSV) is the leading cause of lower respiratory tract infection-related hospitalization in the first year of life. Surfactant dysfunction is central to pathophysiologic mechanisms of various pulmonary diseases including RSV. We hypothesized that RSV severity is associated with single nucleotide polymorphisms (SNPs) of surfactant proteins (SPs). We prospectively enrolled 405 RSV-positive children and divided them into moderate and severe RSV disease. DNA was extracted and genotyped for sixteen specific SP gene SNPs. SP-A1 and A2 haplotypes were assigned. The association of RSV severity with SP gene SNPs was investigated by multivariate logistic regression. A likelihood ratio test was used to test the goodness of fit between two models (one with clinical and demographic data alone and another that included genetic variants). p ≤ 0.05 denotes statistical significance. A molecular dynamics simulation was done to determine the impact of the SFTPA2 rs1965708 on the SP-A behavior under various conditions. Infants with severe disease were more likely to be younger, of lower weight, and exposed to household pets and smoking, as well as having co-infection on admission. A decreased risk of severe RSV was associated with the rs17886395_C of the SFTPA2 and rs2243639_A of the SFTPD, whereas an increased risk was associated with the rs1059047_C of the SFTPA1. RSV severity was not associated with SNPs of SFTPB and SFTPC. An increased risk of severe RSV was associated with the 1A0 genotype of SFTPA2 in its homozygous or heterozygous form with 1A3. A molecular dynamic simulation study of SP-A variants that differ in amino acid 223, an important amino acid change (Q223K) between 1A0 and 1A3, showed no major impact on the behavior of these two variants except for higher thermodynamic stability of the K223 variant. The likelihood ratio test showed that the model with multi-allelic variants along with clinical and demographic data was a better fit to predict RSV severity. In summary, RSV severity was associated with hydrophilic (but not with hydrophobic) SPs gene variants. Collectively, our findings show that SP gene variants may play a key role in RSV infection and have a potential role in prognostication.

Construction and Application of Polygenic Risk Scores in Autoimmune Diseases.

Genome-wide association studies (GWAS) have identified hundreds of genetic variants associated with autoimmune diseases and provided unique mechanistic insights and informed novel treatments. These individual genetic variants on their own typically confer a small effect of disease risk with limited predictive power; however, when aggregated (e.g., via polygenic risk score method), they could provide meaningful risk predictions for a myriad of diseases. In this review, we describe the recent advances in GWAS for autoimmune diseases and the practical application of this knowledge to predict an individual's susceptibility/severity for autoimmune diseases such as systemic lupus erythematosus (SLE) via the polygenic risk score method. We provide an overview of methods for deriving different polygenic risk scores and discuss the strategies to integrate additional information from correlated traits and diverse ancestries. We further advocate for the need to integrate clinical features (e.g., anti-nuclear antibody status) with genetic profiling to better identify patients at high risk of disease susceptibility/severity even before clinical signs or symptoms develop. We conclude by discussing future challenges and opportunities of applying polygenic risk score methods in clinical care.

Reaction processes at step edges on S-decorated Cu(111) and Ag(111) surfaces: MD analysis utilizing machine learning derived potentials.

A variety of complexation, reconstruction, and sulfide formation processes can occur at step edges on the {111} surfaces of coinage metals (M) in the presence of adsorbed S under ultra-high vacuum conditions. Given the cooperative many-atom nature of these reaction processes, Molecular Dynamics (MD) simulation of the associated dynamics is instructive. However, only quite restricted Density Functional Theory (DFT)-level ab initio MD is viable. Thus, for M = Ag and Cu, we instead utilize the DeePMD framework to develop machine-learning derived potentials, retaining near-DFT accuracy for the M-S systems, which should have broad applicability. These potentials are validated by comparison with DFT predictions for various key quantities related to the energetics of S on M(111) surfaces. The potentials are then utilized to perform extensive MD simulations elucidating the above diverse restructuring and reaction processes at step edges. Key observations from MD simulations include the formation of small metal-sulfur complexes, especially MS2; development of a local reconstruction at A-steps featuring an S-decorated {100} motif; and 3D sulfide formation. Additional analysis yields further information on the kinetics for metal-sulfur complex formation, where these complexes can strongly enhance surface mass transport, and on the propensity for sulfide formation.

Integrating 3D genomic and epigenomic data to enhance target gene discovery and drug repurposing in transcriptome-wide association studies.

Transcriptome-wide association studies (TWAS) are popular approaches to test for association between imputed gene expression levels and traits of interest. Here, we propose an integrative method PUMICE (Prediction Using Models Informed by Chromatin conformations and Epigenomics) to integrate 3D genomic and epigenomic data with expression quantitative trait loci (eQTL) to more accurately predict gene expressions. PUMICE helps define and prioritize regions that harbor cis-regulatory variants, which outperforms competing methods. We further describe an extension to our method PUMICE +, which jointly combines TWAS results from single- and multi-tissue models. Across 79 traits, PUMICE + identifies 22% more independent novel genes and increases median chi-square statistics values at known loci by 35% compared to the second-best method, as well as achieves the narrowest credible interval size. Lastly, we perform computational drug repurposing and confirm that PUMICE + outperforms other TWAS methods.