Application of in silico bulked segregant analysis for rapid development of markers linked to Bean common mosaic virus resistance in common bean

Background

Common bean was one of the first crops that benefited from the development and utilization of molecular marker-assisted selection (MAS) for major disease resistance genes. Efficiency of MAS for breeding common bean is still hampered, however, due to the dominance, linkage phase, and loose linkage of previously developed markers. Here we applied in silico bulked segregant analysis (BSA) to the BeanCAP diversity panel, composed of over 500 lines and genotyped with the BARCBEAN_3 6K SNP BeadChip, to develop codominant and tightly linked markers to the I gene controlling resistance to Bean common mosaic virus (BCMV).

Results

We physically mapped the genomic region underlying the I gene. This locus, in the distal arm of chromosome Pv02, contains seven putative NBS-LRR-type disease resistance genes. Two contrasting bulks, containing BCMV host differentials and ten BeanCAP lines with known disease reaction to BCMV, were subjected to in silico BSA for targeting the I gene and flanking sequences. Two distinct haplotypes, containing a cluster of six single nucleotide polymorphisms (SNP), were associated with resistance or susceptibility to BCMV. One-hundred and twenty-two lines, including 115 of the BeanCAP panel, were screened for BCMV resistance in the greenhouse, and all of the resistant or susceptible plants displayed distinct SNP haplotypes as those found in the two bulks. The resistant/susceptible haplotypes were validated in 98 recombinant inbred lines segregating for BCMV resistance. The closest SNP (~25-32 kb) to the distal NBS-LRR gene model for the I gene locus was targeted for conversion to codominant KASP (Kompetitive Allele Specific PCR) and CAPS (Cleaved Amplified Polymorphic Sequence) markers. Both marker systems accurately predicted the disease reaction to BCMV conferred by the I gene in all screened lines of this study.

Conclusions

We demonstrated the utility of the in silico BSA approach using genetically diverse germplasm, genotyped with a high-density SNP chip array, to discover SNP variation at a specific targeted genomic region. In common bean, many disease resistance genes are mapped and their physical genomic position can now be determined, thus the application of this approach will facilitate further development of codominant and tightly linked markers for use in MAS.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-903) contains supplementary material, which is available to authorized users.     

RNAi Experiments in D. melanogaster: Solutions to the Overlooked Problem of Off-Targets Shared by Independent dsRNAs

Background

RNAi technology is widely used to downregulate specific gene products. Investigating the phenotype induced by downregulation of gene products provides essential information about the function of the specific gene of interest. When RNAi is applied in Drosophila melanogaster or Caenorhabditis elegans, often large dsRNAs are used. One of the drawbacks of RNAi technology is that unwanted gene products with sequence similarity to the gene of interest can be down regulated too. To verify the outcome of an RNAi experiment and to avoid these unwanted off-target effects, an additional non-overlapping dsRNA can be used to down-regulate the same gene. However it has never been tested whether this approach is sufficient to reduce the risk of off-targets.

Methodology

We created a novel tool to analyse the occurance of off-target effects in Drosophila and we analyzed 99 randomly chosen genes.

Principal Findings

Here we show that nearly all genes contain non-overlapping internal sequences that do show overlap in a common off-target gene.

Conclusion

Based on our in silico findings, off-target effects should not be ignored and our presented on-line tool enables the identification of two RNA interference constructs, free of overlapping off-targets, from any gene of interest.     

Identification of pneumococcal colonization determinants in the stringent response pathway facilitated by genomic diversity

Background

Understanding genetic determinants of a microbial phenotype generally involves creating and comparing isogenic strains differing at the locus of interest, but the naturally existing genomic and phenotypic diversity of microbial populations has rarely been exploited. Here we report use of a diverse collection of 616 carriage isolates of Streptococcus pneumoniae and their genome sequences to help identify a novel determinant of pneumococcal colonization.

Results

A spontaneously arising laboratory variant (SpnYL101) of a capsule-switched TIGR4 strain (TIGR4:19F) showed reduced ability to establish mouse nasal colonization and lower resistance to non-opsonic neutrophil-mediated killing in vitro, a phenotype correlated with in vivo success. Whole genome sequencing revealed 5 single nucleotide polymorphisms (SNPs) affecting 4 genes in SpnYL101 relative to its ancestor. To evaluate the effect of variation in each gene, we performed an in silico screen of 616 previously published genome sequences to identify pairs of closely-related, serotype-matched isolates that differ at the gene of interest, and compared their resistance to neutrophil-killing. This method allowed rapid examination of multiple candidate genes and found phenotypic differences apparently associated with variation in SP_1645, a RelA/ SpoT homolog (RSH) involved in the stringent response. To establish causality, the alleles corresponding to SP_1645 were switched between the TIGR4:19F and SpnYL101. The wild-type SP_1645 conferred higher resistance to neutrophil-killing and competitiveness in mouse colonization. Using a similar strategy, variation in another RSH gene (TIGR4 locus tag SP_1097) was found to alter resistance to neutrophil-killing.

Conclusions

These results indicate that analysis of naturally existing genomic diversity complements traditional genetics approaches to accelerate genotype-phenotype analysis.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1573-6) contains supplementary material, which is available to authorized users.     

Mapping and analysis of a novel candidate Fusarium wilt resistance gene FOC1 in Brassica oleracea

Background

Cabbage Fusarium wilt is a major disease worldwide that can cause severe yield loss in cabbage (Brassica olerecea). Although markers linked to the resistance gene FOC1 have been identified, no candidate gene for it has been determined so far. In this study, we report the fine mapping and analysis of a candidate gene for FOC1 using a double haploid (DH) population with 160 lines and a F₂ population of 4000 individuals derived from the same parental lines.

Results

We confirmed that the resistance to Fusarium wilt was controlled by a single dominant gene based on the resistance segregation ratio of the two populations. Using InDel primers designed from whole-genome re-sequencing data for the two parental lines (the resistant inbred-line 99–77 and the highly susceptible line 99–91) and the DH population, we mapped the resistance gene to a 382-kb genomic region on chromosome C06. Using the F₂ population, we narrowed the region to an 84-kb interval that harbored ten genes, including four probable resistance genes (R genes): Bol037156, Bol037157, Bol037158 and Bol037161 according to the gene annotations from BRAD, the genomic database for B. oleracea. After correcting the model of the these genes, we re-predicted two R genes in the target region: re-Bol037156 and re-Bol0371578. The latter was excluded after we compared the two genes’ sequences between ten resistant materials and ten susceptible materials. For re-Bol037156, we found high identity among the sequences of the resistant lines, while among the susceptible lines, there were two types of InDels (a 1-bp insertion and a 10-bp deletion), each of which caused a frameshift and terminating mutation in the cDNA sequences. Further sequence analysis of the two InDel loci from 80 lines (40 resistant and 40 susceptible) also showed that all 40 R lines had no InDel mutation while 39 out of 40 S lines matched the two types of loci. Thus re-Bol037156 was identified as a likely candidate gene for FOC1 in cabbage.

Conclusions

This work may lay the foundation for marker-assisted selection as well as for further function analysis of the FOC1 gene.     

Characterization of the MLO gene family in Rosaceae and gene expression analysis in Malus domestica

Background

Powdery mildew (PM) is a major fungal disease of thousands of plant species, including many cultivated Rosaceae. PM pathogenesis is associated with up-regulation of MLO genes during early stages of infection, causing down-regulation of plant defense pathways. Specific members of the MLO gene family act as PM-susceptibility genes, as their loss-of-function mutations grant durable and broad-spectrum resistance.

Results

We carried out a genome-wide characterization of the MLO gene family in apple, peach and strawberry, and we isolated apricot MLO homologs through a PCR-approach. Evolutionary relationships between MLO homologs were studied and syntenic blocks constructed. Homologs that are candidates for being PM susceptibility genes were inferred by phylogenetic relationships with functionally characterized MLO genes and, in apple, by monitoring their expression following inoculation with the PM causal pathogen Podosphaera leucotricha.

Conclusions

Genomic tools available for Rosaceae were exploited in order to characterize the MLO gene family. Candidate MLO susceptibility genes were identified. In follow-up studies it can be investigated whether silencing or a loss-of-function mutations in one or more of these candidate genes leads to PM resistance.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-618) contains supplementary material, which is available to authorized users.     

Analyses of phylogeny, evolution, conserved sequences and genome-wide expression of the ICK/KRP family of plant CDK inhibitors

Background and Aims

The cell cycle is controlled by cyclin-dependent kinases (CDKs), and CDK inhibitors are major regulators of their activities. The ICK/KRP family of CDK inhibitors has been reported in several plants, with seven members in arabidopsis; however, the phylogenetic relationship among members in different species is unknown. Also, there is a need to understand how these genes and proteins are regulated. Furthermore, little information is available on the functional differences among ICK/KRP family members.

Methods

We searched publicly available databases and identified over 120 unique ICK/KRP protein sequences from more than 60 plant species. Phylogenetic analysis was performed using 101 full-length sequences from 40 species and intron–exon organization of ICK/KRP genes in model species. Conserved sequences and motifs were analysed using ICK/KRP protein sequences from arabidopsis (Arabidopsis thaliana), rice (Orysa sativa) and poplar (Populus trichocarpa). In addition, gene expression was examined using microarray data from arabidopsis, rice and poplar, and further analysed by RT-PCR for arabidopsis.

Key Results and Conclusions

Phylogenetic analysis showed that plant ICK/KRP proteins can be grouped into three major classes. Whereas the C-class contains sequences from dicotyledons, monocotyledons and gymnosperms, the A- and B-classes contain only sequences from dicotyledons or monocotyledons, respectively, suggesting that the A- and B-classes might have evolved from the C-class. This classification is also supported by exon–intron organization. Genes in the A- and B- classes have four exons, whereas genes in the C-class have only three exons. Analysis of sequences from arabidopsis, rice and poplar identified conserved sequence motifs, some of which had not been described previously, and putative functional sites. The presence of conserved motifs in different family members is consistent with the classification. In addition, gene expression analysis showed preferential expression of ICK/KRP genes in certain tissues. A model has been proposed for the evolution of this gene family in plants.     

Highly Conserved Non-Coding Sequences and the 18q Critical Region for Short Stature: A Common Mechanism of Disease?

Background

Isolated growth hormone deficiency (IGHD) and multiple pituitary hormone deficiency (MPHD) are heterogeneous disorders with several different etiologies and they are responsible for most cases of short stature. Mutations in different genes have been identified but still many patients did not present mutations in any of the known genes. Chromosomal rearrangements may also be involved in short stature and, among others, deletions of 18q23 defined a critical region for the disorder. No gene was yet identified.

Methodology/Principal Findings

We now report a balanced translocation X;18 in a patient presenting a breakpoint in 18q23 that was surprisingly mapped about 500 Kb distal from the short stature critical region. It separated from the flanking SALL3 gene a region enriched in highly conserved non-coding elements (HCNE) that appeared to be regulatory sequences, active as enhancers or silencers during embryonic development.

Conclusion

We propose that, during pituitary development, the 18q rearrangement may alter expression of 18q genes or of X chromosome genes that are translocated next to the HCNEs. Alteration of expression of developmentally regulated genes by translocation of HCNEs may represent a common mechanism for disorders associated to isolated chromosomal rearrangements.     

Proof by synthesis of Tobacco mosaic virus

Background

Synthetic biology is a discipline that includes making life forms artificially from chemicals. Here, a DNA molecule was enzymatically synthesized in vitro from DNA templates made from oligonucleotides representing the text of the first Tobacco mosaic virus (TMV) sequence elucidated in 1982. No infectious DNA molecule of that seminal reference sequence exists, so the goal was to synthesize it and then build viral chimeras.

Results

RNA was transcribed from synthetic DNA and encapsidated with capsid protein in vitro to make synthetic virions. Plants inoculated with the virions did not develop symptoms. When two nucleotide mutations present in the original sequence, but not present in most other TMV sequences in GenBank, were altered to reflect the consensus, the derivative synthetic virions produced classic TMV symptoms. Chimeras were then made by exchanging TMV capsid protein DNA with Tomato mosaic virus (ToMV) and Barley stripe mosaic virus (BSMV) capsid protein DNA. Virus expressing ToMV capsid protein exhibited altered, ToMV-like symptoms in Nicotiana sylvestris. A hybrid ORF6 protein unknown to nature, created by substituting the capsid protein genes in the virus, was found to be a major symptom determinant in Nicotiana benthamiana. Virus expressing BSMV capsid protein did not have an extended host range to barley, but did produce novel symptoms in N. benthamiana.

Conclusions

This first report of the chemical synthesis and artificial assembly of a plant virus corrects a long-standing error in the TMV reference genome sequence and reveals that unnatural hybrid virus proteins can alter symptoms unexpectedly.     

Relationship between symptoms and gene expression induced by the infection of three strains of Rice dwarf virus

Background

Rice dwarf virus (RDV) is the causal agent of rice dwarf disease, which often results in severe yield losses of rice in East Asian countries. The disease symptoms are stunted growth, chlorotic specks on leaves, and delayed and incomplete panicle exsertion. Three RDV strains, O, D84, and S, were reported. RDV-S causes the most severe symptoms, whereas RDV-O causes the mildest. Twenty amino acid substitutions were found in 10 of 12 virus proteins among three RDV strains.

Methodology/Principal Findings

We analyzed the gene expression of rice in response to infection with the three RDV strains using a 60-mer oligonucleotide microarray to examine the relationship between symptom severity and gene responses. The number of differentially expressed genes (DEGs) upon the infection of RDV-O, -D84, and -S was 1985, 3782, and 6726, respectively, showing a correlation between the number of DEGs and symptom severity. Many DEGs were related to defense, stress response, and development and morphogenesis processes. For defense and stress response processes, gene silencing-related genes were activated by RDV infection and the degree of activation was similar among plants infected with the three RDV strains. Genes for hormone-regulated defense systems were also activated by RDV infection, and the degree of activation seemed to be correlated with the concentration of RDV in plants. Some development and morphogenesis processes were suppressed by RDV infection, but the degree of suppression was not correlated well with the RDV concentration.

Conclusions/Significance

Gene responses to RDV infection were regulated differently depending on the gene groups regulated and the strains infecting. It seems that symptom severity is associated with the degree of gene response in defense-related and development- and morphogenesis-related processes. The titer levels of RDV in plants and the amino acid substitutions in RDV proteins could be involved in regulating such gene responses.