Genetic architecture and evolution of color variation in American black bears.

Color variation is a frequent evolutionary substrate for camouflage in small mammals, but the underlying genetics and evolutionary forces that drive color variation in natural populations of large mammals are mostly unexplained. The American black bear, Ursus americanus (U. americanus), exhibits a range of colors including the cinnamon morph, which has a similar color to the brown bear, U. arctos, and is found at high frequency in the American southwest. Reflectance and chemical melanin measurements showed little distinction between U. arctos and cinnamon U. americanus individuals. We used a genome-wide association for hair color as a quantitative trait in 151 U. americanus individuals and identified a single major locus (p < 10-13). Additional genomic and functional studies identified a missense alteration (R153C) in Tyrosinase-related protein 1 (TYRP1) that likely affects binding of the zinc cofactor, impairs protein localization, and results in decreased pigment production. Population genetic analyses and demographic modeling indicated that the R153C variant arose 9.36 kya in a southwestern population where it likely provided a selective advantage, spreading both northwards and eastwards by gene flow. A different TYRP1 allele, R114C, contributes to the characteristic brown color of U. arctos but is not fixed across the range.

Aberrant regulation of a poison exon caused by a non-coding variant in a mouse model of Scn1a-associated epileptic encephalopathy.

Dravet syndrome (DS) is a developmental and epileptic encephalopathy that results from mutations in the Nav1.1 sodium channel encoded by SCN1A. Most known DS-causing mutations are in coding regions of SCN1A, but we recently identified several disease-associated SCN1A mutations in intron 20 that are within or near to a cryptic and evolutionarily conserved "poison" exon, 20N, whose inclusion is predicted to lead to transcript degradation. However, it is not clear how these intron 20 variants alter SCN1A expression or DS pathophysiology in an organismal context, nor is it clear how exon 20N is regulated in a tissue-specific and developmental context. We address those questions here by generating an animal model of our index case, NM_006920.4(SCN1A):c.3969+2451G>C, using gene editing to create the orthologous mutation in laboratory mice. Scn1a heterozygous knock-in (+/KI) mice exhibited an ~50% reduction in brain Scn1a mRNA and Nav1.1 protein levels, together with characteristics observed in other DS mouse models, including premature mortality, seizures, and hyperactivity. In brain tissue from adult Scn1a +/+ animals, quantitative RT-PCR assays indicated that ~1% of Scn1a mRNA included exon 20N, while brain tissue from Scn1a +/KI mice exhibited an ~5-fold increase in the extent of exon 20N inclusion. We investigated the extent of exon 20N inclusion in brain during normal fetal development in RNA-seq data and discovered that levels of inclusion were ~70% at E14.5, declining progressively to ~10% postnatally. A similar pattern exists for the homologous sodium channel Nav1.6, encoded by Scn8a. For both genes, there is an inverse relationship between the level of functional transcript and the extent of poison exon inclusion. Taken together, our findings suggest that poison exon usage by Scn1a and Scn8a is a strategy to regulate channel expression during normal brain development, and that mutations recapitulating a fetal-like pattern of splicing cause reduced channel expression and epileptic encephalopathy.

A comparison of genotyping-by-sequencing analysis methods on low-coverage crop datasets shows advantages of a new workflow, GB-eaSy.

BACKGROUND: Genotyping-by-sequencing (GBS), a method to identify genetic variants and quickly genotype samples, reduces genome complexity by using restriction enzymes to divide the genome into fragments whose ends are sequenced on short-read sequencing platforms. While cost-effective, this method produces extensive missing data and requires complex bioinformatics analysis. GBS is most commonly used on crop plant genomes, and because crop plants have highly variable ploidy and repeat content, the performance of GBS analysis software can vary by target organism. Here we focus our analysis on soybean, a polyploid crop with a highly duplicated genome, relatively little public GBS data and few dedicated tools. RESULTS: We compared the performance of five GBS pipelines using low-coverage Illumina sequence data from three soybean populations. To address issues identified with existing methods, we developed GB-eaSy, a GBS bioinformatics workflow that incorporates widely used genomics tools, parallelization and automation to increase the accuracy and accessibility of GBS data analysis. Compared to other GBS pipelines, GB-eaSy rapidly and accurately identified the greatest number of SNPs, with SNP calls closely concordant with whole-genome sequencing of selected lines. Across all five GBS analysis platforms, SNP calls showed unexpectedly low convergence but generally high accuracy, indicating that the workflows arrived at largely complementary sets of valid SNP calls on the low-coverage data analyzed. CONCLUSIONS: We show that GB-eaSy is approximately as good as, or better than, other leading software solutions in the accuracy, yield and missing data fraction of variant calling, as tested on low-coverage genomic data from soybean. It also performs well relative to other solutions in terms of the run time and disk space required. In addition, GB-eaSy is built from existing open-source, modular software packages that are regularly updated and commonly used, making it straightforward to install and maintain. While GB-eaSy outperformed other individual methods on the datasets analyzed, our findings suggest that a comprehensive approach integrating the results from multiple GBS bioinformatics pipelines may be the optimal strategy to obtain the largest, most highly accurate SNP yield possible from low-coverage polyploid sequence data.

Distinct roles of the Pumilio and FBF translational repressors during C. elegans vulval development.

The C. elegans PUF and FBF proteins regulate various aspects of germline development by selectively binding to the 3' untranslated region of their target mRNAs and repressing translation. Here, we show that puf-8, fbf-1 and fbf-2 also act in the soma where they negatively regulate vulvaI development. Loss-of-function mutations in puf-8 cause ectopic vulval differentiation when combined with mutations in negative regulators of the EGFR/RAS/MAPK pathway and suppress the vulvaless phenotype caused by mutations that reduce EGFR/RAS/MAPK signalling. PUF-8 acts cell-autonomously in the vulval cells to limit their temporal competence to respond to the extrinsic patterning signals. fbf-1 and fbf-2, however, redundantly inhibit primary vulval cell fate specification in two distinct pathways acting in the soma and in the germline. The FBFs thereby ensure that the inductive signal selects only one vulval precursor cell for the primary cell fate. Thus, translational repressors regulate various aspects of vulval cell fate specification, and they may play a conserved role in modulating signal transduction during animal development.

The C. elegans G-protein-coupled receptor SRA-13 inhibits RAS/MAPK signalling during olfaction and vulval development.

In C. elegans, the RAS/MAPK pathway is used in different tissues to regulate various cell fate decisions. Several positive and negative regulators tightly control the activity of the RAS/MAPK pathway at different steps. We demonstrate a link between a G-protein-coupled receptor signalling pathway and the RAS/MAPK cascade. SRA-13, a member of the SRA family of chemosensory receptors, negatively regulates RAS/MAPK signalling during vulval induction and the olfaction of volatile attractants. Epistasis analysis indicates that SRA-13 inhibits the RAS/MAPK pathway at the level or upstream of MAPK. In both tissues, the vulval precursor cells and the chemosensory neurones, SRA-13 acts through the GPA-5 Galpha protein subunit, suggesting a common mechanism of crosstalk. Moreover, we find that vulval induction is repressed by food withdrawal during larval development and that SRA-13 activity is required for the suppression of vulval induction in response to food starvation. Thus, SRA-13 may serve to adapt the activity of the RAS/MAPK pathway to environmental conditions.