Divergent evolutionary trajectories shape the postmating transcriptional profiles of conspecifically and heterospecifically mated cactophilic Drosophila females.

Postmating-prezygotic (PMPZ) reproductive isolation is hypothesized to result from divergent coevolutionary trajectories of sexual selection and/or sexual conflict in isolated populations. However, the genetic basis of PMPZ incompatibilities between species is poorly understood. Here, we use a comparative framework to compare global gene expression in con- and heterospecifically mated Drosophila mojavensis and D. arizonae female reproductive tracts. We find striking divergence between the species in the female postmating transcriptional response to conspecific mating, including differences in differential expression (DE), alternative splicing (AS), and intron retention (IR). As predicted, heterospecific matings produce disrupted transcriptional profiles, but the overall patterns of misregulation are different between the reciprocal crosses. Moreover, we find a positive correlation between postmating transcriptional divergence between species and levels of transcriptional disruption in heterospecific crosses. This result indicates that mating responsive genes that have diverged more in expression also have more disrupted transcriptional profiles in heterospecifically mated females. Overall, our results provide insights into the evolution of PMPZ isolation and lay the foundation for future studies aimed at identifying specific genes involved in PMPZ incompatibilities and the evolutionary forces that have contributed to their divergence in closely related species.

Contributions of cis- and trans-Regulatory Evolution to Transcriptomic Divergence across Populations in the Drosophila mojavensis Larval Brain.

Natural selection on gene expression was originally predicted to result primarily in cis- rather than trans-regulatory evolution, due to the expectation of reduced pleiotropy. Despite this, numerous studies have ascribed recent evolutionary divergence in gene expression predominantly to trans-regulation. Performing RNA-seq on single isofemale lines from genetically distinct populations of the cactophilic fly Drosophila mojavensis and their F1 hybrids, we recapitulated this pattern in both larval brains and whole bodies. However, we demonstrate that improving the measurement of brain expression divergence between populations by using seven additional genotypes considerably reduces the estimate of trans-regulatory contributions to expression evolution. We argue that the finding of trans-regulatory predominance can result from biases due to environmental variation in expression or other sources of noise, and that cis-regulation is likely a greater contributor to transcriptional evolution across D. mojavensis populations. Lastly, we merge these lines of data to identify several previously hypothesized and intriguing novel candidate genes, and suggest that the integration of regulatory and population-level transcriptomic data can provide useful filters for the identification of potentially adaptive genes.

Assessing the Architecture of Drosophila mojavensis Locomotor Evolution with Bulk Segregant Analysis.

Behavior is frequently predicted to be especially important for evolution in novel environments. If these predictions are accurate, there might be particular patterns of genetic architecture associated with recently diverged behaviors. Specifically, it has been predicted that behaviors linked to population divergence should be underpinned by a few genes of relatively large effect, compared to architectures of intrapopulation behavioral variation, which is considered to be highly polygenic. More mapping studies of behavioral variation between recently diverged populations are needed to continue assessing the generality of these predictions. Here, we used a bulk segregant mapping approach to dissect the genetic architecture of a locomotor trait that has evolved between two populations of the cactophilic fly Drosophila mojavensis We created an F8 mapping population of 1,500 individuals from advanced intercross lines and sequenced the 10% of individuals with the highest and lowest levels of locomotor activity. Using three alternative statistical approaches, we found strong evidence for two relatively large-effect QTL that is localized in a region homologous to a region of densely packed behavior loci in Drosophila melanogaster, suggesting that clustering of behavior genes may display relatively deep evolutionary conservation. Broadly, our data are most consistent with a polygenic architecture, though with several loci explaining a high proportion of variation in comparison to similar behavioral traits. We further note the presence of several antagonistic QTL linked to locomotion and discuss these results in light of theories regarding behavioral evolution and the effect size and direction of QTL for diverging traits in general.

Behavioral evolution accompanying host shifts in cactophilic Drosophila larvae.

For plant utilizing insects, the shift to a novel host is generally accompanied by a complex set of phenotypic adaptations. Many such adaptations arise in response to differences in plant chemistry, competitive environment, or abiotic conditions. One less well-understood factor in the evolution of phytophagous insects is the selective environment provided by plant shape and volume. Does the physical structure of a new plant host favor certain phenotypes? Here, we use cactophilic Drosophila, which have colonized the necrotic tissues of cacti with dramatically different shapes and volumes, to examine this question. Specifically, we analyzed two behavioral traits in larvae, pupation height, and activity that we predicted might be related to the ability to utilize variably shaped hosts. We found that populations of D. mojavensis living on lengthy columnar or barrel cactus hosts have greater activity and pupate higher in a laboratory environment than populations living on small and flat prickly pear cactus cladodes. Crosses between the most phenotypically extreme populations suggest that the genetic architectures of these behaviors are distinct. A comparison of activity in additional cactophilic species that are specialized on small and large cactus hosts shows a consistent trend. Thus, we suggest that greater motility and an associated tendency to pupate higher in the laboratory are potential larval adaptations for life on a large plant where space is more abundant and resources may be more sparsely distributed.

Osteoarthritis-Like Changes in Bardet-Biedl Syndrome Mutant Ciliopathy Mice (Bbs1M390R/M390R): Evidence for a Role of Primary Cilia in Cartilage Homeostasis and Regulation of Inflammation.

Osteoarthritis (OA) is a debilitating inflammation related disease characterized by joint pain and effusion, loss of mobility, and deformity that may result in functional joint failure and significant impact on quality of life. Once thought of as a simple "wear and tear" disease, it is now widely recognized that OA has a considerable metabolic component and is related to chronic inflammation. Defects associated with primary cilia have been shown to be cause OA-like changes in Bardet-Biedl mice. We examined the role of dysfunctional primary cilia in OA in mice through the regulation of the previously identified degradative and pro-inflammatory molecular pathways common to OA. We observed an increase in the presence of pro-inflammatory markers TGFß-1 and HTRA1 as well as cartilage destructive protease MMP-13 but a decrease in DDR-2. We observed a morphological difference in cartilage thickness in Bbs1 M390R/M390R mice compared to wild type (WT). We did not observe any difference in OARSI or Mankin scores between WT and Bbs1M390R/M390R mice. Primary cilia appear to be involved in the upregulation of biomarkers, including pro-inflammatory markers common to OA.