Description: Here, we provide revised gene models for D . ananassae , D. yakuba , and D. simulans , which include untranslated regions and empirically verified intron-exon boundaries, as well as ortholog groups identified using a fuzzy reciprocal-best-hit blast comparison. Using these revised annotations, we perform differential expression testing using the cufflinks suite to provide a broad overview of differential expression between reproductive tissues and the carcass. We identify thousands of genes that are differentially expressed across tissues in D. yakuba and D. simulans , with roughly 60% agreement in expression patterns of orthologs in D. yakuba and D. simulans . We identify several cases of putative polycistronic transcripts, pointing to a combination of transcriptional read-through in the genome as well as putative gene fusion and fission events across taxa. We furthermore identify hundreds of lineage specific genes in each species with no blast hits among transcripts of any other Drosophila species, which are candidates for neofunctionalized proteins and a potential source of genetic novelty.

Description: Oocyte maturation in all species is controlled by a protein complex termed the maturation promoting factor (MPF). MPF comprises a cyclin-dependent kinase (CDK) and its partner cyclin, and it is regulated by dueling regulatory phosphorylation events on the CDK. In Caenorhabditis elegans , the Wee1/Myt1 ortholog WEE-1.3 provides the inhibitory phosphorylations on CDK-1 that keep MPF inactive and halt meiosis. Prior work has shown that depletion of WEE-1.3 in C. elegans results in precocious oocyte maturation in vivo and a highly penetrant infertility phenotype. This study sought to further define the precocious maturation phenotype and to identify novel interactors with WEE-1.3 . We found that WEE-1.3 is expressed throughout the germline and in developing embryos in a perinuclear pattern, and demonstrated that oocytes in WEE-1.3 –depleted germlines have begun to transcribe embryonic genes and exhibit inappropriate expression of proteins normally restricted to fertilized eggs. In addition, we performed an RNAi suppressor screen of the infertile phenotype to identify novel factors that, when co-depleted with WEE-1.3 , restore fertility to these animals. We screened ~1900 essential genes by RNAi feeding and identified 44 (~2% of the tested genes) that are suppressors of the WEE-1.3 depletion phenotype. The suppressors include many previously unidentified players in the meiotic cell cycle and represent a pool of potential WEE-1.3 interacting proteins that function during C. elegans oocyte maturation and zygotic development.

Description: Quantitative trait loci (QTL) mapping of wood properties in conifer species has focused on single time point measurements or on trait means based on heterogeneous wood samples ( e.g. , increment cores), thus ignoring systematic within-tree trends. In this study, functional QTL mapping was performed for a set of important wood properties in increment cores from a 17-yr-old Scots pine ( Pinus sylvestris L.) full-sib family with the aim of detecting wood trait QTL for general intercepts (means) and for linear slopes by increasing cambial age. Two multi-locus functional QTL analysis approaches were proposed and their performances were compared on trait datasets comprising 2 to 9 time points, 91 to 455 individual tree measurements and genotype datasets of amplified length polymorphisms (AFLP), and single nucleotide polymorphism (SNP) markers. The first method was a multilevel LASSO analysis whereby trend parameter estimation and QTL mapping were conducted consecutively; the second method was our Bayesian linear mixed model whereby trends and underlying genetic effects were estimated simultaneously. We also compared several different hypothesis testing methods under either the LASSO or the Bayesian framework to perform QTL inference. In total, five and four significant QTL were observed for the intercepts and slopes, respectively, across wood traits such as earlywood percentage, wood density, radial fiberwidth, and spiral grain angle. Four of these QTL were represented by candidate gene SNPs, thus providing promising targets for future research in QTL mapping and molecular function. Bayesian and LASSO methods both detected similar sets of QTL given datasets that comprised large numbers of individuals.

Description: The resistance gene Pi-ta has been effectively used to control rice blast disease, but some populations of cultivated and wild rice have evolved resistance. Insights into the evolutionary processes that led to this resistance during crop domestication may be inferred from the population history of domesticated and wild rice strains. In this study, we applied a recently developed statistical method, time-dependent Poisson random field model, to examine the evolution of the Pi-ta gene in cultivated and weedy rice. Our study suggests that the Pi-ta gene may have more recently introgressed into cultivated rice, indica and japonica , and U.S. weedy rice from the wild species, O. rufipogon . In addition, the Pi-ta gene is under positive selection in japonica , tropical japonica , U.S. cultivars and U.S. weedy rice . We also found that sequences of two domains of the Pi-ta gene, the nucleotide binding site and leucine-rich repeat domain, are highly conserved among all rice accessions examined. Our results provide a valuable analytical tool for understanding the evolution of disease resistance genes in crop plants.

Description: There has been a renewed interest in investigating the role of stabilizing selection acting on genome-wide traits such as codon usage bias. Codon bias, when synonymous codons are used at unequal frequencies, occurs in a wide variety of taxa. Standard evolutionary models explain the maintenance of codon bias through a balance of genetic drift, mutation and weak purifying selection. The efficacy of selection is expected to be reduced in regions of suppressed recombination. Contrary to observations in Drosophila melanogaster , some recent studies have failed to detect a relationship between the recombination rate, intensity of selection acting at synonymous sites, and the magnitude of codon bias as predicted under these standard models. Here, we examined codon bias in 2798 protein coding loci on the third chromosome of D. pseudoobscura using whole-genome sequences of 47 individuals, representing five common third chromosome gene arrangements. Fine-scale recombination maps were constructed using more than 1 million segregating sites. As expected, recombination was demonstrated to be significantly suppressed between chromosome arrangements, allowing for a direct examination of the relationship between recombination, selection, and codon bias. As with other Drosophila species, we observe a strong mutational bias away from the most frequently used codons. We find the rate of synonymous and nonsynonymous polymorphism is variable between different amino acids. However, we do not observe a reduction in codon bias or the strength of selection in regions of suppressed recombination as expected. Instead, we find that the interaction between weak stabilizing selection and mutational bias likely plays a role in shaping the composition of synonymous codons across the third chromosome in D. pseudoobscura .

Description: Shigella spp. use a repertoire of virulence plasmid-encoded factors to cause shigellosis. These include components of a Type III Secretion Apparatus (T3SA) that is required for invasion of epithelial cells and many genes of unknown function. We constructed an array of 99 deletion mutants comprising all genes encoded by the virulence plasmid (excluding those known to be required for plasmid maintenance) of Shigella flexneri . We screened these mutants for their ability to bind the dye Congo red: an indicator of T3SA function. This screen focused our attention on an operon encoding genes that modify the cell envelope including virK , a gene of partially characterized function. We discovered that virK is required for controlled release of proteins to the culture supernatant. Mutations in virK result in a temperature-dependent overproduction of outer membrane vesicles (OMVs). The periplasmic chaperone/protease DegP, a known regulator of OMV production in Escherichia coli (encoded by a chromosomal gene), was found to similarly control OMV production in S. flexneri . Both virK and degP show genetic interactions with mxiD , a structural component of the T3SA. Our results are consistent with a model in which VirK and DegP relieve the periplasmic stress that accompanies assembly of the T3SA.

Description: Accurate phylogenies are critical to taxonomy as well as studies of speciation processes and other evolutionary patterns. Accurate branch lengths in phylogenies are critical for dating and rate measurements. Such accuracy may be jeopardized by unacknowledged sequencing error. We use simulated data to test a correction for DNA sequencing error in maximum likelihood phylogeny inference. Over a wide range of data polymorphism and true error rate, we found that correcting for sequencing error improves recovery of the branch lengths, even if the assumed error rate is up to twice the true error rate. Low error rates have little effect on recovery of the topology. When error is high, correction improves topological inference; however, when error is extremely high, using an assumed error rate greater than the true error rate leads to poor recovery of both topology and branch lengths. The error correction approach tested here was proposed in 2004 but has not been widely used, perhaps because researchers do not want to commit to an estimate of the error rate. This study shows that correction with an approximate error rate is generally preferable to ignoring the issue.

Description: In wild-type Caenorhabditis elegans , the synapse from motor neuron M4 to pharyngeal terminal bulb (TB) muscles is silent, and the muscles are instead excited by gap junction connections from adjacent muscles. An eat-5 innexin mutant lacking this electrical connection has few TB contractions and is unable to grow well on certain foods. We showed previously that this defect can be overcome by activation of the M4 -〉 TB synapse. To identify genes that negatively regulate synaptic transmission, we isolated new suppressors of eat-5 . To our surprise, these suppressors included null mutations in NPQR-type calcium channel subunit genes unc-2 and unc-36 . Our results are consistent with the hypothesis that Ca 2+ entry through the NPQR-type channel inhibits synaptic transmission by activating the calcium-activated K + channel SLO-1 , thus antagonizing the EGL-19 L-type calcium channel.

Description: Outliers often pose problems in analyses of data in plant breeding, but their influence on the performance of methods for estimating predictive accuracy in genomic prediction studies has not yet been evaluated. Here, we evaluate the influence of outliers on the performance of methods for accuracy estimation in genomic prediction studies using simulation. We simulated 1000 datasets for each of 10 scenarios to evaluate the influence of outliers on the performance of seven methods for estimating accuracy. These scenarios are defined by the number of genotypes, marker effect variance, and magnitude of outliers. To mimic outliers, we added to one observation in each simulated dataset, in turn, 5-, 8-, and 10-times the error SD used to simulate small and large phenotypic datasets. The effect of outliers on accuracy estimation was evaluated by comparing deviations in the estimated and true accuracies for datasets with and without outliers. Outliers adversely influenced accuracy estimation, more so at small values of genetic variance or number of genotypes. A method for estimating heritability and predictive accuracy in plant breeding and another used to estimate accuracy in animal breeding were the most accurate and resistant to outliers across all scenarios and are therefore preferable for accuracy estimation in genomic prediction studies. The performances of the other five methods that use cross-validation were less consistent and varied widely across scenarios. The computing time for the methods increased as the size of outliers and sample size increased and the genetic variance decreased.

Description: Inbred mice exhibit strain-specific variation in susceptibility to atherosclerosis and dyslipidemia that renders them useful in dissecting the genetic architecture of these complex diseases. Traditional quantitative trait locus (QTL) mapping studies using inbred strains often identify large genomic regions, containing many genes, due to limited recombination and/or sample size. This hampers candidate gene identification and translation of these results into possible risk factors and therapeutic targets. An alternative approach is the use of multiparental outbred lines for genetic mapping, such as the Diversity Outbred (DO) mouse panel, which can be more informative than traditional two-parent crosses and can aid in the identification of causal genes and variants associated with QTL. We fed 292 female DO mice either a high-fat, cholesterol-containing (HFCA) diet, to induce atherosclerosis, or a low-fat, high-protein diet for 18 wk and measured plasma lipid levels before and after diet treatment. We measured markers of atherosclerosis in the mice fed the HFCA diet. The mice were genotyped on a medium-density single-nucleotide polymorphism array and founder haplotypes were reconstructed using a hidden Markov model. The reconstructed haplotypes were then used to perform linkage mapping of atherosclerotic lesion size as well as plasma total cholesterol, triglycerides, insulin, and glucose. Among our highly significant QTL we detected a ~100 kb QTL interval for atherosclerosis on Chromosome 6, as well as a 1.4 Mb QTL interval on Chromosome 9 for triglyceride levels at baseline and a coincident 22.2 Mb QTL interval on Chromosome 9 for total cholesterol after dietary treatment. One candidate gene within the Chromosome 6 peak region associated with atherosclerosis is Apobec1 , the apolipoprotein B (ApoB) mRNA-editing enzyme, which plays a role in the regulation of ApoB, a critical component of low-density lipoprotein, by editing ApoB mRNA. This study demonstrates the value of the DO population to improve mapping resolution and to aid in the identification of potential therapeutic targets for cardiovascular disease. Using a DO mouse population fed an HFCA diet, we were able to identify an A/J-specific isoform of Apobec1 that contributes to atherosclerosis.