SAR11 marine bacteria require exogenous reduced sulphur for growth.

Sulphur is a universally required cell nutrient found in two amino acids and other small organic molecules. All aerobic marine bacteria are known to use assimilatory sulphate reduction to supply sulphur for biosynthesis, although many can assimilate sulphur from organic compounds that contain reduced sulphur atoms. An analysis of three complete 'Candidatus Pelagibacter ubique' genomes, and public ocean metagenomic data sets, suggested that members of the ubiquitous and abundant SAR11 alphaproteobacterial clade are deficient in assimilatory sulphate reduction genes. Here we show that SAR11 requires exogenous sources of reduced sulphur, such as methionine or 3-dimethylsulphoniopropionate (DMSP) for growth. Titrations of the algal osmolyte DMSP in seawater medium containing all other macronutrients in excess showed that 1.5 x 10(8) SAR11 cells are produced per nanomole of DMSP. Although it has been shown that other marine alphaproteobacteria use sulphur from DMSP in preference to sulphate, our results indicate that 'Cand. P. ubique' relies exclusively on reduced sulphur compounds that originate from other plankton.

Modulation of neural gene networks by estradiol in old rhesus macaque females.

The postmenopausal decrease in circulating estradiol (E2) levels has been shown to contribute to several adverse physiological and psychiatric effects. To elucidate the molecular effects of E2 on the brain, we examined differential gene expression and DNA methylation (DNAm) patterns in the nonhuman primate brain following ovariectomy (Ov) and subsequent subcutaneous bioidentical E2 chronic treatment. We identified several dysregulated molecular networks, including MAPK signaling and dopaminergic synapse response, that are associated with ovariectomy and shared across two different brain areas, the occipital cortex (OC) and prefrontal cortex (PFC). The finding that hypomethylation (p = 1.6 × 10-51) and upregulation (p = 3.8 × 10-3) of UBE2M across both brain regions provide strong evidence for molecular differences in the brain induced by E2 depletion. Additionally, differential expression (p = 1.9 × 10-4; interaction p = 3.5 × 10-2) of LTBR in the PFC provides further support for the role E2 plays in the brain, by demonstrating that the regulation of some genes that are altered by ovariectomy may also be modulated by Ov followed by hormone replacement therapy (HRT). These results present real opportunities to understand the specific biological mechanisms that are altered with depleted E2. Given E2's potential role in cognitive decline and neuroinflammation, our findings could lead to the discovery of novel therapeutics to slow cognitive decline. Together, this work represents a major step toward understanding molecular changes in the brain that are caused by ovariectomy and how E2 treatment may revert or protect against the negative neuro-related consequences caused by a depletion in estrogen as women approach menopause.

Selective sweeps and parallel mutation in the adaptive recovery from deleterious mutation in Caenorhabditis elegans.

Deleterious mutation poses a serious threat to human health and the persistence of small populations. Although adaptive recovery from deleterious mutation has been well-characterized in prokaryotes, the evolutionary mechanisms by which multicellular eukaryotes recover from deleterious mutation remain unknown. We applied high-throughput DNA sequencing to characterize genomic divergence patterns associated with the adaptive recovery from deleterious mutation using a Caenorhabditis elegans recovery-line system. The C. elegans recovery lines were initiated from a low-fitness mutation-accumulation (MA) line progenitor and allowed to independently evolve in large populations (N ∼ 1000) for 60 generations. All lines rapidly regained levels of fitness similar to the wild-type (N2) MA line progenitor. Although there was a near-zero probability of a single mutation fixing due to genetic drift during the recovery experiment, we observed 28 fixed mutations. Cross-generational analysis showed that all mutations went from undetectable population-level frequencies to a fixed state in 10-20 generations. Many recovery-line mutations fixed at identical timepoints, suggesting that the mutations, if not beneficial, hitchhiked to fixation during selective sweep events observed in the recovery lines. No MA line mutation reversions were detected. Parallel mutation fixation was observed for two sites in two independent recovery lines. Analysis using a C. elegans interactome map revealed many predicted interactions between genes with recovery line-specific mutations and genes with previously accumulated MA line mutations. Our study suggests that recovery-line mutations identified in both coding and noncoding genomic regions might have beneficial effects associated with compensatory epistatic interactions.

Methylome reorganization during in vitro dedifferentiation and regeneration of Populus trichocarpa.

BACKGROUND: Cytosine DNA methylation (5mC) is an epigenetic modification that is important to genome stability and regulation of gene expression. Perturbations of 5mC have been implicated as a cause of phenotypic variation among plants regenerated through in vitro culture systems. However, the pattern of change in 5mC and its functional role with respect to gene expression, are poorly understood at the genome scale. A fuller understanding of how 5mC changes during in vitro manipulation may aid the development of methods for reducing or amplifying the mutagenic and epigenetic effects of in vitro culture and plant transformation. RESULTS: We investigated the in vitro methylome of the model tree species Populus trichocarpa in a system that mimics routine methods for regeneration and plant transformation in the genus Populus (poplar). Using methylated DNA immunoprecipitation followed by high-throughput sequencing (MeDIP-seq), we compared the methylomes of internode stem segments from micropropagated explants, dedifferentiated calli, and internodes from regenerated plants. We found that more than half (56%) of the methylated portion of the genome appeared to be differentially methylated among the three tissue types. Surprisingly, gene promoter methylation varied little among tissues, however, the percentage of body-methylated genes increased from 9% to 14% between explants and callus tissue, then decreased to 8% in regenerated internodes. Forty-five percent of differentially-methylated genes underwent transient methylation, becoming methylated in calli, and demethylated in regenerants. These genes were more frequent in chromosomal regions with higher gene density. Comparisons with an expression microarray dataset showed that genes methylated at both promoters and gene bodies had lower expression than genes that were unmethylated or only promoter-methylated in all three tissues. Four types of abundant transposable elements showed their highest levels of 5mC in regenerated internodes. CONCLUSIONS: DNA methylation varies in a highly gene- and chromosome-differential manner during in vitro differentiation and regeneration. 5mC in redifferentiated tissues was not reset to that in original explants during the study period. Hypermethylation of gene bodies in dedifferentiated cells did not interfere with transcription, and may serve a protective role against activation of abundant transposable elements.

Cross-species epigenetic regulation of nucleus accumbens KCNN3 transcripts by excessive ethanol drinking.

The underlying genetic and epigenetic mechanisms driving functional adaptations in neuronal excitability and excessive alcohol intake are poorly understood. Small-conductance Ca2+-activated K+ (KCa2 or SK) channels encoded by the KCNN family of genes have emerged from preclinical studies as a key contributor to alcohol-induced functional neuroadaptations in alcohol-drinking monkeys and alcohol dependent mice. Here, this cross-species analysis focused on KCNN3 DNA methylation, gene expression, and single nucleotide polymorphisms including alternative promoters in KCNN3 that could influence surface trafficking and function of KCa2 channels. Bisulfite sequencing analysis of the nucleus accumbens tissue from alcohol-drinking monkeys and alcohol dependent mice revealed a differentially methylated region in exon 1A of KCNN3 that overlaps with a predicted promoter sequence. The hypermethylation of KCNN3 in the accumbens paralleled an increase in expression of alternative transcripts that encode apamin-insensitive and dominant-negative KCa2 channel isoforms. A polymorphic repeat in macaque KCNN3 encoded by exon 1 did not correlate with alcohol drinking. At the protein level, KCa2.3 channel expression in the accumbens was significantly reduced in very heavy drinking monkeys. Together, our cross-species findings on epigenetic dysregulation of KCNN3 represent a complex mechanism that utilizes alternative promoters to impact firing of accumbens neurons. Thus, these results provide support for hypermethylation of KCNN3 as a possible key molecular mechanism underlying harmful alcohol intake and alcohol use disorder.

Cross-species epigenetic regulation of nucleus accumbens KCNN3 transcripts by excessive ethanol drinking.

The underlying genetic and epigenetic mechanisms driving functional adaptations in neuronal excitability and excessive alcohol intake are poorly understood. Small-conductance Ca2+-activated K+ (KCa2 or SK) channels encoded by the KCNN family of genes have emerged from preclinical studies as a key contributor to alcohol-induced functional neuroadaptations in alcohol-drinking monkeys and alcohol-dependent mice. Here, this cross-species analysis focused on KCNN3 DNA methylation, gene expression, and single nucleotide polymorphisms, including alternative promoters in KCNN3, that could influence surface trafficking and function of KCa2 channels. Bisulfite sequencing analysis of the nucleus accumbens tissue from alcohol-drinking monkeys and alcohol-dependent mice revealed a differentially methylated region in exon 1A of KCNN3 that overlaps with a predicted promoter sequence. The hypermethylation of KCNN3 in the accumbens paralleled an increase in the expression of alternative transcripts that encode apamin-insensitive and dominant-negative KCa2 channel isoforms. A polymorphic repeat in macaque KCNN3 encoded by exon 1 did not correlate with alcohol drinking. At the protein level, KCa2.3 channel expression in the accumbens was significantly reduced in very heavy-drinking monkeys. Together, our cross-species findings on epigenetic dysregulation of KCNN3 represent a complex mechanism that utilizes alternative promoters to potentially impact the firing of accumbens neurons. Thus, these results provide support for hypermethylation of KCNN3 as a possible key molecular mechanism underlying harmful alcohol intake and alcohol use disorder.

Modulation of neural gene networks by estradiol in old rhesus macaque females.

The postmenopausal decrease in circulating estradiol (E2) levels has been shown to contribute to several adverse physiological and psychiatric effects. To elucidate the molecular effects of E2 on the brain, we examined differential gene expression and DNA methylation (DNAm) patterns in the nonhuman primate brain following ovariectomy (Ov) and subsequent E2 treatment. We identified several dysregulated molecular networks, including MAPK signaling and dopaminergic synapse response, that are associated with ovariectomy and shared across two different brain areas, the occipital cortex (OC) and prefrontal cortex (PFC). The finding that hypomethylation (p=1.6×10-51) and upregulation (p=3.8×10-3) of UBE2M across both brain regions, provide strong evidence for molecular differences in the brain induced by E2 depletion. Additionally, differential expression (p=1.9×10-4; interaction p=3.5×10-2) of LTBR in the PFC, provides further support for the role E2 plays in the brain, by demonstrating that the regulation of some genes that are altered by ovariectomy may also be modulated by Ov followed by hormone replacement therapy (HRT). These results present real opportunities to understand the specific biological mechanisms that are altered with depleted E2. Given E2's potential role in cognitive decline and neuroinflammation, our findings could lead to the discovery of novel therapeutics to slow cognitive decline. Together, this work represents a major step towards understanding molecular changes in the brain that are caused by ovariectomy and how E2 treatment may revert or protect against the negative neuro-related consequences caused by a depletion in estrogen as women approach menopause.

Methylation analysis by targeted bisulfite sequencing in large for gestational age (LGA) newborns: the LARGAN cohort.

BACKGROUND: In 1990, David Barker proposed that prenatal nutrition is directly linked to adult cardiovascular disease. Since then, the relationship between adult cardiovascular risk, metabolic syndrome and birth weight has been widely documented. Here, we used the TruSeq Methyl Capture EPIC platform to compare the methylation patterns in cord blood from large for gestational age (LGA) vs adequate for gestational age (AGA) newborns from the LARGAN cohort. RESULTS: We found 1672 differentially methylated CpGs (DMCs) with a nominal p < 0.05 and 48 differentially methylated regions (DMRs) with a corrected p < 0.05 between the LGA and AGA groups. A systems biology approach identified several biological processes significantly enriched with genes in association with DMCs with FDR < 0.05, including regulation of transcription, regulation of epinephrine secretion, norepinephrine biosynthesis, receptor transactivation, forebrain regionalization and several terms related to kidney and cardiovascular development. Gene ontology analysis of the genes in association with the 48 DMRs identified several significantly enriched biological processes related to kidney development, including mesonephric duct development and nephron tubule development. Furthermore, our dataset identified several DNA methylation markers enriched in gene networks involved in biological pathways and rare diseases of the cardiovascular system, kidneys, and metabolism. CONCLUSIONS: Our study identified several DMCs/DMRs in association with fetal overgrowth. The use of cord blood as a material for the identification of DNA methylation biomarkers gives us the possibility to perform follow-up studies on the same patients as they grow. These studies will not only help us understand how the methylome responds to continuum postnatal growth but also link early alterations of the DNA methylome with later clinical markers of growth and metabolic fitness.

Dynamic DNA cytosine methylation in the Populus trichocarpa genome: tissue-level variation and relationship to gene expression.

BACKGROUND: DNA cytosine methylation is an epigenetic modification that has been implicated in many biological processes. However, large-scale epigenomic studies have been applied to very few plant species, and variability in methylation among specialized tissues and its relationship to gene expression is poorly understood. RESULTS: We surveyed DNA methylation from seven distinct tissue types (vegetative bud, male inflorescence [catkin], female catkin, leaf, root, xylem, phloem) in the reference tree species black cottonwood (Populus trichocarpa). Using 5-methyl-cytosine DNA immunoprecipitation followed by Illumina sequencing (MeDIP-seq), we mapped a total of 129,360,151 36- or 32-mer reads to the P. trichocarpa reference genome. We validated MeDIP-seq results by bisulfite sequencing, and compared methylation and gene expression using published microarray data. Qualitative DNA methylation differences among tissues were obvious on a chromosome scale. Methylated genes had lower expression than unmethylated genes, but genes with methylation in transcribed regions ("gene body methylation") had even lower expression than genes with promoter methylation. Promoter methylation was more frequent than gene body methylation in all tissues except male catkins. Male catkins differed in demethylation of particular transposable element categories, in level of gene body methylation, and in expression range of genes with methylated transcribed regions. Tissue-specific gene expression patterns were correlated with both gene body and promoter methylation. CONCLUSIONS: We found striking differences among tissues in methylation, which were apparent at the chromosomal scale and when genes and transposable elements were examined. In contrast to other studies in plants, gene body methylation had a more repressive effect on transcription than promoter methylation.

Genome resequencing reveals multiscale geographic structure and extensive linkage disequilibrium in the forest tree Populus trichocarpa.

• Plant population genomics informs evolutionary biology, breeding, conservation and bioenergy feedstock development. For example, the detection of reliable phenotype-genotype associations and molecular signatures of selection requires a detailed knowledge about genome-wide patterns of allele frequency variation, linkage disequilibrium and recombination. • We resequenced 16 genomes of the model tree Populus trichocarpa and genotyped 120 trees from 10 subpopulations using 29,213 single-nucleotide polymorphisms. • Significant geographic differentiation was present at multiple spatial scales, and range-wide latitudinal allele frequency gradients were strikingly common across the genome. The decay of linkage disequilibrium with physical distance was slower than expected from previous studies in Populus, with r(2) dropping below 0.2 within 3-6 kb. Consistent with this, estimates of recent effective population size from linkage disequilibrium (N(e) ≈ 4000-6000) were remarkably low relative to the large census sizes of P. trichocarpa stands. Fine-scale rates of recombination varied widely across the genome, but were largely predictable on the basis of DNA sequence and methylation features. • Our results suggest that genetic drift has played a significant role in the recent evolutionary history of P. trichocarpa. Most importantly, the extensive linkage disequilibrium detected suggests that genome-wide association studies and genomic selection in undomesticated populations may be more feasible in Populus than previously assumed.

A genome-wide view of Caenorhabditis elegans base-substitution mutation processes.

Knowledge of mutation processes is central to understanding virtually all evolutionary phenomena and the underlying nature of genetic disorders and cancers. However, the limitations of standard molecular mutation detection methods have historically precluded a genome-wide understanding of mutation rates and spectra in the nuclear genomes of multicellular organisms. We applied two high-throughput DNA sequencing technologies to identify and characterize hundreds of spontaneously arising base-substitution mutations in 10 Caenorhabditis elegans mutation-accumulation (MA)-line nuclear genomes. C. elegans mutation rate estimates were similar to previous calculations based on smaller numbers of mutations. Mutations were distributed uniformly within and among chromosomes and were not associated with recombination rate variation in the MA lines, suggesting that intragenomic variation in genetic hitchhiking and/or background selection are primarily responsible for the chromosomal distribution patterns of polymorphic nucleotides in C. elegans natural populations. A strong mutational bias from G/C to A/T nucleotides was detected in the MA lines, implicating oxidative DNA damage as a major endogenous mutagenic force in C. elegans. The observed mutational bias also suggests that the C. elegans nuclear genome cannot be at equilibrium because of mutation alone. Transversions dominate the spectrum of spontaneous mutations observed here, whereas transitions dominate patterns of allegedly neutral polymorphism in natural populations of C. elegans and many other animal species; this observation challenges the assumption that natural patterns of molecular variation in noncoding regions of the nuclear genome accurately reflect underlying mutation processes.

Genome-wide analysis of the nucleus accumbens identifies DNA methylation signals differentiating low/binge from heavy alcohol drinking.

Alcohol-use disorders encompass a range of drinking levels and behaviors, including low, binge, and heavy drinking. In this regard, investigating the neural state of individuals who chronically self-administer lower doses of alcohol may provide insight into mechanisms that prevent the escalation of alcohol use. DNA methylation is one of the epigenetic mechanisms that stabilizes adaptations in gene expression and has been associated with alcohol use. Thus, we investigated DNA methylation, gene expression, and the predicted neural effects in the nucleus accumbens core (NAcc) of male rhesus macaques categorized as "low" or "binge" drinkers, compared to "alcohol-naïve" and "heavy" drinkers based on drinking patterns during a 12-month alcohol self-administration protocol. Using genome-wide CpG-rich region enrichment and bisulfite sequencing, the methylation levels of 2.6 million CpGs were compared between alcohol-naïve (AN), low/binge (L/BD), and heavy/very heavy (H/VHD) drinking subjects (n = 24). Through regional clustering analysis, we identified nine significant differential methylation regions (DMRs) that specifically distinguished ANs and L/BDs, and then compared those DMRs among H/VHDs. The DMRs mapped to genes encoding ion channels, receptors, cell adhesion molecules, and cAMP, NF-κß and Wnt signaling pathway proteins. Two of the DMRs, linked to PDE10A and PKD2L2, were also differentially methylated in H/VHDs, suggesting an alcohol-dose independent effect. However, two other DMRs, linked to the CCBE1 and FZD5 genes, had L/BD methylation levels that significantly differed from both ANs and H/VHDs. The remaining five DMRs also differentiated L/BDs and ANs. However, H/VHDs methylation levels were not distinguishable from either of the two groups. Functional validation of two DMRs, linked to FZD5 and PDE10A, support their role in regulating gene expression and exon usage, respectively. In summary, the findings demonstrate that L/BD is associated with unique DNA methylation signatures in the primate NAcc, and that the methylation signatures identify synaptic genes that may play a role in preventing the escalation of alcohol use.

Nicotinic receptors in non-human primates: Analysis of genetic and functional conservation with humans.

Nicotinic acetylcholine receptors (nAChRs) are highly conserved between humans and non-human primates. Conservation exists at the level of genomic structure, protein structure and epigenetics. Overall homology of nAChRs at the protein level is 98% in macaques versus 89% in mice, which is highly relevant for evaluating subtype-specific ligands that have different affinities in humans versus rodents. In addition to conservation at the protein level, there is high conservation of genomic structure in terms of intron and exon size and placement of CpG sites that play a key role in epigenetic regulation. Analysis of single nucleotide polymorphisms (SNPs) shows that while the majority of SNPs are not conserved between humans and macaques, some functional polymorphisms are. Most significantly, cynomolgus monkeys express a similar α5 nAChR Asp398Asn polymorphism to the human α5 Asp398Asn polymorphism that has been linked to greater nicotine addiction and smoking related disease. Monkeys can be trained to readily self-administer nicotine, and in an initial study we have demonstrated that cynomolgus monkeys bearing the α5 D398N polymorphism show a reduced behavioral sensitivity to oral nicotine and tend to consume it in a different pattern when compared to wild-type monkeys. Thus the combination of highly homologous nAChR, higher cortical functions and capacity for complex training makes non-human primates a unique model to study in vivo functions of nicotinic receptors. In particular, primate studies on nicotine addiction and evaluation of therapies to prevent or overcome nicotine addiction are likely to be highly predictive of treatment outcomes in humans. This article is part of the Special Issue entitled 'The Nicotinic Acetylcholine Receptor: From Molecular Biology to Cognition'.

Conserved syntenic clusters of protein coding genes are missing in birds.

BACKGROUND: Birds are one of the most highly successful and diverse groups of vertebrates, having evolved a number of distinct characteristics, including feathers and wings, a sturdy lightweight skeleton and unique respiratory and urinary/excretion systems. However, the genetic basis of these traits is poorly understood. RESULTS: Using comparative genomics based on extensive searches of 60 avian genomes, we have found that birds lack approximately 274 protein coding genes that are present in the genomes of most vertebrate lineages and are for the most part organized in conserved syntenic clusters in non-avian sauropsids and in humans. These genes are located in regions associated with chromosomal rearrangements, and are largely present in crocodiles, suggesting that their loss occurred subsequent to the split of dinosaurs/birds from crocodilians. Many of these genes are associated with lethality in rodents, human genetic disorders, or biological functions targeting various tissues. Functional enrichment analysis combined with orthogroup analysis and paralog searches revealed enrichments that were shared by non-avian species, present only in birds, or shared between all species. CONCLUSIONS: Together these results provide a clearer definition of the genetic background of extant birds, extend the findings of previous studies on missing avian genes, and provide clues about molecular events that shaped avian evolution. They also have implications for fields that largely benefit from avian studies, including development, immune system, oncogenesis, and brain function and cognition. With regards to the missing genes, birds can be considered 'natural knockouts' that may become invaluable model organisms for several human diseases.

Using mitogenomic and nuclear ribosomal sequence data to investigate the phylogeny of the Xiphinema americanum species complex.

Nematodes within the Xiphinema americanum species complex are economically important because they vector nepoviruses which cause considerable damage to a variety of agricultural crops. The taxonomy of X. americanum species complex is controversial, with the number of putative species being the subject of debate. Accurate phylogenetic knowledge of this group is highly desirable as it may ultimately reveal genetic differences between species. For this study, nematodes belonging to the X. americanum species complex, including potentially mixed species populations, were collected from 12 geographically disparate locations across the U.S. from different crops and in varying association with nepoviruses. At least four individuals from each population were analyzed. A portion of the 18S nuclear ribosomal DNA (rDNA) gene was sequenced for all individuals while the internal transcribed spacer region 1 (ITS1) of rDNA was cloned and 2 to 6 clones per individual were sequenced. Mitochondrial genomes for numerous individuals were sequenced in parallel using high-throughput DNA sequencing (HTS) technology. Phylogenetic analysis of the 18S rDNA revealed virtually identical sequences across all populations. Analysis of ITS1 rDNA sequences revealed several well-supported clades, with some degree of congruence with geographic location and viral transmission, but also numerous presumably paralogous sequences that failed to form clades with other sequences from the same population. Analysis of mitochondrial DNA (mtDNA) indicated the presence of three distinct monophyletic clades of X. americanum species complex nematodes. Two clades contained nematodes found in association with nepovirus and the third contained divergent mtDNA sequences from three nematode populations from the western U.S. where nepovirus was absent. The inherent heterogeneity in ITS1 rDNA sequence data and lack of informative sites in 18S rDNA analysis suggests that mtDNA may be more useful in sorting out the taxonomic confusion of the X. americanum species complex.

Comparative genomics of a plant-pathogenic fungus, Pyrenophora tritici-repentis, reveals transduplication and the impact of repeat elements on pathogenicity and population divergence.

Pyrenophora tritici-repentis is a necrotrophic fungus causal to the disease tan spot of wheat, whose contribution to crop loss has increased significantly during the last few decades. Pathogenicity by this fungus is attributed to the production of host-selective toxins (HST), which are recognized by their host in a genotype-specific manner. To better understand the mechanisms that have led to the increase in disease incidence related to this pathogen, we sequenced the genomes of three P. tritici-repentis isolates. A pathogenic isolate that produces two known HSTs was used to assemble a reference nuclear genome of approximately 40 Mb composed of 11 chromosomes that encode 12,141 predicted genes. Comparison of the reference genome with those of a pathogenic isolate that produces a third HST, and a nonpathogenic isolate, showed the nonpathogen genome to be more diverged than those of the two pathogens. Examination of gene-coding regions has provided candidate pathogen-specific proteins and revealed gene families that may play a role in a necrotrophic lifestyle. Analysis of transposable elements suggests that their presence in the genome of pathogenic isolates contributes to the creation of novel genes, effector diversification, possible horizontal gene transfer events, identified copy number variation, and the first example of transduplication by DNA transposable elements in fungi. Overall, comparative analysis of these genomes provides evidence that pathogenicity in this species arose through an influx of transposable elements, which created a genetically flexible landscape that can easily respond to environmental changes.

Variation in base-substitution mutation in experimental and natural lineages of Caenorhabditis nematodes.

Variation among lineages in the mutation process has the potential to impact diverse biological processes ranging from susceptibilities to genetic disease to the mode and tempo of molecular evolution. The combination of high-throughput DNA sequencing (HTS) with mutation-accumulation (MA) experiments has provided a powerful approach to genome-wide mutation analysis, though insights into mutational variation have been limited by the vast evolutionary distances among the few species analyzed. We performed a HTS analysis of MA lines derived from four Caenorhabditis nematode natural genotypes: C. elegans N2 and PB306 and C. briggsae HK104 and PB800. Total mutation rates did not differ among the four sets of MA lines. A mutational bias toward G:C→A:T transitions and G:C→T:A transversions was observed in all four sets of MA lines. Chromosome-specific rates were mostly stable, though there was some evidence for a slightly elevated X chromosome mutation rate in PB306. Rates were homogeneous among functional coding sequence types and across autosomal cores, arms, and tips. Mutation spectra were similar among the four MA line sets but differed significantly when compared with patterns of natural base-substitution polymorphism for 13/14 comparisons performed. Our findings show that base-substitution mutation processes in these closely related animal lineages are mostly stable but differ from natural polymorphism patterns in these two species.

Host-selective toxins of Pyrenophora tritici-repentis induce common responses associated with host susceptibility.

Pyrenophora tritici-repentis (Ptr), a necrotrophic fungus and the causal agent of tan spot of wheat, produces one or a combination of host-selective toxins (HSTs) necessary for disease development. The two most studied toxins produced by Ptr, Ptr ToxA (ToxA) and Ptr ToxB (ToxB), are proteins that cause necrotic or chlorotic symptoms respectively. Investigation of host responses induced by HSTs provides better insight into the nature of the host susceptibility. Microarray analysis of ToxA has provided evidence that it can elicit responses similar to those associated with defense. In order to evaluate whether there are consistent host responses associated with susceptibility, a similar analysis of ToxB-induced changes in the same sensitive cultivar was conducted. Comparative analysis of ToxA- and ToxB-induced transcriptional changes showed that similar groups of genes encoding WRKY transcription factors, RLKs, PRs, components of the phenylpropanoid and jasmonic acid pathways are activated. ROS accumulation and photosystem dysfunction proved to be common mechanism-of-action for these toxins. Despite similarities in defense responses, transcriptional and biochemical responses as well as symptom development occur more rapidly for ToxA compared to ToxB, which could be explained by differences in perception as well as by differences in activation of a specific process, for example, ethylene biosynthesis in ToxA treatment. Results of this study suggest that perception of HSTs will result in activation of defense responses as part of a susceptible interaction and further supports the hypothesis that necrotrophic fungi exploit defense responses in order to induce cell death.

Development and evaluation of a genome-wide 6K SNP array for diploid sweet cherry and tetraploid sour cherry.

High-throughput genome scans are important tools for genetic studies and breeding applications. Here, a 6K SNP array for use with the Illumina Infinium® system was developed for diploid sweet cherry (Prunus avium) and allotetraploid sour cherry (P. cerasus). This effort was led by RosBREED, a community initiative to enable marker-assisted breeding for rosaceous crops. Next-generation sequencing in diverse breeding germplasm provided 25 billion basepairs (Gb) of cherry DNA sequence from which were identified genome-wide SNPs for sweet cherry and for the two sour cherry subgenomes derived from sweet cherry (avium subgenome) and P. fruticosa (fruticosa subgenome). Anchoring to the peach genome sequence, recently released by the International Peach Genome Initiative, predicted relative physical locations of the 1.9 million putative SNPs detected, preliminarily filtered to 368,943 SNPs. Further filtering was guided by results of a 144-SNP subset examined with the Illumina GoldenGate® assay on 160 accessions. A 6K Infinium® II array was designed with SNPs evenly spaced genetically across the sweet and sour cherry genomes. SNPs were developed for each sour cherry subgenome by using minor allele frequency in the sour cherry detection panel to enrich for subgenome-specific SNPs followed by targeting to either subgenome according to alleles observed in sweet cherry. The array was evaluated using panels of sweet (n = 269) and sour (n = 330) cherry breeding germplasm. Approximately one third of array SNPs were informative for each crop. A total of 1825 polymorphic SNPs were verified in sweet cherry, 13% of these originally developed for sour cherry. Allele dosage was resolved for 2058 polymorphic SNPs in sour cherry, one third of these being originally developed for sweet cherry. This publicly available genomics resource represents a significant advance in cherry genome-scanning capability that will accelerate marker-locus-trait association discovery, genome structure investigation, and genetic diversity assessment in this diploid-tetraploid crop group.

Genome-wide SNP detection, validation, and development of an 8K SNP array for apple.

As high-throughput genetic marker screening systems are essential for a range of genetics studies and plant breeding applications, the International RosBREED SNP Consortium (IRSC) has utilized the Illumina Infinium® II system to develop a medium- to high-throughput SNP screening tool for genome-wide evaluation of allelic variation in apple (Malus×domestica) breeding germplasm. For genome-wide SNP discovery, 27 apple cultivars were chosen to represent worldwide breeding germplasm and re-sequenced at low coverage with the Illumina Genome Analyzer II. Following alignment of these sequences to the whole genome sequence of 'Golden Delicious', SNPs were identified using SoapSNP. A total of 2,113,120 SNPs were detected, corresponding to one SNP to every 288 bp of the genome. The Illumina GoldenGate® assay was then used to validate a subset of 144 SNPs with a range of characteristics, using a set of 160 apple accessions. This validation assay enabled fine-tuning of the final subset of SNPs for the Illumina Infinium® II system. The set of stringent filtering criteria developed allowed choice of a set of SNPs that not only exhibited an even distribution across the apple genome and a range of minor allele frequencies to ensure utility across germplasm, but also were located in putative exonic regions to maximize genotyping success rate. A total of 7867 apple SNPs was established for the IRSC apple 8K SNP array v1, of which 5554 were polymorphic after evaluation in segregating families and a germplasm collection. This publicly available genomics resource will provide an unprecedented resolution of SNP haplotypes, which will enable marker-locus-trait association discovery, description of the genetic architecture of quantitative traits, investigation of genetic variation (neutral and functional), and genomic selection in apple.