Balancing selection maintains ancient polymorphisms at conserved enhancers for the olfactory receptor genes of a Chinese marine fish.

The study of balancing selection, as a selective force maintaining adaptive genetic variation in gene pools longer than expected by drift, is currently experiencing renewed interest due to the increased availability of new data, methods of analysis, and case studies. In this investigation, evidence of balancing selection operating on conserved enhancers of the olfactory receptor (OR) genes is presented for the Chinese sleeper (Bostrychus sinensis), a coastal marine fish that is emerging as a model species for evolutionary studies in the Northwest Pacific marginal seas. Coupled with tests for Gene Ontology enrichment and transcription factor binding, population genomic data allow for the identification of an OR cluster in the sleeper with a downstream flanking region containing three enhancers that are conserved with human and other fish species. Phylogenetic and population genetic analyses indicate that the enhancers are under balancing selection as evidenced by their translineage polymorphisms, excess common alleles, and increased within-group diversities. Age comparisons between the translineage polymorphisms and most recent common ancestors of neutral genealogies substantiate that the former are old, and thus, due to ancient balancing selection. The survival and reproduction of vertebrates depend on their sense of smell, and thereby, on their ORs. In addition to locus duplication and allelic variation of structural genes, this study highlights a third mechanism by which receptor diversity can be achieved for detecting and responding to the huge variety of environmental odorants (i.e., by balancing selection acting on OR gene expression through their enhancer variability).

Explaining Pathogenicity of Congenital Zika and Guillain-Barré Syndromes: Does Dysregulation of RNA Editing Play a Role?

Previous studies of Zika virus (ZIKV) pathogenesis have focused primarily on virus-driven pathology and neurotoxicity, as well as host-related changes in cell proliferation, autophagy, immunity, and uterine function. It is now hypothesized that ZIKV pathogenesis arises instead as an (unintended) consequence of host innate immunity, specifically, as the side effect of an otherwise well-functioning machine. The hypothesis presented here suggests a new way of thinking about the role of host immune mechanisms in disease pathogenesis, focusing on dysregulation of post-transcriptional RNA editing as a candidate driver of a broad range of observed neurodevelopmental defects and neurodegenerative clinical symptoms in both infants and adults linked with ZIKV infections. The authors collect and synthesize existing evidence of ZIKV-mediated changes in the expression of adenosine deaminases acting on RNA (ADARs), known links between abnormal RNA editing and pathogenesis, as well as ideas for future research directions, including potential treatment strategies.

Strong phenotypic divergence in spite of low genetic structure in the endemic Mangrove Warbler subspecies (Setophaga petechia xanthotera) of Costa Rica.

Despite the enormous advances in genetics, links between phenotypes and genotypes have been made for only a few nonmodel organisms. However, such links can be essential to understand mechanisms of ecological speciation. The Costa Rican endemic Mangrove Warbler subspecies provides an excellent subject to study differentiation with gene flow, as it is distributed along a strong precipitation gradient on the Pacific coast with no strong geographic barriers to isolate populations. Mangrove Warbler populations could be subject to divergent selection driven by precipitation, which influences soil salinity levels, which in turn influences forest structure and food resources. We used single nucleotide polymorphisms (SNPs) and morphological traits to examine the balance between neutral genetic and phenotypic divergence to determine whether selection has acted on traits and genes with functions related to specific environmental variables. We present evidence showing: (a) associations between environmental variables and SNPs, identifying candidate genes related to bill morphology (BMP) and osmoregulation, (b) absence of population genetic structure in neutrally evolving markers, (c) divergence in bill size across the precipitation gradient, and (d) strong phenotypic differentiation (P ST) which largely exceeds neutral genetic differentiation (F ST) in bill size. Our results indicate an important role for salinity, forest structure, and resource availability in maintaining phenotypic divergence of Mangrove Warblers through natural selection. Our findings add to the growing body of literature identifying the processes involved in phenotypic differentiation along environmental gradients in the face of gene flow.

Diversity and Dynamics of "Candidatus Endobugula" and Other Symbiotic Bacteria in Chinese Populations of the Bryozoan, Bugula neritina.

Bugula neritina is a common invasive cosmopolitan bryozoan that harbors (like many sessile marine invertebrates) a symbiotic bacterial (SB) community. Among the SB of B. neritina, "Candidatus Endobugula sertula" continues to receive the greatest attention, because it is the source of bryostatins. The bryostatins are potent bioactive polyketides, which have been investigated for their therapeutic potential to treat various cancers, Alzheimer's disease, and AIDS. In this study, we compare the metagenomics sequences for the 16S ribosomal RNA gene of the SB communities from different geographic and life cycle samples of Chinese B. neritina. Using a variety of approaches for estimating alpha/beta diversity and taxonomic abundance, we find that the SB communities vary geographically with invertebrate and fish mariculture and with latitude and environmental temperature. During the B. neritina life cycle, we find that the diversity and taxonomic abundances of the SB communities change with the onset of host metamorphosis, filter feeding, colony formation, reproduction, and increased bryostatin production. "Ca. Endobugula sertula" is confirmed as the symbiont of the Chinese "Ca. Endobugula"/B. neritina symbiosis. Our study extends our knowledge about B. neritina symbiosis from the New to the Old World and offers new insights into the environmental and life cycle factors that can influence its SB communities, "Ca. Endobugula," and bryostatins more globally.

Correction : Characterization of hybridization within a secondary contact region of the inshore fish, Bostrychus sinensis, in the East China Sea.

Following the publication of this article, the authors have requested that the Acknowledgements section be amended to thank Weidi Yang for his assistance with their Bostrychus sinensis photograph that was chosen for the front cover of the January 2018 issue of the journal. This error has been corrected in both the PDF and HTML versions of the paper. Also, the legends for Supplementary Figures 1 and 2 were not posted online. This error has been corrected in the HTML version of the paper.

Characterization of hybridization within a secondary contact region of the inshore fish, Bostrychus sinensis, in the East China Sea.

The northwest Pacific marginal seas are a primary center of phylogeographic and evolutionary research, because of their dynamic geographic history of falling and rising sea levels during the glaciations and interglaciations of the last one million years. Here we present new molecular and morphological data for geographic samples of the four-eyed sleeper (Bostrychus sinensis), which reinforce the evidence for secondary contact and hybridization between two phylogeographic lineages in the East China Sea. Specifically, we find that the secondary contact region is characterized by a low frequency of hybridization, where mitochondrial DNA introgression is relatively common, whereas F1 hybrids are correspondingly scarce. Furthermore, the adult standard lengths of the two phylogeographic lineages vary geographically in a manner that is consistent with reproductive character displacement. Collectively, the molecular and morphological data document that sleeper hybridization conforms to the classic "tension zone" model, where alleles are lost via reduced hybrid viability and/or positive assortative mating but are then replenished by dispersal from south of the secondary contact region. They also indicate that the two phylogeographic lineages are at an incipient stage of the speciation process. These results and conclusions for the four-eyed sleeper are presented as a case study for future research on the vicariance, secondary contact, and hybridization of marine groups in the northwest Pacific marginal seas.

RNA editing by the host ADAR system affects the molecular evolution of the Zika virus.

Zika virus (ZIKV) is a mosquito-transmitted flavivirus, linked to microcephaly and fetal death in humans. Here, we investigate whether host-mediated RNA editing of adenosines (ADAR) plays a role in the molecular evolution of ZIKV. Using complete coding sequences for the ZIKV polyprotein, we show that potential ADAR substitutions are underrepresented at the ADAR-resistant GA dinucleotides of both the positive and negative strands, that these changes are spatially and temporally clustered (as expected of ADAR editing) for three branches of the viral phylogeny, and that ADAR mutagenesis can be linked to its codon usage. Furthermore, resistant GA dinucleotides are enriched on the positive (but not negative) strand, indicating that the former is under stronger purifying selection than the latter. ADAR editing also affects the evolution of the rhabdovirus sigma. Our study now documents that host ADAR editing is a mutation and evolutionary force of positive- as well as negative-strand RNA viruses.

Role of Host-Driven Mutagenesis in Determining Genome Evolution of Sigma Virus (DMelSV; Rhabdoviridae) in Drosophila melanogaster.

Sigma virus (DMelSV) is ubiquitous in natural populations of Drosophila melanogaster. Host-mediated, selective RNA editing of adenosines to inosines (ADAR) may contribute to control of viral infection by preventing transcripts from being transported into the cytoplasm or being translated accurately; or by increasing the viral genomic mutation rate. Previous PCR-based studies showed that ADAR mutations occur in DMelSV at low frequency. Here we use SOLiDTM deep sequencing of flies from a single host population from Athens, GA, USA to comprehensively evaluate patterns of sequence variation in DMelSV with respect to ADAR. GA dinucleotides, which are weak targets of ADAR, are strongly overrepresented in the positive strand of the virus, consistent with selection to generate ADAR resistance on this complement of the transient, double-stranded RNA intermediate in replication and transcription. Potential ADAR sites in a worldwide sample of viruses are more likely to be "resistant" if the sites do not vary among samples. Either variable sites are less constrained and hence are subject to weaker selection than conserved sites, or the variation is driven by ADAR. We also find evidence of mutations segregating within hosts, hereafter referred to as hypervariable sites. Some of these sites were variable only in one or two flies (i.e., rare); others were shared by four or even all five of the flies (i.e., common). Rare and common hypervariable sites were indistinguishable with respect to susceptibility to ADAR; however, polymorphism in rare sites were more likely to be consistent with the action of ADAR than in common ones, again suggesting that ADAR is deleterious to the virus. Thus, in DMelSV, host mutagenesis is constraining viral evolution both within and between hosts.

Phylogeography of the inshore fish, Bostrychus sinensis, along the Pacific coastline of China.

This study assesses the phylogeography of the Chinese four-eyed sleeper (Bostrychus sinensis) with one mitochondrial and one nuclear genes and two morphological characters. Phylogenetic and population genetic analyses of the sequences reveals two phylogeographic lineages from the East and South China Seas, which are corroborated by the morphological data. The vicariance of the two lineages is attributed to the Pleistocene Ice Age exposure of the Taiwan Strait and consequent connection of Taiwan to the mainland, which thereby introduced an ecological barrier to gene flow between populations in the East and South China Seas. The distributions of the two lineages now overlap in the East China Sea and this secondary contact is attributed to biased northward migration along the two main currents of the Taiwan Strait following its interglacial re-flooding. In conclusion, this study reinforces the importance of "vicariance, then secondary contact" due to Late Pliocene and Pleistocene sea-level changes to the phylogeography of marine species. Specifically, it corroborates the importance of Pleistocene sea-level changes in the Taiwan Strait to the phylogeography of Chinese inshore species.

Transcriptome Changes during the Life Cycle of the Red Sponge, Mycale phyllophila (Porifera, Demospongiae, Poecilosclerida).

Sponges are an ancient metazoan group with broad ecological, evolutionary, and biotechnological importance. As in other marine invertebrates with a biphasic life cycle, the developing sponge undergoes a significant morphological, physiological, and ecological transformation during settlement and metamorphosis. In this study, we compare new transcriptome datasets for three life cycle stages of the red sponge (Mycale phyllophila) to test whether gene expression (as in the model poriferan, Amphimedon queenslandica) also varies more after settlement and metamorphosis. In contrast to A. queenslandica, we find that the transcriptome of M. phyllophila changes more during the earlier pre-competent larva/post-larva transition that spans these defining events. We also find that this transition is marked by a greater frequency of significantly up-regulated Gene Ontology terms including those for morphogenesis, differentiation, and development and that the transcriptomes of its pre-competent larvae and adult are distinct. The life cycle transcriptome variation between M. phyllophila and A. queenslandica may be due to their long separate evolutionary histories and corresponding differences in developmental rates and timing. This study now calls for new transcriptome datasets of M. phyllophila and other sponges, which will allow for tests of the generality of our life cycle expression differences and for the greater exploitation of poriferans in both basic and applied research.

Scombroid fishes provide novel insights into the trait/rate associations of molecular evolution.

The study of which life history traits primarily affect molecular evolutionary rates is often confounded by the covariance of these traits. Scombroid fishes (billfishes, tunas, barracudas, and their relatives) are unusual in that their mass-specific metabolic rate is positively associated with body size. This study exploits this atypical pattern of trait variation, which allows for direct tests of whether mass-specific metabolic rate or body size is the more important factor of molecular evolutionary rates. We inferred a phylogeny for scombroids from a supermatrix of molecular and morphological characters and used new phylogenetic comparative approaches to assess the associations of body size and mass-specific metabolic rate with substitution rate. As predicted by the body size hypothesis, there is a negative correlation between body size and substitution rate. However, unexpectedly, we also find a negative association between mass-specific metabolic and substitution rates. These relationships are supported by analyses of the total molecular data, separate mitochondrial and nuclear genes, and individual loci, and they are robust to phylogenetic uncertainty. The molecular evolutionary rates of scombroids are primarily tied to body size. This study demonstrates that groups with novel patterns of trait variation can be particularly informative for identifying which life history traits are the primary factors of molecular evolutionary rates.

Phylogenetic, expression, and functional analyses of anoctamin homologs in Caenorhabditis elegans.

Ca(2+)-activated Cl(-) channels (CaCCs) are critical to processes such as epithelial transport, membrane excitability, and signal transduction. Anoctamin, or TMEM16, is a family of 10 mammalian transmembrane proteins, 2 of which were recently shown to function as CaCCs. The functions of other family members have not been firmly established, and almost nothing is known about anoctamins in invertebrates. Therefore, we performed a phylogenetic analysis of anoctamins across the animal kingdom and examined the expression and function of anoctamins in the genetically tractable nematode Caenorhabditis elegans. Phylogenetic analyses support five anoctamin clades that are at least as old as the deuterostome/protosome ancestor. This includes a branch containing two Drosophila paralogs that group with mammalian ANO1 and ANO2, the two best characterized CaCCs. We identify two anoctamins in C. elegans (ANOH-1 and ANOH-2) that are also present in basal metazoans. The anoh-1 promoter is active in amphid sensory neurons that detect external chemical and nociceptive cues. Within amphid neurons, ANOH-1::GFP fusion protein is enriched within sensory cilia. RNA interference silencing of anoh-1 reduced avoidance of steep osmotic gradients without disrupting amphid cilia development, chemotaxis, or withdrawal from noxious stimuli, suggesting that ANOH-1 functions in a sensory mode-specific manner. The anoh-2 promoter is active in mechanoreceptive neurons and the spermatheca, but loss of anoh-2 had no effect on motility or brood size. Our study indicates that at least five anoctamin duplicates are evolutionarily ancient and suggests that sensory signaling may be a basal function of the anoctamin protein family.

Microsatellite DNA suggests that group size affects sex-biased dispersal patterns in red colobus monkeys.

Dispersal is a major life history trait of social organisms influencing the behavioral and genetic structure of their groups. Unfortunately, primate dispersal is difficult to quantify, because of the rarity of these events and our inability to ascertain if individuals dispersed or died when they disappear. Socioecological models have been partially developed to understand the ecological causes of different dispersal systems and their social consequences. However, these models have yielded confusing results when applied to folivores. The folivorous red colobus monkey (Procolobus rufomitratus) in Kibale National Park, Uganda is thought to exhibit female-biased dispersal, although both sexes have been observed to disperse and there remains considerable debate over the selective pressures favoring the transfers of males and females and the causes of variation in the proportion of each sex to leave the natal group. We circumvent this problem by using microsatellite DNA data to investigate the prediction that female dispersal will be more frequent in larger groups as compared to smaller ones. The rationale for this prediction is that red colobus exhibit increased within-group competition in bigger groups, which should favor higher female dispersal rates and ultimately lower female relatedness. Genetic data from two unequally sized neighboring groups of red colobus demonstrate increased female relatedness within the smaller group, suggesting females are less likely to disperse when there is less within-group competition. We suggest that the dispersal system is mediated to some degree by scramble competition and group size. Since red colobus group sizes have increased throughout Kibale by over 50% in the last decade, these changes may have major implications for the genetic structure and ultimately the population viability of this endangered primate.

A possible explanation for the population size discrepancy in tuna (genus Thunnus) estimated from mitochondrial DNA and microsatellite data.

A recent study using both mitochondrial DNA (mtDNA) and microsatellite data reported on a population size discrepancy in the eastern tiger salamander where the effective population size (N(e)) estimate of the former exceeded that of the latter. That study suggested, among other hypotheses, that homoplasy of microsatellite alleles is responsible for the discrepancy. In this investigation, we report 10 new cases of a similar discrepancy in five species of tuna. These cases derive from our Bayesian inferences using data from Pacific Bluefin Tuna (Thunnus orientalis) and Yellowfin Tuna (Thunnus albacares), as well as from published estimates of genetic diversity for additional populations of Yellowfin Tuna and three other tuna species. Phylogenetic character analyses of inferred genealogies of Pacific Bluefin and Yellowfin Tuna reveal similar reduced levels of mtDNA and microsatellite homoplasy. Thus, the discrepancy between inferred population sizes from mtDNA and microsatellite data in tuna is most likely not an artifact of the chosen mutation models used in the microsatellite analyses, but may reflect behavioral differences between the sexes such as female-biased philopatry and male-biased dispersal. This explanation now warrants critical testing with more local populations of tuna and with other animal and plant groups that have different life histories.

Primate DNA suggests long-term stability of an African rainforest.

Red colobus monkeys, due to their sensitivity to environmental change, are indicator species of the overall health of their tropical rainforest habitats. As a result of habitat loss and overhunting, they are among the most endangered primates in the world, with very few viable populations remaining. Traditionally, extant indicator species have been used to signify the conditions of their current habitats, but they have also been employed to track past environmental conditions by detecting previous population fluctuations. Kibale National Park (KNP) in Uganda harbors the only remaining unthreatened large population of red colobus. We used microsatellite DNA to evaluate the historical demography of these red colobus and, therefore, the long-term stability of their habitat. We find that the red colobus population throughout KNP has been stable for at least ∼40,000 years. We interpret this result as evidence of long-term forest stability because a change in the available habitat or population movement would have elicited a corresponding change in population size. We conclude that the forest of what is now Kibale National Park may have served as a Late Pleistocene refuge for many East African species.

Unique structure and regulation of the nematode detoxification gene regulator, SKN-1: implications to understanding and controlling drug resistance.

Nematodes parasitize an alarming number of people and agricultural animals globally and cause debilitating morbidity and mortality. Anthelmintics have been the primary tools used to control parasitic nematodes for the past several decades, but drug resistance is becoming a major obstacle. Xenobiotic detoxification pathways defend against drugs and other foreign chemicals in diverse organisms, and evidence is accumulating that they play a role in mediating resistance to anthelmintics in nematodes. Related antioxidation pathways may also provide filarial parasites with protection against host free-radical-mediated immune responses. Upstream regulatory pathways have received almost no attention in nematode parasites, despite their potential to coregulate multiple detoxification and antioxidation genes. The nuclear eurythroid 2-related factor 2 (NRF2) transcription factor mediates inducible detoxification and antioxidation defenses in mammals, and recent studies have demonstrated that it promotes multidrug resistance in some human tumors. Recent studies in the free-living model nematode, Caenorhabditis elegans, have defined the homologous transcription factor, SKN-1, as a master regulator of detoxification and antioxidation genes. Despite similar functions, SKN-1 and NRF2 have important differences in structure and regulatory pathways. Protein alignment and phylogenetic analyses indicate that these differences are shared among many nematodes, making SKN-1 a candidate for specifically targeting nematode detoxification and antioxidation.

Redesigning a library-based genetics class research project through instructional theory and authentic experience.

QUESTION: How can the library-based research project of a genetics course be reinvigorated and made sustainable without sacrificing educational integrity? SETTING: The University of Florida's Health Science Center Library provides the case study. METHODS: Since 1996, the librarian has codeveloped, supported, and graded all components of the project. In 2009, the project evolved from a single-authored paper to a group-work poster, with graded presentations hosted by the library. In 2010, students were surveyed regarding class enhancements. RESULTS: Responses indicated a preference for collaborative work and the poster format and suggested the changes facilitated learning. Instructors reported that the poster format more clearly documented students' understanding of genetics. CONCLUSION: Results suggest project enhancements contributed to greater appreciation, understanding, and application of classroom material and offered a unique and authentic learning experience, without compromising educational integrity. The library benefitted through increased visibility as a partner in the educational mission and development of a sustainable instructional collaboration.

The importance of taxon sampling in genomic studies: an example from the cyclooxygenases of teleost fishes.

Comparative genomic studies must often rely on single model species and exemplars to represent the genetic variation both within and among different major groups, because of technological, financial, and time constraints. This study of the cyclooxygenases from teleost fishes serves as a reminder that caution is required in these cases, since such incomplete taxon sampling can lead to errors in the interpretation and prediction of genome evolution, function, and structure.

A computer simulator for assessing different challenges and strategies of de novo sequence assembly.

This study presents a new computer program for assessing the effects of different factors and sequencing strategies on de novo sequence assembly. The program uses reads from actual sequencing studies or from simulations with a reference genome that may also be real or simulated. The simulated reads can be created with our read simulator. They can be of differing length and coverage, consist of paired reads with varying distance, and include sequencing errors such as color space miscalls to imitate SOLiD data. The simulated or real reads are mapped to their reference genome and our assembly simulator is then used to obtain optimal assemblies that are limited only by the distribution of repeats. By way of this mapping, the assembly simulator determines which contigs are theoretically possible, or conversely (and perhaps more importantly), which are not. We illustrate the application and utility of our new simulation tools with several experiments that test the effects of genome complexity (repeats), read length and coverage, word size in De Bruijn graph assembly, and alternative sequencing strategies (e.g., BAC pooling) on sequence assemblies. These experiments highlight just some of the uses of our simulators in the experimental design of sequencing projects and in the further development of assembly algorithms.