Bonsai: An efficient method for inferring large human pedigrees from genotype data.

Pedigree inference from genotype data is a challenging problem, particularly when pedigrees are sparsely sampled and individuals may be distantly related to their closest genotyped relatives. We present a method that infers small pedigrees of close relatives and then assembles them into larger pedigrees. To assemble large pedigrees, we introduce several formulas and tools including a likelihood for the degree separating two small pedigrees, a generalization of the fast DRUID point estimate of the degree separating two pedigrees, a method for detecting individuals who share background identity-by-descent (IBD) that does not reflect recent common ancestry, and a method for identifying the ancestral branches through which distant relatives are connected. Our method also takes several approaches that help to improve the accuracy and efficiency of pedigree inference. In particular, we incorporate age information directly into the likelihood rather than using ages only for consistency checks and we employ a heuristic branch-and-bound-like approach to more efficiently explore the space of possible pedigrees. Together, these approaches make it possible to construct large pedigrees that are challenging or intractable for current inference methods.

Fast and Robust Identity-by-Descent Inference with the Templated Positional Burrows-Wheeler Transform.

Estimating the genomic location and length of identical-by-descent (IBD) segments among individuals is a crucial step in many genetic analyses. However, the exponential growth in the size of biobank and direct-to-consumer genetic data sets makes accurate IBD inference a significant computational challenge. Here we present the templated positional Burrows-Wheeler transform (TPBWT) to make fast IBD estimates robust to genotype and phasing errors. Using haplotype data simulated over pedigrees with realistic genotyping and phasing errors, we show that the TPBWT outperforms other state-of-the-art IBD inference algorithms in terms of speed and accuracy. For each phase-aware method, we explore the false positive and false negative rates of inferring IBD by segment length and characterize the types of error commonly found. Our results highlight the fragility of most phased IBD inference methods; the accuracy of IBD estimates can be highly sensitive to the quality of haplotype phasing. Additionally, we compare the performance of the TPBWT against a widely used phase-free IBD inference approach that is robust to phasing errors. We introduce both in-sample and out-of-sample TPBWT-based IBD inference algorithms and demonstrate their computational efficiency on massive-scale data sets with millions of samples. Furthermore, we describe the binary file format for TPBWT-compressed haplotypes that results in fast and efficient out-of-sample IBD computes against very large cohort panels. Finally, we demonstrate the utility of the TPBWT in a brief empirical analysis, exploring geographic patterns of haplotype sharing within Mexico. Hierarchical clustering of IBD shared across regions within Mexico reveals geographically structured haplotype sharing and a strong signal of isolation by distance. Our software implementation of the TPBWT is freely available for noncommercial use in the code repository (https://github.com/23andMe/phasedibd, last accessed January 11, 2021).

On the joint distribution of tree height and tree length under the coalescent.

Many statistics that examine genetic variation depend on the underlying shapes of genealogical trees. Under the coalescent model, we investigate the joint distribution of two quantities that describe genealogical tree shape: tree height and tree length. We derive a recursive formula for their exact joint distribution under a demographic model of a constant-sized population. We obtain approximations for the mean and variance of the ratio of tree height to tree length, using them to show that this ratio converges in probability to 0 as the sample size increases. We find that as the sample size increases, the correlation coefficient for tree height and length approaches (π2-6)∕[π2π2-18]≈0.9340. Using simulations, we examine the joint distribution of height and length under demographic models with population growth and population subdivision. We interpret the joint distribution in relation to problems of interest in data analysis, including inference of the time to the most recent common ancestor. The results assist in understanding the influences of demographic histories on two fundamental features of tree shape.

The Effects of Population Size Histories on Estimates of Selection Coefficients from Time-Series Genetic Data.

Many approaches have been developed for inferring selection coefficients from time series data while accounting for genetic drift. These approaches have been motivated by the intuition that properly accounting for the population size history can significantly improve estimates of selective strengths. However, the improvement in inference accuracy that can be attained by modeling drift has not been characterized. Here, by comparing maximum likelihood estimates of selection coefficients that account for the true population size history with estimates that ignore drift by assuming allele frequencies evolve deterministically in a population of infinite size, we address the following questions: how much can modeling the population size history improve estimates of selection coefficients? How much can mis-inferred population sizes hurt inferences of selection coefficients? We conduct our analysis under the discrete Wright-Fisher model by deriving the exact probability of an allele frequency trajectory in a population of time-varying size and we replicate our results under the diffusion model. For both models, we find that ignoring drift leads to estimates of selection coefficients that are nearly as accurate as estimates that account for the true population history, even when population sizes are small and drift is high. This result is of interest because inference methods that ignore drift are widely used in evolutionary studies and can be many orders of magnitude faster than methods that account for population sizes.

SpectralTDF: transition densities of diffusion processes with time-varying selection parameters, mutation rates and effective population sizes.

MOTIVATION: In the Wright-Fisher diffusion, the transition density function describes the time evolution of the population-wide frequency of an allele. This function has several practical applications in population genetics and computing it for biologically realistic scenarios with selection and demography is an important problem. RESULTS: We develop an efficient method for finding a spectral representation of the transition density function for a general model where the effective population size, selection coefficients and mutation parameters vary over time in a piecewise constant manner. AVAILABILITY AND IMPLEMENTATION: The method, called SpectralTDF, is available at https://sourceforge.net/projects/spectraltdf/ CONTACT: yss@berkeley.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Genome-wide association studies in diverse populations.

Genome-wide association (GWA) studies have identified a large number of SNPs associated with disease phenotypes. As most GWA studies have been performed in populations of European descent, this Review examines the issues involved in extending the consideration of GWA studies to diverse worldwide populations. Although challenges exist with issues such as imputation, admixture and replication, investigation of a greater diversity of populations could make substantial contributions to the goal of mapping the genetic determinants of complex diseases for the human population as a whole.

Macrophage migration inhibitory factor acts as a neurotrophin in the developing inner ear.

This study is the first to demonstrate that macrophage migration inhibitory factor (MIF), an immune system 'inflammatory' cytokine that is released by the developing otocyst, plays a role in regulating early innervation of the mouse and chick inner ear. We demonstrate that MIF is a major bioactive component of the previously uncharacterized otocyst-derived factor, which directs initial neurite outgrowth from the statoacoustic ganglion (SAG) to the developing inner ear. Recombinant MIF acts as a neurotrophin in promoting both SAG directional neurite outgrowth and neuronal survival and is expressed in both the developing and mature inner ear of chick and mouse. A MIF receptor, CD74, is found on both embryonic SAG neurons and adult mouse spiral ganglion neurons. Mif knockout mice are hearing impaired and demonstrate altered innervation to the organ of Corti, as well as fewer sensory hair cells. Furthermore, mouse embryonic stem cells become neuron-like when exposed to picomolar levels of MIF, suggesting the general importance of this cytokine in neural development.

Theory and applications of a deterministic approximation to the coalescent model.

Under the coalescent model, the random number nt of lineages ancestral to a sample is nearly deterministic as a function of time when nt is moderate to large in value, and it is well approximated by its expectation E[nt]. In turn, this expectation is well approximated by simple deterministic functions that are easy to compute. Such deterministic functions have been applied to estimate allele age, effective population size, and genetic diversity, and they have been used to study properties of models of infectious disease dynamics. Although a number of simple approximations of E[nt] have been derived and applied to problems of population-genetic inference, the theoretical accuracy of the resulting approximate formulas and the inferences obtained using these approximations is not known, and the range of problems to which they can be applied is not well understood. Here, we demonstrate general procedures by which the approximation nt≈E[nt] can be used to reduce the computational complexity of coalescent formulas, and we show that the resulting approximations converge to their true values under simple assumptions. Such approximations provide alternatives to exact formulas that are computationally intractable or numerically unstable when the number of sampled lineages is moderate or large. We also extend an existing class of approximations of E[nt] to the case of multiple populations of time-varying size with migration among them. Our results facilitate the use of the deterministic approximation nt≈E[nt] for deriving functionally simple, computationally efficient, and numerically stable approximations of coalescent formulas under complicated demographic scenarios.

CORRECTING MODEL MISSPECIFICATION IN RELATIONSHIP ESTIMATES.

The datasets of large genotyping biobanks and direct-to-consumer genetic testing companies contain many related individuals. Until now, it has been widely accepted that the most distant relationships that can be detected are around fifteen degrees (approximately 8 th cousins) and that practical relationship estimates have a ceiling around ten degrees (approximately 5 th cousins). However, we show that these assumptions are incorrect and that they are due to a misapplication of relationship estimators. In particular, relationship estimators are applied almost exclusively to putative relatives who have been identified because they share detectable tracts of DNA identically by descent (IBD). However, no existing relationship estimator conditions on the event that two individuals share at least one detectable segment of IBD anywhere in the genome. As a result, the relationship estimates obtained using existing estimators are dramatically biased for distant relationships, inferring all sufficiently distant relationships to be around ten degrees regardless of the depth of the true relationship. Moreover, existing relationship estimators are derived under a model that assumes that each pair of related individuals shares a single common ancestor (or mating pair of ancestors). This model breaks down for relationships beyond 10 generations in the past because individuals share many thousands of cryptic common ancestors due to pedigree collapse. We first derive a corrected likelihood that conditions on the event that at least one segment is observed between a pair of putative relatives and we demonstrate that the corrected likelihood largely eliminates the bias in estimates of pairwise relationships and provides a more accurate characterization of the uncertainty in these estimates. We then reformulate the relationship inference problem to account for the fact that individuals share many common ancestors, not just one. We demonstrate that the most distant relationship that can be inferred may be forty degrees or more, rather than ten, extending the time-to-common ancestor from approximately 200 years in the past to approximately 600 years in the past or more. This dramatic increase in the range of relationship estimators makes it possible to infer relationships whose common ancestors lived before historical events such as European settlement of the Americas and the Transatlantic Slave Trade, and possibly much earlier.

SIMULATING PEDIGREES ASCERTAINED ON THE BASIS OF OBSERVED IBD SHARING.

In large genotyping datasets, individuals often have thousands of distant cousins with whom they share detectable segments of DNA identically by descent (IBD). The ability to simulate these distant relationships is important for developing and testing methods, carrying out power analyses, and performing population genetic analyses. Because distant relatives are unlikely to share detectable IBD segments by chance, many simulation replicates are needed to sample IBD between any given pair of distant relatives. Exponentially more samples are needed to simulate observable segments of IBD simultaneously among multiple pairs of distant relatives in a single pedigree. Using existing pedigree simulation methods that do not condition on the event that IBD is observed among certain pairs of relatives, the chances of sampling shared IBD patterns that reflect those observed in real data ascertained from large genotyping datasets are vanishingly small, even for pedigrees of modest size. Here, we show how to sample recombination breakpoints on a fixed pedigree while conditioning on the event that specified pairs of individuals share at least one observed segment of IBD. The resulting simulator makes it possible to sample genotypes and IBD segments on pedigrees that reflect those ascertained from biobank scale data.

Reconstructing parent genomes using siblings and other relatives.

Reconstructing the DNA of ancestors from their descendants has the potential to empower phenotypic analyses (including association and genetic nurture studies), improve pedigree reconstruction, and shed light on the ancestral population and phenotypes of ancestors. We developed HAPI-RECAP, a method that reconstructs the DNA of parents from full siblings and their relatives. This tool leverages HAPI2's output, a new phasing approach that applies to siblings (and optionally one or both parents) and reliably infers parent haplotypes but does not link the ungenotyped parents' DNA across chromosomes or between segments flanking ambiguities. By combining IBD between the reconstructed parents and the relatives, HAPI-RECAP resolves the source parent of these segments. Moreover, the method exploits crossovers the children inherited and sex-specific genetic maps to infer the reconstructed parents' sexes. We validated these methods on research participants from both 23andMe, Inc. and the San Antonio Mexican American Family Studies. Given data for one parent, HAPI2 reconstructs large fractions of the missing parent's DNA, between 77.6% and 99.97% among all families, and 90.3% on average in three- and four-child families. When reconstructing both parents, HAPI-RECAP inferred between 33.2% and 96.6% of the parents' genotypes, averaging 70.6% in four-child families. Reconstructed genotypes have average error rates < 10-3, or comparable to those from direct genotyping. HAPI-RECAP inferred the parent sexes 100% correctly given IBD-linked segments and can also reconstruct parents without any IBD. As datasets grow in size, more families will be implicitly collected; HAPI-RECAP holds promise to enable high quality parent genotype reconstruction.

The cytokine macrophage migration inhibitory factor (MIF) acts as a neurotrophin in the developing inner ear of the zebrafish, Danio rerio.

Macrophage migration inhibitory factor (MIF) plays versatile roles in the immune system. MIF is also widely expressed during embryonic development, particularly in the nervous system, although its roles in neural development are only beginning to be understood. Evidence from frogs, mice and zebrafish suggests that MIF has a major role as a neurotrophin in the early development of sensory systems, including the auditory system. Here we show that the zebrafish mif pathway is required for both sensory hair cell (HC) and sensory neuronal cell survival in the ear, for HC differentiation, semicircular canal formation, statoacoustic ganglion (SAG) development, and lateral line HC differentiation. This is consistent with our findings that MIF is expressed in the developing mammalian and avian auditory systems and promotes mouse and chick SAG neurite outgrowth and neuronal survival, demonstrating key instructional roles for MIF in vertebrate otic development.

30 years of [11C]methyl triflate.

Some scientific discoveries are well known only to a core group of researchers working on technical subjects. Nevertheless, they open new research directions, allow existing knowledge to be viewed in entirely new and useful ways, or provide a way to make something that was hard or impossible to make before. Carbon-11 methyl triflate ([11C]MeOTf) is one such advance, facilitating the synthesis of many carbon-11 radio tracers and broadening the range of applications of carbon-11 radiochemistry. The year 2022 marked the 30th anniversary of the original paper in Applied Radiation and Isotopes introducing a simple synthesis of [11C]MeOTf from carbon-11 methyl iodide ([11C]MeI) and it also marked the end of the fruitful career and life of the researcher who developed it, Douglas Jewett. It seems fitting to say a few words on how it came to be and how it has helped advance carbon-11 radiochemistry.

The genetic legacy of African Americans from Catoctin Furnace.

Few African Americans have been able to trace family lineages back to ancestors who died before the 1870 United States Census, the first in which all Black people were listed by name. We analyzed 27 individuals from Maryland's Catoctin Furnace African American Cemetery (1774-1850), identifying 41,799 genetic relatives among consenting research participants in 23andMe, Inc.'s genetic database. One of the highest concentrations of close relatives is in Maryland, suggesting that descendants of the Catoctin individuals remain in the area. We find that many of the Catoctin individuals derived African ancestry from the Wolof or Kongo groups and European ancestry from Great Britain and Ireland. This study demonstrates the power of joint analysis of historical DNA and large datasets generated through direct-to-consumer ancestry testing.

A national survey of 'inactive' physicians in the United States of America: enticements to reentry.

BACKGROUND: Physicians leaving and reentering clinical practice can have significant medical workforce implications. We surveyed inactive physicians younger than typical retirement age to determine their reasons for clinical inactivity and what barriers, real or perceived, there were to reentry into the medical workforce. METHODS: A random sample of 4975 inactive physicians aged under 65 years was drawn from the Physician Masterfile of the American Medical Association in 2008. Physicians were mailed a survey about activity in medicine and perceived barriers to reentry. Chi-square statistics were used for significance tests of the association between categorical variables and t-tests were used to test differences between means. RESULTS: Our adjusted response rate was 36.1%. Respondents were fully retired (37.5%), not currently active in medicine (43.0%) or now active (reentered, 19.4%). Nearly half (49.5%) were in or had practiced primary care. Personal health was the top reason for leaving for fully retired physicians (37.8%) or those not currently active in medicine (37.8%) and the second highest reason for physicians who had reentered (28.8%). For reentered (47.8%) and inactive (51.5%) physicians, the primary reason for returning or considering returning to practice was the availability of part-time work or flexible scheduling. Retired and currently inactive physicians used similar strategies to explore reentry, and 83% of both groups thought it would be difficult; among those who had reentered practice, 35.9% reported it was difficult to reenter. Retraining was uncommon for this group (37.5%). CONCLUSION: Availability of part-time work and flexible scheduling have a strong influence on decisions to leave or reenter clinical practice. Lack of retraining before reentry raises questions about patient safety and the clinical competence of reentered physicians.

Nonsteroidal sulfamate derivatives as new therapeutic approaches for Neurofibromatosis 2 (NF2).

BACKGROUND: Neurofibromatosis 1 and 2, although involving two different tumour suppressor genes (neurofibromin and merlin, respectively), are both cancer predisposition syndromes that disproportionately affect cells of neural crest origin. New therapeutic approaches for both NF1 and NF2 are badly needed. In promising previous work we demonstrated that two non-steroidal analogues of 2-methoxy-oestradiol (2ME2), STX3451(2-(3-bromo-4,5-dimethoxybenzyl)-7-methoxy-6-sulfamoyloxy-1,2,3,4-tetrahydroisoquinoline), and STX2895 (7-Ethyl-6-sulfamoyloxy-2-(3,4,5-trimethoxybenzyl)-1,2,3,4-tetrahydroisoquinoline) reduced tumour cell growth and induced apoptosis in malignant and benign human Neurofibromatosis 1 (NF1) tumour cells. In earlier NF1 mechanism of action studies we found that in addition to their effects on non-classical hormone-sensitive pathways, STX agents acted on the actin- and myosin-cytoskeleton, as well as PI3Kinase and MTOR signaling pathways. Tumour growth in NF2 cells is affected by different inhibitors from those affecting NF1 growth pathways: specifically, NF2 cells are affected by merlin-downstream pathway inhibitors. Because Merlin, the affected tumour suppressor gene in NF2, is also known to be involved in stabilizing membrane-cytoskeletal complexes, as well as in cell proliferation, and apoptosis, we looked for potentially common mechanisms of action in the agents' effects on NF1 and NF2. We set out to determine whether STX agents could therefore also provide a prospective avenue for treatment of NF2. METHODS: STX3451 and STX2895 were tested in dose-dependent studies for their effects on growth parameters of malignant and benign NF2 human tumour cell lines in vitro. The mechanisms of action of STX3451 and STX2895 were also analysed. RESULTS: Although neither of the agents tested affected cell growth or apoptosis in the NF2 tumour cell lines tested through the same mechanisms by which they affect these parameters in NF1 tumour cell lines, both agents disrupted actin- and myosin-based cytoskeletal structures in NF2 cell lines, with subsequent effects on growth and cell death. CONCLUSIONS: Both STX3451 and STX2895 provide new approaches for inducing cell death and lowering tumour burden in NF2 as well as in NF1, which both have limited treatment options.

Parallel Trajectories of Genetic and Linguistic Admixture in a Genetically Admixed Creole Population.

Joint analyses of genes and languages, both of which are transmitted in populations by descent with modification-genes vertically by Mendel's laws, language via combinations of vertical, oblique, and horizontal processes [1-4]-provide an informative approach for human evolutionary studies [5-10]. Although gene-language analyses have employed extensive data on individual genetic variation [11-23], their linguistic data have not considered corresponding long-recognized [24] variability in individual speech patterns, or idiolects. Genetically admixed populations that speak creole languages show high genetic and idiolectal variation-genetic variation owing to heterogeneity in ancestry within admixed groups [25, 26] and idiolectal variation owing to recent language formation from differentiated sources [27-31]. To examine cotransmission of genetic and linguistic variation within populations, we collected genetic markers and speech recordings in the admixed creole-speaking population of Cape Verde, whose Kriolu language traces to West African languages and Portuguese [29, 32-35] and whose genetic ancestry has individual variation in European and continental African contributions [36-39]. In parallel with the combined Portuguese and West African origin of Kriolu, we find that genetic admixture in Cape Verde varies on an axis separating Iberian and Senegambian populations. We observe, analogously to vertical genetic transmission, transmission of idiolect from parents to offspring, as idiolect is predicted by parental birthplace, even after controlling for shared parent-child birthplaces. Further, African genetic admixture correlates with an index tabulating idiolectal features with likely African origins. These results suggest that Cape Verdean genetic and linguistic admixture have followed parallel evolutionary trajectories, with cotransmission of genetic and linguistic variation.

Loci associated with skin pigmentation identified in African populations.

Despite the wide range of skin pigmentation in humans, little is known about its genetic basis in global populations. Examining ethnically diverse African genomes, we identify variants in or near SLC24A5, MFSD12, DDB1, TMEM138, OCA2, and HERC2 that are significantly associated with skin pigmentation. Genetic evidence indicates that the light pigmentation variant at SLC24A5 was introduced into East Africa by gene flow from non-Africans. At all other loci, variants associated with dark pigmentation in Africans are identical by descent in South Asian and Australo-Melanesian populations. Functional analyses indicate that MFSD12 encodes a lysosomal protein that affects melanogenesis in zebrafish and mice, and that mutations in melanocyte-specific regulatory regions near DDB1/TMEM138 correlate with expression of ultraviolet response genes under selection in Eurasians.

Practice of pediatric otolaryngology: results of the future of pediatric education II project.

OBJECTIVES: To define the practice of pediatric otolaryngology compared with general otolaryngology and to estimate pediatric otolaryngology workforce utilization and needs. METHODS: Survey of members of the American Academy of Pediatrics Section on Otolaryngology and Bronchoesophagology and the American Society of Pediatric Otolaryngology and of a random sample of the membership of the American Academy of Otolaryngology-Head and Neck Surgery. RESULTS: Pediatric otolaryngologists were more likely to practice in urban and/or academic settings than were general otolaryngologists. Children (age <18 years) comprised over 88% of the patients of pediatric otolaryngologists and 30% to 35% of the patients of general otolaryngologists. Pediatric otolaryngologists were more likely to see children with complicated diseases such as airway disorders or congenital anomalies than were general otolaryngologists. Pediatric otolaryngologists, unlike general otolaryngologists, reported an increasing volume of pediatric referrals, as well as increased complexity in the patients referred. The surveyed physicians estimated the present number of pediatric otolaryngologists in their communities as approximately 0.2 to 0.3 per 100 000 people. CONCLUSIONS: Most children receiving otolaryngologic care in the United States receive such care from general otolaryngologists. The patient profile and practice setting of the subspecialty of pediatric otolaryngology differ from those of general otolaryngology. The demand for pediatric otolaryngologists appears to be increasing, but many general otolaryngologists do not believe there is an increased need.

Practice of pediatric critical care medicine: results of the Future of Pediatric Education II survey of sections project.

OBJECTIVE: To summarize the demographics and practice patterns of the current pediatric critical care workforce and to identify the key workforce issues that may affect the delivery of pediatric critical care services in the future. DESIGN: A questionnaire designed to analyze current pediatric critical care workforce demographics and future workforce trends. SUBJECTS: Pediatric critical care physicians from the United States were identified from the American Academy of Pediatrics Critical Care Section, from a list of physicians certified in pediatric critical care medicine (PCCM) by the American Board of Pediatrics, and from a list of pediatrician members of the Society for Critical Care Medicine. INTERVENTIONS: None. MEASUREMENTS: PCCM physicians were polled regarding board certification, practice characteristics, professional activities, referral patterns, patient profiles, competition, job satisfaction, and projected retirement age. MAIN RESULTS: A total of 805 PCCM physicians completed the survey. When grouped by age, 40% of the responding PCCM physicians were younger than 40 yrs, 49% were 40 to 49 yrs old, and only 11% were 50 yrs of age or older. The younger group had a higher percentage of female pediatricians than the older groups. For all age groups, the largest proportion of time was devoted to direct patient care time in pediatric critical care. This was especially true for the youngest age group that had the largest amount of patient care time devoted to critical care (43%). Time devoted to research was also significantly higher for the younger age group, although very few respondents reported that they have >50% of their time protected for research. For all age groups, those reporting increases in referral volume and referral complexity over the previous 12 months far outnumbered those reporting decreases. The majority of respondents reported being satisfied with their career choice. In general, respondents were more likely to report that too many rather than too few PCCM physicians were currently being trained. Approximately one third of respondents (34%) planned on leaving the field of critical care medicine before retiring from medicine completely. CONCLUSIONS: PCCM physicians were increasingly women and working for >65 hrs/wk, with a good level of job satisfaction. Competition from a variety of sources seems to affect the work of PCCM physicians. The relatively small percentage of time devoted to research, however, is a finding of great concern.