Multiple Sources of Uncertainty Confound Inference of Historical Human Generation Times.

Wang et al. (2023) recently proposed an approach to infer the history of human generation intervals from changes in mutation profiles over time. As the relative proportions of different mutation types depend on the ages of parents, binning variants by the time they arose allows for the inference of changes in average paternal and maternal generation intervals. Applying this approach to published allele age estimates, Wang et al. (2023) inferred long-lasting sex differences in average generation times and surprisingly found that ancestral generation times of West African populations remained substantially higher than those of Eurasian populations extending tens of thousands of generations into the past. Here, we argue that the results and interpretations in Wang et al. (2023) are primarily driven by noise and biases in input data and a lack of validation using independent approaches for estimating allele ages. With the recent development of methods to reconstruct genome-wide gene genealogies, coalescence times, and allele ages, we caution that downstream analyses may be strongly influenced by uncharacterized biases in their output.

The NR4A agonist, Cytosporone B, attenuates pro-inflammatory mediators in human colorectal cancer tissue ex vivo.

Inflammation is a pivotal pathological factor in colorectal cancer (CRC) initiation and progression, and modulating this inflammatory state has the potential to ameliorate disease progression. NR4A receptors have emerged as key regulators of inflammatory pathways that are important in CRC. Here, we have examined the effect of NR4A agonist, Cytosporone B (CsnB), on colorectal tissue integrity and its effect on the inflammatory profile in CRC tissue ex vivo. Here, we demonstrate concentrations up 100 µM CsnB did not adversely affect tissue integrity as measured using transepithelial electrical resistance, histology and crypt height. Subsequently, we reveal through the use of a cytokine/chemokine array, ELISA and qRT-PCR analysis that multiple pro-inflammatory mediators were significantly increased in CRC tissue compared to control tissue, which were then attenuated with the addition of CsnB (such as IL-1ß, IL-8 and TNFα). Lastly, stratification of the data revealed that CsnB especially alters the inflammatory profile of tumours derived from males who had not undergone chemoradiotherapy. Thus, this study demonstrates that NR4A agonist CsnB does not adversely affect colon tissue structure or functionality and can attenuate the pro-inflammatory state of human CRC tissue ex vivo.

Genetic architecture and selective sweeps after polygenic adaptation to distant trait optima.

Understanding the genetic basis of phenotypic adaptation to changing environments is an essential goal of population and quantitative genetics. While technological advances now allow interrogation of genome-wide genotyping data in large panels, our understanding of the process of polygenic adaptation is still limited. To address this limitation, we use extensive forward-time simulation to explore the impacts of variation in demography, trait genetics, and selection on the rate and mode of adaptation and the resulting genetic architecture. We simulate a population adapting to an optimum shift, modeling sequence variation for 20 QTL for each of 12 different demographies for 100 different traits varying in the effect size distribution of new mutations, the strength of stabilizing selection, and the contribution of the genomic background. We then use random forest regression approaches to learn the relative importance of input parameters in determining a number of aspects of the process of adaptation, including the speed of adaptation, the relative frequency of hard sweeps and sweeps from standing variation, or the final genetic architecture of the trait. We find that selective sweeps occur even for traits under relatively weak selection and where the genetic background explains most of the variation. Though most sweeps occur from variation segregating in the ancestral population, new mutations can be important for traits under strong stabilizing selection that undergo a large optimum shift. We also show that population bottlenecks and expansion impact overall genetic variation as well as the relative importance of sweeps from standing variation and the speed with which adaptation can occur. We then compare our results to two traits under selection during maize domestication, showing that our simulations qualitatively recapitulate differences between them. Overall, our results underscore the complex population genetics of individual loci in even relatively simple quantitative trait models, but provide a glimpse into the factors that drive this complexity and the potential of these approaches for understanding polygenic adaptation.

Evidence of directional and stabilizing selection in contemporary humans.

Modern molecular genetic datasets, primarily collected to study the biology of human health and disease, can be used to directly measure the action of natural selection and reveal important features of contemporary human evolution. Here we leverage the UK Biobank data to test for the presence of linear and nonlinear natural selection in a contemporary population of the United Kingdom. We obtain phenotypic and genetic evidence consistent with the action of linear/directional selection. Phenotypic evidence suggests that stabilizing selection, which acts to reduce variance in the population without necessarily modifying the population mean, is widespread and relatively weak in comparison with estimates from other species.

The Effects of the Marine-Derived Polysaccharides Laminarin and Chitosan on Aspects of Colonic Health in Pigs Challenged with Dextran Sodium Sulphate.

This study examined the effects of dietary supplementation with laminarin or chitosan on colonic health in pigs challenged with dextran sodium sulphate (DSS). Weaned pigs were assigned to: (1) a basal diet (n = 22); (2) a basal diet + laminarin (n = 10); and (3) a basal diet + chitosan (n = 10). On d35, the basal group was split, creating four groups: (1) the basal diet (control); (2) the basal diet + DSS; (3) the basal diet + laminarin + DSS; and (4) the basal diet + chitosan + DSS. From d39-42, the pigs were orally challenged with DSS. On d44, colonic tissue/digesta samples were collected. The basal DSS group had reduced growth, higher pathology score and an increased expression of MMP1, IL13 and IL23 compared with the controls (p < 0.05); these parameters were similar between the DSS-challenged groups (p > 0.05). In the basal DSS group, the relative abundance of beneficial taxa including Prevotella and Roseburia were reduced while Escherichia/Shigella were increased, compared with the controls (p < 0.05). The relative abundance of Escherichia/Shigella was reduced and the molar proportions of acetate were increased in the laminarin DSS group compared with the basal DSS group (p < 0.01), suggesting that laminarin has potential to prevent pathogen proliferation and enhance the volatile fatty acid profile in the colon in a porcine model of colitis.

Tandem duplications lead to novel expression patterns through exon shuffling in Drosophila yakuba.

One common hypothesis to explain the impacts of tandem duplications is that whole gene duplications commonly produce additive changes in gene expression due to copy number changes. Here, we use genome wide RNA-seq data from a population sample of Drosophila yakuba to test this 'gene dosage' hypothesis. We observe little evidence of expression changes in response to whole transcript duplication capturing 5' and 3' UTRs. Among whole gene duplications, we observe evidence that dosage sharing across copies is likely to be common. The lack of expression changes after whole gene duplication suggests that the majority of genes are subject to tight regulatory control and therefore not sensitive to changes in gene copy number. Rather, we observe changes in expression level due to both shuffling of regulatory elements and the creation of chimeric structures via tandem duplication. Additionally, we observe 30 de novo gene structures arising from tandem duplications, 23 of which form with expression in the testes. Thus, the value of tandem duplications is likely to be more intricate than simple changes in gene dosage. The common regulatory effects from chimeric gene formation after tandem duplication may explain their contribution to genome evolution.

A Model of Compound Heterozygous, Loss-of-Function Alleles Is Broadly Consistent with Observations from Complex-Disease GWAS Datasets.

The genetic component of complex disease risk in humans remains largely unexplained. A corollary is that the allelic spectrum of genetic variants contributing to complex disease risk is unknown. Theoretical models that relate population genetic processes to the maintenance of genetic variation for quantitative traits may suggest profitable avenues for future experimental design. Here we use forward simulation to model a genomic region evolving under a balance between recurrent deleterious mutation and Gaussian stabilizing selection. We consider multiple genetic and demographic models, and several different methods for identifying genomic regions harboring variants associated with complex disease risk. We demonstrate that the model of gene action, relating genotype to phenotype, has a qualitative effect on several relevant aspects of the population genetic architecture of a complex trait. In particular, the genetic model impacts genetic variance component partitioning across the allele frequency spectrum and the power of statistical tests. Models with partial recessivity closely match the minor allele frequency distribution of significant hits from empirical genome-wide association studies without requiring homozygous effect sizes to be small. We highlight a particular gene-based model of incomplete recessivity that is appealing from first principles. Under that model, deleterious mutations in a genomic region partially fail to complement one another. This model of gene-based recessivity predicts the empirically observed inconsistency between twin and SNP based estimated of dominance heritability. Furthermore, this model predicts considerable levels of unexplained variance associated with intralocus epistasis. Our results suggest a need for improved statistical tools for region based genetic association and heritability estimation.

Efficient pedigree recording for fast population genetics simulation.

In this paper we describe how to efficiently record the entire genetic history of a population in forwards-time, individual-based population genetics simulations with arbitrary breeding models, population structure and demography. This approach dramatically reduces the computational burden of tracking individual genomes by allowing us to simulate only those loci that may affect reproduction (those having non-neutral variants). The genetic history of the population is recorded as a succinct tree sequence as introduced in the software package msprime, on which neutral mutations can be quickly placed afterwards. Recording the results of each breeding event requires storage that grows linearly with time, but there is a great deal of redundancy in this information. We solve this storage problem by providing an algorithm to quickly 'simplify' a tree sequence by removing this irrelevant history for a given set of genomes. By periodically simplifying the history with respect to the extant population, we show that the total storage space required is modest and overall large efficiency gains can be made over classical forward-time simulations. We implement a general-purpose framework for recording and simplifying genealogical data, which can be used to make simulations of any population model more efficient. We modify two popular forwards-time simulation frameworks to use this new approach and observe efficiency gains in large, whole-genome simulations of one to two orders of magnitude. In addition to speed, our method for recording pedigrees has several advantages: (1) All marginal genealogies of the simulated individuals are recorded, rather than just genotypes. (2) A population of N individuals with M polymorphic sites can be stored in O(N log N + M) space, making it feasible to store a simulation's entire final generation as well as its history. (3) A simulation can easily be initialized with a more efficient coalescent simulation of deep history. The software for recording and processing tree sequences is named tskit.

The Effect of Increasing Inclusion Levels of a Fucoidan-Rich Extract Derived from Ascophyllum nodosum on Growth Performance and Aspects of Intestinal Health of Pigs Post-Weaning.

This study examines the effects of increasing dietary inclusion levels of fucoidan, from a 44% fucoidan extract on the growth performance and intestinal health of pigs post-weaning (PW). Seventy-two newly weaned pigs (8.4 kg (SD 1.06)) were assigned to: (T1) basal diet (BD); (T2) BD + 125 ppm fucoidan; (T3) BD + 250 ppm fucoidan (8 pens/treatment). The appropriate quantity of a 44% fucoidan extract was included to achieve these inclusion levels. Faecal scores were recorded daily. On d15 PW, samples were collected from the intestinal tract from 1 pig/pen from the BD and BD + 250 ppm fucoidan groups. Pigs supplemented with 250 ppm fucoidan had improved faecal scores and increased concentrations of total volatile fatty acids and propionate in the colon (p < 0.05). The fucoidan-rich extract reduced the expression of CLDN5 (duodenum), SCL5A1/SGLT1 and SI (jejunum) and TJP1, FABP2, and SLC5A1 (ileum) (p < 0.05). The extract reduced the relative abundance of Prevotella and Lachnospiraceae (p < 0.05) and increased the abundance of Helicobacter (p < 0.01) in the caecum. However, no negative impact on growth performance or small intestinal morphology was observed. Thus, the inclusion of 250 ppm fucoidan improves faecal consistency without affecting growth performance and therefore warrants further investigation as a supplement for the prevention of PW diarrhoea under more challenging commercial conditions.

Interprofessional Care and Teamwork in the ICU.

OBJECTIVES: We describe the importance of interprofessional care in modern critical care medicine. This review highlights the essential roles played by specific members of the interprofessional care team, including patients and family members, and discusses quality improvement initiatives that require interprofessional collaboration for success. DATA SOURCES: Studies were identified through MEDLINE search using a variety of search phrases related to interprofessional care, critical care provider types, and quality improvement initiatives. Additional articles were identified through a review of the reference lists of identified articles. STUDY SELECTION: Original articles, review articles, and systematic reviews were considered. DATA EXTRACTION: Manuscripts were selected for inclusion based on expert opinion of well-designed or key studies and review articles. DATA SYNTHESIS: "Interprofessional care" refers to care provided by a team of healthcare professionals with overlapping expertise and an appreciation for the unique contribution of other team members as partners in achieving a common goal. A robust body of data supports improvement in patient-level outcomes when care is provided by an interprofessional team. Critical care nurses, advanced practice providers, pharmacists, respiratory care practitioners, rehabilitation specialists, dieticians, social workers, case managers, spiritual care providers, intensivists, and nonintensivist physicians each provide unique expertise and perspectives to patient care, and therefore play an important role in a team that must address the diverse needs of patients and families in the ICU. Engaging patients and families as partners in their healthcare is also critical. Many important ICU quality improvement initiatives require an interprofessional approach, including Awakening and Breathing Coordination, Delirium, Early Exercise/Mobility, and Family Empowerment bundle implementation, interprofessional rounding practices, unit-based quality improvement initiatives, Patient and Family Advisory Councils, end-of-life care, coordinated sedation awakening and spontaneous breathing trials, intrahospital transport, and transitions of care. CONCLUSIONS: A robust body of evidence supports an interprofessional approach as a key component in the provision of high-quality critical care to patients of increasing complexity and with increasingly diverse needs.

The Drosophila melanogaster Genetic Reference Panel.

A major challenge of biology is understanding the relationship between molecular genetic variation and variation in quantitative traits, including fitness. This relationship determines our ability to predict phenotypes from genotypes and to understand how evolutionary forces shape variation within and between species. Previous efforts to dissect the genotype-phenotype map were based on incomplete genotypic information. Here, we describe the Drosophila melanogaster Genetic Reference Panel (DGRP), a community resource for analysis of population genomics and quantitative traits. The DGRP consists of fully sequenced inbred lines derived from a natural population. Population genomic analyses reveal reduced polymorphism in centromeric autosomal regions and the X chromosome, evidence for positive and negative selection, and rapid evolution of the X chromosome. Many variants in novel genes, most at low frequency, are associated with quantitative traits and explain a large fraction of the phenotypic variance. The DGRP facilitates genotype-phenotype mapping using the power of Drosophila genetics.

Preventing Harm in the ICU-Building a Culture of Safety and Engaging Patients and Families.

OBJECTIVE: Preventing harm remains a persistent challenge in the ICU despite evidence-based practices known to reduce the prevalence of adverse events. This review seeks to describe the critical role of safety culture and patient and family engagement in successful quality improvement initiatives in the ICU. We review the evidence supporting the impact of safety culture and provide practical guidance for those wishing to implement initiatives aimed at improving safety culture and more effectively integrate patients and families in such efforts. DATA SOURCES: Literature review using PubMed including evaluation of key studies assessing large-scale quality improvement efforts in the ICU, impact of safety culture on patient outcomes, methodologies for quality improvement commonly used in healthcare, and patient and family engagement. Print and web-based resources from leading patient safety organizations were also searched. STUDY SELECTION: Our group completed a review of original studies, review articles, book chapters, and recommendations from leading patient safety organizations. DATA EXTRACTION: Our group determined by consensus which resources would best inform this review. DATA SYNTHESIS: A strong safety culture is associated with reduced adverse events, lower mortality rates, and lower costs. Quality improvement efforts have been shown to be more effective and sustainable when paired with a strong safety culture. Different methodologies exist for quality improvement in the ICU; a thoughtful approach to implementation that engages frontline providers and administrative leadership is essential for success. Efforts to substantively include patients and families in the processes of quality improvement work in the ICU should be expanded. CONCLUSIONS: Efforts to establish a culture of safety and meaningfully engage patients and families should form the foundation for all safety interventions in the ICU. This review describes an approach that integrates components of several proven quality improvement methodologies to enhance safety culture in the ICU and highlights opportunities to include patients and families.

A second-generation assembly of the Drosophila simulans genome provides new insights into patterns of lineage-specific divergence.

We create a new assembly of the Drosophila simulans genome using 142 million paired short-read sequences and previously published data for strain w(501). Our assembly represents a higher-quality genomic sequence with greater coverage, fewer misassemblies, and, by several indexes, fewer sequence errors. Evolutionary analysis of this genome reference sequence reveals interesting patterns of lineage-specific divergence that are different from those previously reported. Specifically, we find that Drosophila melanogaster evolves faster than D. simulans at all annotated classes of sites, including putatively neutrally evolving sites found in minimal introns. While this may be partly explained by a higher mutation rate in D. melanogaster, we also find significant heterogeneity in rates of evolution across classes of sites, consistent with historical differences in the effective population size for the two species. Also contrary to previous findings, we find that the X chromosome is evolving significantly faster than autosomes for nonsynonymous and most noncoding DNA sites and significantly slower for synonymous sites. The absence of a X/A difference for putatively neutral sites and the robustness of the pattern to Gene Ontology and sex-biased expression suggest that partly recessive beneficial mutations may comprise a substantial fraction of noncoding DNA divergence observed between species. Our results have more general implications for the interpretation of evolutionary analyses of genomes of different quality.

Genome-wide analysis of a long-term evolution experiment with Drosophila.

Experimental evolution systems allow the genomic study of adaptation, and so far this has been done primarily in asexual systems with small genomes, such as bacteria and yeast. Here we present whole-genome resequencing data from Drosophila melanogaster populations that have experienced over 600 generations of laboratory selection for accelerated development. Flies in these selected populations develop from egg to adult ∼20% faster than flies of ancestral control populations, and have evolved a number of other correlated phenotypes. On the basis of 688,520 intermediate-frequency, high-quality single nucleotide polymorphisms, we identify several dozen genomic regions that show strong allele frequency differentiation between a pooled sample of five replicate populations selected for accelerated development and pooled controls. On the basis of resequencing data from a single replicate population with accelerated development, as well as single nucleotide polymorphism data from individual flies from each replicate population, we infer little allele frequency differentiation between replicate populations within a selection treatment. Signatures of selection are qualitatively different than what has been observed in asexual species; in our sexual populations, adaptation is not associated with 'classic' sweeps whereby newly arising, unconditionally advantageous mutations become fixed. More parsimonious explanations include 'incomplete' sweep models, in which mutations have not had enough time to fix, and 'soft' sweep models, in which selection acts on pre-existing, common genetic variants. We conclude that, at least for life history characters such as development time, unconditionally advantageous alleles rarely arise, are associated with small net fitness gains or cannot fix because selection coefficients change over time.

Properties and modeling of GWAS when complex disease risk is due to non-complementing, deleterious mutations in genes of large effect.

Current genome-wide association studies (GWAS) have high power to detect intermediate frequency SNPs making modest contributions to complex disease, but they are underpowered to detect rare alleles of large effect (RALE). This has led to speculation that the bulk of variation for most complex diseases is due to RALE. One concern with existing models of RALE is that they do not make explicit assumptions about the evolution of a phenotype and its molecular basis. Rather, much of the existing literature relies on arbitrary mapping of phenotypes onto genotypes obtained either from standard population-genetic simulation tools or from non-genetic models. We introduce a novel simulation of a 100-kilobase gene region, based on the standard definition of a gene, in which mutations are unconditionally deleterious, are continuously arising, have partially recessive and non-complementing effects on phenotype (analogous to what is widely observed for most Mendelian disorders), and are interspersed with neutral markers that can be genotyped. Genes evolving according to this model exhibit a characteristic GWAS signature consisting of an excess of marginally significant markers. Existing tests for an excess burden of rare alleles in cases have low power while a simple new statistic has high power to identify disease genes evolving under our model. The structure of linkage disequilibrium between causative mutations and significantly associated markers under our model differs fundamentally from that seen when rare causative markers are assumed to be neutral. Rather than tagging single haplotypes bearing a large number of rare causative alleles, we find that significant SNPs in a GWAS tend to tag single causative mutations of small effect relative to other mutations in the same gene. Our results emphasize the importance of evaluating the power to detect associations under models that are genetically and evolutionarily motivated.

The power to detect quantitative trait loci using resequenced, experimentally evolved populations of diploid, sexual organisms.

A novel approach for dissecting complex traits is to experimentally evolve laboratory populations under a controlled environment shift, resequence the resulting populations, and identify single nucleotide polymorphisms (SNPs) and/or genomic regions highly diverged in allele frequency. To better understand the power and localization ability of such an evolve and resequence (E&R) approach, we carried out forward-in-time population genetics simulations of 1 Mb genomic regions under a large combination of experimental conditions, then attempted to detect significantly diverged SNPs. Our analysis indicates that the ability to detect differentiation between populations is primarily affected by selection coefficient, population size, number of replicate populations, and number of founding haplotypes. We estimate that E&R studies can detect and localize causative sites with 80% success or greater when the number of founder haplotypes is over 500, experimental populations are replicated at least 25-fold, population size is at least 1,000 diploid individuals, and the selection coefficient on the locus of interest is at least 0.1. More achievable experimental designs (less replicated, fewer founder haplotypes, smaller effective population size, and smaller selection coefficients) can have power of greater than 50% to identify a handful of SNPs of which one is likely causative. Similarly, in cases where s ≥ 0.2, less demanding experimental designs can yield high power.

Landscape of standing variation for tandem duplications in Drosophila yakuba and Drosophila simulans.

We have used whole genome paired-end Illumina sequence data to identify tandem duplications in 20 isofemale lines of Drosophila yakuba and 20 isofemale lines of D. simulans and performed genome wide validation with PacBio long molecule sequencing. We identify 1,415 tandem duplications that are segregating in D. yakuba as well as 975 duplications in D. simulans, indicating greater variation in D. yakuba. Additionally, we observe high rates of secondary deletions at duplicated sites, with 8% of duplicated sites in D. simulans and 17% of sites in D. yakuba modified with deletions. These secondary deletions are consistent with the action of the large loop mismatch repair system acting to remove polymorphic tandem duplication, resulting in rapid dynamics of gain and loss in duplicated alleles and a richer substrate of genetic novelty than has been previously reported. Most duplications are present in only single strains, suggesting that deleterious impacts are common. Drosophila simulans shows larger numbers of whole gene duplications in comparison to larger proportions of gene fragments in D. yakuba. Drosophila simulans displays an excess of high-frequency variants on the X chromosome, consistent with adaptive evolution through duplications on the D. simulans X or demographic forces driving duplicates to high frequency. We identify 78 chimeric genes in D. yakuba and 38 chimeric genes in D. simulans, as well as 143 cases of recruited noncoding sequence in D. yakuba and 96 in D. simulans, in agreement with rates of chimeric gene origination in D. melanogaster. Together, these results suggest that tandem duplications often result in complex variation beyond whole gene duplications that offers a rich substrate of standing variation that is likely to contribute both to detrimental phenotypes and disease, as well as to adaptive evolutionary change.

Genome sequencing reveals complex speciation in the Drosophila simulans clade.

The three species of the Drosophila simulans clade--the cosmopolitan species, D. simulans, and the two island endemic species, D. mauritiana and D. sechellia--are important models in speciation genetics, but some details of their phylogenetic and speciation history remain unresolved. The order and timing of speciation are disputed, and the existence, magnitude, and timing of gene flow among the three species remain unclear. Here we report on the analysis of a whole-genome four-species sequence alignment that includes all three D. simulans clade species as well as the D. melanogaster reference sequence. The alignment comprises novel, paired short-read sequence data from a single highly inbred line each from D. simulans, D. mauritiana, and D. sechellia. We are unable to reject a species phylogeny with a basal polytomy; the estimated age of the polytomy is 242,000 yr before the present. However, we also find that up to 4.6% of autosomal and 2.2% of X-linked regions have evolutionary histories consistent with recent gene flow between the mainland species (D. simulans) and the two island endemic species (D. mauritiana and D. sechellia). Our findings thus show that gene flow has occurred throughout the genomes of the D. simulans clade species despite considerable geographic, ecological, and intrinsic reproductive isolation. Last, our analysis of lineage-specific changes confirms that the D. sechellia genome has experienced a significant excess of slightly deleterious changes and a dearth of presumed favorable changes. The relatively reduced efficacy of natural selection in D. sechellia is consistent with its derived, persistently reduced historical effective population size.

Breed differences in humoral and cellular responses of lambs to experimental infection with the gastrointestinal nematode Teladorsagia circumcincta.

While Texel lambs have increased resistance to infection with the gastrointestinal nematode Teladorsagia circumcincta compared to Suffolk lambs, the underlying resistance mechanisms are still unknown. The aim of this study was to compare parasitological, humoral and cellular responses of Texel and Suffolk lambs over time following a single experimental infection with T. circumcincta. Gastrointestinal nematode free (but not naïve) lambs received a single oral dose of 3 × 104 infective T. circumcincta larvae. The variables examined included worm burden, mucosal and serum IgA, abomasal mast cells and eosinophils, haematological parameters and plasma pepsinogen. Texel lambs had significantly lower worm burden on day 14 and lower plasma pepsinogen concentration from day 14 onwards than Suffolks and their response in mucosal IgA to infection occurred earlier. The results from the study suggest that an earlier local IgA response in the Texel contributes to the resistant characteristics of the breed, while the increased level of plasma pepsinogen in the Suffolk lambs implies greater abomasal tissue damage arising from the nematode infection.

Abundance and distribution of transposable elements in two Drosophila QTL mapping resources.

Here we present computational machinery to efficiently and accurately identify transposable element (TE) insertions in 146 next-generation sequenced inbred strains of Drosophila melanogaster. The panel of lines we use in our study is composed of strains from a pair of genetic mapping resources: the Drosophila Genetic Reference Panel (DGRP) and the Drosophila Synthetic Population Resource (DSPR). We identified 23,087 TE insertions in these lines, of which 83.3% are found in only one line. There are marked differences in the distribution of elements over the genome, with TEs found at higher densities on the X chromosome, and in regions of low recombination. We also identified many more TEs per base pair of intronic sequence and fewer TEs per base pair of exonic sequence than expected if TEs are located at random locations in the euchromatic genome. There was substantial variation in TE load across genes. For example, the paralogs derailed and derailed-2 show a significant difference in the number of TE insertions, potentially reflecting differences in the selection acting on these loci. When considering TE families, we find a very weak effect of gene family size on TE insertions per gene, indicating that as gene family size increases the number of TE insertions in a given gene within that family also increases. TEs are known to be associated with certain phenotypes, and our data will allow investigators using the DGRP and DSPR to assess the functional role of TE insertions in complex trait variation more generally. Notably, because most TEs are very rare and often private to a single line, causative TEs resulting in phenotypic differences among individuals may typically fail to replicate across mapping panels since individual elements are unlikely to segregate in both panels. Our data suggest that "burden tests" that test for the effect of TEs as a class may be more fruitful.