Allele-specific CAPS marker in a <Emphasis Type="Italic">Ve1</Emphasis> homolog of <Emphasis Type="Italic">Capsicum annuum</Emphasis> for improved selection of <Emphasis Type="Italic">Verticillium dahliae</Emphasis> resistance

Verticillium wilt (Verticillium dahliae) is an economically important disease for many high-value crops. The pathogen is difficult to manage due to the long viability of its resting structures, wide host range, and the inability of fungicides to affect the pathogen once in the plant vascular system. In chile pepper (Capsicum annuum), breeding for resistance to Verticillium wilt is especially challenging due to the limited resistance sources. The dominant Ve locus in tomato (Solanum lycopersicum) contains two closely linked and inversely oriented genes, Ve1 and Ve2. Homologs of Ve1 have been characterized in diverse plant species, and interfamily transfer of Ve1 confers race-specific resistance. Queries in the chile pepper WGS database in NCBI with Ve1 and Ve2 sequences identified one open reading frame (ORF) with homology to the tomato Ve genes. Comparison of the candidate CaVe (Capsicum annuum Ve) gene sequences from susceptible and resistant accessions revealed 16 single nucleotide polymorphisms (SNPs) and several haplotypes. A homozygous haplotype was identified for the susceptible accessions and for resistant accessions. We developed a cleaved amplified polymorphic sequence (CAPS) molecular marker within the coding region of CaVe and screened diverse germplasm that has been previously reported as being resistant to Verticillium wilt in other regions. Based on our phenotyping using the New Mexico V. dahliae isolate, the marker could select resistance accessions with 48% accuracy. This molecular marker is a promising tool towards marker-assisted selection for Verticillium wilt resistance and has the potential to improve the efficacy of chile pepper breeding programs, but does not eliminate the need for a bioassay. Furthermore, this work provides a basis for future research in this important pathosystem.     

High-throughput sequencing of <Emphasis Type="Italic">Prunus ferganensis</Emphasis> indicates that it is a geographical population of <Emphasis Type="Italic">P. persica</Emphasis>

Peach belongs to the genus Prunus, which includes Prunus persica and its relative species, P. mira, P. davidiana, P. kansuensis, and P. ferganensis. Of these, P. ferganensis have been classified as a species, subspecies, or geographical population of P. persica. To explore the genetic difference between P. ferganensis and P. persica, high-throughput sequencing was used in different peach accessions belonging to different species. First, low-depth sequencing data of peach accessions belonging to four categories revealed that similarity between P. ferganensis and P. persica was similar to that between P. persica accessions from different geographical populations. Then, to further detect the genomic variation in P. ferganensis, the P. ferganensis accession “Xinjiang Pan Tao 1” and the P. persica accession “Xia Miao 1” were sequenced with high depth, and sequence reads were assembled. The results showed that the collinearity of “Xinjiang Pan Tao 1” with the reference genome “Lovell” was higher than that of “Xia Miao 1” and “Lovell” peach. Additionally, the number of genetic variants, including single nucleotide polymorphisms (SNPs), structural variations (SVs), and the specific genes annotated from unmapped sequence in “Xia Miao 1” was higher than that in “Xinjiang Pan Tao 1” peach. The data showed that there was a close distance between “Xinjiang Pan Tao 1” (P. ferganensis) and reference genome which belong to P. persica, comparing “Xia Miao 1” (P. persica) and reference ones. The results accompany with phylogenetic tree and structure analysis confirmed that P. ferganensis should be considered as a geographic population of P. persica rather than a subspecies or a distinct species. Furthermore, gene ontology analysis was performed using the gene comprising large-effect variation to understand the phenotypic difference between two accessions. The result revealed that the pathways of gene function affected by SVs but SNPs and insertion-deletions markedly differed between the two peach accessions.     

A new <Emphasis Type="Italic">Em</Emphasis>-like protein from <Emphasis Type="Italic">Lactuca sativa</Emphasis>, <Emphasis Type="Italic">LsEm1</Emphasis>, enhances drought and salt stress tolerance in <Emphasis Type="Italic">Escherichia coli</Emphasis> and rice

Late embryogenesis abundant (LEA) proteins are closely related to abiotic stress tolerance of plants. In the present study, we identified a novel Em-like gene from lettuce, termed LsEm1, which could be classified into group 1 LEA proteins, and shared high homology with Cynara cardunculus Em protein. The LsEm1 protein contained three different 20-mer conserved elements (C-element, N-element, and M-element) in the C-termini, N-termini, and middle-region, respectively. The LsEm1 mRNAs were accumulated in all examined tissues during the flowering and mature stages, with a little accumulation in the roots and leaves during the seedling stage. Furthermore, the LsEm1 gene was also expressed in response to salt, dehydration, abscisic acid (ABA), and cold stresses in young seedlings. The LsEm1 protein could effectively reduce damage to the lactate dehydrogenase (LDH) and protect LDH activity under desiccation and salt treatments. The Escherichia coli cells overexpressing the LsEm1 gene showed a growth advantage over the control under drought and salt stresses. Moreover, LsEm1-overexpressing rice seeds were relatively sensitive to exogenously applied ABA, suggesting that the LsEm1 gene might depend on an ABA signaling pathway in response to environmental stresses. The transgenic rice plants overexpressing the LsEm1 gene showed higher tolerance to drought and salt stresses than did wild-type (WT) plants on the basis of the germination performances, higher survival rates, higher chlorophyll content, more accumulation of soluble sugar, lower relative electrolyte leakage, and higher superoxide dismutase activity under stress conditions. The LsEm1-overexpressing rice lines also showed less yield loss compared with WT rice under stress conditions. Furthermore, the LsEm1 gene had a positive effect on the expression of the OsCDPK9, OsCDPK13, OsCDPK15, OsCDPK25, and rab21 (rab16a) genes in transgenic rice under drought and salt stress conditions, implying that overexpression of these genes may be involved in the enhanced drought and salt tolerance of transgenic rice. Thus, this work paves the way for improvement in tolerance of crops by genetic engineering breeding.     

Development and characterization of <Emphasis Type="Italic">japonica</Emphasis> rice lines carrying the brown planthopper-resistance genes <Emphasis Type="Italic">BPH12</Emphasis> and <Emphasis Type="Italic">BPH6</Emphasis>

The brown planthopper (Nilaparvata lugens Stål; BPH) has become a severe constraint on rice production. Identification and pyramiding BPH-resistance genes is an economical and effective solution to increase the resistance level of rice varieties. All the BPH-resistance genes identified to date have been from indica rice or wild species. The BPH12 gene in the indica rice accession B14 is derived from the wild species Oryza latifolia. Using an F₂ population from a cross between the indica cultivar 93-11 and B14, we mapped the BPH12 gene to a 1.9-cM region on chromosome 4, flanked by the markers RM16459 and RM1305. In this population, BPH12 appeared to be partially dominant and explained 73.8% of the phenotypic variance in BPH resistance. A near-isogenic line (NIL) containing the BPH12 locus in the background of the susceptible japonica variety Nipponbare was developed and crossed with a NIL carrying BPH6 to generate a pyramid line (PYL) with both genes. BPH insects showed significant differences in non-preference in comparisons between the lines harboring resistance genes (NILs and PYL) and Nipponbare. BPH growth and development were inhibited and survival rates were lower on the NIL-BPH12 and NIL-BPH6 plants compared to the recurrent parent Nipponbare. PYL-BPH6 + BPH12 exhibited 46.4, 26.8 and 72.1% reductions in population growth rates (PGR) compared to NIL-BPH12, NIL-BPH6 and Nipponbare, respectively. Furthermore, insect survival rates were the lowest on the PYL-BPH6 + BPH12 plants. These results demonstrated that pyramiding different BPH-resistance genes resulted in stronger antixenotic and antibiotic effects on the BPH insects. This gene pyramiding strategy should be of great benefit for the breeding of BPH-resistant japonica rice varieties.     

Genome-wide association study of <Emphasis Type="Italic">Mycobacterium avium</Emphasis> subspecies <Emphasis Type="Italic">Paratuberculosis</Emphasis> infection in Chinese Holstein

Background

Paratuberculosis is a contagious, chronic and enteric disease in ruminants, which is caused by Mycobacterium avium subspecies paratuberculosis (MAP) infection, resulting in enormous economic losses worldwide. There is currently no effective cure for MAP infection or a vaccine, it is thus important to explore the genetic variants that contribute to host susceptibility to infection by MAP, which may provide a better understanding of the mechanisms of paratuberculosis and benefit animal genetic improvement. Herein we performed a genome-wide association study (GWAS) to identify genomic regions and candidate genes associated with susceptibility to MAP infection in dairy cattle.

Results

Using Illumina Bovine 50?K (54,609 SNPs) and GeneSeek HD (138,893 SNPs) chips, two analytical approaches were performed, GRAMMAR-GC and ROADTRIPS in 937 Chinese Holstein cows, among which individuals genotyped by the 50?K chip were imputed to HD SNPs with Beagle software. Consequently, 15 and 11 significant SNPs (P?<?5?×?10^??5) were identified with GRAMMAR-GC and ROADTDRIPS, respectively. A total of 10 functional genes were in proximity to (i.e., within 1?Mb) these SNPs, including IL4, IL5, IL13, IRF1, MyD88, PACSIN1, DEF6, TDP2, ZAP70 and CSF2. Functional enrichment analysis showed that these genes were involved in immune related pathways, such as interleukin, T cell receptor signaling pathways and inflammatory bowel disease (IBD), implying their potential associations with susceptibility to MAP infection. In addition, by examining the publicly available cattle QTLdb, a previous QTL for MAP was found to be overlapped with one of regions detected currently at 32.5?Mb on BTA23, where the TDP2 gene was anchored.

Conclusions

In conclusion, we identified 26 SNPs located on 15 chromosomes in the Chinese Holstein population using two GWAS strategies with high density SNPs. Integrated analysis of GWAS, biological functions and the reported QTL information helps to detect positional candidate genes and the identification of regions associated with susceptibility to MAP traits in dairy cattle.

     

Association between expression levels and growth trait-related SNPs located in promoters of the <Emphasis Type="Italic">MC4R</Emphasis> and <Emphasis Type="Italic">MSTN</Emphasis> genes in <Emphasis Type="Italic">Spinibarbus hollandi</Emphasis>

Melanocortin 4 receptor: (MC4R) and Myostatin (MSTN) are two important growth trait-related genes in animals. In this study, we showed that two SNPs, MC4R-719A>G and MSTN-519C>T, found in the promoters of the MC4R and MSTN genes, respectively, are both associated with growth traits in Spinibarbus hollandi. Furthermore, we observed that there were significant associations between the expression levels of the MC4R and MSTN genes and these two growth trait-related SNPs. The expression level of MC4R gene in brain was lower in GG genotype fish with extremely high growth performance than that in AA genotype fish with extremely low growth performance. Expression level of the MSTN gene in muscle was lower in TT genotype fish with extremely high growth performance than that in CC and CT genotype fish with lower growth performance. The results indicated that these SNPs located in the promoters of MC4R and MSTN are associated with growth-related traits through modification of gene expression levels. The MSTN and MC4R SNPs may have useful application in effective marker-assisted selection aimed to increase output in S. hollandi.     

Role of Arabidopsis <Emphasis Type="Italic">ABF1</Emphasis>/<Emphasis Type="Italic">3</Emphasis>/<Emphasis Type="Italic">4</Emphasis> during <Emphasis Type="Italic">det1</Emphasis> germination in salt and osmotic stress conditions

Key message

Arabidopsis det1 mutants exhibit salt and osmotic stress resistant germination. This phenotype requires HY5, ABF1, ABF3, and ABF4.

Abstract

While DE-ETIOLATED 1 (DET1) is well known as a negative regulator of light development, here we describe how det1 mutants also exhibit altered responses to salt and osmotic stress, specifically salt and mannitol resistant germination. LONG HYPOCOTYL 5 (HY5) positively regulates both light and abscisic acid (ABA) signalling. We found that hy5 suppressed the det1 salt and mannitol resistant germination phenotype, thus, det1 stress resistant germination requires HY5. We then queried publically available microarray datasets to identify genes downstream of HY5 that were differentially expressed in det1 mutants. Our analysis revealed that ABA regulated genes, including ABA RESPONSIVE ELEMENT BINDING FACTOR 3 (ABF3), are downregulated in det1 seedlings. We found that ABF3 is induced by salt in wildtype seeds, while homologues ABF4 and ABF1 are repressed, and all three genes are underexpressed in det1 seeds. We then investigated the role of ABF3, ABF4, and ABF1 in det1 phenotypes. Double mutant analysis showed that abf3, abf4, and abf1 all suppress the det1 salt/osmotic stress resistant germination phenotype. In addition, abf1 suppressed det1 rapid water loss and open stomata phenotypes. Thus interactions between ABF genes contribute to det1 salt/osmotic stress response phenotypes.

     

Horizontal transfer of the C-termini of <Emphasis Type="Italic">tccC</Emphasis> genes in <Emphasis Type="Italic">Photorhabdus</Emphasis> and <Emphasis Type="Italic">Xenorhabdus</Emphasis>

The high molecular weight insecticidal toxin complexes (Tcs), including four toxin-complex loci (tca, tcb, tcc and tcd), were first identified in Photorhabdus luminescens W14. Each member of tca, tcb or tcc is required for oral toxicity of Tcs. However, the sequence sources of the C-termini of tccC3, tccC4, tccC6 and tccC7 are unknown. Here, we performed a whole genome survey to identify the orthologs of Tc genes, and found 165 such genes in 14 bacterial genomes, including 40 genes homologous to tccC1-7 in P. luminescens TT01. The sequence sources of the C-termini of tccC2-6 were determined by sequence analysis. Further phylogenetic investigations suggested that the C-termini of 6 tccC genes experienced horizontal gene transfer events.     

Tandem duplications of receptor-like kinase genes reshaped the homologous chromosome 1 regions in <Emphasis Type="Italic">Oryza sativa</Emphasis> ssp. <Emphasis Type="Italic">indica</Emphasis> and <Emphasis Type="Italic">japonica</Emphasis> and their possible ancestral genomes

Plant receptor-like kinase (Rlk) genes form a large family, each encoding a protein with a signal motif, a single transmembrane region, and a cytoplasmic kinase domain. Various gene duplications have contributed to the establishment and expansion of the family. Here, we characterized the formation and evolution of the Rlk gene family in cultivated rice and their possible progenitors. Using wheat Rlk gene sequences, we identified orthologs from the genomes of domesticated rice subspecies Oryza sativa ssp. japonica and ssp. indica and their putative progenitors O. glaberrima and O. rufipogon. The four chromosome 1 orthologous regions ranged from 103 to 281 kb comprising 181 syntenic blocks with 75 to 100% sequence identity. These regions contained 11–19 Triticum aestivum kinases (Taks) and 10–15 Lr10 receptor-like kinases (Lrks) organized in clusters and 3–12 transposable elements (TEs). Dot plot analyses showed that the 4 regions had 21–37 conserved catalytic domains, mainly in protein kinases (PKs) and tyrosine kinases (TyrKs) in coupling state. Over 50% of the sequences of glaberrima/rufipogon and japonica/indica pairs were colinear, while japonica/indica displayed a marked sequence expansion with duplicated genes and TEs. A total of 2312 single nucleotide polymorphisms (SNPs) and insertion-deletions (INDELs) were identified between japonica and indica. Duplication of the Rlk genes in O. glaberrima and O. rufipogon occurred after the grass species radiation and before the divergence of O. rufipogon from O. glaberrima; the orthologous Rlk genes from O. japonica and O. indica duplicated after O. sativa separated from O. rufipogon; paralogs, obtained through extensive duplication, happened after the separation of rice from maize. Tandem duplication was the major factor contributing to the gene copy number variation and genome size expansion.     

Genetic complementation analysis of rice sucrose transporter genes in Arabidopsis <Emphasis Type="Italic">SUC2</Emphasis> mutant <Emphasis Type="Italic">atsuc2</Emphasis>

Sucrose transporters (SUTs) play a critical role on the phloem plasma membrane in loading sucrose into the phloem of source leaves for long-distance transport to sink organs. Rice has a small gene family of five SUTs, Oryza sativa SUT1 (OsSUT1) to OsSUT5. To identify rice SUTs that function as phloem loaders, we adopted a growth restoration assay of the severe growth retardation phenotype of atsuc2, a mutant of the best-characterized Arabidopsis phloem loader AtSUC2, by introducing OsSUTs. The rice SUT genes were expressed by two different promoters, the native phloem-specific promoter of AtSUC2 (pAtSUC2) and the constitutive Cauliflower Mosaic Virus 35S (pCaMV35S) promoter. Of all the transgenic atsuc2 plants, only pAtSUC2: OsSUT1 complemented the atsuc2 mutant phenotype in a comparable manner to wild type (WT), and consistent levels of soluble sugars and starch were recovered compared to those of WT. This suggests that OsSUT1 is a functional ortholog of the Arabidopsis AtSUC2 and functions as an apoplastic phloem loader. In addition, ossut1 mutants were produced via anther culture and their primary carbohydrate levels and growth phenotypes were indistinguishable from those of WT. This suggests that the rice phloem loader OsSUT1 function may not be essential for rice vegetative growth under normal conditions.     

Genome-wide characterization and stress-responsive expression profiling of <Emphasis Type="Italic">MCM</Emphasis> genes in <Emphasis Type="Italic">Brassica oleracea</Emphasis> and <Emphasis Type="Italic">Brassica rapa</Emphasis>

The Minichromosome maintenance protein [MCM (2-7)] complex is associated with helicase activity for replication fork formation during DNA replication. We identified and characterized each 12 putative MCM genes from Brassica oleracea and Brassica rapa. MCM genes were classified into nine groups according to their evolutionary relationships. A high number of syntenic regions were present on chromosomes C03 and A03 in B. oleracea and B. rapa, respectively, compared to the other chromosomes. Expression analysis showed that most of the MCM(2-7) helicase-subunit genes and their coregulating MCM genes were upregulated during hydroxyurea (HU) induced stress in B. oleracea. In B. rapa, MCM(2-7) helicase genes BrMCM2_2, BrMCM7_1, BrMCM7_2 and their co-regulating genes were upregulated during replication stress. During cold stress, BoMCM6 in B. oleracea and BrMCM5 in B. rapa were remarkably upregulated. During salt stress, BoMCM6_2, BoMCM7_1, BoMCM8, BoMCM9, and BoMCM10 were markedly upregulated in B. oleracea. Hence, our study identified the candidate MCM family genes those possess abiotic stress-responsive behavior and DNA replication stress tolerance. As the first genome-wide analysis of MCM genes in B. oleracea and B. rapa, this work provides a foundation to develop stress responsive plants. Further functional and molecular studies on MCM genes will be helpful to enhance stress tolerance in plants.     

Association between <Emphasis Type="Italic">cagA</Emphasis>, <Emphasis Type="Italic">vacAi</Emphasis>, and <Emphasis Type="Italic">dupA</Emphasis> genes of <Emphasis Type="Italic">Helicobacter pylori</Emphasis> and gastroduodenal pathologies in Chilean patients

In addition to the already known cagA gene, novel genetic markers have been associated with Helicobacter pylori (H. pylori) virulence: the dupA and vacAi genes. These genes might play an important role as specific markers to determine the clinical outcome of the disease, especially the vacAi gene, which has been expected to be a good marker of severe pathologies like gastric adenocarcinoma. In the present study, the association of cagA, dupA, and vacAi genes with gastroduodenal pathologies in Chilean patients was studied. One hundred and thirty-two patients positive for H. pylori were divided into two groups—non-severe and severe gastric pathologies—and investigated for the presence of cagA, dupA, and vacAi H. pylori virulence genes by PCR. The cagA gene was detected in 20/132 patients (15.2%), the vacAi1 gene was detected in 54/132 patients (40.9%), the vacAi2 gene was detected in 26/132 patients (19.7%), and the dupA gene was detected in 50/132 (37.9%) patients. Logistic regression model analysis showed that the vacAi1 isoform gene in the infected strains and the severity of the diseases outcome were highly associated, causing severe gastric damage that may lead to gastric cancer (p < 0.0001; OR = 8.75; 95% CI 3.54–21.64). Conversely, cagA (p = 0.3507; OR = 1.62; 95% CI 0.59–4.45) and vacAi2 (p = 0.0114; OR = 3.09; 95% CI 1.26–7.60) genes were not associated with damage, while the dupA gene was associated significantly with non-severe clinical outcome (p = 0.0032; OR = 0.25; 95% CI 0.09–0.65). In addition, dupA gene exerts protection against severe gastric pathologies induced by vacAi1 by delaying the outcome of the disease by approximately 20 years.     

Development of species diagnostic SNP markers for quality control genotyping in four rice (<Emphasis Type="Italic">Oryza</Emphasis> L.) species

Species misclassification (misidentification) and handling errors have been frequently reported in various plant species conserved at diverse gene banks, which could restrict use of germplasm for correct purpose. The objectives of the present study were to (i) determine the extent of genotyping error (reproducibility) on DArTseq-based single-nucleotide polymorphisms (SNPs); (ii) determine the proportion of misclassified accessions across 3134 samples representing three African rice species complex (Oryza glaberrima, O. barthii, and O. longistaminata) and an Asian rice (O. sativa), which are conserved at the AfricaRice gene bank; and (iii) develop species- and sub-species (ecotype)-specific diagnostic SNP markers for rapid and low-cost quality control (QC) analysis. Genotyping error estimated from 15 accessions, each replicated from 2 to 16 times, varied from 0.2 to 3.1%, with an overall average of 0.8%. Using a total of 3134 accessions genotyped with 31,739 SNPs, the proportion of misclassified samples was 3.1% (97 of the 3134 accessions). Excluding the 97 misclassified accessions, we identified a total of 332 diagnostic SNPs that clearly discriminated the three indigenous African species complex from Asian rice (156 SNPs), O. longistaminata accessions from both O. barthii and O. glaberrima (131 SNPs), and O. sativa spp. indica from O. sativa spp. japonica (45 SNPs). Using chromosomal position, minor allele frequency, and polymorphic information content as selection criteria, we recommended a subset of 24 to 36 of the 332 diagnostic SNPs for routine QC genotyping, which would be highly useful in determining the genetic identity of each species and correct human errors during routine gene bank operations.     

Improvement of seedling and panicle blast resistance in <Emphasis Type="Italic">Xian</Emphasis> rice varieties following <Emphasis Type="Italic">Pish</Emphasis> introgression

Rice blast is a serious disease caused by the filamentous ascomycetous fungus Magnaporthe oryzae. Incorporating disease resistance genes in rice varieties and characterizing the distribution of M. oryzae isolates form the foundation for enhancing rice blast resistance. In this study, the blast resistance gene Pish was observed to be differentially distributed in the genomes of rice sub-species. Specifically, Pish was present in 80.5% of Geng varieties, but in only 2.3% of Xian varieties. Moreover, Pish conferred resistance against only 23.5% of the M. oryzae isolates from the Geng-planting regions, but against up to 63.2% of the isolates from the Xian-planting regions. Thus, Pish may be an elite resistance gene for improving rice blast resistance in Xian varieties. Therefore, near-isogenic lines (NILs) with Pish and the polygene pyramid lines (PPLs) PPL^Pish/Pi1, PPL^Pish/Pi54, and PPL^Pish/Pi33 in the Xian background Yangdao 6 were generated using a molecular marker-assisted selection method. The results suggested that (1) Pish significantly improved rice blast resistance in Xian varieties, which exhibited considerably improved seedling and panicle blast resistance after Pish was introduced; (2) PPLs with Pish were more effective than the NILs with Pish regarding seedling and panicle blast resistance; (3) the PPL seedling and panicle blast resistance was improved by the complementary and overlapping effects of different resistance genes; and (4) the stability of NIL and PPL resistance varied under different environmental conditions, with only PPL^Pish/Pi54 exhibiting highly stable resistance in three natural disease nurseries (Jianyang, Jinggangshan, and Huangshan). This study provides new blast resistance germplasm resources and describes a novel molecular strategy for enhancing rice blast resistance.