Next‐generation sequencing (NGS)‐based identification of induced mutations in a doubly mutagenized tomato (Solanum lycopersicum) population

The identification of mutations in targeted genes has been significantly simplified by the advent of TILLING (Targeting Induced Local Lesions In Genomes), speeding up the functional genomic analysis of animals and plants. Next‐generation sequencing (NGS) is gradually replacing classical TILLING for mutation detection, as it allows the analysis of a large number of amplicons in short durations. The NGS approach was used to identify mutations in a population of Solanum lycopersicum (tomato) that was doubly mutagenized by ethylmethane sulphonate (EMS). Twenty‐five genes belonging to carotenoids and folate metabolism were PCR‐amplified and screened to identify potentially beneficial alleles. To augment efficiency, the 600‐bp amplicons were directly sequenced in a non‐overlapping manner in Illumina MiSeq, obviating the need for a fragmentation step before library preparation. A comparison of the different pooling depths revealed that heterozygous mutations could be identified up to 128‐fold pooling. An evaluation of six different software programs (camba , crisp , gatk unified genotyper , lofreq , snver and vipr ) revealed that no software program was robust enough to predict mutations with high fidelity. Among these, crisp and camba predicted mutations with lower false discovery rates. The false positives were largely eliminated by considering only mutations commonly predicted by two different software programs. The screening of 23.47 Mb of tomato genome yielded 75 predicted mutations, 64 of which were confirmed by Sanger sequencing with an average mutation density of 1/367 Kb. Our results indicate that NGS combined with multiple variant detection tools can reduce false positives and significantly speed up the mutation discovery rate.     

The within‐host dynamics of infection in trans‐generationally primed flour beetles

Many taxa exhibit plastic immune responses initiated after primary microbial exposure that provide increased protection against disease‐induced mortality and the fitness costs of infection. In several arthropod species, this protection can even be passed from parents to offspring through a phenomenon called trans‐generational immune priming. Here, we first demonstrate that trans‐generational priming is a repeatable phenomenon in flour beetles (Tribolium castaneum) primed and infected with Bacillus thuringiensis (Bt). We then quantify the within‐host dynamics of microbes and host physiological responses in infected offspring from primed and unprimed mothers by monitoring bacterial density and using mRNA‐seq to profile host gene expression, respectively, over the acute infection period. We find that priming increases inducible resistance against Bt around a critical temporal juncture where host septicaemic trajectories, and consequently survival, may be determined in unprimed individuals. Our results identify a highly differentially expressed biomarker of priming, containing an EIF4‐e domain, in uninfected individuals, as well as several other candidate genes. Moreover, the induction and decay dynamics of gene expression over time suggest a metabolic shift in primed individuals. The identified bacterial and gene expression dynamics are likely to influence patterns of bacterial fitness and disease transmission in natural populations.     

A modified universal fast walking method for single-tube transposon mapping

An enhanced universal fast walking (UFW) method adapted for the mapping of transposons is described. This protocol combines the original UFW method with the use of agarase to unravel composite nucleotide sequence, thereby forgoing molecular cloning steps and the use of restriction enzymes and ligases necessary in other available genome walking methods such as the prominent inverse PCR. The minuscule automatable chemistry of UFW is completed within one reaction vessel using a constant enzyme buffer, and the intrinsic DNA fingerprints, from which amplicons may be quantitatively recovered, offer quality assurance. The core steps of the protocol, spanning half a day or less, comprise first-strand synthesis, primer destruction, random-ended-primer annealing, distal branched-end repair, second-primer destruction, lariat formation and final amplification. Distinctively, no starting or intermediate templates are wasted during the reaction series, thus achieving yields comparable to direct PCR. Ultimate per-reaction walk-lengths are schematically illimitable and sequence-ready amplicons can be produced immediately from prevalent single-copy genomic walk origins. The core UFW protocol may be applied, as described here, to expedited transposon boundary retrieval, but is also applicable to general genome walking and cDNA walking, as well as viral and other insertional element mapping.     

Tag jumps illuminated – reducing sequence‐to‐sample misidentifications in metabarcoding studies

Metabarcoding of environmental samples on second‐generation sequencing platforms has rapidly become a valuable tool for ecological studies. A fundamental assumption of this approach is the reliance on being able to track tagged amplicons back to the samples from which they originated. In this study, we address the problem of sequences in metabarcoding sequencing outputs with false combinations of used tags (tag jumps). Unless these sequences can be identified and excluded from downstream analyses, tag jumps creating sequences with false, but already used tag combinations, can cause incorrect assignment of sequences to samples and artificially inflate diversity. In this study, we document and investigate tag jumping in metabarcoding studies on Illumina sequencing platforms by amplifying mixed‐template extracts obtained from bat droppings and leech gut contents with tagged generic arthropod and mammal primers, respectively. We found that an average of 2.6% and 2.1% of sequences had tag combinations, which could be explained by tag jumping in the leech and bat diet study, respectively. We suggest that tag jumping can happen during blunt‐ending of pools of tagged amplicons during library build and as a consequence of chimera formation during bulk amplification of tagged amplicons during library index PCR. We argue that tag jumping and contamination between libraries represents a considerable challenge for Illumina‐based metabarcoding studies, and suggest measures to avoid false assignment of tag jumping‐derived sequences to samples.     

Hexanoic acid protects tomato plants against Botrytis cinerea by priming defence responses and reducing oxidative stress

Treatment with the resistance priming inducer hexanoic acid (Hx) protects tomato plants from Botrytis cinerea by activating defence responses. To investigate the molecular mechanisms underlying hexanoic acid‐induced resistance (Hx‐IR), we compared the expression profiles of three different conditions: Botrytis‐infected plants (Inf), Hx‐treated plants (Hx) and Hx‐treated + infected plants (Hx+Inf). The microarray analysis at 24 h post‐inoculation showed that Hx and Hx+Inf plants exhibited the differential expression and priming of many Botrytis‐induced genes. Interestingly, we found that the activation by Hx of other genes was not altered by the fungus at this time point. These genes may be considered to be specific targets of the Hx priming effect and may help to elucidate its mechanisms of action. It is noteworthy that, in Hx and Hx+Inf plants, there was up‐regulation of proteinase inhibitor genes, DNA‐binding factors, enzymes involved in plant hormone signalling and synthesis, and, remarkably, the genes involved in oxidative stress. Given the relevance of the oxidative burst occurring in plant–pathogen interactions, the effect of Hx on this process was studied in depth. We showed by specific staining that reactive oxygen species (ROS) accumulation in Hx+Inf plants was reduced and more restricted around infection sites. In addition, these plants showed higher ratios of reduced to oxidized glutathione and ascorbate, and normal levels of antioxidant activities. The results obtained indicate that Hx protects tomato plants from B. cinerea by regulating and priming Botrytis‐specific and non‐specific genes, preventing the harmful effects of oxidative stress produced by infection.     

Phylogenetic and proteomic analysis of an anaerobic toluene-degrading community

Aims: This study intended to unravel the physiological interplay in an anaerobic microbial community that degrades toluene under sulfate‐reducing conditions combining proteomic and genetic techniques. Methods and Results: An enriched toluene‐degrading community (Zz5‐7) growing in batch cultures was investigated by DNA‐ and protein‐based analyses. The affiliation and diversity of the community were analysed using 16S ribosomal RNA (rRNA) genes as a phylogenetic marker as well as bssA and dsrAB genes as functional markers. Metaproteome analysis was carried out by a global protein extraction and a subsequent protein separation by two‐dimensional gel electrophoresis (2‐DE). About 85% of the proteins in the spots were identified by nano‐liquid chromatography coupled with electrospray mass spectrometry (nano‐LC–ESI‐MS/MS) analysis. DNA sequencing of bssA and the most abundant dsrAB amplicons revealed high similarities to a member of the Desulfobulbaceae, which was also predominant according to 16S rRNA gene amplicons. Metaproteome analysis provided 202 unambiguous protein identifications derived from 236 unique protein spots. The proteins involved in anaerobic toluene activation, dissimilatory sulfate reduction, hydrogen production/consumption and autotrophic carbon fixation were mainly affiliated to members of the Desulfobulbaceae and several other Deltaproteobacteria. Conclusion: Phylogenetic and metaproteomic analyses revealed a member of the Desulfobulbaceae as the key player of anaerobic toluene degradation in a sulfate‐reducing consortium. Significance and Impact of the Study: This is the first study that combines genetic and proteomic analyses to indicate the interactions in an anaerobic toluene‐degrading microbial consortium.     

Divergent immune priming responses across flour beetle life stages and populations

Growing evidence shows that low doses of pathogens may prime the immune response in many insects, conferring subsequent protection against infection in the same developmental stage (within‐life stage priming), across life stages (ontogenic priming), or to offspring (transgenerational priming). Recent work also suggests that immune priming is a costly response. Thus, depending on host and pathogen ecology and evolutionary history, tradeoffs with other fitness components may constrain the evolution of priming. However, the relative impacts of priming at different life stages and across natural populations remain unknown. We quantified immune priming responses of 10 natural populations of the red flour beetle Tribolium castaneum, primed and infected with the natural insect pathogen Bacillus thuringiensis. We found that priming responses were highly variable both across life stages and populations, ranging from no detectable response to a 13‐fold survival benefit. Comparing across stages, we found that ontogenic immune priming at the larval stage conferred maximum protection against infection. Finally, we found that various forms of priming showed sex‐specific associations that may represent tradeoffs or shared mechanisms. These results indicate the importance of sex‐, life stage‐, and population‐specific selective pressures that can cause substantial divergence in priming responses even within a species. Our work highlights the necessity of further work to understand the mechanistic basis of this variability.     

Testing genotyping strategies for ultra‐deep sequencing of a co‐amplifying gene family: MHC class I in a passerine bird

Characterization of highly duplicated genes, such as genes of the major histocompatibility complex (MHC), where multiple loci often co‐amplify, has until recently been hindered by insufficient read depths per amplicon. Here, we used ultra‐deep Illumina sequencing to resolve genotypes at exon 3 of MHC class I genes in the sedge warbler (Acrocephalus schoenobaenus). We sequenced 24 individuals in two replicates and used this data, as well as a simulated data set, to test the effect of amplicon coverage (range: 500–20 000 reads per amplicon) on the repeatability of genotyping using four different genotyping approaches. A third replicate employed unique barcoding to assess the extent of tag jumping, that is swapping of individual tag identifiers, which may confound genotyping. The reliability of MHC genotyping increased with coverage and approached or exceeded 90% within‐method repeatability of allele calling at coverages of >5000 reads per amplicon. We found generally high agreement between genotyping methods, especially at high coverages. High reliability of the tested genotyping approaches was further supported by our analysis of the simulated data set, although the genotyping approach relying primarily on replication of variants in independent amplicons proved sensitive to repeatable errors. According to the most repeatable genotyping method, the number of co‐amplifying variants per individual ranged from 19 to 42. Tag jumping was detectable, but at such low frequencies that it did not affect the reliability of genotyping. We thus demonstrate that gene families with many co‐amplifying genes can be reliably genotyped using HTS, provided that there is sufficient per amplicon coverage.     

Single‐nucleotide polymorphism discovery in Leptographium longiclavatum,a mountain pine beetle‐associated symbiotic fungus,using whole‐genome resequencing

Single‐nucleotide polymorphisms (SNPs) are rapidly becoming the standard markers in population genomics studies; however, their use in nonmodel organisms is limited due to the lack of cost‐effective approaches to uncover genome‐wide variation, and the large number of individuals needed in the screening process to reduce ascertainment bias. To discover SNPs for population genomics studies in the fungal symbionts of the mountain pine beetle (MPB), we developed a road map to discover SNPs and to produce a genotyping platform. We undertook a whole‐genome sequencing approach of Leptographium longiclavatum in combination with available genomics resources of another MPB symbiont, Grosmannia clavigera. We sequenced 71 individuals pooled into four groups using the Illumina sequencing technology. We generated between 27 and 30 million reads of 75 bp that resulted in a total of 1, 181 contigs longer than 2 kb and an assembled genome size of 28.9 Mb (N₅₀ = 48 kb, average depth = 125x). A total of 9052 proteins were annotated, and between 9531 and 17 266 SNPs were identified in the four pools. A subset of 206 genes (containing 574 SNPs, 11% false positives) was used to develop a genotyping platform for this species. Using this roadmap, we developed a genotyping assay with a total of 147 SNPs located in 121 genes using the Illumina^® Sequenom iPLEX Gold. Our preliminary genotyping (success rate = 85%) of 304 individuals from 36 populations supports the utility of this approach for population genomics studies in other MPB fungal symbionts and other fungal nonmodel species.     

Microdevice‐based solid‐phase polymerase chain reaction for rapid detection of pathogenic microorganisms

We demonstrate the integration of DNA amplification and detection functionalities developed on a lab‐on‐a‐chip microdevice utilizing solid‐phase polymerase chain reaction (SP‐PCR) for point‐of‐need (PON) DNA analyses. First, the polycarbonate microdevice was fabricated by thermal bonding to contain microchambers as reservoirs for performing SP‐PCR. Next, the microchambers were subsequently modified with polyethyleneimine and glutaraldehyde for immobilizing amine‐modified forward primers. During SP‐PCR, the immobilized forward primers and freely diffusing fluorescence‐labeled reverse primers cooperated to generate target amplicons, which remained covalently attached to the microchambers for the fluorescence detection. The SP‐PCR microdevice was used for the direct identifications of two widely detected foodborne pathogens, namely Salmonella spp. and Staphylococcus aureus, and an alga causing harmful algal blooms annually in South Korea, Cochlodinium polykrikoides. The SP‐PCR microdevice would be versatilely applied in PON testing as a universal platform for the fast identification of foodborne pathogens and environmentally threatening biogenic targets.     

β‐Aminobutyric Acid Primed Expression of WRKY and Defence Genes in Brassica carinata Against Alternaria Blight

Foliar spray with BABA led to a significant reduction of lesion development in Brassica carinata caused by Alternaria brassicae. To get better insight into molecular mechanisms underlying priming of defence responses by BABA, expression pattern of BcWRKY genes and marker genes for the SA and JA pathway namely PR‐1 and PDF 1.2 was examined. Q‐RT‐PCR analysis revealed priming of BcWRKY70, BcWRKY11 and BcWRKY53 gene expression in BABA‐pretreated Brassica plants challenged with pathogen. However, the expression of BcWRKY72 and BcWRKY18 remained unchanged. Furthermore, BcWRKY7 gene was found to be upregulated in water‐treated plants in response to pathogen indicating its role in susceptibility. In addition, BABA application potentiated expression of defence genes PR‐1, PDF1.2 and PAL in response to the pathogen. In conclusion, BABA‐primed expression of BcWRKY70, BcWRKY11 and BcWRKY53 genes is strongly correlated with enhanced expression of PR‐1, PDF1.2 and PAL hence suggesting their role in BABA‐induced resistance.     

I kappa B kinase alpha (IKKα) activity is required for functional maturation of dendritic cells and acquired immunity to infection

Dendritic cells (DC) are required for priming antigen‐specific T cells and acquired immunity to many important human pathogens, including Mycobacteriuim tuberculosis (TB) and influenza. However, inappropriate priming of auto‐reactive T cells is linked with autoimmune disease. Understanding the molecular mechanisms that regulate the priming and activation of naïve T cells is critical for development of new improved vaccines and understanding the pathogenesis of autoimmune diseases. The serine/threonine kinase IKKα (CHUK) has previously been shown to have anti‐inflammatory activity and inhibit innate immunity. Here, we show that IKKα is required in DC for priming antigen‐specific T cells and acquired immunity to the human pathogen Listeria monocytogenes. We describe a new role for IKKα in regulation of IRF3 activity and the functional maturation of DC. This presents a unique role for IKKα in dampening inflammation while simultaneously promoting adaptive immunity that could have important implications for the development of new vaccine adjuvants and treatment of autoimmune diseases.     

Linear Mixed Model Selection for False Discovery Rate Control in Microarray Data Analysis

Summary In a microarray experiment, one experimental design is used to obtain expression measures for all genes. One popular analysis method involves fitting the same linear mixed model for each gene, obtaining gene‐specific p‐values for tests of interest involving fixed effects, and then choosing a threshold for significance that is intended to control false discovery rate (FDR) at a desired level. When one or more random factors have zero variance components for some genes, the standard practice of fitting the same full linear mixed model for all genes can result in failure to control FDR. We propose a new method that combines results from the fit of full and selected linear mixed models to identify differentially expressed genes and provide FDR control at target levels when the true underlying random effects structure varies across genes.     

Bacteria‐type‐specific biparental immune priming in the pipefish Syngnathus typhle

The transfer of acquired and specific immunity against previously encountered bacteria from mothers to offspring boosts the immune response of the next generation and supports the development of a successful pathogen defense. While most studies claim that the transfer of immunity is a maternal trait, in the sex‐role‐reversed pipefish Syngnathus typhle, fathers nurse the embryos over a placenta‐like structure, which opens the door for additional paternal immune priming. We examined the potential and persistence of bacteria‐type‐specific parental immune priming in the pipefish S. typhle over maturation time using a fully reciprocal design with two different bacteria species (Vibrio spp. and Tenacibaculum maritimum). Our results suggest that S. typhle is able to specifically prime the next generation against prevalent local bacteria and to a limited extent even also against newly introduced bacteria species. Long‐term protection was thereby maintained only against prevailing Vibrio bacteria. Maternal and paternal transgenerational immune priming can complement each other, as they affect different pathways of the offspring immune system and come with distinct degree of specificity. The differential regulation of DNA‐methylation genes upon parental bacteria exposure in premature pipefish offspring indicates that epigenetic regulation processes are involved in transferring immune‐related information across generations. The identified trade‐offs between immune priming and reproduction determine TGIP as a costly trait, which might constrain the evolution of long‐lasting TGIP, if parental and offspring generations do not share the same parasite assembly.     

Universal fast walking applied to cDNA

The elucidation of cDNA sequence remains problematic in cases such as genes possessing long coding regions, low expression levels, or poor library coverage. The recently described Universal Fast Walk (UFW) procedure offers a means of determining DNA sequence adjacent to characterised regions. To date, however, the approach has been applied only to genomic DNA. We demonstrate the first successful application of the UFW procedure to the elucidation of cDNA sequence, a previously unknown region of the large tammar wallaby ATRX gene in the theoretically more challenging 3' direction. To do this, we modified the previously published method by including an initial linear amplification and a final, fully nested PCR. We also exchanged buffers between preparative enzyme reactions to ensure optimal conditions for successive steps. These additional steps ensured a product not observed in their absence. UFW, therefore, represents a powerful alternative mechanism for the cloning and sequencing of cDNA, harnessing the exquisite sensitivity and specificity of fully nested PCR in challenging cloning scenarios where conventional 5' or 3' RACE may fail.     

The Good,the Bad,and the Ugly: in search of gold standards for assessing functional genetic screen quality

Variable screen quality, off‐target effects, and unclear false discovery rates often hamper large‐scale functional genomic screens in mammalian cells. Hart et al (2014) introduce gold standard reference sets of essential and non‐essential genes, aiming at standardizing the analysis of genome‐wide screens. This work provides a framework to compare both the quality and analysis methods of functional genetic screens.