Recombination affects allele-specific expression of deleterious variants in human populations.

How the genetic composition of a population changes through stochastic processes, such as genetic drift, in combination with deterministic processes, such as selection, is critical to understanding how phenotypes vary in space and time. Here, we show how evolutionary forces affecting selection, including recombination and effective population size, drive genomic patterns of allele-specific expression (ASE). Integrating tissue-specific genotypic and transcriptomic data from 1500 individuals from two different cohorts, we demonstrate that ASE is less often observed in regions of low recombination, and loci in high or normal recombination regions are more efficient at using ASE to underexpress harmful mutations. By tracking genetic ancestry, we discriminate between ASE variability due to past demographic effects, including subsequent bottlenecks, versus local environment. We observe that ASE is not randomly distributed along the genome and that population parameters influencing the efficacy of natural selection alter ASE levels genome wide.

Soybean aphids adapted to host-plant resistance by down regulating putative effectors and up regulating transposable elements.

In agricultural systems, crops equipped with host-plant resistance (HPR) have enhanced protection against pests, and are used as a safe and sustainable tool in pest management. In soybean, HPR can control the soybean aphid (Aphis glycines), but certain aphid populations have overcome this resistance (i.e., virulence). The molecular mechanisms underlying aphid virulence to HPR are unknown, but likely involve effector proteins that are secreted by aphids to modulate plant defenses. Another mechanism to facilitate adaptation is through the activity of transposable elements, which can become activated by stress. In this study, we performed RNA sequencing of virulent and avirulent soybean aphids fed susceptible or resistant (Rag1 + Rag2) soybean. Our goal was to better understand the molecular mechanisms underlying soybean aphid virulence. Our data showed that virulent aphids mostly down regulate putative effector genes relative to avirulent aphids, especially when aphids were fed susceptible soybean. Decreased expression of effectors may help evade HPR plant defenses. Virulent aphids also transcriptionally up regulate a diverse set of transposable elements and nearby genes, which is consistent with stress adaptation. Our work demonstrates two mechanisms of pest adaptation to resistance, and identifies effector gene targets for future functional testing.

Ecological contributions to body shape evolution in salamanders of the genus Eurycea (Plethodontidae).

BACKGROUND: Body shape can be both a consequence and cause of a species' evolution and ecology. There are many examples of phenotypes associated with specific ecological niches, likely as a result of specific selective regimes. A classic example of this is the phenotypic change associated with colonization of caves, including body and limb elongation. However, studies explicitly testing for differences in body shape between cave-dwelling and non-cave-dwelling lineages have been limited and so the role of the cave environment in determining morphological characteristics is still not completely understood. Here we examine variation in body shape among 405 individuals representing 20 species in the salamander genus Eurycea (Plethodontidae) and select outgroups exhibiting great diversity in morphology, ecological niche, and life history. RESULTS: After analyzing morphometric data in a phylogenetic context using phylogenetic MANOVA and examination of the phylomorphospace, we found significant differences in body shape among cave-dwelling and non-cave-dwelling species and between aquatic and terrestrial species. Notably, limb elongation and reduced body and tail size characterized cave-dwelling species. Terrestrial species also exhibited elongation of the limbs and digits. We also observed differences in shape variance among paedomorphic and biphasic species. Our results suggest that the functional limitations imposed by habitat and life history played a key role in the evolution of body shape in this group in the context of their phylogenetic history.

Aberrant PRDM9 expression impacts the pan-cancer genomic landscape.

The binding of PRDM9 to chromatin is a key step in the induction of DNA double-strand breaks associated with meiotic recombination hotspots; it is normally expressed solely in germ cells. We interrogated 1879 cancer samples in 39 different cancer types and found that PRDM9 is unexpectedly expressed in 20% of these tumors even after stringent gene homology correction. The expression levels of PRDM9 in tumors are significantly higher than those found in healthy neighboring tissues and in healthy nongerm tissue databases. Recurrently mutated regions located within 5 Mb of the PRDM9 loci, as well as differentially expressed genes in meiotic pathways, correlate with PRDM9 expression. In samples with aberrant PRDM9 expression, structural variant breakpoints frequently neighbor the DNA motif recognized by PRDM9, and there is an enrichment of structural variants at sites of known meiotic PRDM9 activity. This study is the first to provide evidence of an association between aberrant expression of the meiosis-specific gene PRDM9 with genomic instability in cancer.

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.

Cyto-nuclear discordance suggests complex evolutionary history in the cave-dwelling salamander, Eurycea lucifuga.

Our understanding of the evolutionary history and ecology of cave-associated species has been driven historically by studies of morphologically adapted cave-restricted species. Our understanding of the evolutionary history and ecology of nonrestricted cave species, troglophiles, is limited to a few studies, which present differing accounts of troglophiles' relationship with the cave habitat, and its impact on population dynamics. Here, we used phylogenetics, demographic statistics, and population genetic methods to study lineage divergence, dates of divergence, and population structure in the Cave Salamander, Eurycea lucifuga, across its range. In order to perform these analyses, we sampled 233 individuals from 49 populations, using sequence data from three gene loci as well as genotyping data from 19 newly designed microsatellite markers. We find, as in many other species studied in a phylogeographic context, discordance between patterns inferred from mitochondrial relationships and those inferred by nuclear markers indicating a complicated evolutionary history in this species. Our results suggest Pleistocene-based divergence among three main lineages within E. lucifuga corresponding to the western, central, and eastern regions of the range, similar to patterns seen in species separated in multiple refugia during climatic shifts. The conflict between mitochondrial and nuclear patterns is consistent with what we would expect from secondary contact between regional populations following expansion from multiple refugia.