Search

The impact of population bottlenecks is an important factor to consider when assessing species survival. Population declines can considerably limit the evolutionary potential of species and make them more susceptible to stochastic events. New Zealand has a well documented history of decline of endemic avifauna related to human colonization. Here, we investigate the genetic effects of a recent population decline in the endangered kea (Nestor notabilis). Kea have undergone a long-lasting persecution between the late 1800s to 1970s where an estimated 150,000 kea were culled under a governmental bounty scheme. Kea now number 1,000–5,000 individuals in the wild and it is likely that the recent population decline may have reduced the genetic diversity of the species. Comparison of contemporary (n = 410), historical (n = 15) and fossil samples (n = 4) showed a loss of mitochondrial diversity since the end of the last glaciation (Otiran Glacial) but no loss of overall genetic diversity associated with the cull. Microsatellite data indicated a recent bottleneck for only one population and a range-wide decline in Ne dating back some 300 – 6,000 years ago, a period predating European arrival in NZ. These results suggest that despite a recent human persecution, kea might have experienced a large population decline before stabilizing in numbers prior to human settlement of New Zealand in response to Holocene changes in habitat distribution. Our study therefore highlights the need to understand the respective effects of climate change and human activities on endangered species dynamics when proposing conservation guidelines.

Set-specific contingent attentional capture is a particularly strong form of capture that occurs when multiple attentional sets guide visual search (e.g., "search for green letters" and "search for orange letters"). In this type of capture, a potential target that matches one attentional set (e.g. a green stimulus) impairs the ability to identify a temporally proximal target that matches another attentional set (e.g. an orange stimulus). In the present study, we investigated whether set-specific capture stems from a bottleneck in working memory or from a depletion of limited resources that are distributed across multiple attentional sets. In each trial, participants searched a rapid serial visual presentation (RSVP) stream for up to three target letters (T1-T3) that could appear in any of three target colors (orange, green, or lavender). The most revealing findings came from trials in which T1 and T2 matched different attentional sets and were both identified. In these trials, T3 accuracy was lower when it did not match T1's set than when it did match, but only when participants failed to identify T2. These findings support a bottleneck model of set-specific capture in which a limited-capacity mechanism in working memory enhances only one attentional set at a time, rather than a resource model in which processing capacity is simultaneously distributed across multiple attentional sets.

There has been extensive debate over whether certain classes of genes are more likely than others to contain the causal variants responsible for phenotypic differences in complex traits between individuals. One hypothesis states that input/output genes positioned in signal transduction bottlenecks are more likely than other genes to contain causal natural variation. The IME1 gene resides at such a signaling bottleneck in the yeast sporulation pathway, suggesting that it may be more likely to contain causal variation than other genes in the sporulation pathway. Through crosses between natural isolates of yeast, we demonstrate that the specific causal nucleotides responsible for differences in sporulation efficiencies reside not only in IME1 but also in the genes that surround IME1 in the signaling pathway, including RME1, RSF1, RIM15, and RIM101. Our results support the hypothesis that genes at the critical decision making points in signaling cascades will be enriched for causal variants responsible for phenotypic differences., Author Summary Distinguishing the small number of genetic variants that impact phenotypes from the huge number of innocuous variants within an individual's genome is a difficult problem. Several hypotheses concerning the location of causal variants have been put forward based on the fact that genes are often organized into signaling cascades where the activation of a gene at the top of a pathway in turn activates large numbers of downstream genes. One hypothesis states that causal variations are more likely to reside in the genes at the top of these pathways because their effects are amplified by the signaling cascade. Here we provide support for this hypothesis by showing that causal genetic variants in yeast sporulation cluster around a gene at the top of the sporulation signaling cascade. Our result suggests a way to focus the search for causal genetic variants, including those that cause disease, on a smaller number of genes that are more likely to harbor important variations.

The HIV-1 epidemic in sub-Saharan Africa is driven largely by heterosexual transmission of non-subtype B viruses, of which subtypes C and A are predominant. Previous studies of subtype B and subtype C transmission pairs have suggested that a single variant from the chronically infected partner can establish infection in their newly infected partner. However, in subtype A infected individuals from a sex worker cohort and subtype B individuals from STD clinics, infection was frequently established by multiple variants. This study examined over 1750 single-genome amplified viral sequences derived from epidemiologically linked subtype C and subtype A transmission pairs very early after infection. In 90% (18/20) of the pairs, HIV-1 infection is initiated by a single viral variant that is derived from the quasispecies of the transmitting partner. In addition, the virus initiating infection in individuals who were infected by someone other than their spouse was characterized to determine if genital infections mitigated the severe genetic bottleneck observed in a majority of epidemiologically linked heterosexual HIV-1 transmission events. In nearly 50% (3/7) of individuals infected by someone other than their spouse, multiple genetic variants from a single individual established infection. A statistically significant association was observed between infection by multiple genetic variants and an inflammatory genital infection in the newly infected individual. Thus, in the vast majority of HIV-1 transmission events in cohabiting heterosexual couples, a single genetic variant establishes infection. Nevertheless, this severe genetic bottleneck can be mitigated by the presence of inflammatory genital infections in the at risk partner, suggesting that this restriction on genetic diversity is imposed in large part by the mucosal barrier., Author Summary Previous studies of HIV transmission have yielded conflicting results regarding the genetic heterogeneity of the virus establishing infection in the newly infected individual. In this study of populations from Zambia and Rwanda that are infected by two distinct viral genetic subtypes, we compared viral sequences that encode the entry-mediating envelope glycoproteins from newly infected individuals (recipients) and their spouses (donors) very early after infection, as well as newly infected individuals infected by someone other than their spouse. In spite of the genetically diverse virus population in the donor, approximately 90% of newly infected individuals were infected by a single viral variant, while the rest were infected by multiple viral variants. The homogeneity of the virus population in the newly infected recipient, as well as the presence, in some cases, of identical virus variants in the donor, allowed us to precisely identify the transmitted variant. We were able to examine the clinical history of each newly infected individual and observed that all individuals infected by multiple variants also showed evidence of inflammatory genital infections. Our results suggest that the genital mucosa provides a natural barrier to infection by multiple genetic variants of HIV-1, but that this barrier can be lowered by inflammatory genital infections.

Mathematical models (MMs) have been used to study the kinetics of influenza A virus infections under antiviral therapy, and to characterize the efficacy of antivirals such as neuraminidase inhibitors (NAIs). NAIs prevent viral neuraminidase from cleaving sialic acid receptors that bind virus progeny to the surface of infected cells, thereby inhibiting their release, suppressing infection spread. When used to study treatment with NAIs, MMs represent viral release implicitly as part of viral replication. Consequently, NAIs in such MMs do not act specifically and exclusively on virus release. We compared a MM with an explicit representation of viral release (i.e., distinct from virus production) to a simple MM without explicit release, and investigated whether parameter estimation and the estimation of NAI efficacy were affected by the use of a simple MM. Since the release rate of influenza A virus is not well-known, a broad range of release rates were considered. If the virus release rate is greater than ∼0.1 h-1, the simple MM provides accurate estimates of infection parameters, but underestimates NAI efficacy, which could lead to underdosing and the emergence of NAI resistance. In contrast, when release is slower than ∼0.1 h-1, the simple MM accurately estimates NAI efficacy, but it can significantly overestimate the infectious lifespan (i.e., the time a cell remains infectious and producing free virus), and it will significantly underestimate the total virus yield and thus the likelihood of resistance emergence. We discuss the properties of, and a possible lower bound for, the influenza A virus release rate.

RNA viruses typically occur in genetically diverse populations due to their error-prone genome replication. Genetic diversity is thought to be important in allowing RNA viruses to explore sequence space, facilitating adaptation to changing environments and hosts. Some arboviruses that infect both a mosquito vector and a mammalian host are known to experience population bottlenecks in their vectors, which may constrain their genetic diversity and could potentially lead to extinction events via Muller's ratchet. To examine this potential challenge of bottlenecks for arbovirus perpetuation, we studied Venezuelan equine encephalitis virus (VEEV) enzootic subtype IE and its natural vector Culex (Melanoconion) taeniopus, as an example of a virus-vector interaction with a long evolutionary history. Using a mixture of marked VEEV clones to infect C. taeniopus and real-time RT-PCR to track these clones during mosquito infection and dissemination, we observed severe bottleneck events that resulted in a significant drop in the number of clones present. At higher initial doses, the midgut was readily infected and there was a severe bottleneck at the midgut escape. Following a lower initial dose, the major bottleneck occurred at initial midgut infection. A second, less severe bottleneck was identified at the salivary gland infection stage following intrathoracic inoculation. Our results suggest that VEEV consistently encounters bottlenecks during infection, dissemination and transmission by its natural enzootic vector. The potential impacts of these bottlenecks on viral fitness and transmission, and the viral mechanisms that prevent genetic drift leading to extinction, deserve further study., Author Summary The ability of arboviruses to perpetuate in nature given that they must infect two disparate hosts (the mosquito vector and the vertebrate host) remains a mystery. We studied how viral genetic diversity is impacted by the dual host transmission cycle. Our studies of an enzootic cycle using Venezuelan equine encephalitis virus (VEEV) and its natural mosquito, Culex taeniopus, revealed the stages of infection that result in a viral population bottleneck. Using a set of marked VEEV clones and repeated sampling at various time points following C. taeniopus infection, we determined the number of clones in various mosquito tissues culminating in transmission. Bottlenecks were identified but the stage of occurrence was dependent on the dose that initiated infection. Understanding the points at which mosquito-borne viruses are constrained will shed light on the ways in which virus diversity varies, leading to selection of mutants that may result in host range changes or alterations in virulence.

After endocytic uptake, influenza viruses transit early endosomal compartments and eventually reach late endosomes. There, the viral glycoprotein hemagglutinin (HA) triggers fusion between endosomal and viral membrane, a critical step that leads to release of the viral segmented genome destined to reach the cell nucleus. Endosomal maturation is a complex process involving acidification of the endosomal lumen as well as endosome motility along microtubules. While the pH drop is clearly critical for the conformational change and membrane fusion activity of HA, the effect of intracellular transport dynamics on the progress of infection remains largely unclear. In this study, we developed a comprehensive mathematical model accounting for the first steps of influenza virus infection. We calibrated our model with experimental data and challenged its predictions using recombinant viruses with altered pH sensitivity of HA. We identified the time point of virus-endosome fusion and thereby the diffusion distance of the released viral genome to the nucleus as a critical bottleneck for efficient virus infection. Further, we concluded and supported experimentally that the viral RNA is subjected to cytosolic degradation strongly limiting the probability of a successful genome import into the nucleus., Author Summary Influenza A virus carries its segmented genome inside a lipid envelope. Since genome replication occurs inside the nucleus, the main goal of virus infection is to deliver all genome segments through the cytoplasm into the nucleus. After endocytic uptake, influenza viruses transit early endosomal compartments and eventually reach late endosomes. Within a complex maturation process, the endosomal lumen acidifies while the vesicles are transported trough the cytosol. If and how these early processes affect virus infection remained mostly speculative. To reach a better understanding and to quantify the physical interplay between membrane fusion, genome diffusion and infection, we developed a mathematical model that comprises all initial steps of virus infection until genome delivery. We calibrated our model using experimental data and challenged its predictions using recombinant viruses to introduce perturbations. Our results provide a theoretical framework to understand the importance of the endosomal virus passage before membrane fusion and genome release. We further unraveled RNA degradation as a previously unknown limiting factor for virus infection. Our work will help to make predictions and evaluate newly occurring virus strains, regarding their infection efficiency in a given host cell, by simply considering their pH sensitivity.

The 3′ untranslated genome region (UTR) of arthropod-borne viruses is characterized by enriched direct repeats (DRs) and stem-loop structures. Despite many years of theoretical and experimental study, on-going positive selection on the 3′UTR had never been observed in ‘real-time,’ and the role of the arbovirus 3′UTR remains poorly understood. We observed a lineage-specific 3′UTR sequence pattern in all available Asian lineage of the mosquito-borne alphavirus, chikungunya virus (CHIKV) (1958–2009), including complicated mutation and duplication patterns of the long DRs. Given that a longer genome is usually associated with less efficient replication, we hypothesized that the fixation of these genetic changes in the Asian lineage 3′UTR was due to their beneficial effects on adaptation to vectors or hosts. Using reverse genetic methods, we examined the functional importance of each direct repeat. Our results suggest that adaptation to mosquitoes, rather than to mammalian hosts, is a major evolutionary force on the CHIKV 3′UTR. Surprisingly, the Asian 3′UTR appeared to be inferior to its predicted ancestral sequence for replication in both mammals and mosquitoes, suggesting that its fixation in Asia was not a result of directional selection. Rather, it may have resulted from a population bottleneck during its introduction from Africa to Asia. We propose that this introduction of a 3′UTR with deletions led to genetic drift and compensatory mutations associated with the loss of structural/functional constraints, followed by two independent beneficial duplications and fixation due to positive selection. Our results provide further evidence that the limited epidemic potential of the Asian CHIKV strains resulted from founder effects that reduced its fitness for efficient transmission by mosquitoes there., Author Summary The 3′ untranslated genome region (UTR) of arbovirus (arthropod-borne virus) RNA genomes is characterized by enriched direct repeats (DRs) and stem-loop structures, which are thought to play a specific role in maintaining efficient transmission in multiple hosts. However, this hypothesis is vague and has little experimental support. Based on our observation of a distinct, lineage-specific sequence structure in the chikungunya virus (CHIKV) Asian lineage 3′UTR, we tried to understand the underlying driving forces on the arboviral 3′UTR evolution. Specifically, we sought to determine whether the dramatic genetic changes in the Asian lineage 3′UTR were the result of adaptation into the new host or vector populations. Using reverse genetic methods, we obtained results that suggested that the DRs have a significant effect on the viral fitness level in mosquitoes but not in mammals. Interestingly, instead of a directional selection from its ancestral state, our results suggest that the evolution of the CHIKV Asian lineage 3′UTR involved a population bottleneck with a deleterious deletion of DRs, followed by accumulation of mutations due to the loss of structural/functional constraints. Later, the duplication of two regions that are beneficial in mosquitoes led to the fixation of this 3′UTR sequence in Asian lineage, possibly facilitating the 1958 Thailand outbreak.

ELISA is the main approach for the sensitive quantification of protein biomarkers in body fluids and is currently employed in clinical laboratories for the measurement of clinical markers. As such, it also constitutes the main methodological approach for biomarker validation and further qualification. For the latter, specific assay performance requirements have to be met, as described in respective guidelines of regulatory agencies. Even though many clinical ELISA assays in serum are regularly used, ELISA clinical applications in urine are significantly less. The scope of our study was to evaluate ELISA assay analytical performance in urine for a series of potential biomarkers for bladder cancer, as a first step towards their large scale clinical validation. Seven biomarkers (Secreted protein acidic and rich in cysteine, Survivin, Slit homolog 2 protein, NRC-Interacting Factor 1, Histone 2B, Proteinase-3 and Profilin-1) previously described in the literature as having differential expression in bladder cancer were included in the study. A total of 11 commercially available ELISA tests for these markers were tested by standard curve analysis, assay reproducibility, linearity and spiking experiments. The results show disappointing performance with coefficients of variation>20% for the vast majority of the tests performed. Only 3 assays (for Secreted protein acidic and rich in cysteine, Survivin and Slit homolog 2 protein) passed the accuracy thresholds and were found suitable for further application in marker quantification. These results collectively reflect the difficulties in developing urine-based ELISA assays of sufficient analytical performance for clinical application, presumably attributed to the urine matrix itself and/or presence of markers in various isoforms.

In mammals, observations of rapid shifts in mitochondrial DNA (mtDNA) variants between generations have led to the creation of the bottleneck theory for the transmission of mtDNA. The bottleneck could be attributed to a marked decline of mtDNA content in germ cells giving rise to the next generation, to a small effective number of mtDNA segregation units resulting from homoplasmic nucleoids rather than the single mtDNA molecule serving as the units of segregation, or to the selective transmission of a subgroup of the mtDNA population to the progeny. We have previously determined mtDNA copy number in single germ cells and shown that the bottleneck occurs without the reduction in germline mtDNA content. Recently one study suggested that the bottleneck is driven by a remarkable decline of mtDNA copies in early primordial germ cells (PGCs), while another study reported that the mtDNA genetic bottleneck results from replication of a subpopulation of the mtDNA genome during postnatal oocyte maturation and not during embryonic oogenesis, despite a detected a reduction in mtDNA content in early PGCs. To clarify these contradictory results, we examined the mtDNA copy number in PGCs isolated from transgenic mice expressing fluorescent proteins specifically in PGCs as in the aforementioned two other studies. We provide clear evidence to confirm that no remarkable reduction in mtDNA content occurs in PGCs and reinforce that the bottleneck is generated without reduction of mtDNA content in germ cells., Author Summary Mutations of mtDNA are responsible for many types of mitochondrial diseases in humans, including myopathy and neurological disorders. Females carrying a mixture of mutant and wild-type mtDNA variants transmit a variable amount of mutant mtDNA to each offspring. The proportion of mutated mtDNA inherited from the mother determines the onset and severity of diseases. Studies have suggested that the mtDNA genome is transmitted through a bottleneck, but the underlying mechanism remains controversial. By detecting mtDNA copy number in single cells, we previously showed that the bottleneck occurs without reduction of mtDNA content in germline cells. However, recently a study reported a marked decline of mtDNA copies in embryonic germ cells and attributed this reduction to the creation of the bottleneck. Yet another study concluded that the bottleneck occurs during postnatal oocyte maturation and not during embryonic oogenesis. To resolve these controversies, we examined mtDNA copies in embryonic germ cells identified using the same methodology as in the other two studies. We show solid evidence to confirm our previous findings. This confirmation is important because the understanding of mtDNA content in female germ cells will facilitate the development of therapeutic strategies preventing the transmission of mitochondrial diseases from mother to offspring.

The western flower thrips, Frankliniella occidentalis (Pergande), is an invasive species and the most economically important pest within the insect order Thysanoptera. F. occidentalis, which is endemic to North America, was initially detected in Kunming in southwestern China in 2000 and since then it has rapidly invaded several other localities in China where it has greatly damaged greenhouse vegetables and ornamental crops. Controlling this invasive pest in China requires an understanding of its genetic makeup and migration patterns. Using the mitochondrial COI gene and 10 microsatellites, eight of which were newly isolated and are highly polymorphic, we investigated the genetic structure and the routes of range expansion of 14 F. occidentalis populations in China. Both the mitochondrial and microsatellite data revealed that the genetic diversity of F. occidentalis of the Chinese populations is lower than that in its native range. Two previously reported cryptic species (or ecotypes) were found in the study. The divergence in the mitochondrial COI of two Chinese cryptic species (or ecotypes) was about 3.3% but they cannot be distinguished by nuclear markers. Hybridization might produce such substantial mitochondrial-nuclear discordance. Furthermore, we found low genetic differentiation (global F(ST) = 0.043, P

Why is the human brain fundamentally limited when attempting to execute two tasks at the same time or in close succession? Two classical paradigms, psychological refractory period (PRP) and task switching, have independently approached this issue, making significant advances in our understanding of the architecture of cognition. Yet, there is an apparent contradiction between the conclusions derived from these two paradigms. The PRP paradigm, on the one hand, suggests that the simultaneous execution of two tasks is limited solely by a passive structural bottleneck in which the tasks are executed on a first-come, first-served basis. The task-switching paradigm, on the other hand, argues that switching back and forth between task configurations must be actively controlled by a central executive system (the system controlling voluntary, planned, and flexible action). Here we have explicitly designed an experiment mixing the essential ingredients of both paradigms: task uncertainty and task simultaneity. In addition to a central bottleneck, we obtain evidence for active processes of task setting (planning of the appropriate sequence of actions) and task disengaging (suppression of the plan set for the first task in order to proceed with the next one). Our results clarify the chronometric relations between these central components of dual-task processing, and in particular whether they operate serially or in parallel. On this basis, we propose a hierarchical model of cognitive architecture that provides a synthesis of task-switching and PRP paradigms., The authors investigate the cognitive processes that constrain our ability to perform multiple mental tasks by clarifying the relationship between task chronometry and the dynamics of task selection.

Carbon steels are widely used in the oil and gas industry from downhole tubing to transport trunk lines. Microbes form biofilms, some of which cause the so-called microbiologically influenced corrosion (MIC) of carbon steels. MIC by sulfate reducing bacteria (SRB) is often a leading cause in MIC failures. Electrogenic SRB sessile cells harvest extracellular electrons from elemental iron oxidation for energy production in their metabolism. A previous study suggested that electron mediators riboflavin and flavin adenine dinucleotide (FAD) both accelerated the MIC of 304 stainless steel by the Desulfovibrio vulgaris biofilm that is a corrosive SRB biofilm. Compared with stainless steels, carbon steels are usually far more prone to SRB attacks because SRB biofilms form much denser biofilms on carbon steel surfaces with a sessile cell density that is two orders of magnitude higher. In this work, C1018 carbon steel coupons were used in tests of MIC by D. vulgaris with and without an electron mediator. Experimental weight loss and pit depth data conclusively confirmed that both riboflavin and FAD were able to accelerate D. vulgaris attack against the carbon steel considerably. It has important implications in MIC failure analysis and MIC mitigation in the oil and gas industry.

As wild populations decline, ex situ propagation provides a potential bank of genetic diversity and a hedge against extinction. These programs are unlikely to succeed if captive populations do not recover from the severe bottleneck imposed when they are founded with a limited number of individuals from remnant populations. In small captive populations allelic richness may be lost due to genetic drift, leading to a decline in fitness. Wild populations of the Hawaiian tree snail Achatinella lila, a hermaphroditic snail with a long life history, have declined precipitously due to introduced predators and other human impacts. A captive population initially thrived after its founding with seven snails, exceeding 600 captive individuals in 2009, but drastically declined in the last five years. Measures of fitness were examined from 2,018 captive snails that died between 1998 and 2012, and compared with genotypic data for six microsatellite loci from a subset of these deceased snails (N = 335), as well as live captive snails (N = 198) and wild snails (N = 92). Surprisingly, the inbreeding coefficient (Fis) declined over time in the captive population, and is now approaching values observed in the 2013 wild population, despite a significant decrease in allelic richness. However, adult annual survival and fecundity significantly declined in the second generation. These measures of fitness were positively correlated with heterozygosity. Snails with higher measures of heterozygosity had more offspring, and third generation offspring with higher measures of heterozygosity were more likely to reach maturity. These results highlight the importance of maintaining genetic diversity in captive populations, particularly those initiated with a small number of individuals from wild remnant populations. Genetic rescue may allow for an increase in genetic diversity in the captive population, as measures of heterozygosity and rarified allelic richness were higher in wild tree snails.

Parsing a mental operation into components, characterizing the parallel or serial nature of this flow, and understanding what each process ultimately contributes to response time are fundamental questions in cognitive neuroscience. Here we show how a simple theoretical model leads to an extended set of predictions concerning the distribution of response time and its alteration by simultaneous performance of another task. The model provides a synthesis of psychological refractory period and random-walk models of response time. It merely assumes that a task consists of three consecutive stages-perception, decision based on noisy integration of evidence, and response-and that the perceptual and motor stages can operate simultaneously with stages of another task, while the central decision process constitutes a bottleneck. We designed a number-comparison task that provided a thorough test of the model by allowing independent variations in number notation, numerical distance, response complexity, and temporal asynchrony relative to an interfering probe task of tone discrimination. The results revealed a parsing of the comparison task in which each variable affects only one stage. Numerical distance affects the integration process, which is the only step that cannot proceed in parallel and has a major contribution to response time variability. The other stages, mapping the numeral to an internal quantity and executing the motor response, can be carried out in parallel with another task. Changing the duration of these processes has no significant effect on the variance.

Background: Although there are different strains of HIV-1 in a chronically infected individual, only one or limited virus strains are successfully transmitted to a new individual. The reason for this "transmission bottleneck" is as yet unknown. Methodology/Principal Findings: A human cervical explant model was used to measure HIV-1 transmission efficiency of viral strains from chronic infections, and transmitter/founder variants. We also evaluated the genetic characteristics of HIV-1 variants in the inoculums compared to those transmitted across the cervical mucosa. Eight different HIV-1 isolates were used in this study, six chronic isolates and two transmitter/founder viruses. The transmission efficiency of the chronic and transmitter/founder virus isolates and the viral diversity of chronic isolates before and after viral transmission were assessed. The results indicate that transmitter/founder viruses did not display higher transmission efficiency than chronic HIV-1 isolates. Furthermore, no evidence for a difference in diversity was found between the inoculums and transmitted virus strains. Phylogenetic analysis indicated that the sequences of variants in the inoculums and those present in transmitted virus intermingled irrespective of co-receptor usage. In addition, the inoculum and transmitted variants had a similar pairwise distance distribution. Conclusion: There was no selection of a single or limited number of viral variants during HIV-1 transmission across the cervical mucosa in the organ culture model, indicating that the cervical mucosa alone may not produce the transmission bottleneck of HIV-1 infection observed in vivo. © 2012 Shen et al.

Background Breast milk transmission of HIV-1 remains a major route of pediatric infection. Defining the characteristics of viral variants to which breastfeeding infants are exposed is important for understanding the genetic bottleneck that occurs in the majority of mother-to-child transmissions. The blood-milk epithelial barrier markedly restricts the quantity of HIV-1 in breast milk, even in the absence of antiretroviral drugs. The basis of this restriction and the genetic relationship between breast milk and blood variants are not well established. Methodology/Principal Findings We compared 356 HIV-1 subtype C gp160 envelope (env) gene sequences from the plasma and breast milk of 13 breastfeeding women. A trend towards lower viral population diversity and divergence in breast milk was observed, potentially indicative of clonal expansion within the breast. No differences in potential N-linked glycosylation site numbers or in gp160 variable loop amino acid lengths were identified. Genetic compartmentalization was evident in only one out of six subjects in whom contemporaneously obtained samples were studied. However, in samples that were collected 10 or more days apart, six of seven subjects were classified as having compartmentalized viral populations, highlighting the necessity of contemporaneous sampling for genetic compartmentalization studies. We found evidence of CXCR4 co-receptor using viruses in breast milk and blood in nine out of the thirteen subjects, but no evidence of preferential localization of these variants in either tissue. Conclusions/Significance Despite marked restriction of HIV-1 quantities in milk, our data indicate intermixing of virus between blood and breast milk. Thus, we found no evidence that a restriction in viral genotype diversity in breast milk accounts for the genetic bottleneck observed following transmission. In addition, our results highlight the rapidity of HIV-1 env evolution and the importance of sample timing in analyses of gene flow.

Forward genetics screens with N-ethyl-N-nitrosourea (ENU) provide a powerful way to illuminate gene function and generate mouse models of human disease; however, the identification of causative mutations remains a limiting step. Current strategies depend on conventional mapping, so the propagation of affected mice requires non-lethal screens; accurate tracking of phenotypes through pedigrees is complex and uncertain; out-crossing can introduce unexpected modifiers; and Sanger sequencing of candidate genes is inefficient. Here we show how these problems can be efficiently overcome using whole-genome sequencing (WGS) to detect the ENU mutations and then identify regions that are identical by descent (IBD) in multiple affected mice. In this strategy, we use a modification of the Lander-Green algorithm to isolate causative recessive and dominant mutations, even at low coverage, on a pure strain background. Analysis of the IBD regions also allows us to calculate the ENU mutation rate (1.54 mutations per Mb) and to model future strategies for genetic screens in mice. The introduction of this approach will accelerate the discovery of causal variants, permit broader and more informative lethal screens to be used, reduce animal costs, and herald a new era for ENU mutagenesis., Author Summary Damaging mutations in single genes are an important source of information about the causes of disease, including more complex genetic disease; but these single gene disorders are typically rare in humans. An important strategy for identifying new disease mechanisms is to introduce multiple random mutations in mice and test the mice for biological differences; these mice act as models of human disease. However, discovering the disease-causing mutation is time-consuming and complex, requiring further generations of breeding. In this study we demonstrate a method that overcomes these problems by sequencing the entire genomes of multiple mice that have inherited a disease-causing mutation from a common ancestor. We use an algorithm that uses knowledge of all the mutations carried by the sequenced mice to identify the regions of the genome and mutations that are common to all the mice. Using this method we can rapidly link biological traits to genetic mutations. In contrast to current approaches, our strategy does not require large amounts of breeding, and it permits more accurate measurement of a wider range of traits; consequently its introduction will significantly reduce the number of mice required in the future, increase the number of traits that can be detected, and accelerate the discovery of new pathways and gene functions relevant to human diseases.

The innate immune system plays important roles in a number of disparate processes. Foremost, innate immunity is a first responder to invasion by pathogens and triggers early defensive responses and recruits the adaptive immune system. The innate immune system also responds to endogenous damage signals that arise from tissue injury. Recently it has been found that innate immunity plays an important role in neuroprotection against ischemic stroke through the activation of the primary innate immune receptors, Toll-like receptors (TLRs). Using several large-scale transcriptomic data sets from mouse and mouse macrophage studies we identified targets predicted to be important in controlling innate immune processes initiated by TLR activation. Targets were identified as genes with high betweenness centrality, so-called bottlenecks, in networks inferred from statistical associations between gene expression patterns. A small set of putative bottlenecks were identified in each of the data sets investigated including interferon-stimulated genes (Ifit1, Ifi47, Tgtp and Oasl2) as well as genes uncharacterized in immune responses (Axud1 and Ppp1r15a). We further validated one of these targets, Ifit1, in mouse macrophages by showing that silencing it suppresses induction of predicted downstream genes by lipopolysaccharide (LPS)-mediated TLR4 activation through an unknown direct or indirect mechanism. Our study demonstrates the utility of network analysis for identification of interesting targets related to innate immune function, and highlights that Ifit1 can exert a positive regulatory effect on downstream genes.

The population of the globally endangered hog deer (Axis porcinus) has declined severely across its geographic range. Intensive monitoring of its demographic and genetic status is necessary. Northeast India is a stronghold of the species; however, in recent years the population has been getting fragmented, and it is vulnerable to extinction. We examined the demographic and genetic structure of a small hog deer population in the floating meadows of Keibul Lamjao National Park (KLNP), located on the western fringe of the Indo-Burma biodiversity hotspot for conservation planning. We used a double-observer distance sampling method to derive the hog deer abundance and population structure. We also derived the genetic diversity of the population through microsatellite screening and bottleneck detection. Our study revealed that the abundance of the deer in the park was 1.82–4.32 individuals/km2. The adult male to female ratio and fawn to doe ratio were 36.2 ± 1.9 males/100 females and 16.5 ± 0.4 fawns/100 females, respectively. The mean number of alleles at 23 loci was 2.70 ± 0.18, the observed heterozygosity (HO) ranged from 0.26 to 0.63 (mean 0.42 ± 0.02), the expected heterozygosity (HE) ranged from 0.23 to 0.73 (χ = 0.51 ± 0.03), and the polymorphic information content (PIC) ranged from 0.2 to 0.67 (χ = 0.43 ± 0.03). The observed allelic distribution reveals that the population has not encountered any genetic bottleneck in the recent past. Although the population is declining, it still retains some rare alleles, and the genetic diversity is 50%. This diversity will probably not affect the short-term population growth but may affect the evolutionary potential by limiting the selection flexibility. Conservation measures coupled with a scientifically sound management regime may help the persistence of the population in the region at a time when the population still retains rare alleles and maintains reproductive fitness.

CRISPR-Staphylococcus aureus Cas9 (CRISPR-SaCas9) has been harnessed as an effective in vivo genome-editing tool to manipulate genomes. However, off-target effects remain a major bottleneck that precludes safe and reliable applications in genome editing. Here, we characterize the off-target effects of wild-type (WT) SaCas9 at single-nucleotide (single-nt) resolution and describe a directional screening system to identify novel SaCas9 variants with desired properties in human cells. Using this system, we identified enhanced-fidelity SaCas9 (efSaCas9) (variant Mut268 harboring the single mutation of N260D), which could effectively distinguish and reject single base-pair mismatches. We demonstrate dramatically reduced off-target effects (approximately 2- to 93-fold improvements) of Mut268 compared to WT using targeted deep-sequencing analyses. To understand the structural origin of the fidelity enhancement, we find that N260, located in the REC3 domain, orchestrates an extensive network of contacts between REC3 and the guide RNA-DNA heteroduplex. efSaCas9 can be broadly used in genome-editing applications that require high fidelity. Furthermore, this study provides a general strategy to rapidly evolve other desired CRISPR-Cas9 traits besides enhanced fidelity, to expand the utility of the CRISPR toolkit., CRISPR-SaCas9 technology offers promising solutions to treat genetic disorders by correcting pathogenic mutations, but off-target effects remain a major issue. This study characterizes the off-target activity of wild-type SaCas9 and uses a human cell-based screening system to identify high-fidelity SaCas9 variants, showing that they can dramatically reduce off-target effects.

We approach the C. elegans connectome as an information processing network that receives input from about 90 sensory neurons, processes that information through a highly recurrent network of about 80 interneurons, and it produces a coordinated output from about 120 motor neurons that control the nematode’s muscles. We focus on the feedforward flow of information from sensory neurons to motor neurons, and apply a recently developed network analysis framework referred to as the “hourglass effect”. The analysis reveals that this feedforward flow traverses a small core (“hourglass waist”) that consists of 10-15 interneurons. These are mostly the same interneurons that were previously shown (using a different analytical approach) to constitute the “rich-club” of the C. elegans connectome. This result is robust to the methodology that separates the feedforward from the feedback flow of information. The set of core interneurons remains mostly the same when we consider only chemical synapses or the combination of chemical synapses and gap junctions. The hourglass organization of the connectome suggests that C. elegans has some similarities with encoder-decoder artificial neural networks in which the input is first compressed and integrated in a low-dimensional latent space that encodes the given data in a more efficient manner, followed by a decoding network through which intermediate-level sub-functions are combined in different ways to compute the correlated outputs of the network. The core neurons at the hourglass waist represent the information bottleneck of the system, balancing the representation accuracy and compactness (complexity) of the given sensory information., Author summary The C. elegans nematode is the only species for which the complete wiring diagram (“connectome”) of its neural system has been mapped. The connectome provides architectural constraints that limit the scope of possible functions of a neural system. In this work, we identify one such architectural constraint: the C. elegans connectome includes a small set (10-15) of neurons that compress and integrate the information provided by the much larger set of sensory neurons. These intermediate-level neurons encode few sub-functions that are combined and re-used in different ways to activate the circuits of motor neurons, which drive all higher-level complex functions of the organism such as feeding or locomotion. We refer to this encoding-decoding structure as “hourglass architecture” and identify the core neurons at the “waist” of the hourglass. We also discuss the similarities between this property of the C. elegans connectome and artificial neural networks. The hourglass architecture opens a new way to think about, and experiment with, intermediate-level neurons between input and output neural circuits.

Hopea hainanensisMerrill & Chun (Dipterocarpaceae) is an endangered tree species restricted to Hainan Island, China and a small part of Northern Vietnam. On Hainan Island, it is an important indicator species for tropical forests. However, because of its highly valued timber,H.hainanensishas suffered from overexploitation, leading to a sharp population decline. To facilitate the conservation of this species, genetic diversity and population structure were assessed using 12 SSR markers for 10 populations sampled across Hainan Island. Compared to non-threatenedHopeaspecies,H.hainanensisexhibited reduced overall genetic diversity and increased population differentiation (AMOVA:FST= 0.23). Bayesian model-based clustering and principal coordinate analysis consistently assignedH.hainanensisindividuals into three genetic groups, which were found to be widespread and overlapping geographically. A Mantel test found no correlation between genetic and geographical distances (r= 0.040,p= 0.418). The observed genetic structure suggests that long-distance gene flow occurred amongH.hainanensispopulations prior to habitat fragmentation. A recent population bottleneck was revealed, which may cause rapid loss of genetic diversity and increased differentiation across populations. Based on these findings, appropriate strategies for the long-term conservation of the endangered speciesH.hainanensisare proposed.

The beet armyworm, Spodoptera exigua, is a significant agricultural pest of numerous crops and has caused serious economic losses in China. To effectively control this pest, we analyzed its genetic variation, population genetic structure and demographic history. We used mitochondrial DNA (mtDNA) fragments of the cytochrome oxidase subunit I (COI) and eight nuclear microsatellite loci to investigate genetic diversity and population genetic structure of S. exigua populations at 14 sampling sites in western China. Both mtDNA and microsatellite data indicated low levels of genetic diversity among all populations. A moderate genetic differentiation among some S. exigua populations was detected. Neighbor-joining dendrograms, STRUCTURE, and principal coordinate analysis (PCoA) revealed two genetically distinct groups: the KEL group and the remaining population group. Isolation by distance (IBD) results showed a weak significant correlation between geographic distance and genetic differentiation. Haplotype networks, neutrality testing, and mismatch distribution analysis indicated that the beet armyworm experienced a recent rapid expansion without a recent genetic bottleneck in western China. Thus, the results of this population genetic study can help with the development of strategies for managing this highly migratory pest.

Transmission between hosts is a critical part of the viral lifecycle. Recent studies of viral transmission have used genome sequence data to evaluate the number of particles transmitted between hosts, and the role of selection as it operates during the transmission process. However, the interpretation of sequence data describing transmission events is a challenging task. We here present a novel and comprehensive framework for using short-read sequence data to understand viral transmission events, designed for influenza virus, but adaptable to other viral species. Our approach solves multiple shortcomings of previous methods for this purpose; for example, we consider transmission as an event involving whole viruses, rather than sets of independent alleles. We demonstrate how selection during transmission and noisy sequence data may each affect naive inferences of the population bottleneck, accounting for these in our framework so as to achieve a correct inference. We identify circumstances in which selection for increased viral transmission may or may not be identified from data. Applying our method to experimental data in which transmission occurs in the presence of strong selection, we show that our framework grants a more quantitative insight into transmission events than previous approaches, inferring the bottleneck in a manner that accounts for selection, both for within-host virulence, and for inherent viral transmissibility. Our work provides new opportunities for studying transmission processes in influenza, and by extension, in other infectious diseases., Author summary In order to spread, pathogens must not only be able to grow within an infected host, but also transmit to found new infections. Population genetics can exploit genome sequence data to provide a great deal of insight into transmission processes. For example, the number of particles which found a new infection determine the extent to which genetic diversity is passed from host to host. The identification of genetic variants which increase the propensity of a pathogen to transmit from host to host is a valuable step in understanding how an infection might spread. Here we set out a new population genetic framework for understanding transmission events from genome sequence data collected before and after transmission. Our approach corrects for the shortcomings of existing methods for this purpose, setting out a new baseline for the statistical analysis of transmission events. We demonstrate the ability of our method to draw novel quantitative insights by application to data from simulated and real transmission events.

Like other crop species, barley, the fourth most important crop worldwide, suffers from the genetic bottleneck effect, where further improvements in performance through classical breeding methods become difficult. Therefore, indirect selection methods are of great interest. Here, genomic prediction (GP) based on 33,005 SNP markers and, alternatively, metabolic prediction (MP) based on 128 metabolites with sampling at two different time points in one year, were applied to predict multi-year agronomic traits in the nested association mapping (NAM) population HEB-25. We found prediction abilities of up to 0.93 for plant height with SNP markers and of up to 0.61 for flowering time with metabolites. Interestingly, prediction abilities in GP increased after reducing the number of incorporated SNP markers. The estimated effects of GP and MP were highly concordant, indicating MP as an interesting alternative to GP, being able to reflect a stable genotype-specific metabolite profile. In MP, sampling at an early developmental stage outperformed sampling at a later stage. The results confirm the value of GP for future breeding. With MP, an interesting alternative was also applied successfully. However, based on our results, usage of MP alone cannot be recommended in barley. Nevertheless, MP can assist in unravelling physiological pathways for the expression of agronomically important traits.

The genetic variability and population structure of three Bhutanese traditional horse breeds were assessed through genotyping of 74 horses (Boeta 25, Sharta 14 and Yuta 35) for 29 microsatellite DNA loci. Altogether, 282 alleles were detected across 29 polymorphic loci. The allelic diversity (NE) (Boeta 4.94; Sharta 4.65; Yuta 5.30) and gene diversities (HE) (Boeta 0.78; Sharta 0.77; Yuta 0.79) were high. None of the breeds deviated significantly from the Hardy-Weinberg equilibrium. There was no sign of significant population bottleneck for all the breeds. The inbreeding estimates (FIS) of the breeds were low (Boeta 0.023; Sharta 0.001; Yuta 0.021). Analysis of molecular variance showed 0.6% of the total genetic variation among breeds, 1.9% among individuals and 97.5% within individuals. The global FIT, FST, and FIS estimates for the population were 0.025, 0.006 and 0.019 respectively. The analysis of population structure failed to distinguish subpopulations in traditional horses and this was supported by a high genetic exchange among the breeds. Overall, the results of this study suggest a rich genetic diversity in the traditional horse despite a very low genetic differentiation among the breeds in Bhutan.

With fewer than 200 tigers (Panthera tigris tigris) left in Nepal, that are generally confined to five protected areas across the Terai Arc Landscape, genetic studies are needed to provide crucial information on diversity and connectivity for devising an effective country-wide tiger conservation strategy. As part of the Nepal Tiger Genome Project, we studied landscape change, genetic variation, population structure, and gene flow of tigers across the Terai Arc Landscape by conducting Nepal's first comprehensive and systematic scat-based, non-invasive genetic survey. Of the 770 scat samples collected opportunistically from five protected areas and six presumed corridors, 412 were tiger (57%). Out of ten microsatellite loci, we retain eight markers that were used in identifying 78 individual tigers. We used this data set to examine population structure, genetic variation, contemporary gene flow, and potential population bottlenecks of tigers in Nepal. We detected three genetic clusters consistent with three demographic sub-populations and found moderate levels of genetic variation (H-e = 0.61, A(R) = 3.51) and genetic differentiation (F-ST = 0.14) across the landscape. We detected 3-7 migrants, confirming the potential for dispersal-mediated gene flow across the landscape. We found evidence of a bottleneck signature likely caused by large-scale land-use change documented in the last two centuries in the Terai forest. Securing tiger habitat including functional forest corridors is essential to enhance gene flow across the landscape and ensure long-term tiger survival. This requires cooperation among multiple stakeholders and careful conservation planning to prevent detrimental effects of anthropogenic activities on tigers. USAID [FOG AID-367-G-11-00001]; Katheryn Fuller Fellowship (WWF US); Wildlife Conservation Society (WCS); Virginia Tech (VT) This work was supported by USAID FOG AID-367-G-11-00001. Netra Sharma and Bronwyn Llewellyn were involved with the project and helped in conceptualization. The Katheryn Fuller Fellowship (WWF US), Wildlife Conservation Society (WCS), and Virginia Tech (VT) supported Kanchan Thapa.

Diamondback terrapins (Malaclemys terrapin) were a popular food item in early twentieth century America, and were consumed in soup with sherry. Intense market demand for terrapin meat resulted in population declines, notably along the Atlantic seaboard. Efforts to supply terrapins to markets resulted in translocation events, as individuals were moved about to stock terrapin farms. However, in 1920 the market for turtle soup buckled with the enactment of the eighteenth amendment to the United States’ Constitution—which initiated the prohibition of alcoholic drinks—and many terrapin fisheries dumped their stocks into local waters. We used microsatellite data to show that patterns of genetic diversity along the terrapin’s coastal range are consistent with historical accounts of translocation and cultivation activities. We identified possible instances of human-mediated dispersal by estimating gene flow over historical and contemporary timescales, Bayesian model testing, and bottleneck tests. We recovered six genotypic clusters along the Gulf and Atlantic coasts with varying degrees of admixture, including increased contemporary gene flow from Texas to South Carolina, from North Carolina to Maryland, and from North Carolina to New York. In addition, Bayesian models incorporating translocation events outperformed stepping-stone models. Finally, we were unable to detect population bottlenecks, possibly due to translocation reintroducing genetic diversity into bottlenecked populations. Our data suggest that current patterns of genetic diversity in the terrapin were altered by the demand for turtle soup followed by the enactment of alcohol prohibition. In addition, our study shows that population genetic tools can elucidate metapopulation dynamics in taxa with complex genetic histories impacted by anthropogenic activities.

The study of the demographic history of species can help to understand the negative impact of recent population declines in organisms of conservation concern. Here, we use neutral molecular markers to explore the genetic consequences of the recent population decline and posterior recovery of the Eurasian eagle owl (Bubo bubo) in the Iberian Peninsula. During the last century, the species was the object of extermination programs, suffering direct persecution by hunters until the 70’s. Moreover, during the last decades the eagle owl was severely impacted by increased mortality due to electrocution and the decline of its main prey species, the European rabbit (Oryctolagus cuniculus). In recent times, the decrease of direct persecution and the implementation of some conservation schemes have allowed the species’ demographic recovery. Yet, it remains unknown to which extent the past population decline and the later expansion have influenced the current species’ pattern of genetic diversity. We used eight microsatellite markers to genotype 235 eagle owls from ten Spanish subpopulations and analyse the presence of genetic signatures attributable to the recent population fluctuations experienced by the species. We found moderate levels of differentiation among the studied subpopulations and Bayesian analyses revealed the existence of three genetic clusters that grouped subpopulations from central, south-western and southeastern Spain. The observed genetic structure could have resulted from recent humaninduced population fragmentation, a patchy distribution of prey populations and/or the philopatric behaviour and habitat selection of the species. We detected an old population bottleneck, which occurred approximately 10,000 years ago, and significant signatures of recent demographic expansions. However, we did not find genetic signatures for a recent bottleneck, which may indicate that population declines were not severe enough to leave detectable signals on the species genetic makeup or that such signals have been eroded by the rapid demographic recovery experienced by the species in recent years

Cornus florida (flowering dogwood) and C. nuttallii (Pacific dogwood) are North American native tree species that belong to the big-bracted group of dogwoods. Cornus species are highly valued for their ornamental characteristics, and have fruits that contain high fat content for animals. Also, they are an important understory tree in natural forests. Dogwood anthracnose, caused by Discula destructiva, was observed in the late 1970s on the east and west coasts of the United States and by 1991 had quickly spread throughout most of the native ranges of C. florida and C. nuttalli. We investigated the genetic diversity and population structure of 93 D. destructiva isolates using 47 microsatellite loci developed from the sequenced genome of the type strain of D. destructiva. Clone-corrected data indicated low genetic diversity and the presence of four genetic clusters that corresponded to two major geographic areas, the eastern United States and the Pacific Northwest, and to the two collection time periods when the isolates were collected (pre- and post-1993). Linkage disequilibrium was present in five out of six subpopulations, suggesting that the fungus only reproduced asexually. Evidence of population bottlenecks was indicated across four identified genetic clusters, and was probably the result of the limited number of founding individuals on both coasts. These results support the hypothesis that D. destructiva is an exotic pathogen with independent introductions on the east and west coasts of North America. We also tested the cross-amplification of these microsatellite primers to other Discula species. Genomic DNA from 17 isolates of four other Discula species and two isolates of Juglanconis species (formerly Melanconis species) were amplified by 17 of 47 primer pairs. These primers may be useful for investigating the genetic diversity and population structure of these Discula species.

Genetic spaces are often described in terms of fitness landscapes or genotype-to-phenotype maps, where each genetic sequence is associated with phenotypic properties and linked to other genotypes that are a single mutational step away. The positions close to a genotype make up its “mutational landscape” and, in aggregate, determine the short-term evolutionary potential of a population. Populations with wider ranges of phenotypes in their mutational neighborhood are known to be more evolvable. Likewise, those with fewer phenotypic changes available in their local neighborhoods are more mutationally robust. Here, we examine whether forces that change the distribution of phenotypes available by mutation profoundly alter subsequent evolutionary dynamics.We compare evolved populations of digital organisms that were subject to either static or cyclically-changing environments. For each of these, we examine diversity of the phenotypes that are produced through mutations in order to characterize the local genotype-phenotype map. We demonstrate that environmental change can push populations toward more evolvable mutational landscapes where many alternate phenotypes are available, though purely deleterious mutations remain suppressed. Further, we show that populations in environments with harsh changes switch phenotypes more readily than those in environments with more benign changes. We trace this effect to repeated population bottlenecks in the harsh environments, which result in shorter coalescence times and keep populations in regions of the mutational landscape where the phenotypic shifts in question are more likely to occur. Typically, static environments select solely for immediate optimization, at the expensive of long-term evolvability. In contrast, we show that with changing environments, short-term pressures to deal with immediate challenges can align with long-term pressures to explore a more productive portion of the mutational landscape.

Imagine if we could compute across phenotype data as easily as genomic data; this article calls for efforts to realize this vision and discusses the potential benefits., Despite a large and multifaceted effort to understand the vast landscape of phenotypic data, their current form inhibits productive data analysis. The lack of a community-wide, consensus-based, human- and machine-interpretable language for describing phenotypes and their genomic and environmental contexts is perhaps the most pressing scientific bottleneck to integration across many key fields in biology, including genomics, systems biology, development, medicine, evolution, ecology, and systematics. Here we survey the current phenomics landscape, including data resources and handling, and the progress that has been made to accurately capture relevant data descriptions for phenotypes. We present an example of the kind of integration across domains that computable phenotypes would enable, and we call upon the broader biology community, publishers, and relevant funding agencies to support efforts to surmount today's data barriers and facilitate analytical reproducibility.

The lake minnow (Eupallasella percnurus) is critically endangered. In this paper we characterize the genetic properties of this fish over its range of occurrence in Poland and propose the use of this knowledge in its active protection. Twelve populations of lake minnow from across its range in Poland were investigated. 13 microsatellite loci were investigated to evaluate genetic variation and distance among populations. The magnitude of the genetic bottleneck or founder effects was investigated. In the studied populations, the allelic diversity and heterozygosity showed that genetic variation in this species is low. At most loci, only 2–3 alleles per population were detected. The average number of alleles detected across all loci was 35, and ranged from 24 to 53. The average observed heterozygosity (Ho) across all investigated loci was 0.38 (range 0.21–0.59); the average expected heterozygosity (He) was 0.36 (range 0.18–0.55). The populations remained in Hardy-Weinberg equilibrium. The average Garza-Williamson M index value for all populations was low (0.47), suggesting a reduction in genetic variation due to a founder effect or a genetic bottleneck. Genetic distance among populations was high or very high (FST range: 0.20–0.64; δμ2 range: 1.32–16.98); this was likely a consequence of low gene flow among isolated populations, a founder effect or other genetic bottleneck, and strong genetic drift. The large genetic differences among the investigated lake minnow populations are likely to also exist among other populations of this species, and knowledge of these differences should inform active protection programs based on translocation of wild or cultivated fish of this species. The method presented here can potentially be applied to any population of lake minnows or closely related species.

The 21st century has thus far been marked by a heroic effort in genomic science and technology. If not yet upon us, the age of personalized genomic medicine appears imminent. Deriving medically relevant insight from these advances requires the ability to interpret the genetic variation and similarity observed in the population. One approach for interpreting genetic variation is the use of bioinformatics methods; simply put, these approaches are unified by a reliance on molecular/biological information (DNA, RNA, and protein sequence and annotation, protein structure, etc.) [1–4]. Many bioinformatics classifiers, primarily focused on amino acid substitution variants (cSNVs), have been and continue to be developed, typically achieving classification accuracies much better than random and thus supporting the use of molecular information [5–9]. These methods appear, however, to have reached a performance bottleneck; currently realized limits in performance all but forbid the consultation of bioinformatics cSNV classifiers in clinical settings. This bottleneck might be, in part, the result of simplifying assumptions inherent to the prediction of qualitative, dichotomous, or categorical phenotypes from individual missense variants. A complementary and beneficial strategy could be greater use of endophenotypes, quantitative measurements that are related to phenotypes via shared underlying genetics [10,11]. Endophenotypes can include phenomena at diverse biological scales; some examples include protein catalytic rate or melting temperature (stability), cell growth rate, and blood pressure. In this perspective we make a case for the increased use of endophenotypes, beginning with a brief overview of in silico bioinformatics methods for assessing phenotypic impact of cSNVs and the performance reported in recently published comparative studies. We compare the utility of endophenotypes and phenotypes in the context of CFTR cSNVs in cystic fibrosis, and of LDLR cSNV impact on cardiovascular-related diseases such as familial hypercholesterolemia. We also provide examples of bioinformatics methods that have been used to predict endophenotypes from cSNVs.

Color vision in primates is variable across species, and it represents a rare trait in which the genetic mechanisms underlying phenotypic variation are fairly well-understood. Research on primate color vision has largely focused on adaptive explanations for observed variation, but it remains unclear why some species have trichromatic or polymorphic color vision while others are red-green color blind. Lemurs, in particular, are highly variable. While some species are polymorphic, many closely-related species are strictly dichromatic. We provide the first characterization of color vision in a wild population of red-bellied lemurs (Eulemur rubriventer, Ranomafana National Park, Madagascar) with a sample size (87 individuals; NX chromosomes = 134) large enough to detect even rare variants (0.95 probability of detection at ≥ 3% frequency). By sequencing exon 5 of the X-linked opsin gene we identified opsin spectral sensitivity based on known diagnostic sites and found this population to be dichromatic and monomorphic for a long wavelength allele. Apparent fixation of this long allele is in contrast to previously published accounts of Eulemur species, which exhibit either polymorphic color vision or only the medium wavelength opsin. This unexpected result may represent loss of color vision variation, which could occur through selective processes and/or genetic drift (e.g., genetic bottleneck). To indirectly assess the latter scenario, we genotyped 55 adult red-bellied lemurs at seven variable microsatellite loci and used heterozygosity excess and M-ratio tests to assess if this population may have experienced a recent genetic bottleneck. Results of heterozygosity excess but not M-ratio tests suggest a bottleneck might have occurred in this red-bellied lemur population. Therefore, while selection may also play a role, the unique color vision observed in this population might have been influenced by a recent genetic bottleneck. These results emphasize the need to consider adaptive and nonadaptive mechanisms of color vision evolution in primates.

The Siamese crocodile (Crocodylus siamensis) and Saltwater crocodile (C. porosus) are two of the most endangered animals in Thailand. Their numbers have been reduced severely by hunting and habitat fragmentation. A reintroduction plan involving captive-bred populations that are used commercially is important and necessary as a conservation strategy to aid in the recovery of wild populations. Here, the genetic diversity and population structure of 69 individual crocodiles, mostly members of captive populations, were analyzed using both mitochondrial D-loop DNA and microsatellite markers. The overall haplotype diversity was 0.924-0.971 and the mean expected heterozygosity across 22 microsatellite loci was 0.578-0.701 for the two species. This agreed with the star-like shaped topology of the haplotype network, which suggests a high level of genetic diversity. The mean ratio of the number of alleles to the allelic range (M ratio) for the populations of both species was considerably lower than the threshold of 0.68, which was interpreted as indicative of a historical genetic bottleneck. Microsatellite markers provided evidence of introgression for three individual crocodiles, which suggest that hybridization might have occurred between C. siamensis and C. porosus. D-loop sequence analysis detected bi-directional hybridization between male and female individuals of the parent species. Therefore, identification of genetically non-hybrid and hybrid individuals is important for long-term conservation management. Relatedness values were low within the captive populations, which supported their genetic integrity and the viability of a breeding and reintroduction management plan. This work constitutes the first step in establishing an appropriate source population from a scientifically managed perspective for an in situ/ex situ conservation program and reintroduction of crocodile individuals to the wild in Thailand.

With the constant improvement in cost-efficiency and quality of Next Generation Sequencing technologies, shotgun-sequencing approaches -such as metagenomics- have nowadays become the methods of choice for studying and classifying microorganisms from various habitats. The production of data has dramatically increased over the past years and processing and analysis steps are becoming more and more of a bottleneck. Limiting factors are partly the availability of computational resources, but mainly the bioinformatics expertise in establishing and applying appropriate processing and analysis pipelines. Fortunately, a large diversity of specialized software tools is nowadays available. Nevertheless, choosing the most appropriate methods for answering specific biological questions can be rather challenging, especially for non-bioinformaticians. In order to provide a comprehensive overview and guide for the microbiological scientific community, we assessed the most common and freely available metagenome assembly tools with respect to their output statistics, their sensitivity for low abundant community members and variability in resulting community profiles as well as their ease-of-use. In contrast to the highly anticipated "Critical Assessment of Metagenomic Interpretation" (CAMI) challenge, which uses general mock community-based assembler comparison we here tested assemblers on real Illumina metagenome sequencing data from natural communities of varying complexity sampled from forest soil and algal biofilms. Our observations clearly demonstrate that different assembly tools can prove optimal, depending on the sample type, available computational resources and, most importantly, the specific research goal. In addition, we present detailed descriptions of the underlying principles and pitfalls of publically available assembly tools from a microbiologist's perspective, and provide guidance regarding the user-friendliness, sensitivity and reliability of the resulting phylogenetic profiles.

Whole-exome or gene targeted resequencing in hundreds to thousands of individuals has shown that the majority of genetic variants are at low frequency in human populations. Rare variants are enriched for functional mutations and are expected to explain an important fraction of the genetic etiology of human disease, therefore having a potential medical interest. In this work, we analyze the whole-exome sequences of French-Canadian individuals, a founder population with a unique demographic history that includes an original population bottleneck less than 20 generations ago, followed by a demographic explosion, and the whole exomes of French individuals sampled from France. We show that in less than 20 generations of genetic isolation from the French population, the genetic pool of French-Canadians shows reduced levels of diversity, higher homozygosity, and an excess of rare variants with low variant sharing with Europeans. Furthermore, the French-Canadian population contains a larger proportion of putatively damaging functional variants, which could partially explain the increased incidence of genetic disease in the province. Our results highlight the impact of population demography on genetic fitness and the contribution of rare variants to the human genetic variation landscape, emphasizing the need for deep cataloguing of genetic variants by resequencing worldwide human populations in order to truly assess disease risk., Author Summary Recent resequencing of the whole genome or the coding part of the genome (the exome) in thousands of individuals has described a large excess of low frequency variants in humans, probably arising as a consequence of recent rapid growth in human population sizes. Most rare variants are private to specific populations and are enriched for functional mutations, thus potentially having some medical relevance. In this study, we analyze whole-exome sequences from over a hundred individuals from the French-Canadian population, which was founded less than 400 years ago by about 8,500 French settlers who colonized the province between the 17th and 18th centuries. We show that in a remarkably short period of time this population has accumulated substantial differences, including an excess of rare, functional and potentially damaging variants, when compared to the original European population. Our results show the effects of population history on genetic variation that may have an impact on genetic fitness and disease, and have implications in the design of genetic studies, highlighting the importance of extending deep resequencing to worldwide human populations.

The amount of genetic diversity in a finite biological population mostly depends on the interactions among evolutionary forces and the effective population size (N(e)) as well as the time since population establishment. Because the N(e) estimation helps to explore population demographic history, and allows one to predict the behavior of genetic diversity through time, N(e) is a key parameter for the genetic management of small and isolated populations. Here, we explored an N(e)-based approach using a bighorn sheep population on Tiburon Island, Mexico (TI) as a model. We estimated the current (N(crnt)) and ancestral stable (N(stbl)) inbreeding effective population sizes as well as summary statistics to assess genetic diversity and the demographic scenarios that could explain such diversity. Then, we evaluated the feasibility of using TI as a source population for reintroduction programs. We also included data from other bighorn sheep and artiodactyl populations in the analysis to compare their inbreeding effective size estimates. The TI population showed high levels of genetic diversity with respect to other managed populations. However, our analysis suggested that TI has been under a genetic bottleneck, indicating that using individuals from this population as the only source for reintroduction could lead to a severe genetic diversity reduction. Analyses of the published data did not show a strict correlation between H(E) and N(crnt) estimates. Moreover, we detected that ancient anthropogenic and climatic pressures affected all studied populations. We conclude that the estimation of N(crnt) and N(stbl) are informative genetic diversity estimators and should be used in addition to summary statistics for conservation and population management planning.

Background: Viral outbreaks, such as the 2014 ebolavirus, can spread rapidly and have complex evolutionary dynamics, including coinfection and bulk transmission of multiple viral populations. Genomic surveillance can be hindered when the spread of the outbreak exceeds the evolutionary rate, in which case consensus approaches will have limited resolution. Deep sequencing of infected patients can identify genomic variants present in intrahost populations at subclonal frequencies (i.e.

The Gillespie Stochastic Simulation Algorithm (GSSA) and its variants are cornerstone techniques to simulate reaction kinetics in situations where the concentration of the reactant is too low to allow deterministic techniques such as differential equations. The inherent limitations of the GSSA include the time required for executing a single run and the need for multiple runs for parameter sweep exercises due to the stochastic nature of the simulation. Even very efficient variants of GSSA are prohibitively expensive to compute and perform parameter sweeps. Here we present a novel variant of the exact GSSA that is amenable to acceleration by using graphics processing units (GPUs). We parallelize the execution of a single realization across threads in a warp (fine-grained parallelism). A warp is a collection of threads that are executed synchronously on a single multi-processor. Warps executing in parallel on different multi-processors (coarse-grained parallelism) simultaneously generate multiple trajectories. Novel data-structures and algorithms reduce memory traffic, which is the bottleneck in computing the GSSA. Our benchmarks show an 8×−120× performance gain over various state-of-the-art serial algorithms when simulating different types of models.

All existing sea otter, Enhydra lutris, populations have suffered at least one historic population bottleneck stemming from the fur trade extirpations of the eighteenth and nineteenth centuries. We examined genetic variation, gene flow, and population structure at five microsatellite loci in samples from five pre-fur trade populations throughout the sea otter's historical range: California, Oregon, Washington, Alaska, and Russia. We then compared those values to genetic diversity and population structure found within five modern sea otter populations throughout their current range: California, Prince William Sound, Amchitka Island, Southeast Alaska and Washington. We found twice the genetic diversity in the pre-fur trade populations when compared to modern sea otters, a level of diversity that was similar to levels that are found in other mammal populations that have not experienced population bottlenecks. Even with the significant loss in genetic diversity modern sea otters have retained historical structure. There was greater gene flow before extirpation than that found among modern sea otter populations but the difference was not statistically significant. The most dramatic effect of pre fur trade population extirpation was the loss of genetic diversity. For long term conservation of these populations increasing gene flow and the maintenance of remnant genetic diversity should be encouraged.

Experiments in systems neuroscience can be seen as consisting of three steps: (1) selecting the signals we are interested in, (2) probing the system with carefully chosen stimuli, and (3) getting data out of the brain. Here I discuss how emerging techniques in molecular biology are starting to improve these three steps. To estimate its future impact on experimental neuroscience, I will stress the analogy of ongoing progress with that of microprocessor production techniques. These techniques have allowed computers to simplify countless problems; because they are easier to use than mechanical timers, they are even built into toasters. Molecular biology may advance even faster than computer speeds and has made immense progress in understanding and designing molecules. These advancements may in turn produce impressive improvements to each of the three steps, ultimately shifting the bottleneck from obtaining data to interpreting it.

BackgroundThe widespread popularity of genomic applications is threatened by the "bioinformatics bottleneck" resulting from uncertainty about the cost and infrastructure needed to meet increasing demands for next-generation sequence analysis. Cloud computing services have been discussed as potential new bioinformatics support systems but have not been evaluated thoroughly.ResultsWe present benchmark costs and runtimes for common microbial genomics applications, including 16S rRNA analysis, microbial whole-genome shotgun (WGS) sequence assembly and annotation, WGS metagenomics and large-scale BLAST. Sequence dataset types and sizes were selected to correspond to outputs typically generated by small- to midsize facilities equipped with 454 and Illumina platforms, except for WGS metagenomics where sampling of Illumina data was used. Automated analysis pipelines, as implemented in the CloVR virtual machine, were used in order to guarantee transparency, reproducibility and portability across different operating systems, including the commercial Amazon Elastic Compute Cloud (EC2), which was used to attach real dollar costs to each analysis type. We found considerable differences in computational requirements, runtimes and costs associated with different microbial genomics applications. While all 16S analyses completed on a single-CPU desktop in under three hours, microbial genome and metagenome analyses utilized multi-CPU support of up to 120 CPUs on Amazon EC2, where each analysis completed in under 24 hours for less than $60. Representative datasets were used to estimate maximum data throughput on different cluster sizes and to compare costs between EC2 and comparable local grid servers.ConclusionsAlthough bioinformatics requirements for microbial genomics depend on dataset characteristics and the analysis protocols applied, our results suggests that smaller sequencing facilities (up to three Roche/454 or one Illumina GAIIx sequencer) invested in 16S rRNA amplicon sequencing, microbial single-genome and metagenomics WGS projects can achieve cost-efficient bioinformatics support using CloVR in combination with Amazon EC2 as an alternative to local computing centers.

Recent findings suggest a convergence of time scales between ecological and evolutionary processes which is usually explained in terms of rapid micro evolution resulting in evolution on ecological time scales. A similar convergence, however, can also emerge when slow ecological processes take place on evolutionary time scales. A good example of such a slow ecological process is the colonization of remote aquatic habitats by passively dispersed zooplankton. Using variation at the protein coding mitochondrial COI gene, we investigated the balance between mutation and migration as drivers of genetic diversity in two Branchipodopsis fairy shrimp species (Crustacea, Anostraca) endemic to remote temporary rock pool clusters at the summit of isolated mountaintops in central South Africa. We showed that both species colonized the region almost simultaneously c. 0.8 My ago, but exhibit contrasting patterns of regional genetic diversity and demographic history. The haplotype network of the common B. cf. wolfi showed clear evidence of 11 long distance dispersal events (up to 140 km) with five haplotypes that are shared among distant inselbergs, as well as some more spatially isolated derivates. Similar patterns were not observed for B. drakensbergensis presumably since this rarer species experienced a genetic bottleneck. We conclude that the observed genetic patterns reflect rare historic colonization events rather than frequent ongoing gene flow. Moreover, the high regional haplotype diversity combined with a high degree of haplotype endemicity indicates that evolutionary- (mutation) and ecological (migration) processes in this system operate on similar time scales. ispartof: PLoS One vol:6 issue:11 ispartof: location:United States status: published

Phylogenetic tree reconstruction is traditionally based on multiple sequence alignments (MSAs) and heavily depends on the validity of this information bottleneck. With increasing sequence divergence, the quality of MSAs decays quickly. Alignment-free methods, on the other hand, are based on abstract string comparisons and avoid potential alignment problems. However, in general they are not biologically motivated and ignore our knowledge about the evolution of sequences. Thus, it is still a major open question how to define an evolutionary distance metric between divergent sequences that makes use of indel information and known substitution models without the need for a multiple alignment. Here we propose a new evolutionary distance metric to close this gap. It uses finite-state transducers to create a biologically motivated similarity score which models substitutions and indels, and does not depend on a multiple sequence alignment. The sequence similarity score is defined in analogy to pairwise alignments and additionally has the positive semi-definite property. We describe its derivation and show in simulation studies and real-world examples that it is more accurate in reconstructing phylogenies than competing methods. The result is a new and accurate way of determining evolutionary distances in and beyond the twilight zone of sequence alignments that is suitable for large datasets., Comment: to appear in PLoS ONE

Background Mother-to-child transmission of HIV-1 remains a significant problem in the resource-constrained settings where anti-retroviral therapy is still not widely available. Understanding the earliest events during HIV-1 transmission and characterizing the newly transmitted or founder virus is central to intervention efforts. In this study, we analyzed the viral env quasispecies of six mother-infant transmission pairs (MIPs) and characterized the genetic features of envelope glycoprotein that could influence HIV-1 subtype C perinatal transmission. Methodology and Findings The V1-V5 region of env was amplified from 6 MIPs baseline samples and 334 DNA sequences in total were analyzed. A comparison of the viral population derived from the mother and infant revealed a severe genetic bottleneck occurring during perinatal transmission, which was characterized by low sequence diversity in the infant. Phylogenetic analysis indicates that most likely in all our infant subjects a single founder virus was responsible for establishing infection. Furthermore, the newly transmitted viruses from the infant had significantly fewer potential N-linked glycosylation sites in Env V1-V5 region and showed a propensity to encode shorter variable loops compared to the nontransmitted viruses. In addition, a similar intensity of selection was seen between mothers and infants with a higher rate of synonymous (dS) compared to nonsynonymous (dN) substitutions evident (dN/dS

Uganda is the only country where the chronic and acute forms of human African Trypanosomiasis (HAT) or sleeping sickness both occur and are separated by < 100 km in areas north of Lake Kyoga. In Uganda, Glossina fuscipes fuscipes is the main vector of the Trypanosoma parasites responsible for these diseases as well for the animal African Trypanosomiasis (AAT), or Nagana. We used highly polymorphic microsatellite loci and a mitochondrial DNA (mtDNA) marker to provide fine scale spatial resolution of genetic structure of G. f. fuscipes from 42 sampling sites from the northern region of Uganda where a merger of the two disease belts is feared. Based on microsatellite analyses, we found that G. f. fuscipes in northern Uganda are structured into three distinct genetic clusters with varying degrees of interconnectivity among them. Based on genetic assignment and spatial location, we grouped the sampling sites into four genetic units corresponding to northwestern Uganda in the Albert Nile drainage, northeastern Uganda in the Lake Kyoga drainage, western Uganda in the Victoria Nile drainage, and a transition zone between the two northern genetic clusters characterized by high level of genetic admixture. An analysis using HYBRIDLAB supported a hybrid swarm model as most consistent with tsetse genotypes in these admixed samples. Results of mtDNA analyses revealed the presence of 30 haplotypes representing three main haplogroups, whose location broadly overlaps with the microsatellite defined clusters. Migration analyses based on microsatellites point to moderate migration among the northern units located in the Albert Nile, Achwa River, Okole River, and Lake Kyoga drainages, but not between the northern units and the Victoria Nile drainage in the west. Effective population size estimates were variable with low to moderate sizes in most populations and with evidence of recent population bottlenecks, especially in the northeast unit of the Lake Kyoga drainage. Our microsatellite and mtDNA based analyses indicate that G. f. fuscipes movement along the Achwa and Okole rivers may facilitate northwest expansion of the Rhodesiense disease belt in Uganda. We identified tsetse migration corridors and recommend a rolling carpet approach from south of Lake Kyoga northward to minimize disease dispersal and prevent vector re-colonization. Additionally, our findings highlight the need for continuing tsetse monitoring efforts during and after control., Author summary Northern Uganda is an epidemiologically important region affected by human African trypanosomiasis (HAT) because it harbors both forms of the HAT disease (T. b. gambiense and T. b. rhodesiense). The geographic location of this region creates the risk that these distinct foci could merge, which would complicate diagnosis and treatment, and may result in recombination between the two parasite strains with as yet unknown consequences. Both strains require a tsetse fly vector for transmission, and in Uganda, G. f. fuscipes is the major vector of HAT. Controlling the vector remains one of the most effective strategies for controlling trypanosome parasites. However, vector control efforts may not be sustainable in terms of long term reduction in G. f. fuscipes populations due to population rebounds. Population genetics data can allow us to determine the likely source of population rebounds and to establish a more robust control strategy. In this study, we build on a previous broad spatial survey of G. f. fuscipes genetic structure in Uganda by adding more than 30 novel sampling sites that are strategically spaced across a region of northern Uganda that, for historical and political reasons, was severely understudied and faces particularly high disease risk. We identify natural population breaks, migration corridors and a hybrid zone with evidence of free interbreeding of G. f. fuscipes across the geographic region that spans the two HAT disease foci. We also find evidence of low effective population sizes and population bottlenecks in some areas that have been subjects of past control but remain regions of high tsetse density, which stresses the risk of population rebounds if monitoring is not explicitly incorporated into the control strategy. We use these results to make suggestions that will enhance the design and implementation of control activities in northern Uganda.