Genome-wide functional screen of 3'UTR variants uncovers causal variants for human disease and evolution.

3' untranslated region (3'UTR) variants are strongly associated with human traits and diseases, yet few have been causally identified. We developed the massively parallel reporter assay for 3'UTRs (MPRAu) to sensitively assay 12,173 3'UTR variants. We applied MPRAu to six human cell lines, focusing on genetic variants associated with genome-wide association studies (GWAS) and human evolutionary adaptation. MPRAu expands our understanding of 3'UTR function, suggesting that simple sequences predominately explain 3'UTR regulatory activity. We adapt MPRAu to uncover diverse molecular mechanisms at base pair resolution, including an adenylate-uridylate (AU)-rich element of LEPR linked to potential metabolic evolutionary adaptations in East Asians. We nominate hundreds of 3'UTR causal variants with genetically fine-mapped phenotype associations. Using endogenous allelic replacements, we characterize one variant that disrupts a miRNA site regulating the viral defense gene TRIM14 and one that alters PILRB abundance, nominating a causal variant underlying transcriptional changes in age-related macular degeneration.

Prediction by statistical overlap theory of fraction of baseline occupied by chromatographic peaks.

The average minimum resolution required for separating adjacent single-component peaks (SCPs) in one-dimensional chromatograms is an important metric in statistical overlap theory (SOT). However, its value changes with changing chromatographic conditions in non-intuitive ways, when SOT predicts the average number of peaks (maxima). A more stable and easily understood value of resolution is obtained on making a different prediction. A general equation is derived for the sum of all separated and superposed widths of SCPs in a chromatogram. The equation is a function of the saturation α, a metric of chromatographic crowdedness, and is expressed in dimensionless form by dividing by the duration of the chromatogram. This dimensionless function, f(α), is also the cumulative distribution function of the probability of separating adjacent SCPs. Simulations based on the clustering of line segments representing SCPs verify expressions for f(α) calculated from five functions for the distribution of intervals between adjacent SCPs. Synthetic chromatograms are computed with different saturations, distributions of intervals, and distribution of SCP amplitudes. The chromatograms are analyzed by calculating the sum of the widths of peaks at different relative responses, dividing the sum by the duration of the chromatograms, and graphing the reduced sum against relative response. For small values of relative response, the reduced sum approaches the fraction of baseline that is occupied by chromatographic peaks. This fraction can be identified with f(α), if the saturation α is defined with the average minimum resolution equaling 1.5. The identification is general and is independent of the saturation, the interval distribution, or the amplitude distribution. This constant value of resolution corresponds to baseline resolution, which simplifies the interpretation of SOT.

Transcriptomic signatures across human tissues identify functional rare genetic variation.

Rare genetic variants are abundant across the human genome, and identifying their function and phenotypic impact is a major challenge. Measuring aberrant gene expression has aided in identifying functional, large-effect rare variants (RVs). Here, we expanded detection of genetically driven transcriptome abnormalities by analyzing gene expression, allele-specific expression, and alternative splicing from multitissue RNA-sequencing data, and demonstrate that each signal informs unique classes of RVs. We developed Watershed, a probabilistic model that integrates multiple genomic and transcriptomic signals to predict variant function, validated these predictions in additional cohorts and through experimental assays, and used them to assess RVs in the UK Biobank, the Million Veterans Program, and the Jackson Heart Study. Our results link thousands of RVs to diverse molecular effects and provide evidence to associate RVs affecting the transcriptome with human traits.

Estimation of low-level components lost through chromatographic separations with finite detection limits.

The search for biomarkers allowing the assessment of disease by early diagnosis is facilitated by liquid chromatography. However, it is not clear how many components are lost due to being present in concentrations below the detection limit and/or being obscured by chromatographic peak overlap. First, we extend the study of missing components undertaken by Enke and Nagels, who employed the log-normal probability density function (pdf) for the distribution of signal intensities (and concentrations) of three mixtures. The Weibull and exponential pdfs, which have a higher probability of small-concentration components than the log-normal pdf, are also investigated. Results show that assessments of the loss of low-intensity signals by curve fitting are ambiguous. Next, we simulate synthetic chromatograms to compare the loss of peaks from superposition (overlap) with neighboring peaks to the loss arising from lying below the limit of detection (LOD) imposed by a finite signal-to-noise ratio (SNR). The simulations are made using amplitude pdfs based on the Enke-Nagels data as functions of relative column efficiency, i.e., saturation, and SNR. Results show that at the highest efficiencies, the lowest-amplitude peaks are lost below the LOD. However, at small and medium efficiencies, peak overlap is the dominant loss mechanism, suggesting that low-level components will not be found easily in liquid chromatography with single channel detectors regardless of SNR. A simple treatment shows that a multichannel detector, e.g., a mass spectrometer, is necessary to expose more low-level components.

Identification of rare-disease genes using blood transcriptome sequencing and large control cohorts.

It is estimated that 350 million individuals worldwide suffer from rare diseases, which are predominantly caused by mutation in a single gene1. The current molecular diagnostic rate is estimated at 50%, with whole-exome sequencing (WES) among the most successful approaches2-5. For patients in whom WES is uninformative, RNA sequencing (RNA-seq) has shown diagnostic utility in specific tissues and diseases6-8. This includes muscle biopsies from patients with undiagnosed rare muscle disorders6,9, and cultured fibroblasts from patients with mitochondrial disorders7. However, for many individuals, biopsies are not performed for clinical care, and tissues are difficult to access. We sought to assess the utility of RNA-seq from blood as a diagnostic tool for rare diseases of different pathophysiologies. We generated whole-blood RNA-seq from 94 individuals with undiagnosed rare diseases spanning 16 diverse disease categories. We developed a robust approach to compare data from these individuals with large sets of RNA-seq data for controls (n = 1,594 unrelated controls and n = 49 family members) and demonstrated the impacts of expression, splicing, gene and variant filtering strategies on disease gene identification. Across our cohort, we observed that RNA-seq yields a 7.5% diagnostic rate, and an additional 16.7% with improved candidate gene resolution.

Dominance of the tiger: The displacement of Aedes aegypti by Aedes albopictus in parts of the Torres Strait, Australia

Most of the inhabited islands in the Torres Strait region of Australia have experienced dengue outbreaks transmitted by Aedes aegypti at various times since at least the 1890s. However, another potential dengue vector, Aedes albopictus, the Asian tiger mosquito, was detected for the first time in 2005 and it expanded across most of the Torres Strait within a few years. In 2016, a survey of container-inhabiting mosquitoes was conducted in all island communities and Ae. aegypti was undetectable on most of the islands which the species had previously occupied, and had been replaced by Ae. albopictus. It is suspected that competitive displacement was responsible for the changes in species distribution. Aedes aegypti was only detected on Boigu Island and Thursday Island. Recent dengue outbreaks in the Torres Strait have apparently been driven by both Ae. albopictus and Ae. aegypti. The findings have major implications on management of dengue outbreaks in the region.

Theory of the probability of total resolution in chromatograms with systematic variation of average peak spacing and peak width.

An equation is proposed for the probability that all mixture constituents are separated, when the density (i.e., average number of eluting constituents per time) and width of single-component peaks (SCPs) vary systematically. The probability Pr that m SCPs are separated is modeled as the product of the m - 1 probabilities that adjacent pairs of SCPs are separated. Pr is then expressed as the geometric mean of the probability product raised to the power of m - 1. This geometric mean is approximated by an arithmetic mean equaling the probability that adjacent SCPs are separated, as calculated from previously developed statistical overlap theory (SOT) for variable SCP density and width. The theory is tested using previously reported and current in-house simulations of isocratic chromatograms of SCPs with random differences in standard chemical potential. In such chromatograms, more SCPs elute at short times than long times, and their widths are less at short times than long times. The average difference between 179 previously reported and currently predicted values of 100 x Pr is about 0.6, when 100 x Pr > 5. The theory requires numerical computation of one integral but can be approximated by an analytic equation for SOT probabilities close to one. For SCPs having retention times exceeding twice the void time, this equation simplifies to a previous SOT expression, with the gradient peak capacity replaced by the isocratic peak capacity. The versatility of the Pr theory is tested using three other models of chromatograms, in which the density and width of SCPs vary. The Pr predictions agree with simulation for all three models.

Likelihood of total resolution in selective comprehensive two-dimensional liquid chromatography with parallel processing: Simulation and theory.

The probability Pr(sLC×LC) that all peaks are separated by a resolution of 1.5 or more in selective comprehensive two-dimensional liquid chromatography (sLC × LC) is computed for simple model systems of 5 to 60 peaks and first-dimension (1D) gradient times of 100 to 2000 s. The computations include mimics of a commercial instrument, whose fixed second-dimension (2D) gradient time and use of one cycle time for initialization reduces Pr(sLC×LC) relative to an earlier report. For serial sLC × LC, in which a single device collects and transfers 1D multiplets to the second dimension, Pr(sLC×LC) under practical conditions is predicted to be only slightly larger than the probability of total resolution in LC × LC for separations of the same duration in each case. To increase Pr(sLC×LC), two model systems are proposed based on parallel processing, in which one device collects multiplets from the first separation while a second device simultaneously transfers fractions from previously collected multiplets to the second dimension for further separation. A sum of probabilities guideline is proposed by which optimal fixed 2D gradient times, ranging from 9.5 to 12 s, are found for both serial and parallel models. The increases of Pr(sLC×LC) based on parallel processing are modest; the largest is only 0.062 for one system and 0.106 for the other, relative to the serial model. A theory is derived that rationalizes the modesty of the increase, which was unexpected. It shows that Pr(sLC×LC) equals the probability of total resolution in the first dimension, plus the product of the probability that all 1D multiplets are transferred to the second dimension and the probability that all multiplets are separated in the second dimension. The theory shows that, although parallel processing is better than serial processing for multiplet transfer, the ability to leverage this gain is offset by the limited probability that all multiplets are then actually separated in the second dimension, which is only about 0.55 for conditions where the change from serial to parallel processing is most beneficial. With these findings in hand, two scenarios are examined for future consideration: one in which the 2D peak capacity is doubled, and another in which multiplets are always transferred to the second dimension. The latter shows considerable promise for increasing Pr(sLC×LC) substantially beyond its counterpart in LC × LC. For example, a 50% probability of separating all peaks in a 15-component mixture can be reached in 1150 s using LC × LC. The same probability can be reached in the same time for a sample with nearly twice as many components (27) in the case of sLC × LC, assuming transfer of all multiplets to the second dimension. These findings will be useful to those considering systematic approaches to developing 2D-LC methods for moderately complex mixtures, and to those interested in instrument development for 2D-LC.

Incorporation of Biological Knowledge Into the Study of Gene-Environment Interactions.

A growing knowledge base of genetic and environmental information has greatly enabled the study of disease risk factors. However, the computational complexity and statistical burden of testing all variants by all environments has required novel study designs and hypothesis-driven approaches. We discuss how incorporating biological knowledge from model organisms, functional genomics, and integrative approaches can empower the discovery of novel gene-environment interactions and discuss specific methodological considerations with each approach. We consider specific examples where the application of these approaches has uncovered effects of gene-environment interactions relevant to drug response and immunity, and we highlight how such improvements enable a greater understanding of the pathogenesis of disease and the realization of precision medicine.

The impact of rare variation on gene expression across tissues.

Rare genetic variants are abundant in humans and are expected to contribute to individual disease risk. While genetic association studies have successfully identified common genetic variants associated with susceptibility, these studies are not practical for identifying rare variants. Efforts to distinguish pathogenic variants from benign rare variants have leveraged the genetic code to identify deleterious protein-coding alleles, but no analogous code exists for non-coding variants. Therefore, ascertaining which rare variants have phenotypic effects remains a major challenge. Rare non-coding variants have been associated with extreme gene expression in studies using single tissues, but their effects across tissues are unknown. Here we identify gene expression outliers, or individuals showing extreme expression levels for a particular gene, across 44 human tissues by using combined analyses of whole genomes and multi-tissue RNA-sequencing data from the Genotype-Tissue Expression (GTEx) project v6p release. We find that 58% of underexpression and 28% of overexpression outliers have nearby conserved rare variants compared to 8% of non-outliers. Additionally, we developed RIVER (RNA-informed variant effect on regulation), a Bayesian statistical model that incorporates expression data to predict a regulatory effect for rare variants with higher accuracy than models using genomic annotations alone. Overall, we demonstrate that rare variants contribute to large gene expression changes across tissues and provide an integrative method for interpretation of rare variants in individual genomes.

FIRE: functional inference of genetic variants that regulate gene expression.

MOTIVATION: Interpreting genetic variation in noncoding regions of the genome is an important challenge for personal genome analysis. One mechanism by which noncoding single nucleotide variants (SNVs) influence downstream phenotypes is through the regulation of gene expression. Methods to predict whether or not individual SNVs are likely to regulate gene expression would aid interpretation of variants of unknown significance identified in whole-genome sequencing studies. RESULTS: We developed FIRE (Functional Inference of Regulators of Expression), a tool to score both noncoding and coding SNVs based on their potential to regulate the expression levels of nearby genes. FIRE consists of 23 random forests trained to recognize SNVs in cis-expression quantitative trait loci (cis-eQTLs) using a set of 92 genomic annotations as predictive features. FIRE scores discriminate cis-eQTL SNVs from non-eQTL SNVs in the training set with a cross-validated area under the receiver operating characteristic curve (AUC) of 0.807, and discriminate cis-eQTL SNVs shared across six populations of different ancestry from non-eQTL SNVs with an AUC of 0.939. FIRE scores are also predictive of cis-eQTL SNVs across a variety of tissue types. AVAILABILITY AND IMPLEMENTATION: FIRE scores for genome-wide SNVs in hg19/GRCh37 are available for download at https://sites.google.com/site/fireregulatoryvariation/. CONTACT: nilah@stanford.edu. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

A symbiotic-like biologically-driven regenerating fabric.

Living organisms constantly maintain their structural and biochemical integrity by the critical means of response, healing, and regeneration. Inanimate objects, on the other hand, are axiomatically considered incapable of responding to damage and healing it, leading to the profound negative environmental impact of their continuous manufacturing and trashing. Objects with such biological properties would be a significant step towards sustainable technology. In this work we present a feasible strategy for driving regeneration in fabric by means of integration with a bacterial biofilm to obtain a symbiotic-like hybrid - the fabric provides structural framework to the biofilm and supports its growth, whereas the biofilm responds to mechanical tear by synthesizing a silk protein engineered to self-assemble upon secretion from the cells. We propose the term crossbiosis to describe this and other hybrid systems combining organism and object. Our strategy could be implemented in other systems and drive sensing of integrity and response by regeneration in other materials as well.

Orthogonality measurements for multidimensional chromatography in three and higher dimensional separations.

Orthogonality metrics (OMs) for three and higher dimensional separations are proposed as extensions of previously developed OMs, which were used to evaluate the zone utilization of two-dimensional (2D) separations. These OMs include correlation coefficients, dimensionality, information theory metrics and convex-hull metrics. In a number of these cases, lower dimensional subspace metrics exist and can be readily calculated. The metrics are used to interpret previously generated experimental data. The experimental datasets are derived from Gilar's peptide data, now modified to be three dimensional (3D), and a comprehensive 3D chromatogram from Moore and Jorgenson. The Moore and Jorgenson chromatogram, which has 25 identifiable 3D volume elements or peaks, displayed good orthogonality values over all dimensions. However, OMs based on discretization of the 3D space changed substantially with changes in binning parameters. This example highlights the importance in higher dimensions of having an abundant number of retention times as data points, especially for methods that use discretization. The Gilar data, which in a previous study produced 21 2D datasets by the pairing of 7 one-dimensional separations, was reinterpreted to produce 35 3D datasets. These datasets show a number of interesting properties, one of which is that geometric and harmonic means of lower dimensional subspace (i.e., 2D) OMs correlate well with the higher dimensional (i.e., 3D) OMs. The space utilization of the Gilar 3D datasets was ranked using OMs, with the retention times of the datasets having the largest and smallest OMs presented as graphs. A discussion concerning the orthogonality of higher dimensional techniques is given with emphasis on molecular diversity in chromatographic separations. In the information theory work, an inconsistency is found in previous studies of orthogonality using the 2D metric often identified as %O. A new choice of metric is proposed, extended to higher dimensions, characterized by mixes of ordered and random retention times, and applied to the experimental datasets. In 2D, the new metric always equals or exceeds the original one. However, results from both the original and new methods are given.

Allele-specific expression reveals interactions between genetic variation and environment.

Identifying interactions between genetics and the environment (GxE) remains challenging. We have developed EAGLE, a hierarchical Bayesian model for identifying GxE interactions based on associations between environmental variables and allele-specific expression. Combining whole-blood RNA-seq with extensive environmental annotations collected from 922 human individuals, we identified 35 GxE interactions, compared with only four using standard GxE interaction testing. EAGLE provides new opportunities for researchers to identify GxE interactions using functional genomic data.

The impact of structural variation on human gene expression.

Structural variants (SVs) are an important source of human genetic diversity, but their contribution to traits, disease and gene regulation remains unclear. We mapped cis expression quantitative trait loci (eQTLs) in 13 tissues via joint analysis of SVs, single-nucleotide variants (SNVs) and short insertion/deletion (indel) variants from deep whole-genome sequencing (WGS). We estimated that SVs are causal at 3.5-6.8% of eQTLs-a substantially higher fraction than prior estimates-and that expression-altering SVs have larger effect sizes than do SNVs and indels. We identified 789 putative causal SVs predicted to directly alter gene expression: most (88.3%) were noncoding variants enriched at enhancers and other regulatory elements, and 52 were linked to genome-wide association study loci. We observed a notable abundance of rare high-impact SVs associated with aberrant expression of nearby genes. These results suggest that comprehensive WGS-based SV analyses will increase the power of common- and rare-variant association studies.

Population- and individual-specific regulatory variation in Sardinia.

Genetic studies of complex traits have mainly identified associations with noncoding variants. To further determine the contribution of regulatory variation, we combined whole-genome and transcriptome data for 624 individuals from Sardinia to identify common and rare variants that influence gene expression and splicing. We identified 21,183 expression quantitative trait loci (eQTLs) and 6,768 splicing quantitative trait loci (sQTLs), including 619 new QTLs. We identified high-frequency QTLs and found evidence of selection near genes involved in malarial resistance and increased multiple sclerosis risk, reflecting the epidemiological history of Sardinia. Using family relationships, we identified 809 segregating expression outliers (median z score of 2.97), averaging 13.3 genes per individual. Outlier genes were enriched for proximal rare variants, providing a new approach to study large-effect regulatory variants and their relevance to traits. Our results provide insight into the effects of regulatory variants and their relationship to population history and individual genetic risk.

Genome-wide association study of behavioral, physiological and gene expression traits in outbred CFW mice.

Although mice are the most widely used mammalian model organism, genetic studies have suffered from limited mapping resolution due to extensive linkage disequilibrium (LD) that is characteristic of crosses among inbred strains. Carworth Farms White (CFW) mice are a commercially available outbred mouse population that exhibit rapid LD decay in comparison to other available mouse populations. We performed a genome-wide association study (GWAS) of behavioral, physiological and gene expression phenotypes using 1,200 male CFW mice. We used genotyping by sequencing (GBS) to obtain genotypes at 92,734 SNPs. We also measured gene expression using RNA sequencing in three brain regions. Our study identified numerous behavioral, physiological and expression quantitative trait loci (QTLs). We integrated the behavioral QTL and eQTL results to implicate specific genes, including Azi2 in sensitivity to methamphetamine and Zmynd11 in anxiety-like behavior. The combination of CFW mice, GBS and RNA sequencing constitutes a powerful approach to GWAS in mice.

Impact of the X Chromosome and sex on regulatory variation.

The X Chromosome, with its unique mode of inheritance, contributes to differences between the sexes at a molecular level, including sex-specific gene expression and sex-specific impact of genetic variation. Improving our understanding of these differences offers to elucidate the molecular mechanisms underlying sex-specific traits and diseases. However, to date, most studies have either ignored the X Chromosome or had insufficient power to test for the sex-specific impact of genetic variation. By analyzing whole blood transcriptomes of 922 individuals, we have conducted the first large-scale, genome-wide analysis of the impact of both sex and genetic variation on patterns of gene expression, including comparison between the X Chromosome and autosomes. We identified a depletion of expression quantitative trait loci (eQTL) on the X Chromosome, especially among genes under high selective constraint. In contrast, we discovered an enrichment of sex-specific regulatory variants on the X Chromosome. To resolve the molecular mechanisms underlying such effects, we generated chromatin accessibility data through ATAC-sequencing to connect sex-specific chromatin accessibility to sex-specific patterns of expression and regulatory variation. As sex-specific regulatory variants discovered in our study can inform sex differences in heritable disease prevalence, we integrated our data with genome-wide association study data for multiple immune traits identifying several traits with significant sex biases in genetic susceptibilities. Together, our study provides genome-wide insight into how genetic variation, the X Chromosome, and sex shape human gene regulation and disease.

An Efficient Multiple-Testing Adjustment for eQTL Studies that Accounts for Linkage Disequilibrium between Variants.

Methods for multiple-testing correction in local expression quantitative trait locus (cis-eQTL) studies are a trade-off between statistical power and computational efficiency. Bonferroni correction, though computationally trivial, is overly conservative and fails to account for linkage disequilibrium between variants. Permutation-based methods are more powerful, though computationally far more intensive. We present an alternative correction method called eigenMT, which runs over 500 times faster than permutations and has adjusted p values that closely approximate empirical ones. To achieve this speed while also maintaining the accuracy of permutation-based methods, we estimate the effective number of independent variants tested for association with a particular gene, termed Meff, by using the eigenvalue decomposition of the genotype correlation matrix. We employ a regularized estimator of the correlation matrix to ensure Meff is robust and yields adjusted p values that closely approximate p values from permutations. Finally, using a common genotype matrix, we show that eigenMT can be applied with even greater efficiency to studies across tissues or conditions. Our method provides a simpler, more efficient approach to multiple-testing correction than existing methods and fits within existing pipelines for eQTL discovery.