Enabling precision medicine via standard communication of HTS provenance, analysis, and results.

A personalized approach based on a patient's or pathogen's unique genomic sequence is the foundation of precision medicine. Genomic findings must be robust and reproducible, and experimental data capture should adhere to findable, accessible, interoperable, and reusable (FAIR) guiding principles. Moreover, effective precision medicine requires standardized reporting that extends beyond wet-lab procedures to computational methods. The BioCompute framework (https://w3id.org/biocompute/1.3.0) enables standardized reporting of genomic sequence data provenance, including provenance domain, usability domain, execution domain, verification kit, and error domain. This framework facilitates communication and promotes interoperability. Bioinformatics computation instances that employ the BioCompute framework are easily relayed, repeated if needed, and compared by scientists, regulators, test developers, and clinicians. Easing the burden of performing the aforementioned tasks greatly extends the range of practical application. Large clinical trials, precision medicine, and regulatory submissions require a set of agreed upon standards that ensures efficient communication and documentation of genomic analyses. The BioCompute paradigm and the resulting BioCompute Objects (BCOs) offer that standard and are freely accessible as a GitHub organization (https://github.com/biocompute-objects) following the "Open-Stand.org principles for collaborative open standards development." With high-throughput sequencing (HTS) studies communicated using a BCO, regulatory agencies (e.g., Food and Drug Administration [FDA]), diagnostic test developers, researchers, and clinicians can expand collaboration to drive innovation in precision medicine, potentially decreasing the time and cost associated with next-generation sequencing workflow exchange, reporting, and regulatory reviews.

Evidence for the fixation of gene duplications by positive selection in Drosophila.

Gene duplications play a key role in the emergence of novel traits and in adaptation. But despite their centrality to evolutionary processes, it is still largely unknown how new gene duplicates are initially fixed within populations and later maintained in genomes. Long-standing debates on the evolution of gene duplications could be settled by determining the relative importance of genetic drift vs. positive selection in the fixation of new gene duplicates. Using the Drosophila Global Diversity Lines (GDL), we have combined genome-wide SNP polymorphism data with a novel set of copy number variant calls and gene expression profiles to characterize the polymorphic phase of new genes. We found that approximately half of the roughly 500 new complete gene duplications segregating in the GDL lead to significant increases in the expression levels of the duplicated genes and that these duplications are more likely to be found at lower frequencies, suggesting a negative impact on fitness. However, we also found that six of the nine gene duplications that are fixed or close to fixation in at least one of the five populations in our study show signs of being under positive selection, and that these duplications are likely beneficial because of dosage effects, with a possible role for additional mutations in two duplications. Our work suggests that in Drosophila, theoretical models that posit that gene duplications are immediately beneficial and fixed by positive selection are most relevant to explain the long-term evolution of gene duplications in this species.

Contrasting X-linked and autosomal diversity across 14 human populations.

Contrasting the genetic diversity of the human X chromosome (X) and autosomes has facilitated understanding historical differences between males and females and the influence of natural selection. Previous studies based on smaller data sets have left questions regarding how empirical patterns extend to additional populations and which forces can explain them. Here, we address these questions by analyzing the ratio of X-to-autosomal (X/A) nucleotide diversity with the complete genomes of 569 females from 14 populations. Results show that X/A diversity is similar within each continental group but notably lower in European (EUR) and East Asian (ASN) populations than in African (AFR) populations. X/A diversity increases in all populations with increasing distance from genes, highlighting the stronger impact of diversity-reducing selection on X than on the autosomes. However, relative X/A diversity (between two populations) is invariant with distance from genes, suggesting that selection does not drive the relative reduction in X/A diversity in non-Africans (0.842 ± 0.012 for EUR-to-AFR and 0.820 ± 0.032 for ASN-to-AFR comparisons). Finally, an array of models with varying population bottlenecks, expansions, and migration from the latest studies of human demographic history account for about half of the observed reduction in relative X/A diversity from the expected value of 1. They predict values between 0.91 and 0.94 for EUR-to-AFR comparisons and between 0.91 and 0.92 for ASN-to-AFR comparisons. Further reductions can be predicted by more extreme demographic events in excess of those captured by the latest studies but, in the absence of these, also by historical sex-biased demographic events or other processes.

Phenotype fingerprinting of bipolar disorder prodrome.

BACKGROUND: Detecting prodromal symptoms of bipolar disorder (BD) has garnered significant attention in recent research, as early intervention could potentially improve therapeutic efficacy and improve patient outcomes. The heterogeneous nature of the prodromal phase in BD, however, poses considerable challenges for investigators. Our study aimed to identify distinct prodromal phenotypes or "fingerprints" in patients diagnosed with BD and subsequently examine correlations between these fingerprints and relevant clinical outcomes. METHODS: 20,000 veterans diagnosed with BD were randomly selected for this study. K-means clustering analysis was performed on temporal graphs of the clinical features of each patient. We applied what we call "temporal blurring" to each patient image in order to allow clustering to focus on the clinical features, and not cluster patients based upon their varying temporal patterns in diagnosis, which lead to the desired types of clusters. We evaluated several outcomes including mortality rate, hospitalization rate, mean number of hospitalizations, mean length of stay, and the occurrence of a psychosis diagnosis within one year following the initial BD diagnosis. To determine the statistical significance of the observed differences for each outcome, we conducted appropriate tests, such as ANOVA or Chi-square. RESULTS: Our analysis yielded 8 clusters which appear to represent distinct phenotypes with differing clinical attributes. Each of these clusters also has statistically significant differences across all outcomes (p < 0.0001). The clinical features in many of the clusters were consistent with findings in the literature concerning prodromal symptoms in patients with BD. One cluster, notably characterized by patients lacking discernible prodromal symptoms, exhibited the most favorable results across all measured outcomes. CONCLUSION: Our study successfully identified distinct prodromal phenotypes in patients diagnosed with BD. We also found that these distinct prodromal phenotypes are associated with different clinical outcomes.

Estradiol-mediated protection against high-fat diet induced anxiety and obesity is associated with changes in the gut microbiota in female mice.

Decreased estrogens during menopause are associated with increased risk of anxiety, depression, type 2 diabetes and obesity. Similarly, depleting estrogens in rodents by ovariectomy, combined with a high-fat diet (HFD), increases anxiety and adiposity. How estrogens and diet interact to affect anxiety and metabolism is poorly understood. Mounting evidence indicates that gut microbiota influence anxiety and metabolism. Here, we investigated the effects of estradiol (E) and HFD on anxiety, metabolism, and their correlation with changes in gut microbiota in female mice. Adult C57BL/6J mice were ovariectomized, implanted with E or vehicle-containing capsules and fed a standard diet or HFD. Anxiety-like behavior was assessed and neuronal activation was measured by c-fos immunoreactivity throughout the brain using iDISCO. HFD increased anxiety-like behavior, while E reduced this HFD-dependent anxiogenic effect. Interestingly, E decreased neuronal activation in brain regions involved in anxiety and metabolism. E treatment also altered gut microbes, a subset of which were associated with anxiety-like behavior. These findings provide insight into gut microbiota-based therapies for anxiety and metabolic disorders associated with declining estrogens in menopausal women.

Survey of Global Genetic Diversity Within the Drosophila Immune System.

Numerous studies across a wide range of taxa have demonstrated that immune genes are routinely among the most rapidly evolving genes in the genome. This observation, however, does not address what proportion of immune genes undergo strong selection during adaptation to novel environments. Here, we determine the extent of very recent divergence in genes with immune function across five populations of Drosophila melanogaster and find that immune genes do not show an overall trend of recent rapid adaptation. Our population-based approach uses a set of carefully matched control genes to account for the effects of demography and local recombination rate, allowing us to identify whether specific immune functions are putative targets of strong selection. We find evidence that viral-defense genes are rapidly evolving in Drosophila at multiple timescales. Local adaptation to bacteria and fungi is less extreme and primarily occurs through changes in recognition and effector genes rather than large-scale changes to the regulation of the immune response. Surprisingly, genes in the Toll pathway, which show a high rate of adaptive substitution between the D. melanogaster and D. simulans lineages, show little population differentiation. Quantifying the flies for resistance to a generalist Gram-positive bacterial pathogen, we found that this genetic pattern of low population differentiation was recapitulated at the phenotypic level. In sum, our results highlight the complexity of immune evolution and suggest that Drosophila immune genes do not follow a uniform trajectory of strong directional selection as flies encounter new environments.

Extensive local adaptation within the chemosensory system following Drosophila melanogaster's global expansion.

How organisms adapt to new environments is of fundamental biological interest, but poorly understood at the genetic level. Chemosensory systems provide attractive models to address this problem, because they lie between external environmental signals and internal physiological responses. To investigate how selection has shaped the well-characterized chemosensory system of Drosophila melanogaster, we have analysed genome-wide data from five diverse populations. By couching population genomic analyses of chemosensory protein families within parallel analyses of other large families, we demonstrate that chemosensory proteins are not outliers for adaptive divergence between species. However, chemosensory families often display the strongest genome-wide signals of recent selection within D. melanogaster. We show that recent adaptation has operated almost exclusively on standing variation, and that patterns of adaptive mutations predict diverse effects on protein function. Finally, we provide evidence that chemosensory proteins have experienced relaxed constraint, and argue that this has been important for their rapid adaptation over short timescales.

Global diversity lines - a five-continent reference panel of sequenced Drosophila melanogaster strains.

Reference collections of multiple Drosophila lines with accumulating collections of "omics" data have proven especially valuable for the study of population genetics and complex trait genetics. Here we present a description of a resource collection of 84 strains of Drosophila melanogaster whose genome sequences were obtained after 12 generations of full-sib inbreeding. The initial rationale for this resource was to foster development of a systems biology platform for modeling metabolic regulation by the use of natural polymorphisms as perturbations. As reference lines, they are amenable to repeated phenotypic measurements, and already a large collection of metabolic traits have been assayed. Another key feature of these strains is their widespread geographic origin, coming from Beijing, Ithaca, Netherlands, Tasmania, and Zimbabwe. After obtaining 12.5× coverage of paired-end Illumina sequence reads, SNP and indel calls were made with the GATK platform. Thorough quality control was enabled by deep sequencing one line to >100×, and single-nucleotide polymorphisms and indels were validated using ddRAD-sequencing as an orthogonal platform. In addition, a series of preliminary population genetic tests were performed with these single-nucleotide polymorphism data for assessment of data quality. We found 83 segregating inversions among the lines, and as expected these were especially abundant in the African sample. We anticipate that this will make a useful addition to the set of reference D. melanogaster strains, thanks to its geographic structuring and unusually high level of genetic diversity.

Using Mendelian inheritance to improve high-throughput SNP discovery.

Restriction site-associated DNA sequencing or genotyping-by-sequencing (GBS) approaches allow for rapid and cost-effective discovery and genotyping of thousands of single-nucleotide polymorphisms (SNPs) in multiple individuals. However, rigorous quality control practices are needed to avoid high levels of error and bias with these reduced representation methods. We developed a formal statistical framework for filtering spurious loci, using Mendelian inheritance patterns in nuclear families, that accommodates variable-quality genotype calls and missing data--both rampant issues with GBS data--and for identifying sex-linked SNPs. Simulations predict excellent performance of both the Mendelian filter and the sex-linkage assignment under a variety of conditions. We further evaluate our method by applying it to real GBS data and validating a subset of high-quality SNPs. These results demonstrate that our metric of Mendelian inheritance is a powerful quality filter for GBS loci that is complementary to standard coverage and Hardy-Weinberg filters. The described method, implemented in the software MendelChecker, will improve quality control during SNP discovery in nonmodel as well as model organisms.

Analyses of X-linked and autosomal genetic variation in population-scale whole genome sequencing.

The ratio of genetic diversity on chromosome X to that on the autosomes is sensitive to both natural selection and demography. On the basis of whole-genome sequences of 69 females, we report that whereas this ratio increases with genetic distance from genes across populations, it is lower in Europeans than in West Africans independent of proximity to genes. This relative reduction is most parsimoniously explained by differences in demographic history without the need to invoke natural selection.