Dissection and integration of bursty transcriptional dynamics for complex systems.

RNA velocity estimation is a potentially powerful tool to reveal the directionality of transcriptional changes in single-cell RNA-sequencing data, but it lacks accuracy, absent advanced metabolic labeling techniques. We developed an approach, TopicVelo, that disentangles simultaneous, yet distinct, dynamics by using a probabilistic topic model, a highly interpretable form of latent space factorization, to infer cells and genes associated with individual processes, thereby capturing cellular pluripotency or multifaceted functionality. Focusing on process-associated cells and genes enables accurate estimation of process-specific velocities via a master equation for a transcriptional burst model accounting for intrinsic stochasticity. The method obtains a global transition matrix by leveraging cell topic weights to integrate process-specific signals. In challenging systems, this method accurately recovers complex transitions and terminal states, while our use of first-passage time analysis provides insights into transient transitions. These results expand the limits of RNA velocity, empowering future studies of cell fate and functional responses.

Coexistence of many species under a random competition-colonization trade-off.

The competition-colonization (CC) trade-off is a well-studied coexistence mechanism for metacommunities. In this setting, it is believed that the coexistence of all species requires their traits to satisfy restrictive conditions limiting their similarity. To investigate whether diverse metacommunities can assemble in a CC trade-off model, we study their assembly from a probabilistic perspective. From a pool of species with parameters (corresponding to traits) sampled at random, we compute the probability that any number of species coexist and characterize the set of species that emerges through assembly. Remarkably, almost exactly half of the species in a large pool typically coexist, with no saturation as the size of the pool grows, and with little dependence on the underlying distribution of traits. Through a mix of analytical results and simulations, we show that this unlimited niche packing emerges as assembly actively moves communities toward overdispersed configurations in niche space. Our findings also apply to a realistic assembly scenario where species invade one at a time from a fixed regional pool. When diversity arises de novo in the metacommunity, richness still grows without bound, but more slowly. Together, our results suggest that the CC trade-off can support the robust emergence of diverse communities, even when coexistence of the full species pool is exceedingly unlikely.

Evaluating evasion strategies in zebrafish larvae.

An effective evasion strategy allows prey to survive encounters with predators. Prey are generally thought to escape in a direction that is either random or serves to maximize the minimum distance from the predator. Here, we introduce a comprehensive approach to determine the most likely evasion strategy among multiple hypotheses and the role of biomechanical constraints on the escape response of prey fish. Through a consideration of six strategies with sensorimotor noise and previous kinematic measurements, our analysis shows that zebrafish larvae generally escape in a direction orthogonal to the predator's heading. By sensing only the predator's heading, this orthogonal strategy maximizes the distance from fast-moving predators, and, when operating within the biomechanical constraints of the escape response, it provides the best predictions of prey behavior among all alternatives. This work demonstrates a framework for resolving the strategic basis of evasion in predator-prey interactions, which could be applied to a broad diversity of animals.

Hyperedge prediction and the statistical mechanisms of higher-order and lower-order interactions in complex networks.

Complex networked systems often exhibit higher-order interactions, beyond dyadic interactions, which can dramatically alter their observed behavior. Consequently, understanding hypergraphs from a structural perspective has become increasingly important. Statistical, group-based inference approaches are well suited for unveiling the underlying community structure and predicting unobserved interactions. However, these approaches often rely on two key assumptions: that the same groups can explain hyperedges of any order and that interactions are assortative, meaning that edges are formed by nodes with the same group memberships. To test these assumptions, we propose a group-based generative model for hypergraphs that does not impose an assortative mechanism to explain observed higher-order interactions, unlike current approaches. Our model allows us to explore the validity of the assumptions. Our results indicate that the first assumption appears to hold true for real networks. However, the second assumption is not necessarily accurate; we find that a combination of general statistical mechanisms can explain observed hyperedges. Finally, with our approach, we are also able to determine the importance of lower and high-order interactions for predicting unobserved interactions. Our research challenges the conventional assumptions of group-based inference methodologies and broadens our understanding of the underlying structure of hypergraphs.

Dissociable Contributions of Basolateral Amygdala and Ventrolateral Orbitofrontal Cortex to Flexible Learning Under Uncertainty.

Reversal learning measures the ability to form flexible associations between choice outcomes with stimuli and actions that precede them. This type of learning is thought to rely on several cortical and subcortical areas, including the highly interconnected orbitofrontal cortex (OFC) and basolateral amygdala (BLA), and is often impaired in various neuropsychiatric and substance use disorders. However, the unique contributions of these regions to stimulus- and action-based reversal learning have not been systematically compared using a chemogenetic approach particularly before and after the first reversal that introduces new uncertainty. Here, we examined the roles of ventrolateral OFC (vlOFC) and BLA during reversal learning. Male and female rats were prepared with inhibitory designer receptors exclusively activated by designer drugs targeting projection neurons in these regions and tested on a series of deterministic and probabilistic reversals during which they learned about stimulus identity or side (left or right) associated with different reward probabilities. Using a counterbalanced within-subject design, we inhibited these regions prior to reversal sessions. We assessed initial and pre-/post-reversal changes in performance to measure learning and adjustments to reversals, respectively. We found that inhibition of the ventrolateral orbitofrontal cortex (vlOFC), but not BLA, eliminated adjustments to stimulus-based reversals. Inhibition of BLA, but not vlOFC, selectively impaired action-based probabilistic reversal learning, leaving deterministic reversal learning intact. vlOFC exhibited a sex-dependent role in early adjustment to action-based reversals, but not in overall learning. These results reveal dissociable roles for BLA and vlOFC in flexible learning and highlight a more crucial role for BLA in learning meaningful changes in the reward environment.

Deep autoregressive generative models capture the intrinsics embedded in T-cell receptor repertoires.

T-cell receptors (TCRs) play an essential role in the adaptive immune system. Probabilistic models for TCR repertoires can help decipher the underlying complex sequence patterns and provide novel insights into understanding the adaptive immune system. In this work, we develop TCRpeg, a deep autoregressive generative model to unravel the sequence patterns of TCR repertoires. TCRpeg largely outperforms state-of-the-art methods in estimating the probability distribution of a TCR repertoire, boosting the average accuracy from 0.672 to 0.906 measured by the Pearson correlation coefficient. Furthermore, with promising performance in probability inference, TCRpeg improves on a range of TCR-related tasks: profiling TCR repertoire probabilistically, classifying antigen-specific TCRs, validating previously discovered TCR motifs, generating novel TCRs and augmenting TCR data. Our results and analysis highlight the flexibility and capacity of TCRpeg to extract TCR sequence information, providing a novel approach for deciphering complex immunogenomic repertoires.

Probabilistic forecasts of international bilateral migration flows.

We propose a method for forecasting global human migration flows. A Bayesian hierarchical model is used to make probabilistic projections of the 39,800 bilateral migration flows among the 200 most populous countries. We generate out-of-sample forecasts for all bilateral flows for the 2015 to 2020 period, using models fitted to bilateral migration flows for five 5-y periods from 1990 to 1995 through 2010 to 2015. We find that the model produces well-calibrated out-of-sample forecasts of bilateral flows, as well as total country-level inflows, outflows, and net flows. The mean absolute error decreased by 61% using our method, compared to a leading model of international migration. Out-of-sample analysis indicated that simple methods for forecasting migration flows offered accurate projections of bilateral migration flows in the near term. Our method matched or improved on the out-of-sample performance using these simple deterministic alternatives, while also accurately assessing uncertainty. We integrate the migration flow forecasting model into a fully probabilistic population projection model to generate bilateral migration flow forecasts by age and sex for all flows from 2020 to 2025 through 2040 to 2045.

Assessing transcriptomic heterogeneity of single-cell RNASeq data by bulk-level gene expression data.

BACKGROUND: Single-cell RNA sequencing (sc-RNASeq) data illuminate transcriptomic heterogeneity but also possess a high level of noise, abundant missing entries and sometimes inadequate or no cell type annotations at all. Bulk-level gene expression data lack direct information of cell population composition but are more robust and complete and often better annotated. We propose a modeling framework to integrate bulk-level and single-cell RNASeq data to address the deficiencies and leverage the mutual strengths of each type of data and enable a more comprehensive inference of their transcriptomic heterogeneity. Contrary to the standard approaches of factorizing the bulk-level data with one algorithm and (for some methods) treating single-cell RNASeq data as references to decompose bulk-level data, we employed multiple deconvolution algorithms to factorize the bulk-level data, constructed the probabilistic graphical models of cell-level gene expressions from the decomposition outcomes, and compared the log-likelihood scores of these models in single-cell data. We term this framework backward deconvolution as inference operates from coarse-grained bulk-level data to fine-grained single-cell data. As the abundant missing entries in sc-RNASeq data have a significant effect on log-likelihood scores, we also developed a criterion for inclusion or exclusion of zero entries in log-likelihood score computation. RESULTS: We selected nine deconvolution algorithms and validated backward deconvolution in five datasets. In the in-silico mixtures of mouse sc-RNASeq data, the log-likelihood scores of the deconvolution algorithms were strongly anticorrelated with their errors of mixture coefficients and cell type specific gene expression signatures. In the true bulk-level mouse data, the sample mixture coefficients were unknown but the log-likelihood scores were strongly correlated with accuracy rates of inferred cell types. In the data of autism spectrum disorder (ASD) and normal controls, we found that ASD brains possessed higher fractions of astrocytes and lower fractions of NRGN-expressing neurons than normal controls. In datasets of breast cancer and low-grade gliomas (LGG), we compared the log-likelihood scores of three simple hypotheses about the gene expression patterns of the cell types underlying the tumor subtypes. The model that tumors of each subtype were dominated by one cell type persistently outperformed an alternative model that each cell type had elevated expression in one gene group and tumors were mixtures of those cell types. Superiority of the former model is also supported by comparing the real breast cancer sc-RNASeq clusters with those generated by simulated sc-RNASeq data. CONCLUSIONS: The results indicate that backward deconvolution serves as a sensible model selection tool for deconvolution algorithms and facilitates discerning hypotheses about cell type compositions underlying heterogeneous specimens such as tumors.

kalis: a modern implementation of the Li & Stephens model for local ancestry inference in R.

BACKGROUND: Approximating the recent phylogeny of N phased haplotypes at a set of variants along the genome is a core problem in modern population genomics and central to performing genome-wide screens for association, selection, introgression, and other signals. The Li & Stephens (LS) model provides a simple yet powerful hidden Markov model for inferring the recent ancestry at a given variant, represented as an N × N distance matrix based on posterior decodings. RESULTS: We provide a high-performance engine to make these posterior decodings readily accessible with minimal pre-processing via an easy to use package kalis, in the statistical programming language R. kalis enables investigators to rapidly resolve the ancestry at loci of interest and developers to build a range of variant-specific ancestral inference pipelines on top. kalis exploits both multi-core parallelism and modern CPU vector instruction sets to enable scaling to hundreds of thousands of genomes. CONCLUSIONS: The resulting distance matrices accessible via kalis enable local ancestry, selection, and association studies in modern large scale genomic datasets.

Connectome alterations following perinatal deafness in the cat.

Following sensory deprivation, areas and networks in the brain may adapt and reorganize to compensate for the loss of input. These adaptations are manifestations of compensatory crossmodal plasticity, which has been documented in both human and animal models of deafness-including the domestic cat. Although there are abundant examples of structural plasticity in deaf felines from retrograde tracer-based studies, there is a lack of diffusion-based knowledge involving this model compared to the current breadth of human research. The purpose of this study was to explore white matter structural adaptations in the perinatally-deafened cat via tractography, increasing the methodological overlap between species. Plasticity was examined by identifying unique group connections and assessing altered connectional strength throughout the entirety of the brain. Results revealed a largely preserved connectome containing a limited number of group-specific or altered connections focused within and between sensory networks, which is generally corroborated by deaf feline anatomical tracer literature. Furthermore, five hubs of cortical plasticity and altered communication following perinatal deafness were observed. The limited differences found in the present study suggest that deafness-induced crossmodal plasticity is largely built upon intrinsic structural connections, with limited remodeling of underlying white matter.

Characterizing multimorbidity in ALIVE: comparing single and ensemble clustering methods.

Multimorbidity, defined as having 2 or more chronic conditions, is a growing public health concern, but research in this area is complicated by the fact that multimorbidity is a highly heterogenous outcome. Individuals in a sample may have a differing number and varied combinations of conditions. Clustering methods, such as unsupervised machine learning algorithms, may allow us to tease out the unique multimorbidity phenotypes. However, many clustering methods exist, and choosing which to use is challenging because we do not know the true underlying clusters. Here, we demonstrate the use of 3 individual algorithms (partition around medoids, hierarchical clustering, and probabilistic clustering) and a clustering ensemble approach (which pools different clustering approaches) to identify multimorbidity clusters in the AIDS Linked to the Intravenous Experience cohort study. We show how the clusters can be compared based on cluster quality, interpretability, and predictive ability. In practice, it is critical to compare the clustering results from multiple algorithms and to choose the approach that performs best in the domain(s) that aligns with plans to use the clusters in future analyses.

Could differential under-reporting of loneliness between men and women bias the gender-specific association between loneliness duration and rate of memory decline? A probabilistic bias analysis of effect modification.

Gender is an observed effect modifier of the association between loneliness and memory aging. However, this effect modification may be a result of information bias due to differential loneliness under-reporting by gender. We applied probabilistic bias analyses to examine whether effect modification of the loneliness-memory decline relationship by gender is retained under three simulation scenarios with various magnitudes of differential loneliness under-reporting between men and women. Data were from biennial interviews with adults aged 50+ in the US Health and Retirement Study from 1996-2016 (5,646 women and 3,386 men). Loneliness status (yes vs. no) was measured from 1996-2004 using the CES-D loneliness item and memory was measured from 2004-2016. Simulated sensitivity and specificity of the loneliness measure were informed by a validation study using the UCLA Loneliness Scale as a gold standard. The likelihood of observing effect modification by gender was higher than 90% in all simulations, although the likelihood reduced with an increasing difference in magnitude of the loneliness under-reporting between men and women. The gender difference in loneliness under-reporting did not meaningfully affect the observed effect modification by gender in our simulations. Our simulation approach may be promising to quantify potential information bias in effect modification analyses.

Probabilistic colocalization of genetic variants from complex and molecular traits: promise and limitations.

Colocalization analysis has emerged as a powerful tool to uncover the overlapping of causal variants responsible for both molecular and complex disease phenotypes. The findings from colocalization analysis yield insights into the molecular pathways of complex diseases. In this paper, we conduct an in-depth investigation of the promise and limitations of the available colocalization analysis approaches. Focusing on variant-level colocalization approaches, we first establish the connections between various existing methods. We proceed to discuss the impacts of various controllable analytical factors and uncontrollable practical factors on outcomes of colocalization analysis through realistic simulations and real data examples. We identify a single analytical factor, the specification of prior enrichment levels, which can lead to severe inflation of false-positive colocalization findings. Meanwhile, the combination of many other analytical and practical factors all lead to diminished power. Consequently, we recommend the following strategies for the best practice of colocalization analysis: (1) estimating prior enrichment level from the observed data and (2) separating fine-mapping and colocalization analysis. Our analysis of 4,091 complex traits and the multi-tissue expression quantitative trait loci (eQTL) data from the GTEx (v.8) suggests that colocalizations of molecular QTLs and causal complex trait associations are widespread. However, only a small proportion can be confidently identified from currently available data due to a lack of power. Our findings set a benchmark for current and future integrative genetic association analysis applications.

Multi-trait transcriptome-wide association studies with probabilistic Mendelian randomization.

A transcriptome-wide association study (TWAS) integrates data from genome-wide association studies and gene expression mapping studies for investigating the gene regulatory mechanisms underlying diseases. Existing TWAS methods are primarily univariate in nature, focusing on analyzing one outcome trait at a time. However, many complex traits are correlated with each other and share a common genetic basis. Consequently, analyzing multiple traits jointly through multivariate analysis can potentially improve the power of TWASs. Here, we develop a method, moPMR-Egger (multiple outcome probabilistic Mendelian randomization with Egger assumption), for analyzing multiple outcome traits in TWAS applications. moPMR-Egger examines one gene at a time, relies on its cis-SNPs that are in potential linkage disequilibrium with each other to serve as instrumental variables, and tests its causal effects on multiple traits jointly. A key feature of moPMR-Egger is its ability to test and control for potential horizontal pleiotropic effects from instruments, thus maximizing power while minimizing false associations for TWASs. In simulations, moPMR-Egger provides calibrated type I error control for both causal effects testing and horizontal pleiotropic effects testing and is more powerful than existing univariate TWAS approaches in detecting causal associations. We apply moPMR-Egger to analyze 11 traits from 5 trait categories in the UK Biobank. In the analysis, moPMR-Egger identified 13.15% more gene associations than univariate approaches across trait categories and revealed distinct regulatory mechanisms underlying systolic and diastolic blood pressures.

A comprehensive comparison on cell-type composition inference for spatial transcriptomics data.

Spatial transcriptomics (ST) technologies allow researchers to examine transcriptional profiles along with maintained positional information. Such spatially resolved transcriptional characterization of intact tissue samples provides an integrated view of gene expression in its natural spatial and functional context. However, high-throughput sequencing-based ST technologies cannot yet reach single cell resolution. Thus, similar to bulk RNA-seq data, gene expression data at ST spot-level reflect transcriptional profiles of multiple cells and entail the inference of cell-type composition within each ST spot for valid and powerful subsequent analyses. Realizing the critical importance of cell-type decomposition, multiple groups have developed ST deconvolution methods. The aim of this work is to review state-of-the-art methods for ST deconvolution, comparing their strengths and weaknesses. In particular, we construct ST spots from single-cell level ST data to assess the performance of 10 methods, with either ideal reference or non-ideal reference. Furthermore, we examine the performance of these methods on spot- and bead-level ST data by comparing estimated cell-type proportions to carefully matched single-cell ST data. In comparing the performance on various tissues and technological platforms, we concluded that RCTD and stereoscope achieve more robust and accurate inferences.

RdRp-based sensitive taxonomic classification of RNA viruses for metagenomic data.

With advances in library construction protocols and next-generation sequencing technologies, viral metagenomic sequencing has become the major source for novel virus discovery. Conducting taxonomic classification for metagenomic data is an important means to characterize the viral composition in the underlying samples. However, RNA viruses are abundant and highly diverse, jeopardizing the sensitivity of comparison-based classification methods. To improve the sensitivity of read-level taxonomic classification, we developed an RNA-dependent RNA polymerase (RdRp) gene-based read classification tool RdRpBin. It combines alignment-based strategy with machine learning models in order to fully exploit the sequence properties of RdRp. We tested our method and compared its performance with the state-of-the-art tools on the simulated and real sequencing data. RdRpBin competes favorably with all. In particular, when the query RNA viruses share low sequence similarity with the known viruses ($\sim 0.4$), our tool can still maintain a higher F-score than the state-of-the-art tools. The experimental results on real data also showed that RdRpBin can classify more RNA viral reads with a relatively low false-positive rate. Thus, RdRpBin can be utilized to classify novel and diverged RNA viruses.

Sensitivity analyses of probabilistic and deterministic DTI tractography methodologies for studying arm muscle architecture.

PURPOSE: To determine the sensitivity profiles of probabilistic and deterministic DTI tractography methods in estimating geometric properties in arm muscle anatomy. METHODS: Spin-echo diffusion-weighted MR images were acquired in the dominant arm of 10 participants. Both deterministic and probabilistic tractography were performed in two different muscle architectures of the parallel-structured biceps brachii (and the pennate-structured flexor carpi ulnaris. Muscle fascicle geometry estimates and number of fascicles were evaluated with respect to tractography turning angle, polynomial fitting order, and SNR. The DTI tractography estimated fascicle lengths were compared with measurements obtained from conventional cadaveric dissection and ultrasound modalities. RESULTS: The probabilistic method generally estimated fascicle lengths closer to ranges reported by conventional methods than the deterministic method, most evident in the biceps brachii (p > 0.05), consisting of longer, arc-like fascicles. For both methods, a wide turning angle (50º-90°) generated fascicle lengths that were in close agreement with conventional methods, most evident in the flexor carpi ulnaris (p > 0.05), consisting of shorter, feather-like fascicles. The probabilistic approach produced at least two times more fascicles than the deterministic approach. For both approaches, second-order fitting yielded about double the complete tracts as third-order fitting. In both muscles, as SNR decreased, deterministic tractography produced less fascicles but consistent geometry (p > 0.05), whereas probabilistic tractography produced a consistent number but altered geometry of fascicles (p < 0.001). CONCLUSION: Findings from this study provide best practice recommendations for implementing DTI tractography in skeletal muscle and will inform future in vivo studies of healthy and pathological muscle structure.

PET image denoising based on denoising diffusion probabilistic model.

PURPOSE: Due to various physical degradation factors and limited counts received, PET image quality needs further improvements. The denoising diffusion probabilistic model (DDPM) was a distribution learning-based model, which tried to transform a normal distribution into a specific data distribution based on iterative refinements. In this work, we proposed and evaluated different DDPM-based methods for PET image denoising. METHODS: Under the DDPM framework, one way to perform PET image denoising was to provide the PET image and/or the prior image as the input. Another way was to supply the prior image as the network input with the PET image included in the refinement steps, which could fit for scenarios of different noise levels. 150 brain [[Formula: see text]F]FDG datasets and 140 brain [[Formula: see text]F]MK-6240 (imaging neurofibrillary tangles deposition) datasets were utilized to evaluate the proposed DDPM-based methods. RESULTS: Quantification showed that the DDPM-based frameworks with PET information included generated better results than the nonlocal mean, Unet and generative adversarial network (GAN)-based denoising methods. Adding additional MR prior in the model helped achieved better performance and further reduced the uncertainty during image denoising. Solely relying on MR prior while ignoring the PET information resulted in large bias. Regional and surface quantification showed that employing MR prior as the network input while embedding PET image as a data-consistency constraint during inference achieved the best performance. CONCLUSION: DDPM-based PET image denoising is a flexible framework, which can efficiently utilize prior information and achieve better performance than the nonlocal mean, Unet and GAN-based denoising methods.

Probabilistic Refinement of Focused Ultrasound Thalamotomy Targeting for Parkinson's Disease Tremor.

BACKGROUND: There remains high variability in clinical outcomes when the same magnetic resonance image-guided focused ultrasound (MRgFUS) thalamotomy target is used for both essential tremor (ET) and tremor-dominant Parkinson's disease (TDPD). OBJECTIVE: Our goal is to refine the MRgFUS thalamotomy target for TDPD versus ET. METHODS: We retrospectively performed voxel-wise efficacy and structural connectivity mapping using 3-12-month post-procedure hand tremor scores for a multicenter cohort of 32 TDPD patients and a previously published cohort of 79 ET patients, and 24-hour T1-weighted post-MRgFUS brain images. We validated our findings using Unified Parkinson's Disease Rating Scale part III scores for an independent cohort of nine TDPD patients. RESULTS: The post-MRgFUS clinical improvements were 45.9% ± 35.9%, 55.5% ± 36%, and 46.1% ± 18.6% for ET, multicenter TDPD and validation TDPD cohorts, respectively. The TDPD and ET efficacy maps differed significantly (ppermute < 0.05), with peak TDPD improvement (87%) at x = -13.5; y = -15.0; z = 1.5, ~3.5 mm anterior and 3 mm dorsal to the ET target. Discriminative connectivity projections were to the motor and premotor regions in TDPD, and to the motor and somatosensory regions in ET. The disorder-specific voxel-wise efficacy map could be used to estimate outcome in TDPD patients with high accuracy (R = 0.8; R2 = 0.64; P < 0.0001). The model was validated using the independent cohort of nine TDPD patients (R = 0.73; R2 = 0.53; P = 0.025-voxel analysis). CONCLUSION: We demonstrated that the most effective MRgFUS thalamotomy target in TDPD is in the ventral intermediate nucleus/ventralis oralis posterior border region. This finding offers new insights into the thalamic regions instrumental in tremor control, with pivotal implications for improving treatment outcomes. © 2024 The Author(s). Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

Sex estimation of the adult Neandertal Regourdou 1 (Montignac, France): Implications for sexing human fossil remains.

Sex is a biological trait fundamental to the study of hominin fossils. Among the many questions that can be addressed are those related to taxonomy, biological variability, sexual dimorphism, paleoobstetrics, funerary selection, and paleodemography. While new methodologies such as paleogenomics or paleoproteomics can be used to determine sex, they have not been systematically applied to Pleistocene human remains due to their destructive nature. Therefore, we estimated sex from the coxal bone of the newly discovered pelvic remains of the Regourdou 1 Neandertal (Southwest France, MIS 5) based on morphological and metric data employing two methods that have been recently revised and shown to be reliable in multiple studies. Both methods calculate posterior probabilities of the estimate. The right coxal bone of Regourdou 1 was partially reconstructed providing additional traits for sex estimation. These methods were cross validated on 14 sufficiently preserved coxal bones of specimens from the Neandertal lineage. Our results show that the Regourdou 1 individual, whose postcranial skeleton is not robust, is a male, and that previous sex attributions of comparative Neandertal specimens are largely in agreement with those obtained here. Our results encourage additional morphological research of fossil hominins in order to develop a set of methods that are applicable, reliable, and reproducible.