Transcriptional effects of a positive feedback circuit in Drosophila melanogaster.

BACKGROUND: Synthetic systems that use positive feedback have been developed to control human disease vectors and crop pests. The tTAV system, which has been deployed in several insect species, relies on a positive feedback circuit that can be inhibited via dietary tetracycline. Although insects carrying tTAV fail to survive until adulthood in the absence of tetracycline, the exact reason for its lethality, as well as the transcriptomic effects of an active positive feedback circuit, remain unknown. RESULTS: We engineered the tTAV system in Drosophila melanogaster and investigated the effects of tTAV genome integration locus on the whole fly transcriptome during larval and adult life stages in four transgenic fly strains using gene expression microarrays. We found that while there were widespread effects on the transcriptome, the gene expression differences after removal of tetracycline were not consistent between integration sites. No specific region of the genome was affected, no common set of genes or pathways, nor did the integration site affect the transcripts in cis. CONCLUSION: Although the positive feedback tTAV system is effective at killing insect larvae regardless of where it is inserted in the genome, it does not exhibit a specific, consistent transcriptional signature. Instead, each insertion site is associated with broad, but different, transcriptional effects. Our results suggest that lethality may not be caused by a direct effect on transcription of a set of key genes or pathways. Instead, we propose that rather than a specific action of a tTAV protein, it is the stochastic transcriptional effects specific to each insertion site that contribute to the tTAV-induced mortality.

Stability of underdominant genetic polymorphisms in population networks.

Heterozygote disadvantage is potentially a potent driver of population genetic divergence. Also referred to as underdominance, this phenomena describes a situation where a genetic heterozygote has a lower overall fitness than either homozygote. Attention so far has mostly been given to underdominance within a single population and the maintenance of genetic differences between two populations exchanging migrants. Here we explore the dynamics of an underdominant system in a network of multiple discrete, yet interconnected, populations. Stability of genetic differences in response to increases in migration in various topological networks is assessed. The network topology can have a dominant and occasionally non-intuitive influence on the genetic stability of the system.

Convergent adaptation of human lactase persistence in Africa and Europe.

A SNP in the gene encoding lactase (LCT) (C/T-13910) is associated with the ability to digest milk as adults (lactase persistence) in Europeans, but the genetic basis of lactase persistence in Africans was previously unknown. We conducted a genotype-phenotype association study in 470 Tanzanians, Kenyans and Sudanese and identified three SNPs (G/C-14010, T/G-13915 and C/G-13907) that are associated with lactase persistence and that have derived alleles that significantly enhance transcription from the LCT promoter in vitro. These SNPs originated on different haplotype backgrounds from the European C/T-13910 SNP and from each other. Genotyping across a 3-Mb region demonstrated haplotype homozygosity extending >2.0 Mb on chromosomes carrying C-14010, consistent with a selective sweep over the past approximately 7,000 years. These data provide a marked example of convergent evolution due to strong selective pressure resulting from shared cultural traits-animal domestication and adult milk consumption.

Dynamics of a combined Medea-underdominant population transformation system.

BACKGROUND: Transgenic constructs intended to be stably established at high frequencies in wild populations have been demonstrated to "drive" from low frequencies in experimental insect populations. Linking such population transformation constructs to genes which render them unable to transmit pathogens could eventually be used to stop the spread of vector-borne diseases like malaria and dengue. RESULTS: Generally, population transformation constructs with only a single transgenic drive mechanism have been envisioned. Using a theoretical modelling approach we describe the predicted properties of a construct combining autosomal Medea and underdominant population transformation systems. We show that when combined they can exhibit synergistic properties which in broad circumstances surpass those of the single systems. CONCLUSION: With combined systems, intentional population transformation and its reversal can be achieved readily. Combined constructs also enhance the capacity to geographically restrict transgenic constructs to targeted populations. It is anticipated that these properties are likely to be of particular value in attracting regulatory approval and public acceptance of this novel technology.

Alignment-free estimation of nucleotide diversity.

MOTIVATION: Sequencing capacity is currently growing more rapidly than CPU speed, leading to an analysis bottleneck in many genome projects. Alignment-free sequence analysis methods tend to be more efficient than their alignment-based counterparts. They may, therefore, be important in the long run for keeping sequence analysis abreast with sequencing. RESULTS: We derive and implement an alignment-free estimator of the number of pairwise mismatches, . Our implementation of , pim, is based on an enhanced suffix array and inherits the superior time and memory efficiency of this data structure. Simulations demonstrate that is accurate if mutations are distributed randomly along the chromosome. While real data often deviates from this ideal, remains useful for identifying regions of low genetic diversity using a sliding window approach. We demonstrate this by applying it to the complete genomes of 37 strains of Drosophila melanogaster, and to the genomes of two closely related Drosophila species, D.simulans and D.sechellia. In both cases, we detect the diversity minimum and discuss its biological implications.

From genes to games: cooperation and cyclic dominance in meiotic drive.

Evolutionary change can be described on a genotypic level or a phenotypic level. Evolutionary game theory is typically thought of as a phenotypic approach, although it is frequently argued that it can also be used to describe population genetic evolution. Interpreting the interaction between alleles in a diploid genome as a two player game leads to interesting alternative perspectives on genetic evolution. Here we focus on the case of meiotic drive and illustrate how meiotic drive can be directly and precisely interpreted as a social dilemma, such as the prisoners dilemma or the snowdrift game, in which the drive allele takes more than its fair share. Resistance to meiotic drive can lead to the well understood cyclic dominance found in the rock-paper-scissors game. This perspective is well established for the replicator dynamics, but there is still considerable ground for mutual inspiration between the two fields. For example, evolutionary game theorists can benefit from considering the stochastic evolutionary dynamics arising from finite population size. Population geneticists can benefit from game theoretic tools and perspectives on genetic evolution.

Stability properties of underdominance in finite subdivided populations.

IN ISOLATED populations underdominance leads to bistable evolutionary dynamics: below a certain mutant allele frequency the wildtype succeeds. Above this point, the potentially underdominant mutant allele fixes. In subdivided populations with gene flow there can be stable states with coexistence of wildtypes and mutants: polymorphism can be maintained because of a migration-selection equilibrium, i.e., selection against rare recent immigrant alleles that tend to be heterozygous. We focus on the stochastic evolutionary dynamics of systems where demographic fluctuations in the coupled populations are the main source of internal noise. We discuss the influence of fitness, migration rate, and the relative sizes of two interacting populations on the mean extinction times of a group of potentially underdominant mutant alleles. We classify realistic initial conditions according to their impact on the stochastic extinction process. Even in small populations, where demographic fluctuations are large, stability properties predicted from deterministic dynamics show remarkable robustness. Fixation of the mutant allele becomes unlikely but the time to its extinction can be long.

Genome Sequence of Arthrobacter globiformis Phage KeAlii from Hawai'i.

Here, we report the genome sequence of bacteriophage KeAlii, a Siphoviridae that infects Arthrobacter globiformis strain B-2979, from Honolulu, Hawai'i. The 41,850-bp genome contains 66 predicted protein-coding genes and 1 gene that encodes a tRNA for tryptophan. Genome comparisons suggest KeAlii is closely related to actinobacteriophage Adolin.

Rapid formation of distinct hybrid lineages after secondary contact of two fish species (Cottus sp.).

Homoploid hybridization after secondary contact between related species can lead to mixtures of genotypes which have the potential for rapid adaptation to new environmental conditions. Here, we focus on a case where anthropogenic changes within the past 200 years have allowed the hybridization between two fish species (Cottus rhenanus and Cottus perifretum) in the Netherlands. Specifically, we address the question of the dynamics of the emergence of these hybrids and invasion of the river systems. Using a set of 81 mostly ancestry-informative SNP markers, as well as broad sample coverage in and around the area of the initial contact, we find a structured hybrid swarm with at least three distinct hybrid lineages that have emerged out of this secondary contact situation. We show that genetically coherent groups can occur at geographically distant locations, while geographically adjacent groups can be genetically different, indicating that some form of reproductive isolation between the lineages is already effective. Using a newly developed modelling approach, we test the relative influence of founding admixture, drift and migration on the allele compositions of the sampling sites. We find that the allele frequency distributions can best be explained if continued gene flow between the parental species and the hybrid lineages is invoked. Genome mapping of the invasive lineage in the Rhine shows that major chromosomal rearrangements were not involved in creating this distinct lineage. Our results show that hybridization after secondary contact can quickly lead to multiple independent new lineages that have the capacity to form hybrid species.

The collective-risk social dilemma and the prevention of simulated dangerous climate change.

Will a group of people reach a collective target through individual contributions when everyone suffers individually if the target is missed? This "collective-risk social dilemma" exists in various social scenarios, the globally most challenging one being the prevention of dangerous climate change. Reaching the collective target requires individual sacrifice, with benefits to all but no guarantee that others will also contribute. It even seems tempting to contribute less and save money to induce others to contribute more, hence the dilemma and the risk of failure. Here, we introduce the collective-risk social dilemma and simulate it in a controlled experiment: Will a group of people reach a fixed target sum through successive monetary contributions, when they know they will lose all their remaining money with a certain probability if they fail to reach the target sum? We find that, under high risk of simulated dangerous climate change, half of the groups succeed in reaching the target sum, whereas the others only marginally fail. When the risk of loss is only as high as the necessary average investment or even lower, the groups generally fail to reach the target sum. We conclude that one possible strategy to relieve the collective-risk dilemma in high-risk situations is to convince people that failure to invest enough is very likely to cause grave financial loss to the individual. Our analysis describes the social window humankind has to prevent dangerous climate change.

The genetic structure of Pacific Islanders.

Human genetic diversity in the Pacific has not been adequately sampled, particularly in Melanesia. As a result, population relationships there have been open to debate. A genome scan of autosomal markers (687 microsatellites and 203 insertions/deletions) on 952 individuals from 41 Pacific populations now provides the basis for understanding the remarkable nature of Melanesian variation, and for a more accurate comparison of these Pacific populations with previously studied groups from other regions. It also shows how textured human population variation can be in particular circumstances. Genetic diversity within individual Pacific populations is shown to be very low, while differentiation among Melanesian groups is high. Melanesian differentiation varies not only between islands, but also by island size and topographical complexity. The greatest distinctions are among the isolated groups in large island interiors, which are also the most internally homogeneous. The pattern loosely tracks language distinctions. Papuan-speaking groups are the most differentiated, and Austronesian or Oceanic-speaking groups, which tend to live along the coastlines, are more intermixed. A small "Austronesian" genetic signature (always <20%) was detected in less than half the Melanesian groups that speak Austronesian languages, and is entirely lacking in Papuan-speaking groups. Although the Polynesians are also distinctive, they tend to cluster with Micronesians, Taiwan Aborigines, and East Asians, and not Melanesians. These findings contribute to a resolution to the debates over Polynesian origins and their past interactions with Melanesians. With regard to genetics, the earlier studies had heavily relied on the evidence from single locus mitochondrial DNA or Y chromosome variation. Neither of these provided an unequivocal signal of phylogenetic relations or population intermixture proportions in the Pacific. Our analysis indicates the ancestors of Polynesians moved through Melanesia relatively rapidly and only intermixed to a very modest degree with the indigenous populations there.

Wild gut microbiomes reveal individuals, species, and location as drivers of variation in two critically endangered Hawaiian honeycreepers.

BACKGROUND: The gut microbiome of animals is an important component that has strong influence on the health, fitness, and behavior of its host. Most research in the microbiome field has focused on human populations and commercially important species. However, researchers are now considering the link between endangered species conservation and the microbiome. In Hawai'i, several threats (e.g., avian malaria and habitat loss) have caused widespread population declines of Hawaiian honeycreepers (subfamily: Carduelinae). These threats can have a significant effect on the avian gut microbiome and may even lead to disruption of microbial function. However, the gut microbiome of honeycreeper in the wild has yet to be explored. METHODS: We collected 13 and 42 fecal samples, respectively, from two critically endangered honeycreeper species, the 'akikiki (Oreomystis bairdi) and the 'akeke'e (Loxops caeruleirostris). The 16S rRNA gene was sequenced and processed though a MOTHUR-based bioinformatics pipeline. Bacterial ASVs were identified using the DADA2 program and bacterial community analyses, including alpha and beta diversity measures, were conducted using R packages Phyloseq and vegan. RESULTS: A total of 8,958 bacterial ASVs were identified from the fecal samples. Intraspecific differences in the gut microbiome among individual birds explained most of the variation present in the dataset, however differences between species did exist. Both species had distinct microbiomes with minimal overlap in beta diversity. 'Akikiki had a more diverse microbiome compared to 'akeke'e. Additionally, small but stastically significant differences in beta diversity also exist between sampling location and sexes in 'akikiki. CONCLUSION: 'Akikiki and 'akeke'e are currently the focus of captive breeding efforts and plans to translocate the two species to other islands are underway. This baseline knowledge will help inform management decisions for these honeycreeper species in their native habitats, on other islands, and in captivity.

Evaluation of Gene Knockouts by CRISPR as Potential Targets for the Genetic Engineering of the Mosquito Culex quinquefasciatus.

Culex quinquefasciatus mosquitoes are a globally widespread vector of several human and animal pathogens. Their biology and behavior allow them to thrive in proximity to urban areas, rendering them a constant public health threat. Their mixed bird/mammal feeding behavior further offers a vehicle for zoonotic pathogens transmission to people and, separately, poses a threat to the conservation of insular birds. The advent of CRISPR has led to the development of novel technologies for the genetic engineering of wild mosquito populations. Yet, research into Cx. quinquefasciatus has been lagging compared to other disease vectors. Here, we use this tool to disrupt a set of five pigmentation genes in Cx. quinquefasciatus that, when altered, lead to visible, homozygous-viable phenotypes. We further validate this approach in separate laboratories and in two distinct strains of Cx. quinquefasciatus that are relevant to potential future public health and bird conservation applications. We generate a double-mutant line, demonstrating the possibility of sequentially combining multiple such mutations in a single individual. Lastly, we target two loci, doublesex in the sex-determination pathway and proboscipedia, a hox gene, demonstrating the flexibility of these methods applied to novel targets. Our work provides a platform of seven validated loci that could be used for targeted mutagenesis in Cx. quinquefasciatus and the future development of genetic suppression strategies for this species. Furthermore, the mutant lines generated here could have widespread utility to the research community using this model organism, as they could be used as targets for transgene delivery, where a copy of the disrupted gene could be included as an easily scored transgenesis marker.

African human diversity, origins and migrations.

The continent of Africa is thought to be the site of origin of all modern humans and is the more recent origin of millions of African Americans. Although Africa has the highest levels of human genetic diversity both within and between populations, it is under-represented in studies of human genetics. Recent advances have been made in understanding the origins of modern humans within Africa, the rate of adaptations due to positive selection, the routes taken in the first migrations of modern humans out of Africa, and the degree of admixture with archaic populations. Africa is also in dire need of effective medical interventions, and studies of genetic variation in Africans will shed light on the genetic basis of diseases and resistance to infectious diseases. Thus, we have tremendous potential to learn about human variation and evolutionary history and to positively impact human health care from studies of genetic diversity in Africa.

Using underdominance to bi-stably transform local populations.

Underdominance refers to natural selection against individuals with a heterozygous genotype. Here, we analyze a single-locus underdominant system of two large local populations that exchange individuals at a certain migration rate. The system can be characterized by fixed points in the joint allele frequency space. We address the conditions under which underdominance can be applied to transform a local population that is receiving wildtype immigrants from another population. In a single population, underdominance has the benefit of complete removal of genetically modified alleles (reversibility) and coexistence is not stable. The two population system that exchanges migrants can result in internal stable states, where coexistence is maintained, but with additional release of wildtype individuals the system can be reversed to a fully wildtype state. This property is critically controlled by the migration rate. We approximate the critical minimum frequency required to result in a stable population transformation. We also concentrate on the destabilizing effects of fitness and migration rate asymmetry. Practical implications of our results are discussed in the context of utilizing underdominance to genetically modify wild populations. This is of importance especially for genetic pest management strategies, where locally stable and potentially reversible transformations of populations of disease vector species are of interest.

Mutation, selection and the future of human evolution.

Several recent analyses provide growing evidence of the influence of positive selection acting in the ancestors of modern humans. Additionally, the best way to explain current fluctuations in neutral variation across the genome is by including negative selection against a high rate of deleterious mutants. We suggest that explaining these predicted high deleterious mutation rates in humans could require the inclusion of additional factors, such as inbreeding and prezygotic selection, in addition to rank-order selection and fitness interactions among mutations. We also suggest that some forms of selection, rather than being relaxed in modern humans, are probably still acting and might intensify in the near future, and make some predictions about the next several millennia of human evolution.

Two-locus epistasis with sexually antagonistic selection: a genetic Parrondo's paradox.

An example is provided where, with antagonistic selection and epistatic interaction of alleles at two loci, an autosomal allele can rise in frequency, persist in the population, and even continue to fixation, despite having an apparently lower average fitness than the alternative allele, in a process similar to Parrondo's paradox.

Gene expression across tissues, sex, and life stages in the sea urchin Tripneustes gratilla [Toxopneustidae, Odontophora, Camarodonta].

The pan-tropical sea urchin Tripneustes gratilla is an ecologically and economically important shallow water algal grazer. The aquaculture of T. gratilla has spurred growing interest in the population biology of the species, and by extension the generation of more molecular resources. To this purpose, de novo transcriptomes of T. gratilla were generated for two adults, a male and a female, as well as for a cohort of approximately 1000 plutei larvae. Gene expression profiles of three adult tissue samples were quantified and compared. These samples were of gonadal tissue, the neural ring, and pooled tube feet and pedicellariae. Levels of shared and different gene expression between sexes, as well as across functional categories of interest, including the immune system, toxins, genes involved in fertilization, and sensory genes are highlighted. Differences in expression of isoforms between the sexes and Sex determining Region Y-related High Mobility Group box groups is observed. Additionally an expansion of the tumor suppressor DMBT1 is observed in T. gratilla when compared to the annotated genome of the sea urchin Strongylocentrotus purpuratus. The draft transcriptome of T. gratilla is presented here in order to facilitate more genomic level analysis of emerging model sea urchin systems.

Positive selection can create false hotspots of recombination.

Simulations of positive directional selection, under parameter values appropriate for approximating human genetic diversity and rates of recombination, reveal that the effects of strong selective sweeps on patterns of linkage disequilibrium (LD) mimic the pattern expected with recombinant hotspots.

Evaluation of real-time PCR amplification efficiencies to detect PCR inhibitors.

Real-time PCR analysis is a sensitive template DNA quantitation strategy that has recently gained considerable attention in the forensic community. However, the utility of real-time PCR methods extends beyond quantitation and allows for simultaneous evaluation of template DNA extraction quality. This study presents a computational method that allows analysts to identify problematic samples with statistical reliability by comparing the amplification efficiencies of unknown template DNA samples with clean standards. In this study, assays with varying concentrations of tannic acid are used to evaluate and adjust sample-specific amplification efficiency calculation methods in order to optimize their inhibitor detection capabilities. Kinetic outlier detection and prediction boundaries are calculated to identify amplification efficiency outliers. Sample-specific amplification efficiencies calculated over a four-cycle interval starting at the threshold cycle can be used to detect reliably the presence of 0.4 ng of tannic acid in a 25 microL PCR reaction. This approach provides analysts with a precise measure of inhibition severity when template samples are compromised. Early detection of problematic samples allows analysts the opportunity to consider inhibitor mitigation strategies prior to genotype or DNA sequence analysis, thereby facilitating sample processing in high-throughput forensic operations.