py_ped_sim - A flexible forward genetic simulator for complex family pedigree analysis.

Background: Large-scale family pedigrees are commonly used across medical, evolutionary, and forensic genetics. These pedigrees are tools for identifying genetic disorders, tracking evolutionary patterns, and establishing familial relationships via forensic genetic identification. However, there is a lack of software to accurately simulate different pedigree structures along with genomes corresponding to those individuals in a family pedigree. This limits simulation-based evaluations of methods that use pedigrees. Results: We have developed a python command-line-based tool called py_ped_sim that facilitates the simulation of pedigree structures and the genomes of individuals in a pedigree. py_ped_sim represents pedigrees as directed acyclic graphs, enabling conversion between standard pedigree formats and integration with the forward population genetic simulator, SLiM. Notably, py_ped_sim allows the simulation of varying numbers of offspring for a set of parents, with the capacity to shift the distribution of sibship sizes over generations. We additionally add simulations for events of misattributed paternity, which offers a way to simulate half-sibling relationships. We validated the accuracy of our software by simulating genomes onto diverse family pedigree structures, showing that the estimated kinship coefficients closely approximated expected values. Conclusions: py_ped_sim is a user-friendly and open-source solution for simulating pedigree structures and conducting pedigree genome simulations. It empowers medical, forensic, and evolutionary genetics researchers to gain deeper insights into the dynamics of genetic inheritance and relatedness within families.

Decreased accuracy of forensic DNA mixture analysis for groups with lower genetic diversity.

Forensic investigation of DNA samples from multiple contributors has become commonplace. These complex analyses use statistical frameworks accounting for multiple levels of uncertainty in allelic contributions from different individuals, particularly for samples containing few molecules of DNA. These methods have been thoroughly tested along some axes of variation, but less attention has been paid to accuracy across human genetic variation. Here, we quantify the accuracy of DNA mixture analysis over 244 human groups. We find higher false inclusion rates for mixtures with more contributors, and for groups with lower genetic diversity. Even for two-contributor mixtures where one contributor is known and the reference group is correctly specified, false inclusion rates are 1e-5 or higher for 56 out of 244 groups. This means that, depending on multiple testing, some false inclusions may be expected. These false positives could be lessened with more selective and conservative use of DNA mixture analysis.

The Urinary Metabolomic Fingerprint in Extremely Preterm Infants on Total Parenteral Nutrition vs. Enteral Feeds.

Total Parenteral Nutrition (TPN), which uses intravenous administration of nutrients, minerals and vitamins, is essential for sustaining premature infants until they transition to enteral feeds, but there is limited information on metabolomic differences between infants on TPN and enteral feeds. We performed untargeted global metabolomics on urine samples collected between 23-30 days of life from 314 infants born <29 weeks gestational age from the TOLSURF and PROP cohorts. Principal component analysis across all metabolites showed a separation of infants solely on TPN compared to infants who had transitioned to enteral feeds, indicating global metabolomic differences between infants based on feeding status. Among 913 metabolites that passed quality control filters, 609 varied in abundance between infants on TPN vs. enteral feeds at p < 0.05. Of these, 88% were in the direction of higher abundance in the urine of infants on enteral feeds. In a subset of infants in a longitudinal analysis, both concurrent and delayed changes in metabolite levels were observed with the initiation of enteral feeds. These infants had higher concentrations of essential amino acids, lipids, and vitamins, which are necessary for growth and development, suggesting the nutritional benefit of an enteral feeding regimen.

The hydrocortisone-responsive urinary metabolome of premature infants.

BACKGROUND: Extremely premature infants are at risk for circulatory collapse or respiratory failure that are often treated with hydrocortisone (HC); however, there is no information on the metabolic consequences of this therapy. METHODS: Longitudinal urine samples from infants <28 weeks gestation in the Trial of Late Surfactant were analyzed by untargeted UHPLC:MS/MS. Fourteen infants who received a tapering course of HC beginning at 3 mg/kg/day for ≥9 days were compared to 14 matched control infants. A secondary cross-sectional analysis by logistic regression used urines from 314 infants. RESULTS: Of 1145 urinary metabolites detected, abundance of 219, representing all the major biochemical pathways, changed at p < 0.05 in the HC-treated group with 90% decreasing; 3 cortisol derivatives increased ~2-fold with HC therapy. Only 11% of regulated metabolites remained responsive at the lowest HC dose. Regulated metabolites included two steroids and thiamin that are associated with lung inflammation in infants. HC responsiveness was confirmed in 57% of metabolites by cross-sectional analysis. CONCLUSIONS: HC treatment of premature infants influenced in a dose-dependent manner abundance of 19% of identified urinary metabolites of diverse biochemical systems, primarily reducing concentrations. These findings indicate that exposure to HC reversibly impacts the nutritional status of premature infants. IMPACT: Hydrocortisone treatment of premature infants with respiratory failure or circulatory collapse alters levels of a subset of urinary metabolites representing all major biochemical pathways. This is the first description of the scope, magnitude, timing and reversibility of metabolomic changes in infants in response to hydrocortisone, and it confirms corticosteroid regulation of three biochemicals that are associated with lung inflammatory status. The findings indicate a dose-dependency of hydrocortisone for metabolomic and anti-inflammatory effects, that prolonged therapy may lower the supply of many nutrients, and that monitoring concentrations of cortisol and inflammation markers may be a useful clinical approach during corticosteroid therapy.

Associations between forensic loci and expression levels of neighboring genes may compromise medical privacy.

A set of 20 short tandem repeats (STRs) is used by the US criminal justice system to identify suspects and to maintain a database of genetic profiles for individuals who have been previously convicted or arrested. Some of these STRs were identified in the 1990s, with a preference for markers in putative gene deserts to avoid forensic profiles revealing protected medical information. We revisit that assumption, investigating whether forensic genetic profiles reveal information about gene-expression variation or potential medical information. We find six significant correlations (false discovery rate = 0.23) between the forensic STRs and the expression levels of neighboring genes in lymphoblastoid cell lines. We explore possible mechanisms for these associations, showing evidence compatible with forensic STRs causing expression variation or being in linkage disequilibrium with a causal locus in three cases and weaker or potentially spurious associations in the other three cases. Together, these results suggest that forensic genetic loci may reveal expression levels and, perhaps, medical information.