Nucleotide diversity inflation as a genome-wide response to experimental lifespan extension in Drosophila melanogaster

Background

Evolutionary theory predicts that antagonistically selected alleles, such as those with divergent pleiotropic effects in early and late life, may often reach intermediate population frequencies due to balancing selection, an elusive process when sought out empirically. Alternatively, genetic diversity may increase as a result of positive frequency-dependent selection and genetic purging in bottlenecked populations.

Results

While experimental evolution systems with directional phenotypic selection typically result in at least local heterozygosity loss, we report that selection for increased lifespan in Drosophila melanogaster leads to an extensive genome-wide increase of nucleotide diversity in the selected lines compared to replicate control lines, pronounced in regions with no or low recombination, such as chromosome 4 and centromere neighborhoods. These changes, particularly in coding sequences, are most consistent with the operation of balancing selection and the antagonistic pleiotropy theory of aging and life history traits that tend to be intercorrelated. Genes involved in antioxidant defenses, along with multiple lncRNAs, were among those most affected by balancing selection. Despite the overwhelming genetic diversification and the paucity of selective sweep regions, two genes with functions important for central nervous system and memory, Ptp10D and Ank2, evolved under positive selection in the longevity lines.

Conclusions

Overall, the ‘evolve-and-resequence’ experimental approach proves successful in providing unique insights into the complex evolutionary dynamics of genomic regions responsible for longevity.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3485-0) contains supplementary material, which is available to authorized users.

Nucleolar Dominance and Repression of 45S Ribosomal RNA Genes in Hybrids between Xenopus borealis and X. muelleri (2n = 36)

Nucleolar dominance is a dramatic disruption in the formation of nucleoli and the expression of ribosomal RNA (rRNA) genes, characteristic of some plant and animal hybrids. Here, we report that F1 hybrids produced from reciprocal crosses between 2 sister species of Xenopus clawed frogs, X. muelleri and X. borealis, undergo nucleolar dominance somewhat distinct from a pattern previously reported in hybrids between phylogenetically more distant Xenopus species. Patterns of nucleolar development, 45S rRNA expression, and gene copy inheritance were investigated using a combination of immunostaining, pyrosequencing, droplet digital PCR, flow cytometry, and epigenetic inhibition. In X. muelleri × X. borealis hybrids, typically only 1 nucleolus is formed, and 45S rRNA genes are predominantly expressed from 1 progenitor's alleles, X. muelleri, regardless of the cross-direction. These changes are accompanied by an extensive (∼80%) loss of rRNA gene copies in the hybrids relative to their parents, with the transcriptionally underdominant variant (X. borealis) being preferentially lost. Chemical treatment of hybrid larvae with a histone deacetylase inhibitor resulted in a partial derepression of the underdominant variant. Together, these observations shed light on the genetic and epigenetic basis of nucleolar dominance as an underappreciated manifestation of genetic conflicts within a hybrid genome.

Genomic Signatures of Speciation in Sympatric and Allopatric Hawaiian Picture-Winged Drosophila

The Hawaiian archipelago provides a natural arena for understanding adaptive radiation and speciation. The Hawaiian Drosophila are one of the most diverse endemic groups in Hawaiì with up to 1,000 species. We sequenced and analyzed entire genomes of recently diverged species of Hawaiian picture-winged Drosophila, Drosophila silvestris and Drosophila heteroneura from Hawaiì Island, in comparison with Drosophila planitibia, their sister species from Maui, a neighboring island where a common ancestor of all three had likely occurred. Genome-wide single nucleotide polymorphism patterns suggest the more recent origin of D. silvestris and D. heteroneura, as well as a pervasive influence of positive selection on divergence of the three species, with the signatures of positive selection more prominent in sympatry than allopatry. Positively selected genes were significantly enriched for functional terms related to sensory detection and mating, suggesting that sexual selection played an important role in speciation of these species. In particular, sequence variation in Olfactory receptor and Gustatory receptor genes seems to play a major role in adaptive radiation in Hawaiian pictured-winged Drosophila.

Nucleolar dominance and maternal control of 45S rDNA expression

Using a system of interspecies hybrids, trihybrids, and recombinants with varying proportions of genomes from three distinct Xenopus species, we provide evidence for de novo epigenetic silencing of paternal 45 S ribosomal ribonucleic acid (rRNA) genes and their species-dependent expression dominance that escapes transcriptional inactivation after homologous recombination. The same pattern of imprinting is maintained in the offspring from mothers being genetic males (ZZ) sex-reversed to females, indicating that maternal control of ribosomal deoxyribonucleic acid (rDNA) expression is not sex-chromosome linked. Nucleolar dominance (nucleolus underdevelopment) in Xenopus hybrids appears to be associated with a major non-Mendelian reduction in the number of 45 S rDNA gene copies rather than a specific pattern of their expression. The loss of rRNA gene copies in F1 hybrids was non-random with respect to the parental species, with the transcriptionally dominant variant preferentially removed from hybrid zygotes. This dramatic disruption in the structure and function of 45 S rDNA impacts transcriptome patterns of small nucleolar RNAs and messenger RNAs, with genes from the ribosome and oxidative stress pathways being among the most affected. Unorthodoxies of rDNA inheritance and expression may be interpreted as hallmarks of genetic conflicts between parental genomes, as well as defensive epigenetic mechanisms employed to restore genome integrity.

Experimental evolution of recombination and crossover interference in Drosophila caused by directional selection for stress-related traits

BACKGROUND

Population genetics predicts that tight linkage between new and/or pre-existing beneficial and deleterious alleles should decrease the efficiency of natural selection in finite populations. By decoupling beneficial and deleterious alleles and facilitating the combination of beneficial alleles, recombination accelerates the formation of high-fitness genotypes. This may impose indirect selection for increased recombination. Despite the progress in theoretical understanding, interplay between recombination and selection remains a controversial issue in evolutionary biology. Even less satisfactory is the situation with crossover interference, which is a deviation of double-crossover frequency in a pair of adjacent intervals from the product of recombination rates in the two intervals expected on the assumption of crossover independence. Here, we report substantial changes in recombination and interference in three long-term directional selection experiments with Drosophila melanogaster: for desiccation (~50 generations), hypoxia, and hyperoxia tolerance (>200 generations each).

RESULTS

For all three experiments, we found a high interval-specific increase of recombination frequencies in selection lines (up to 40-50% per interval) compared to the control lines. We also discovered a profound effect of selection on interference as expressed by an increased frequency of double crossovers in selection lines. Our results show that changes in interference are not necessarily coupled with increased recombination.

CONCLUSIONS

Our results support the theoretical predictions that adaptation to a new environment can promote evolution toward higher recombination. Moreover, this is the first evidence of selection for different recombination-unrelated traits potentially leading, not only to evolution toward increased crossover rates, but also to changes in crossover interference, one of the fundamental features of recombination.

Insights into the evolution of longevity from the bowhead whale genome

The bowhead whale (Balaena mysticetus) is estimated to live over 200 years and is possibly the longest-living mammal. These animals should possess protective molecular adaptations relevant to age-related diseases, particularly cancer. Here, we report the sequencing and comparative analysis of the bowhead whale genome and two transcriptomes from different populations. Our analysis identifies genes under positive selection and bowhead-specific mutations in genes linked to cancer and aging. In addition, we identify gene gain and loss involving genes associated with DNA repair, cell-cycle regulation, cancer, and aging. Our results expand our understanding of the evolution of mammalian longevity and suggest possible players involved in adaptive genetic changes conferring cancer resistance. We also found potentially relevant changes in genes related to additional processes, including thermoregulation, sensory perception, dietary adaptations, and immune response. Our data are made available online (http://www.bowhead-whale.org) to facilitate research in this long-lived species.

Modeling the Regulatory Mechanisms by Which NLRX1 Modulates Innate Immune Responses to Helicobacter pylori Infection

Helicobacter pylori colonizes half of the world’s population as the dominant member of the gastric microbiota resulting in a lifelong chronic infection. Host responses toward the bacterium can result in asymptomatic, pathogenic or even favorable health outcomes; however, mechanisms underlying the dual role of H. pylori as a commensal versus pathogenic organism are not well characterized. Recent evidence suggests mononuclear phagocytes are largely involved in shaping dominant immunity during infection mediating the balance between host tolerance and succumbing to overt disease. We combined computational modeling, bioinformatics and experimental validation in order to investigate interactions between macrophages and intracellular H. pylori. Global transcriptomic analysis on bone marrow-derived macrophages (BMDM) in a gentamycin protection assay at six time points unveiled the presence of three sequential host response waves: an early transient regulatory gene module followed by sustained and late effector responses. Kinetic behaviors of pattern recognition receptors (PRRs) are linked to differential expression of spatiotemporal response waves and function to induce effector immunity through extracellular and intracellular detection of H. pylori. We report that bacterial interaction with the host intracellular environment caused significant suppression of regulatory NLRC3 and NLRX1 in a pattern inverse to early regulatory responses. To further delineate complex immune responses and pathway crosstalk between effector and regulatory PRRs, we built a computational model calibrated using time-series RNAseq data. Our validated computational hypotheses are that: 1) NLRX1 expression regulates bacterial burden in macrophages; and 2) early host response cytokines down-regulate NLRX1 expression through a negative feedback circuit. This paper applies modeling approaches to characterize the regulatory role of NLRX1 in mechanisms of host tolerance employed by macrophages to respond to and/or to co-exist with intracellular H. pylori.

Macrophages modulate bacterial persistence and gastric pathology during Helicobacter pylori infection (INC4P.346)

Based on preliminary data suggesting that macrophages are critical regulators of Helicobacter pylori colonization and gastric pathology, we investigated how macrophage phenotypes shape the outcome of infection and persistence. The loss of PPARγ in the myeloid compartment favors a 5 to 10-fold decrease in bacterial loads at the expense of more severe gastric lesions. PPARγ-deficient macrophages present altered control of transcription factors such as NFκB, which results in a pro-inflammatory phenotype. Stomach lamina propria analyses at weeks 0 to 4 post-infection revealed increased levels of CD11b+ F4/80hi CD64+ CR3CR1+ macrophages in WT mice, while PPARγKO mice failed to expand and maintain such population. Macrophage depletion using clodronate liposomes resulted in a significant reduction of gastric H. pylori in WT mice, thus abrogating the differences in bacterial loads observed between WT and PPARγKO mice. In vitro co-culture of H. pylori with bone marrow derived macrophages showed that either PPARγ loss or pharmacological blockade enhances bactericidal activity, which is associated to significant differential expression of chil1, etv5, iigp1, ptger4, sqle, osm, hspa2 and rptoros levels as revealed by global transcriptome analyses. Hence, macrophages facilitate H. pylori infection by 1) serving as bacterial reservoirs and allowing intracellular replication, and 2) favoring a gastric regulatory response that favors persistence.

Modulation of immune responses toward Helicobacter pylori infection by NLRs (INM6P.339)

Helicobacter pylori (HP) colonizes 50% of the world’s population resulting in a decades-long gastric infection. Bacterial interaction with host intracellular environment occurs via injection of bacterial components through a TIVSS or intracellular replication. HP has been recognized for its ability to modulate intracellular NOD-like receptors (NLR). Host responses toward the bacterium can result in asymptomatic, pathogenic or even favorable immunity. Mechanisms underlying the dual role of HP as a commensal versus pathogen are not completely understood. We combined computational modeling, bioinformatics and experimental validation to investigate intracellular host-HP interactions. Global transcriptomic analysis on bone marrow-derived macrophages (BMDM) in a gentamycin protection assay unveiled that intracellular colonization of HP upregulated NOD1, NOD2, NLRP3, NLRC5 and inflammasome components (Caspase-1 and -11) but suppressed regulatory NLRX1 which was inversely correlated to TRAF6, NF-B, proinflammatory cytokines and reactive oxygen species. Loss of NLRX1 facilitates bacterial clearance in BMDM and infected mice. Lastly, we constructed a computational model to shed light on complex immune responses and pathway crosstalk regulated by NLRX1 during infection. In conclusion, NLRX1 is associated with chronic bacterial persistence during H. pylori infection and it may represent an immune evasion mechanism employed by the bacterium to facilitate long-term host colonization.

Cell size and cancer: a new solution to Peto's paradox?

Cancer, one of the leading health concerns for humans, is by no means a human-unique malady. Accumulating evidence shows that cancer kills domestic and wild animals at a similar rate to humans and can even pose a conservation threat to certain species. Assuming that each physiologically active and proliferating cell is at risk of malignant transformation, any evolutionary increase in the number of cells (and thus body mass) will lead to a higher cancer frequency, all else being equal. However, available data fail to support the prediction that bigger animals are affected by cancer more than smaller ones. The unexpected lack of correlation between body size (and life span) and cancer risk across taxa was dubbed Peto's paradox. In this perspective, several plausible explanations of Peto's paradox are presented, with the emphasis on a largely underappreciated relation of cell size to both metabolism and cell division rates across species, which we believe are key factors underlying the paradox. We conclude that larger organisms have bigger and slowly dividing cells with lower energy turnover, all significantly reducing the risk of cancer initiation. Solving Peto's paradox will enhance our understanding the evolution of cancer and may provide new implications for cancer prevention and treatment.

Host-defense peptides from skin secretions of the octoploid frogs Xenopus vestitus and Xenopus wittei (Pipidae): insights into evolutionary relationships

The primary structures of host-defense peptides have proved useful in elucidating the evolution history of frogs. Peptidomic analysis was used to compare the diversity of host-defense peptides in norepinephrine-stimulated skin secretions from the octoploid frogs, Xenopus vestitus (Kivu clawed frog) and Xenopus wittei (De Witte's clawed frog) in the family Pipidae. Structural characterization demonstrated that the X. vestitus peptides belong to the magainin (3 peptides), peptide glycine-leucine-amide (PGLa; 4 peptides), xenopsin-precursor fragment (XPF; 1 peptide), and caerulein-precursor fragment (CPF; 5 peptides) families. The X. wittei peptides comprise magainin (4 peptides), PGLa (1 peptide), XPF (2 peptides), and CPF (7 peptides). In addition, secretions from both species contain caerulein, identical to the peptide from Xenopus laevis, but X. wittei secretions contains the novel peptide [R4K]xenopsin. The variability in the numbers of paralogs in each peptide family indicates a selective silencing of the host-defense peptide genes following the polyploidization events. The primary structures of the peptides provide insight into phylogenetic relationships among the octoploid Xenopus frogs. The data support a sister-group relationship between X. vestitus and Xenopus lenduensis, suggestive of bifurcating speciation after allopolyploidization, whereas X. wittei is more closely related to the Xenopus amieti-Xenopus andrei group suggesting a common tetraploid ancestor. Consistent with previous data, the CPF peptides showed the highest growth inhibitory activity against bacteria with CPF-W6 (GIGSLLAKAAKLAAGLV.NH2) combining high antimicrobial potency against Staphylococcus aureus (MIC=4 μM) with relatively low hemolytic activity (LC50=190 μM).

EvoCor: a platform for predicting functionally related genes using phylogenetic and expression profiles

The wealth of publicly available gene expression and genomic data provides unique opportunities for computational inference to discover groups of genes that function to control specific cellular processes. Such genes are likely to have co-evolved and be expressed in the same tissues and cells. Unfortunately, the expertise and computational resources required to compare tens of genomes and gene expression data sets make this type of analysis difficult for the average end-user. Here, we describe the implementation of a web server that predicts genes involved in affecting specific cellular processes together with a gene of interest. We termed the server 'EvoCor', to denote that it detects functional relationships among genes through evolutionary analysis and gene expression correlation. This web server integrates profiles of sequence divergence derived by a Hidden Markov Model (HMM) and tissue-wide gene expression patterns to determine putative functional linkages between pairs of genes. This server is easy to use and freely available at http://pilot-hmm.vbi.vt.edu/.

Immunoregulatory mechanisms of microRNAs during Clostridium difficult infection (MUC9P.819)

MicroRNAs are small non-coding RNA molecules that play pivotal roles in the development and functionality of innate and adaptive immune cells, making them essential during the regulation of bacterial, viral and other immune infections. The role of miRNAs during host-pathogen interactions has recently been highlighted. However, the mechanisms by which miRNAs exert immunoregulatory actions are poorly understood. Our aim is to investigate the impact of Clostidium difficile infection (CDI) on host regulatory and effector pathways by focusing on the disruption of mucosal homeostasis by miRNAs. Sequencing studies on colonic samples and mucosal CD4+ T cells revealed a consistent upregulation of miR146b during CDI in mice, which correlated with worsened disease severity and upregulated MCP-1, IL-1β, IL-6 and IL-17. Moreover, upregulation of miR146b correlated with increased levels of Th17 cells in the spleen, MLN and colonic lamina propria of CDI mice, while abrogating anti-inflammatory responses characterized by IL-10 production. In vivo miR146b inhibition by using locked nucleic acid resulted in accelerated recovery as well as a significant induction of IL-10 production by CD4+, CD8+ and NK T cells and reduced expression of IL-17 by CD4+RORγt+ T cells. C. difficile infection induces upregulation of miR146b at the gut mucosa that contributes to pathogenic Th17 responses and impared immunoregulation.

Evidence for maternal imprinting of 45S ribosomal RNA genes in Xenopus hybrids

We discovered that gene clusters of 45S ribosomal RNA in Xenopus hybrid frogs are maternally imprinted, similar to X chromosome inactivation in marsupial females. Paternal expression was partly restored after chemical inhibition of histone deacetylation during larval stages. This provides a new spectacular example of epigenetic silencing and first evidence of genomic imprinting in amphibians.

Genome differentiation of Drosophila melanogaster from a microclimate contrast in Evolution Canyon, Israel

The opposite slopes of "Evolution Canyon" in Israel have served as a natural model system of adaptation to a microclimate contrast. Long-term studies of Drosophila melanogaster populations inhabiting the canyon have exhibited significant interslope divergence in thermal and drought stress resistance, candidate genes, mobile elements, habitat choice, mating discrimination, and wing-shape variation, all despite close physical proximity of the contrasting habitats, as well as substantial interslope migration. To examine patterns of genetic differentiation at the genome-wide level, we used high coverage sequencing of the flies' genomes. A total of 572 genes were significantly different in allele frequency between the slopes, 106 out of which were associated with 74 significantly overrepresented gene ontology (GO) terms, particularly so with response to stimulus and developmental and reproductive processes, thus corroborating previous observations of interslope divergence in stress response, life history, and mating functions. There were at least 37 chromosomal "islands" of interslope divergence and low sequence polymorphism, plausible signatures of selective sweeps, more abundant in flies derived from one (north-facing) of the slopes. Positive correlation between local recombination rate and the level of nucleotide polymorphism was also found.

Vitamin D receptor regulates mucosal immune responses to Clostridium difficile infection (P3076)

Clostridium difficile is an anaerobic bacterium that has re-emerged as a facultative pathogen and can cause diarrhea, colitis or even death. Vitamin D deficiency contributes to a 5-fold increase in healthcare costs during C. difficile infection (CDI). We investigated the role of vitamin D signaling during CDI by using vitamin D receptor (VDR) -/- and wild-type (WT) mice. The loss of VDR increased disease activity, weight loss, and the severity of colonic inflammatory lesions following CDI. RNAseq analyses showed upregulation of sentrin/SUMO-specific proteases (SENPs), downregulation of phosphatidylinositol 3-kinase regulatory subunit and cAMP-specific 3’5’-cyclic phophodiesterase in colons of VDR-/- mice. SENPs mediate deSUMOylation reactions that can affect both the transcriptional activity of nuclear receptors such as VDR and peroxisome proliferator-activated receptor (PPAR) γ as well as their ability to antagonize NF-kB activity. These results provide a mechanistic basis for the upregulation of IFNγ, IL-6, IL-17 and IL-1β in colons of VDR-/- mice infected with C. difficile. In addition, miR-760-5p was down-regulated in VDR-/- mice and upregulated in WT mice after CDI. Of note, miR-760 has 128 mRNA targets, including lanthionine synthetase component C-like 2 (LANCL2), which is upstream of PPARγ. We provide novel molecular evidence in vivo suggesting novel mechanisms involving LANCL2 and PPARγ by which vitamin D regulates C. difficile-associated disease and immunopathology.

Modeling the role of peroxisome proliferator-activated receptor γ and microRNA-146 in mucosal immune responses to Clostridium difficile

Clostridium difficile is an anaerobic bacterium that has re-emerged as a facultative pathogen and can cause nosocomial diarrhea, colitis or even death. Peroxisome proliferator-activated receptor (PPAR) γ has been implicated in the prevention of inflammation in autoimmune and infectious diseases; however, its role in the immunoregulatory mechanisms modulating host responses to C. difficile and its toxins remains largely unknown. To characterize the role of PPARγ in C. difficile-associated disease (CDAD), immunity and gut pathology, we used a mouse model of C. difficile infection in wild-type and T cell-specific PPARγ null mice. The loss of PPARγ in T cells increased disease activity and colonic inflammatory lesions following C. difficile infection. Colonic expression of IL-17 was upregulated and IL-10 downregulated in colons of T cell-specific PPARγ null mice. Also, both the loss of PPARγ in T cells and C. difficile infection favored Th17 responses in spleen and colonic lamina propria of mice with CDAD. MicroRNA (miRNA)-sequencing analysis and RT-PCR validation indicated that miR-146b was significantly overexpressed and nuclear receptor co-activator 4 (NCOA4) suppressed in colons of C. difficile-infected mice. We next developed a computational model that predicts the upregulation of miR-146b, downregulation of the PPARγ co-activator NCOA4, and PPARγ, leading to upregulation of IL-17. Oral treatment of C. difficile-infected mice with the PPARγ agonist pioglitazone ameliorated colitis and suppressed pro-inflammatory gene expression. In conclusion, our data indicates that miRNA-146b and PPARγ activation may be implicated in the regulation of Th17 responses and colitis in C. difficile-infected mice.

Hybridization between the African clawed frogs Xenopus laevis and Xenopus muelleri (Pipidae) increases the multiplicity of antimicrobial peptides in skin secretions of female offspring

Peptidomic analysis was used to compare the distribution of host-defense peptides in norepinephrine-stimulated skin secretions from laboratory-generated female F1 hybrids of the common clawed frog Xenopus laevis (Daudin, 1802) and Mueller's clawed frog Xenopus muelleri (Peters, 1844) with the corresponding distribution in skin secretions from the parent species. A total of 18 peptides were identified in secretions from the hybrid frogs. Eleven peptides (magainin-1, magainin-2, CPF-1, CPF-3, CPF-4, CPF-5, CPF-6, CPF-7, XPF-1, XPF-2, and PGLa) were identified in secretions of both the hybrids and X. laevis. Four peptides (magainin-M1, XPF-M1, CPF-M1, and tigerinin-M1) were previously found in skin secretions of X. muelleri but magainin-M2 and CPF-M2 from X. muelleri were not detected. Three previously undescribed peptides (magainin-LM1, PGLa-LM1, and CPF-LM1) were purified from the secretions of the hybrid frogs that were not detected in secretions from either X. laevis or X. muelleri. Magainin-LM1 differs from magainin-2 from X. laevis by a single amino acid substitution (Gly(13)→Ala) but PGLa-LM1 and CPF-LM1 differ appreciably in structure from orthologs in the parent species. CPF-LM1 shows potent, broad-spectrum antimicrobial activity and is hemolytic. The data indicate that hybridization increases the multiplicity of skin host-defense peptides in skin secretions. As the female F1 hybrids are fertile, hybridization may represent an adaptive strategy among Xenopus species to increase protection against pathogenic microorganisms in the environment.

Modulation of immune responses to Clostridium difficile by peroxisome proliferator-activated receptor γ and miRNA-146b (164.14)

Clostridium difficile is typically a harmless anaerobic bacterium but has recently re-emerged as a pathogen that can cause nosocomial diarrhea, colitis and death. To investigate the role of peroxisome proliferator-activated receptor (PPAR) γ in modulating immune responses to C. difficile we have performed several studies using a mouse model of antibiotic-induced C. difficile-associated disease (CDAD). The loss of PPARγ in T cells increased disease activity, body weight loss, and colonic inflammatory lesions after infection. It also resulted in upregulated IL-17 and MCP-1, and downregulated IL-10 expression, suggesting that a Th17 phenotype predominates during CDAD in mice lacking T cell PPARγ. Treatment of C. difficile-infected mice with PPAR γ agonists ameliorates disease severity and colonic lesions. RNA-seq results in colonic specimens indicated that three miRNAs were significantly overexpressed in infected mice: mmu-miR-146b, mmu-miR-1940, and mmu-miR-1298 (FDR P < 0.05). Real-time PCR results validated that infection results in overexpression of miR-146b, a molecule involved in regulating immunity and inflammation. Furthermore, NCOA4, a co-activator of PPARγ and target of miR-146b, was down-regulated in colons of infected mice. We provide novel evidence supporting a role for PPARγ in regulating C. difficile-related inflammation and immunopathology and explore the potential role of miRNA in modulating host responses to C. difficile.

Genetic divergence between freshwater and marine morphs of alewife (Alosa pseudoharengus): a 'next-generation' sequencing analysis

Alewife Alosa pseudoharengus, a small clupeid fish native to Atlantic Ocean, has recently (∼150 years ago) invaded the North American Great Lakes and despite challenges of freshwater environment its populations exploded and disrupted local food web structures. This range expansion has been accompanied by dramatic changes at all levels of organization. Growth rates, size at maturation, or fecundity are only a few of the most distinct morphological and life history traits that contrast the two alewife morphs. A question arises to what extent these rapidly evolving differences between marine and freshwater varieties result from regulatory (including phenotypic plasticity) or structural mutations. To gain insights into expression changes and sequence divergence between marine and freshwater alewives, we sequenced transcriptomes of individuals from Lake Michigan and Atlantic Ocean. Population specific single nucleotide polymorphisms were rare but interestingly occurred in sequences of genes that also tended to show large differences in expression. Our results show that the striking phenotypic divergence between anadromous and lake alewives can be attributed to massive regulatory modifications rather than coding changes.

Misexpression of testicular microRNA in sterile Xenopus hybrids points to tetrapod-specific microRNAs associated with male fertility

Spermatogenesis is one of the most complex biological processes undergone by any organism, making it susceptible to perturbations that result in male sterility. Research has demonstrated that mutant phenotypes can be obtained from the disruption of epigenetic modifications, which are commonly microRNA guided. Employing the Xenopus system, whereby homogametic interspecies males are always sterile, thus violating Haldane's Rule, we deep-sequenced testes-specific small-RNAs to identify microRNAs most frequently misexpressed between sterile hybrids and their fertile parental taxa. Using these data, we cross-referenced our expression information with previously published mouse (Mus musculus) data and identified a subset of seven microRNAs common to both (miR-338, miR-222, miR-18, miR-30, miR-10, miR-196, and miR-365). We propose that these microRNAs are likely critical for spermatogenesis in all tetrapods, having retained testicular expression across ~350 million years of evolution (Amphibian-Mammal split). Gene targets of six of these microRNAs are known, and all the six associate with zinc and zinc finger proteins (both previously found critical in male fertility), and three with Hox genes (some of which have also previously been deemed critical for testicular development and male fertility). Expression information for these targets revealed that all those associated with zinc have previously been found to express in mammalian testes. One Hox target has known mammalian testicular expression, two have close relatives with known mammalian testicular expression, and two more are associated with proteins known to have mammalian testicular expression. In addition, miR-222 has prior association with spermatogenesis, and miR-30 has been found to be abundantly expressed in both mouse and human testes.

Snapshot of DNA methylation changes associated with hybridization in Xenopus

Hybridization often results in dramatic genome reconfigurations including epigenetic changes that control gene expression. Here we survey methylation patterns of interspecific Xenopus F1 hybrids relative to parental species X. laevis and X. muelleri, using methyl-sensitive amplification polymorphisms (MSAPs). Out of a total of 546 MSAP markers, 364 were effective in elucidating the difference in methylation patterns between the hybrids and the parental species. Principal coordinate analysis of methylated fragments revealed four distinct clusters with the two parental species separate from hybrid males and females. On average, hybrids were characterized by a higher proportion (70.6%) of methylated fragments compared with the parental species (64.5%), and this difference was consistent with previously observed disruptions of hybrid transcriptomes. The proportion of methylated fragments did not correlate with variation in genome size, as measured with flow cytometry. The levels of methylation in sterile hybrid males (73.8%) were higher than in fertile hybrid females (68.6%), but this difference was not statistically significant. A total of 76 methylated fragments (20.9%) were hybrid-unique, presumably originating from methylation alterations in hybrid genomes.