Identification and validation of supervariants reveal novel loci associated with human white matter microstructure.

As an essential part of the central nervous system, white matter coordinates communications between different brain regions and is related to a wide range of neurodegenerative and neuropsychiatric disorders. Previous genome-wide association studies (GWASs) have uncovered loci associated with white matter microstructure. However, GWASs suffer from limited reproducibility and difficulties in detecting multi-single-nucleotide polymorphism (multi-SNP) and epistatic effects. In this study, we adopt the concept of supervariants, a combination of alleles in multiple loci, to account for potential multi-SNP effects. We perform supervariant identification and validation to identify loci associated with 22 white matter fractional anisotropy phenotypes derived from diffusion tensor imaging. To increase reproducibility, we use United Kingdom (UK) Biobank White British (n = 30,842) data for discovery and internal validation, and UK Biobank White but non-British (n = 1927) data, Europeans from the Adolescent Brain Cognitive Development study (n = 4399) data, and Europeans from the Human Connectome Project (n = 319) data for external validation. We identify 23 novel loci on the discovery set that have not been reported in the previous GWASs on white matter microstructure. Among them, three supervariants on genomic regions 5q35.1, 8p21.2, and 19q13.32 have P-values lower than 0.05 in the meta-analysis of the three independent validation data sets. These supervariants contain genetic variants located in genes that have been related to brain structures, cognitive functions, and neuropsychiatric diseases. Our findings provide a better understanding of the genetic architecture underlying white matter microstructure.

A comprehensive patient-specific prediction model for temporomandibular joint osteoarthritis progression.

Temporomandibular joint osteoarthritis (TMJ OA) is a prevalent degenerative disease characterized by chronic pain and impaired jaw function. The complexity of TMJ OA has hindered the development of prognostic tools, posing a significant challenge in timely, patient-specific management. Addressing this gap, our research employs a comprehensive, multidimensional approach to advance TMJ OA prognostication. We conducted a prospective study with 106 subjects, 74 of whom were followed up after 2 to 3 y of conservative treatment. Central to our methodology is the development of an innovative, open-source predictive modeling framework, the Ensemble via Hierarchical Predictions through Nested cross-validation tool (EHPN). This framework synergistically integrates 18 feature selection, statistical, and machine learning methods to yield an accuracy of 0.87, with an area under the ROC curve of 0.72 and an F1 score of 0.82. Our study, beyond technical advancements, emphasizes the global impact of TMJ OA, recognizing its unique demographic occurrence. We highlight key factors influencing TMJ OA progression. Using SHAP analysis, we identified personalized prognostic predictors: lower values of headache, lower back pain, restless sleep, condyle high gray level-GL-run emphasis, articular fossa GL nonuniformity, and long-run low GL emphasis; and higher values of superior joint space, mouth opening, saliva Vascular-endothelium-growth-factor, Matrix-metalloproteinase-7, serum Epithelial-neutrophil-activating-peptide, and age indicate recovery likelihood. Our multidimensional and multimodal EHPN tool enhances clinicians' decision-making, offering a transformative translational infrastructure. The EHPN model stands as a significant contribution to precision medicine, offering a paradigm shift in the management of temporomandibular disorders and potentially influencing broader applications in personalized healthcare.

FPLS-DC: functional partial least squares through distance covariance for imaging genetics.

MOTIVATION: Imaging genetics integrates imaging and genetic techniques to examine how genetic variations influence the function and structure of organs like the brain or heart, providing insights into their impact on behavior and disease phenotypes. The use of organ-wide imaging endophenotypes has increasingly been used to identify potential genes associated with complex disorders. However, analyzing organ-wide imaging data alongside genetic data presents two significant challenges: high dimensionality and complex relationships. To address these challenges, we propose a novel, nonlinear inference framework designed to partially mitigate these issues. RESULTS: We propose a functional partial least squares through distance covariance (FPLS-DC) framework for efficient genome wide analyses of imaging phenotypes. It consists of two components. The first component utilizes the FPLS-derived base functions to reduce image dimensionality while screening genetic markers. The second component maximizes the distance correlation between genetic markers and projected imaging data, which is a linear combination of the FPLS-basis functions, using simulated annealing algorithm. In addition, we proposed an iterative FPLS-DC method based on FPLS-DC framework, which effectively overcomes the influence of inter-gene correlation on inference analysis. We efficiently approximate the null distribution of test statistics using a gamma approximation. Compared to existing methods, FPLS-DC offers computational and statistical efficiency for handling large-scale imaging genetics. In real-world applications, our method successfully detected genetic variants associated with the hippocampus, demonstrating its value as a statistical toolbox for imaging genetic studies. AVAILABILITY AND IMPLEMENTATION: The FPLS-DC method we propose opens up new research avenues and offers valuable insights for analyzing functional and high-dimensional data. In addition, it serves as a useful tool for scientific analysis in practical applications within the field of imaging genetics research. The R package FPLS-DC is available in Github: https://github.com/BIG-S2/FPLSDC.

The interaction effects of age, APOE and common environmental risk factors on human brain structure.

Mounting evidence suggests considerable diversity in brain aging trajectories, primarily arising from the complex interplay between age, genetic, and environmental risk factors, leading to distinct patterns of micro- and macro-cerebral aging. The underlying mechanisms of such effects still remain unclear. We conducted a comprehensive association analysis between cerebral structural measures and prevalent risk factors, using data from 36,969 UK Biobank subjects aged 44-81. Participants were assessed for brain volume, white matter diffusivity, Apolipoprotein E (APOE) genotypes, polygenic risk scores, lifestyles, and socioeconomic status. We examined genetic and environmental effects and their interactions with age and sex, and identified 726 signals, with education, alcohol, and smoking affecting most brain regions. Our analysis revealed negative age-APOE-ε4 and positive age-APOE-ε2 interaction effects, respectively, especially in females on the volume of amygdala, positive age-sex-APOE-ε4 interaction on the cerebellar volume, positive age-excessive-alcohol interaction effect on the mean diffusivity of the splenium of the corpus callosum, positive age-healthy-diet interaction effect on the paracentral volume, and negative APOE-ε4-moderate-alcohol interaction effects on the axial diffusivity of the superior fronto-occipital fasciculus. These findings highlight the need of considering age, sex, genetic, and environmental joint effects in elucidating normal or abnormal brain aging.

Merging or ensembling: integrative analysis in multiple neuroimaging studies.

The aim of this paper is to systematically investigate merging and ensembling methods for spatially varying coefficient mixed effects models (SVCMEM) in order to carry out integrative learning of neuroimaging data obtained from multiple biomedical studies. The "merged" approach involves training a single learning model using a comprehensive dataset that encompasses information from all the studies. Conversely, the "ensemble" approach involves creating a weighted average of distinct learning models, each developed from an individual study. We systematically investigate the prediction accuracy of the merged and ensemble learners under the presence of different degrees of interstudy heterogeneity. Additionally, we establish asymptotic guidelines for making strategic decisions about when to employ either of these models in different scenarios, along with deriving optimal weights for the ensemble learner. To validate our theoretical results, we perform extensive simulation studies. The proposed methodology is also applied to 3 large-scale neuroimaging studies.

Urinary metabolite concentrations of phthalate and plasticizers in infancy and childhood in the UNC baby connectome project.

INTRODUCTION: Existing evidence suggests that exposure to phthalates is higher among younger age groups. However, limited knowledge exists on how phthalate exposure, as well as exposure to replacement plasticizers, di(isononyl) cyclohexane-1,2-dicarboxylate (DINCH) and di-2-ethylhexyl terephthalate (DEHTP), change from infancy through early childhood. METHODS: Urine samples were collected across the first 5 years of life from typically developing infants and young children enrolled between 2017 and 2020 in the longitudinal UNC Baby Connectome Project. From 438 urine samples among 187 participants, we quantified concentrations of monobutyl phthalate (MnBP), mono-3-carboxypropyl phthalate (MCPP), monoisobutyl phthalate (MiBP), monoethyl phthalate (MEP), monobenzyl phthalate (MBzP), and metabolites of di(2-ethylhexyl) phthalate (DEHP), diisonoyl phthalate (DiNP), DINCH and DEHTP. Specific gravity (SG) adjusted metabolite and molar sum concentrations were compared across age groups. Intraclass correlation coefficients (ICCs) were calculated among 122 participants with multiple urine specimens (373 samples). RESULTS: Most phthalate metabolites showed high detection frequencies (>80% of samples). Replacement plasticizers DINCH (58-60%) and DEHTP (>97%) were also commonly found. DiNP metabolites were less frequently detected (<10%). For some metabolites, SG-adjusted concentrations were inversely associated with age, with the highest concentrations found in the first year of life. ICCs revealed low to moderate reliability in metabolite measurements (ρ = 0.10-0.48) suggesting a high degree of within-individual variation in exposure among this age group. The first 6 months (compared to remaining age groups) showed an increased ratio of carboxylated metabolites of DEHP and DEHTP, compared to other common metabolites, but no clear age trends for DINCH metabolite ratios were observed. CONCLUSION: Metabolites of phthalates and replacements plasticizers were widely detected in infancy and early childhood, with the highest concentrations observed in the first year of life for several metabolites. Higher proportions of carboxylated metabolites of DEHP and DEHTP in younger age groups indicate potential differences in metabolism during infancy.

Source-free unsupervised domain adaptation: A survey.

Unsupervised domain adaptation (UDA) via deep learning has attracted appealing attention for tackling domain-shift problems caused by distribution discrepancy across different domains. Existing UDA approaches highly depend on the accessibility of source domain data, which is usually limited in practical scenarios due to privacy protection, data storage and transmission cost, and computation burden. To tackle this issue, many source-free unsupervised domain adaptation (SFUDA) methods have been proposed recently, which perform knowledge transfer from a pre-trained source model to the unlabeled target domain with source data inaccessible. A comprehensive review of these works on SFUDA is of great significance. In this paper, we provide a timely and systematic literature review of existing SFUDA approaches from a technical perspective. Specifically, we categorize current SFUDA studies into two groups, i.e., white-box SFUDA and black-box SFUDA, and further divide them into finer subcategories based on different learning strategies they use. We also investigate the challenges of methods in each subcategory, discuss the advantages/disadvantages of white-box and black-box SFUDA methods, conclude the commonly used benchmark datasets, and summarize the popular techniques for improved generalizability of models learned without using source data. We finally discuss several promising future directions in this field.

Disease severity drives risk of venous thrombotic events in women with sickle cell disease in a single-center retrospective study.

Background: Estrogen-containing hormonal contraception (HC) is a well-established risk factor for venous thromboembolism (VTE). Women with sickle cell disease (SCD) also have an increased risk of VTE. However, it is unknown if exposure to HC exacerbates the risk of VTE in women with SCD. Objectives: Assess the impact of HC on VTE risk in women with SCD and explore additional risk factors contributing to VTE development. Methods: We analyzed a retrospective cohort of women of reproductive age (15-49 years) with SCD at the University of North Carolina from 2010 to 2022. Results: We identified 370 women with SCD, and 93 (25.1%) had a history of VTE. Among 219 women exposed to HC, 38 of 184 (20.6%) had a VTE while actively using HC, whereas 20 of 151 (13.2%) women never exposed to HC had a VTE. Of the patients exposed to HC, 64 of 184 (34.7%) were on estrogen-containing HC, with 120 of 184 (65.3%) using progestin-only formulations. Cox regression analysis found that progestin-only formulations increased VTE risk (hazard ratio: 2.03; 95% CI: 1.107-3.726, P < .05). However, when accounting for disease severity, the association between progestin-only treatment and VTE risk was not significant. Indeed, a nuanced analysis revealed that both severe (odds ratio: 11.79; 95% CI: 5.14-27.06; P < .001) and moderate (odds ratio: 4.37; 95% CI: 1.77-10.76; P = .001) disease increased risk compared with mild disease. Neither genotype nor hydroxyurea use influenced VTE risk. Conclusion: Overall, we found that increased thrombotic risk is more likely influenced by disease status than HC exposure and should play a role in shared decision-making with patients.

The Interaction Effects of Age, Sex, APOE and Common Health Risk Factors on Human Brain Functions.

Recent studies have shed light on the complex nonlinear changes in brain functions across the lifespan, demonstrating the variability in the individual cognitive and neural development during aging. This variability is influenced by factors such as sex, age, genetics, and modifiable health risk factors (MHRFs), which collectively shape unique patterns of brain functional connectivities (FCs) across different regions. However, their joint effects and underlying mechanisms remain unclear. We conduct a comprehensive analysis to jointly examine the association of common risk factors with brain functional measures, using data from 36,630 UK Biobank participants aged 44-81. Participants were assessed for age, sex, Apolipoprotein E ( APOE ) genotypes, ten common MHRFs, and brain FCs measured via resting-state functional magnetic resonance imaging. Using the fine-grained HCP-MMP parcellation and Ji-12 network atlases, we identified 91 associations with network functional connectivity (NFC) and 102 associations with network edge strength (NES) measures. Hypertension, BMI, and education emerged as the top three influential factors across networks. Notably, a negative interaction between sex and APOE-ε4 ( APOE4 ) genotype was observed, with male APOE4 carriers showing greater reductions in NFC between the cingulo- opercular (CON) and posterior multimodal (PMN) networks. Additionally, a negative age-BMI interaction on NES between the visual and dorsal attention (DAN) networks suggested that higher BMI accelerates the decline in visual-DAN connectivity. A positive age-hypertension interaction between the frontoparietal (FPN) and default mode (DMN) networks indicated a more rapid decrease in functional segregation associated with hypertension. We also identified sex- education interactions, showing more pronounced positive effects on CON-FPN networks in females and PMN-DMN networks in males. Further interactions involving sex and other MHRFs, such as smoking, alcohol consumption, diabetes, and BMI, revealed that smoking, alcohol, and BMI had more detrimental effects in males, while diabetes had a more pronounced negative impact in females within specific networks. These findings underscore the necessity of jointly considering sex, age, genetic factors, and MHRFs to accurately delineate the multifactorial alterations in the FCs during brain aging.

Alzheimer's Disease Knowledge Graph Enhances Knowledge Discovery and Disease Prediction.

Background: Alzheimer's disease (AD), a progressive neurodegenerative disorder, continues to increase in prevalence without any effective treatments to date. In this context, knowledge graphs (KGs) have emerged as a pivotal tool in biomedical research, offering new perspectives on drug repurposing and biomarker discovery by analyzing intricate network structures. Our study seeks to build an AD-specific knowledge graph, highlighting interactions among AD, genes, variants, chemicals, drugs, and other diseases. The goal is to shed light on existing treatments, potential targets, and diagnostic methods for AD, thereby aiding in drug repurposing and the identification of biomarkers. Results: We annotated 800 PubMed abstracts and leveraged GPT-4 for text augmentation to enrich our training data for named entity recognition (NER) and relation classification. A comprehensive data mining model, integrating NER and relationship classification, was trained on the annotated corpus. This model was subsequently applied to extract relation triplets from unannotated abstracts. To enhance entity linking, we utilized a suite of reference biomedical databases and refine the linking accuracy through abbreviation resolution. As a result, we successfully identified 3,199,276 entity mentions and 633,733 triplets, elucidating connections between 5,000 unique entities. These connections were pivotal in constructing a comprehensive Alzheimer's Disease Knowledge Graph (ADKG). We also integrated the ADKG constructed after entity linking with other biomedical databases. The ADKG served as a training ground for Knowledge Graph Embedding models with the high-ranking predicted triplets supported by evidence, underscoring the utility of ADKG in generating testable scientific hypotheses. Further application of ADKG in predictive modeling using the UK Biobank data revealed models based on ADKG outperforming others, as evidenced by higher values in the areas under the receiver operating characteristic (ROC) curves. Conclusion: The ADKG is a valuable resource for generating hypotheses and enhancing predictive models, highlighting its potential to advance AD's disease research and treatment strategies.

Leveraging Brain Modularity Prior for Interpretable Representation Learning of fMRI.

Resting-state functional magnetic resonance imaging (rs-fMRI) can reflect spontaneous neural activities in the brain and is widely used for brain disorder analysis. Previous studies focus on extracting fMRI representations using machine/deep learning methods, but these features typically lack biological interpretability. The human brain exhibits a remarkable modular structure in spontaneous brain functional networks, with each module comprised of functionally interconnected brain regions-of-interest (ROIs). However, existing learning-based methods cannot adequately utilize such brain modularity prior. In this paper, we propose a brain modularity-constrained dynamic representation learning framework for interpretable fMRI analysis, consisting of dynamic graph construction, dynamic graph learning via a novel modularity-constrained graph neural network (MGNN), and prediction and biomarker detection. The designed MGNN is constrained by three core neurocognitive modules (i.e., salience network, central executive network, and default mode network), encouraging ROIs within the same module to share similar representations. To further enhance discriminative ability of learned features, we encourage the MGNN to preserve network topology of input graphs via a graph topology reconstruction constraint. Experimental results on 534 subjects with rs-fMRI scans from two datasets validate the effectiveness of the proposed method. The identified discriminative brain ROIs and functional connectivities can be regarded as potential fMRI biomarkers to aid in clinical diagnosis.

Multi-organ imaging-derived polygenic indexes for brain and body health.

The UK Biobank (UKB) imaging project is a crucial resource for biomedical research, but is limited to 100,000 participants due to cost and accessibility barriers. Here we used genetic data to predict heritable imaging-derived phenotypes (IDPs) for a larger cohort. We developed and evaluated 4,375 IDP genetic scores (IGS) derived from UKB brain and body images. When applied to UKB participants who were not imaged, IGS revealed links to numerous phenotypes and stratified participants at increased risk for both brain and somatic diseases. For example, IGS identified individuals at higher risk for Alzheimer's disease and multiple sclerosis, offering additional insights beyond traditional polygenic risk scores of these diseases. When applied to independent external cohorts, IGS also stratified those at high disease risk in the All of Us Research Program and the Alzheimer's Disease Neuroimaging Initiative study. Our results demonstrate that, while the UKB imaging cohort is largely healthy and may not be the most enriched for disease risk management, it holds immense potential for stratifying the risk of various brain and body diseases in broader external genetic cohorts.

Advancing microbial production through artificial intelligence-aided biology.

Microbial cell factories (MCFs) have been leveraged to construct sustainable platforms for value-added compound production. To optimize metabolism and reach optimal productivity, synthetic biology has developed various genetic devices to engineer microbial systems by gene editing, high-throughput protein engineering, and dynamic regulation. However, current synthetic biology methodologies still rely heavily on manual design, laborious testing, and exhaustive analysis. The emerging interdisciplinary field of artificial intelligence (AI) and biology has become pivotal in addressing the remaining challenges. AI-aided microbial production harnesses the power of processing, learning, and predicting vast amounts of biological data within seconds, providing outputs with high probability. With well-trained AI models, the conventional Design-Build-Test (DBT) cycle has been transformed into a multidimensional Design-Build-Test-Learn-Predict (DBTLP) workflow, leading to significantly improved operational efficiency and reduced labor consumption. Here, we comprehensively review the main components and recent advances in AI-aided microbial production, focusing on genome annotation, AI-aided protein engineering, artificial functional protein design, and AI-enabled pathway prediction. Finally, we discuss the challenges of integrating novel AI techniques into biology and propose the potential of large language models (LLMs) in advancing microbial production.

Eye-brain connections revealed by multimodal retinal and brain imaging genetics.

The retina, an anatomical extension of the brain, forms physiological connections with the visual cortex of the brain. Although retinal structures offer a unique opportunity to assess brain disorders, their relationship to brain structure and function is not well understood. In this study, we conducted a systematic cross-organ genetic architecture analysis of eye-brain connections using retinal and brain imaging endophenotypes. We identified novel phenotypic and genetic links between retinal imaging biomarkers and brain structure and function measures from multimodal magnetic resonance imaging (MRI), with many associations involving the primary visual cortex and visual pathways. Retinal imaging biomarkers shared genetic influences with brain diseases and complex traits in 65 genomic regions, with 18 showing genetic overlap with brain MRI traits. Mendelian randomization suggests bidirectional genetic causal links between retinal structures and neurological and neuropsychiatric disorders, such as Alzheimer's disease. Overall, our findings reveal the genetic basis for eye-brain connections, suggesting that retinal images can help uncover genetic risk factors for brain disorders and disease-related changes in intracranial structure and function.

Pan-cancer subclonal mutation analysis of 7,827 tumors predicts clinical outcome.

Intra-tumor heterogeneity is an important driver of tumor evolution and therapy response. Advances in precision cancer treatment will require understanding of mutation clonality and subclonal architecture. Currently the slow computational speed of subclonal reconstruction hinders large cohort studies. To overcome this bottleneck, we developed Clonal structure identification through Pairwise Penalization, or CliPP, which clusters subclonal mutations using a regularized likelihood model. CliPP reliably processed whole-genome and whole-exome sequencing data from over 12,000 tumor samples within 24 hours, thus enabling large-scale downstream association analyses between subclonal structures and clinical outcomes. Through a pan-cancer investigation of 7,827 tumors from 32 cancer types, we found that high subclonal mutational load (sML), a measure of latency time in tumor evolution, was significantly associated with better patient outcomes in 16 cancer types with low to moderate tumor mutation burden (TMB). In a cohort of prostate cancer patients participating in an immunotherapy clinical trial, high sML was indicative of favorable response to immune checkpoint blockade. This comprehensive study using CliPP underscores sML as a key feature of cancer. sML may be essential for linking mutation dynamics with immunotherapy response in the large population of non-high TMB cancers.

Structural MRI Harmonization via Disentangled Latent Energy-Based Style Translation.

Multi-site brain magnetic resonance imaging (MRI) has been widely used in clinical and research domains, but usually is sensitive to non-biological variations caused by site effects (e.g., field strengths and scanning protocols). Several retrospective data harmonization methods have shown promising results in removing these non-biological variations at feature or whole-image level. Most existing image-level harmonization methods are implemented through generative adversarial networks, which are generally computationally expensive and generalize poorly on independent data. To this end, this paper proposes a disentangled latent energy-based style translation (DLEST) framework for image-level structural MRI harmonization. Specifically, DLEST disentangles site-invariant image generation and site-specific style translation via a latent autoencoder and an energy-based model. The autoencoder learns to encode images into low-dimensional latent space, and generates faithful images from latent codes. The energy-based model is placed in between the encoding and generation steps, facilitating style translation from a source domain to a target domain implicitly. This allows highly generalizable image generation and efficient style translation through the latent space. We train our model on 4,092 T1-weighted MRIs in 3 tasks: histogram comparison, acquisition site classification, and brain tissue segmentation. Qualitative and quantitative results demonstrate the superiority of our approach, which generally outperforms several state-of-the-art methods.

An atlas of trait associations with resting-state and task-evoked human brain functional organizations in the UK Biobank.

Functional magnetic resonance imaging (fMRI) has been widely used to identify brain regions linked to critical functions, such as language and vision, and to detect tumors, strokes, brain injuries, and diseases. It is now known that large sample sizes are necessary for fMRI studies to detect small effect sizes and produce reproducible results. Here we report a systematic association analysis of 647 traits with imaging features extracted from resting-state and task-evoked fMRI data of more than 40,000 UK Biobank participants. We used a parcellation-based approach to generate 64,620 functional connectivity measures to reveal fine-grained details about cerebral cortex functional organizations. The difference between functional organizations at rest and during task was examined, and we have prioritized important brain regions and networks associated with a variety of human traits and clinical outcomes. For example, depression was most strongly associated with decreased connectivity in the somatomotor network. We have made our results publicly available and developed a browser framework to facilitate the exploration of brain function-trait association results (http://fmriatlas.org/).