Characterization of the TyrR Regulon in the Rhizobacterium Enterobacter ludwigii UW5 Reveals Overlap with the CpxR Envelope Stress Response.

The TyrR transcription factor controls the expression of genes for the uptake and biosynthesis of aromatic amino acids in Escherichia coli In the plant-associated and clinically significant proteobacterium Enterobacter ludwigii UW5, the TyrR orthologue was previously shown to regulate genes that encode enzymes for synthesis of the plant hormone indole-3-acetic acid and for gluconeogenesis, indicating a broader function for the transcription factor. This study aimed to delineate the TyrR regulon of E. ludwigii by comparing the transcriptomes of the wild type and a tyrR deletion strain. In E. ludwigii, TyrR positively or negatively regulates the expression of over 150 genes. TyrR downregulated expression of envelope stress response regulators CpxR and CpxP through interaction with a DNA binding site in the intergenic region between divergently transcribed cpxP and cpxR Repression of cpxP was alleviated by tyrosine. Methyltransferase gene dmpM, which is possibly involved in antibiotic synthesis, was strongly activated in the presence of tyrosine and phenylalanine by TyrR binding to its promoter region. TyrR also regulated expression of genes for aromatic catabolism and anaerobic respiration. Our findings suggest that the E. ludwigii TyrR regulon has diverged from that of E. coli to include genes for survival in the diverse environments that this bacterium inhabits and illustrate the expansion and plasticity of transcription factor regulons.IMPORTANCE Genome-wide RNA sequencing revealed a broader regulatory role for the TyrR transcription factor in the ecologically versatile bacterium Enterobacter ludwigii beyond that of aromatic amino acid synthesis and transport that constitute the role of the TyrR regulon of E. coli In E. ludwigii, a plant symbiont and human gut commensal, the TyrR regulon is expanded to include genes that are beneficial for plant interactions and response to stresses. Identification of the genes regulated by TyrR provides insight into the mechanisms by which the bacterium adapts to its environment.

Confidently identifying the correct K value using the ΔK method: When does K = 2?

Populations delineated based on genetic data are commonly used for wildlife conservation and management. Many studies use the program structure combined with the ΔK method to identify the most probable number of populations (K). We recently found K = 2 was identified more often when studies used ΔK compared to studies that did not. We suggested two reasons for this: hierarchical population structure leads to underestimation, or the ΔK method does not evaluate K = 1 causing an overestimation. The present contribution aims to develop a better understanding of the limits of the method using one, two and three population simulations across migration scenarios. From these simulations we identified the "best K" using model likelihood and ΔK. Our findings show that mean probability plots and ΔK are unable to resolve the correct number of populations once migration rate exceeds 0.005. We also found a strong bias towards selecting K = 2 using the ΔK method. We used these data to identify the range of values where the ΔK statistic identifies a value of K that is not well supported. Finally, using the simulations and a review of empirical data, we found that the magnitude of ΔK corresponds to the level of divergence between populations. Based on our findings, we suggest researchers should use the ΔK method cautiously; they need to report all relevant data, including the magnitude of ΔK, and an estimate of connectivity for the research community to assess whether meaningful genetic structure exists within the context of management and conservation.

Heritability of body size in the polar bears of Western Hudson Bay.

Among polar bears (Ursus maritimus), fitness is dependent on body size through males' abilities to win mates, females' abilities to provide for their young and all bears' abilities to survive increasingly longer fasting periods caused by climate change. In the Western Hudson Bay subpopulation (near Churchill, Manitoba, Canada), polar bears have declined in body size and condition, but nothing is known about the genetic underpinnings of body size variation, which may be subject to natural selection. Here, we combine a 4449-individual pedigree and an array of 5,433 single nucleotide polymorphisms (SNPs) to provide the first quantitative genetic study of polar bears. We used animal models to estimate heritability (h2 ) among polar bears handled between 1966 and 2011, obtaining h2 estimates of 0.34-0.48 for strictly skeletal traits and 0.18 for axillary girth (which is also dependent on fatness). We genotyped 859 individuals with the SNP array to test for marker-trait association and combined p-values over genetic pathways using gene-set analysis. Variation in all traits appeared to be polygenic, but we detected one region of moderately large effect size in body length near a putative noncoding RNA in an unannotated region of the genome. Gene-set analysis suggested that variation in body length was associated with genes in the regulatory cascade of cyclin expression, which has previously been associated with body size in mice. A greater understanding of the genetic architecture of body size variation will be valuable in understanding the potential for adaptation in polar bear populations challenged by climate change.

The K = 2 conundrum.

Assessments of population genetic structure have become an increasing focus as they can provide valuable insight into patterns of migration and gene flow. structure, the most highly cited of several clustering-based methods, was developed to provide robust estimates without the need for populations to be determined a priori. structure introduces the problem of selecting the optimal number of clusters, and as a result, the ΔK method was proposed to assist in the identification of the "true" number of clusters. In our review of 1,264 studies using structure to explore population subdivision, studies that used ΔK were more likely to identify K = 2 (54%, 443/822) than studies that did not use ΔK (21%, 82/386). A troubling finding was that very few studies performed the hierarchical analysis recommended by the authors of both ΔK and structure to fully explore population subdivision. Furthermore, extensions of earlier simulations indicate that, with a representative number of markers, ΔK frequently identifies K = 2 as the top level of hierarchical structure, even when more subpopulations are present. This review suggests that many studies may have been over- or underestimating population genetic structure; both scenarios have serious consequences, particularly with respect to conservation and management. We recommend publication standards for population structure results so that readers can assess the implications of the results given their own understanding of the species biology.

Circumpolar Genetic Structure and Recent Gene Flow of Polar Bears: A Reanalysis.

Recently, an extensive study of 2,748 polar bears (Ursus maritimus) from across their circumpolar range was published in PLOS ONE, which used microsatellites and mitochondrial haplotypes to apparently show altered population structure and a dramatic change in directional gene flow towards the Canadian Archipelago-an area believed to be a future refugium for polar bears as their southernmost habitats decline under climate change. Although this study represents a major international collaborative effort and promised to be a baseline for future genetics work, methodological shortcomings and errors of interpretation undermine some of the study's main conclusions. Here, we present a reanalysis of this data in which we address some of these issues, including: (1) highly unbalanced sample sizes and large amounts of systematically missing data; (2) incorrect calculation of FST and of significance levels; (3) misleading estimates of recent gene flow resulting from non-convergence of the program BayesAss. In contrast to the original findings, in our reanalysis we find six genetic clusters of polar bears worldwide: the Hudson Bay Complex, the Western and Eastern Canadian Arctic Archipelago, the Western and Eastern Polar Basin, and-importantly-we reconfirm the presence of a unique and possibly endangered cluster of bears in Norwegian Bay near Canada's expected last sea-ice refugium. Although polar bears' abundance, distribution, and population structure will certainly be negatively affected by ongoing-and increasingly rapid-loss of Arctic sea ice, these genetic data provide no evidence of strong directional gene flow in response to recent climate change.

Temporal dynamics of linkage disequilibrium in two populations of bighorn sheep.

Linkage disequilibrium (LD) is the nonrandom association of alleles at two markers. Patterns of LD have biological implications as well as practical ones when designing association studies or conservation programs aimed at identifying the genetic basis of fitness differences within and among populations. However, the temporal dynamics of LD in wild populations has received little empirical attention. In this study, we examined the overall extent of LD, the effect of sample size on the accuracy and precision of LD estimates, and the temporal dynamics of LD in two populations of bighorn sheep (Ovis canadensis) with different demographic histories. Using over 200 microsatellite loci, we assessed two metrics of multi-allelic LD, D', and χ ('2). We found that both populations exhibited high levels of LD, although the extent was much shorter in a native population than one that was founded via translocation, experienced a prolonged bottleneck post founding, followed by recent admixture. In addition, we observed significant variation in LD in relation to the sample size used, with small sample sizes leading to depressed estimates of the extent of LD but inflated estimates of background levels of LD. In contrast, there was not much variation in LD among yearly cross-sections within either population once sample size was accounted for. Lack of pronounced interannual variability suggests that researchers may not have to worry about interannual variation when estimating LD in a population and can instead focus on obtaining the largest sample size possible.

Design of a 9K illumina BeadChip for polar bears (Ursus maritimus) from RAD and transcriptome sequencing.

Single-nucleotide polymorphisms (SNPs) offer numerous advantages over anonymous markers such as microsatellites, including improved estimation of population parameters, finer-scale resolution of population structure and more precise genomic dissection of quantitative traits. However, many SNPs are needed to equal the resolution of a single microsatellite, and reliable large-scale genotyping of SNPs remains a challenge in nonmodel species. Here, we document the creation of a 9K Illumina Infinium BeadChip for polar bears (Ursus maritimus), which will be used to investigate: (i) the fine-scale population structure among Canadian polar bears and (ii) the genomic architecture of phenotypic traits in the Western Hudson Bay subpopulation. To this end, we used restriction-site associated DNA (RAD) sequencing from 38 bears across their circumpolar range, as well as blood/fat transcriptome sequencing of 10 individuals from Western Hudson Bay. Six-thousand RAD SNPs and 3000 transcriptomic SNPs were selected for the chip, based primarily on genomic spacing and gene function respectively. Of the 9000 SNPs ordered from Illumina, 8042 were successfully printed, and - after genotyping 1450 polar bears - 5441 of these SNPs were found to be well clustered and polymorphic. Using this array, we show rapid linkage disequilibrium decay among polar bears, we demonstrate that in a subsample of 78 individuals, our SNPs detect known genetic structure more clearly than 24 microsatellites genotyped for the same individuals and that these results are not driven by the SNP ascertainment scheme. Here, we present one of the first large-scale genotyping resources designed for a threatened species.

Genomic resources notes accepted 1 June 2013-31 July 2013.

This article documents the public availability of (i) raw transcriptome sequence data, assembled contigs and UniProt BLAST hits from common crossbill (Loxia curvirostra) and Plasmodium relictum (lineage SGS1) obtained from a controlled infection experiment; and (ii) raw transcriptome sequence data and 66 596 SNPs for the white-tailed deer (Odocoileus virginianus).

Genomic resources notes accepted 1 August 2013-30 September 2013.

This article documents the public availability of raw transcriptome sequence data and 63,020 SNPs for the polar bear (Ursus maritimus).

Short reads, circular genome: skimming solid sequence to construct the bighorn sheep mitochondrial genome.

As sequencing technology improves, an increasing number of projects aim to generate full genome sequence, even for nonmodel taxa. These projects may be feasibly conducted at lower read depths if the alignment can be aided by previously developed genomic resources from a closely related species. We investigated the feasibility of constructing a complete mitochondrial (mt) genome without preamplification or other targeting of the sequence. Here we present a full mt genome sequence (16,463 nucleotides) for the bighorn sheep (Ovis canadensis) generated though alignment of SOLiD short-read sequences to a reference genome. Average read depth was 1240, and each base was covered by at least 36 reads. We then conducted a phylogenomic analysis with 27 other bovid mitogenomes, which placed bighorn sheep firmly in the Ovis clade. These results show that it is possible to generate a complete mitogenome by skimming a low-coverage genomic sequencing library. This technique will become increasingly applicable as the number of taxa with some level of genome sequence rises.