Genetic mapping and development of co-segregating markers of <Emphasis Type="Italic">RpsQ</Emphasis>, which provides resistance to <Emphasis Type="Italic">Phytophthora sojae</Emphasis> in soybean

Key message

The RpsQ Phytophthora resistance locus was finely mapped to a 118-kb region on soybean chromosome 3. A best candidate gene was predicted and three co-segregating gene markers were developed.

Abstract

Phytophthora root rot (PRR), caused by Phytophthora sojae, is a major threat to sustainable soybean production. The use of genetically resistant cultivars is considered the most effective way to control this disease. The Chinese soybean cultivar Qichadou 1 exhibited a broad spectrum resistance, with a distinct resistance phenotype, following inoculation with 36 Chinese P. sojae isolates. Genetic analyses indicated that the disease resistance in Qichadou 1 is controlled by a single dominant gene. This gene locus was designated as RpsQ and mapped to a 118-kb region between BARCSOYSSR_03_0165 and InDel281 on soybean chromosome 3, and co-segregated with Insert11, Insert144 and SNP276. Within this region, there was only one gene Glyma.03g27200 encoding a protein with a typical serine/threonine protein kinase structure, and the expression pattern analysis showed that this gene induced by P. sojae infection, which was suggested as a best candidate gene of RpsQ. Candidate gene specific marker Insert144 was used to distinguish RpsQ from the other known Rps genes on chromosome 3. Identical polymerase chain reaction amplification products were produced for cultivars Qichadou 1 (RpsQ) and Ludou 4 (Rps9). All other cultivars carrying Rps genes on chromosome 3 produced different PCR products, which all lacked a 144-bp fragment present in Qichadou 1 and Ludou 4. The phenotypes of the analyzed cultivars combined with the physical position of the PRR resistance locus, candidate gene analyses, and the candidate gene marker test revealed RpsQ and Rps9 are likely the same gene, and confer resistance to P. sojae.

     

Next-generation sequencing to identify candidate genes and develop diagnostic markers for a novel Phytophthora resistance gene, <Emphasis Type="Italic">RpsHC18</Emphasis>, in soybean

Key message

A novel Phytophthora sojae resistance gene RpsHC18 was identified and finely mapped on soybean chromosome 3. Two NBS–LRR candidate genes were identified and two diagnostic markers of RpsHC18 were developed.

Abstract

Phytophthora root rot caused by Phytophthora sojae is a destructive disease of soybean. The most effective disease-control strategy is to deploy resistant cultivars carrying Phytophthora-resistant Rps genes. The soybean cultivar Huachun 18 has a broad and distinct resistance spectrum to 12 P. sojae isolates. Quantitative trait loci sequencing (QTL-seq), based on the whole-genome resequencing (WGRS) of two extreme resistant and susceptible phenotype bulks from an F_2:3 population, was performed, and one 767-kb genomic region with ΔSNP-index ≥ 0.9 on chromosome 3 was identified as the RpsHC18 candidate region in Huachun 18. The candidate region was reduced to a 146-kb region by fine mapping. Nonsynonymous SNP and haplotype analyses were carried out in the 146-kb region among ten soybean genotypes using WGRS. Four specific nonsynonymous SNPs were identified in two nucleotide-binding sites–leucine-rich repeat (NBS–LRR) genes, RpsHC18-NBL1 and RpsHC18-NBL2, which were considered to be the candidate genes. Finally, one specific SNP marker in each candidate gene was successfully developed using a tetra-primer ARMS-PCR assay, and the two markers were verified to be specific for RpsHC18 and to effectively distinguish other known Rps genes. In this study, we applied an integrated genomic-based strategy combining WGRS with traditional genetic mapping to identify RpsHC18 candidate genes and develop diagnostic markers. These results suggest that next-generation sequencing is a precise, rapid and cost-effective way to identify candidate genes and develop diagnostic markers, and it can accelerate Rps gene cloning and marker-assisted selection for breeding of P. sojae-resistant soybean cultivars.

     

Over-expression of the <Emphasis Type="Italic">Pseudomonas syringae</Emphasis> harpin-encoding gene <Emphasis Type="Italic">hrpZm</Emphasis> confers enhanced tolerance to Phytophthora root and stem rot in transgenic soybean

Phytophthora root and stem rot (PRR) caused by Phytophthora sojae is one of the most devastating diseases reducing soybean (Glycine max) production all over the world. Harpin proteins in many plant pathogenic bacteria were confirmed to enhance disease and insect resistance in crop plants. Here, a harpin protein-encoding gene hrpZpsta from the P. syringae pv. tabaci strain Psta218 was codon-optimized (renamed hrpZm) and introduced into soybean cultivars Williams 82 and Shennong 9 by Agrobacterium-mediated transformation. Three independent transgenic lines over-expressing hrpZm were obtained and exhibited stable and enhanced tolerance to P. sojae infection in T₂–T₄ generations compared to the non-transformed (NT) and empty vector (EV)-transformed plants. Quantitative real-time PCR (qRT-PCR) analysis revealed that the expression of salicylic acid-dependent genes PR1, PR12, and PAL, jasmonic acid-dependent gene PPO, and hypersensitive response (HR)-related genes GmNPR1 and RAR was significantly up-regulated after P. sojae inoculation. Moreover, the activities of defense-related enzymes such as phenylalanine ammonia lyase (PAL), polyphenoloxidase (PPO), peroxidase, and superoxide dismutase also increased significantly in the transgenic lines compared to the NT and EV-transformed plants after inoculation. Our results suggest that over-expression of the hrpZm gene significantly enhances PRR tolerance in soybean by eliciting resistance responses mediated by multiple defense signaling pathways, thus providing an alternative approach for development of soybean varieties with improved tolerance against the soil-borne pathogen PRR.     

Genome-wide analysis and characterization of molecular evolution of the <Emphasis Type="Italic">HCT</Emphasis> gene family in pear (<Emphasis Type="Italic">Pyrus bretschneideri</Emphasis>)

Lignin is a major component of stone cells in pear fruit, which significantly affects fruit quality. Hydroxycinnamoyl CoA: shikimate hydroxycinnamoyl transferase (HCT), a recently discovered enzyme in plants, is an important gene that participates in the formation of lignin. Although HCT gene cloning and expression patterns have been studied in several species, including pear, there is still no extensive genome-wide bioinformatics analysis on the whole gene family, and the evolutionary history of HCT gene family is still unknown. A total of 82 HCT genes were identified in pear, most of which have one or two exons, and all with the conserved HXXXD motif and transferase domains. Based on the structural characteristics and phylogenetic analysis of these sequences, the HCT gene family genes could be classified into four main groups. Structural analysis also revealed that 25 % of HCT genes share a MYB binding site. Expansion of the HCT gene family mostly occurred before the divergence between Arabidopsis and Rosaceae, with whole-genome duplication or segmental duplication events playing the most important role in the expansion of the HCT gene family in pear. At the same time, purifying selection also played a critical role in the evolution of HCT genes. Five of the 82 HCT genes were verified by qRT-PCR to correspond to the pattern of stone cell formation during pear fruit development. The genome-wide identification, chromosome localization, gene structures, synteny, and expression analyses of pear HCT genes provide an overall insight into HCT gene family and their potential involvement in growth and development of stone cells.     

Genome-Wide Analysis of the ZRT,IRT-Like Protein (ZIP) Family and Their Responses to Metal Stress in <Emphasis Type="Italic">Populus trichocarpa</Emphasis>

The ZRT, IRT-like protein (ZIP) family plays an important role in the transport of zinc (Zn) and iron (Fe) across the cell membrane in many different species. However, studies on ZIP family are mainly limited in herbaceous species; hence, we investigated functional divergence of ZIP family in Populus trichocarpa. We identified 21 ZIP genes in P. trichocarpa and classified them into four groups based on phylogenetic analysis. Structural analyses revealed that most of the PtrZIP transporters have eight transmembrane domains (TMDs). PtrZIP members were unequally positioned in 19 P. trichocarpa linkage groups (LGs), with six tandem duplications and four segmental duplications. The promoter regions of PtrZIP genes contain Zn, Fe, copper (Cu), and other metal stress-related cis-elements. Additionally, tissue-specific expression of PtrZIP genes showed that most of them had relatively high expression levels in the root. Quantitative real-time RT-PCR (qRT-PCR) analysis revealed that the expression of PtrZIP genes were induced not only under deficiency or excess condition of Zn, Fe, Cu and manganese (Mn) but also under excess condition of cadmium (Cd) and lead (Pb) stress. These findings indicated that PtrZIP genes may have played potential roles in metal transporters. Genome-wide analysis of PtrZIP genes in P. trichocarpa provided more comprehensive insights on the structure and function of this gene family.     

A betaine aldehyde dehydrogenase gene from <Emphasis Type="Italic">Ammopiptanthus nanus</Emphasis> enhances tolerance of <Emphasis Type="Italic">Arabidopsis</Emphasis> to high salt and drought stresses

YUCCA is an important enzyme which catalyzes a key rate-limiting step in the tryptophan-dependent pathway for auxin biosynthesis and implicated in several processes during plant growth and development. Genome wide analyses of YUCCA genes have been performed in Arabidopsis, rice, tomato, and Populus, but have never been characterized in soybean, one of the most important oil crops in the world. In this study, 22 GmYUCCA genes (GmYUCCA1-22) were identified and named based on soybean whole-genome sequence. Phylogenetic analysis of YUCCA proteins from Glycine max, Arabidopsis, Oryza sativa, tomato, and Populus euphratica revealed that GmYUCCA proteins could be divided into four subfamilies. Quantitative real-time RT-PCR (qRT-PCR) analysis showed that GmYUCCA genes have diverse expression patterns in different tissues and under various stress treatments. Compared to the wild type (WT), the transgenic GmYUCCA5 Arabidopsis plants displayed downward curling of the leaf blade margin, evident apical dominance, higher plant height, and shorter length of siliques. Our results provide a comprehensive analysis of the soybean YUCCA gene family and lay a solid foundation for further experiments in order to functionally characterize these gene members during soybean growth and development.     

Identification,characterization, and expression of the <Emphasis Type="Italic">SWEET</Emphasis> gene family in <Emphasis Type="Italic">Phalaenopsis equestris</Emphasis> and <Emphasis Type="Italic">Dendrobium officinale</Emphasis>

Sugars are important molecules that function not only as primary metabolites, but also as nutrients and signal molecules in plants. The sugar transport protein genes family SWEET has been recently identified. The availability of the Dendrobium officinale and Phalaenopsis equestris genome sequences offered the opportunity to study the SWEET gene family in this two orchid species. We identified 22 and 16 putative SWEET genes, respectively, in the genomes of D. officinale and P. equestris using comprehensive bioinformatics analysis. Based on phylogenetic comparisons with SWEET proteins from Arabidopsis and rice, the DoSWEET and PeSWEET proteins could be divided into four clades; among these, clade II specifically lacked PeSWEETs and clade IV specifically lacked DoSWEETs, and there were orthologs present between D. officinale and P. equestris. Protein sequence alignments suggest that there is a predicted serine phosphorylation site in each of the highly conserved MtN3/saliva domain regions. Gene expression analysis in four tissues showed that three PeSWEET genes were most highly expressed in the flower, leaf, stem, and root, suggesting that these genes might play important roles in growth and development in P. equestris. Analysis of gene expression in different floral organs showed that five PeSWEET genes were highly expressed in the column (gynostemium), implying their possible involvement in reproductive development in this species. The expression patterns of seven PeSWEETs in response to different abiotic stresses showed that three genes were upregulated significantly in response to high temperature and two genes were differently expressed at low temperature. The results of this study lay the foundation for further functional analysis of SWEET genes in orchids.     

Identification and polymorphism of pectinase genes <Emphasis Type="Italic">PGU</Emphasis> in the <Emphasis Type="Italic">Saccharomyces bayanus</Emphasis> complex

Pectinase (endo-polygalacturonase) is the key enzyme splitting plant pectin. The corresponding single gene PGU1 is documented for the yeast S. cerevisiae. On the basis of phylogenetic analysis of the PGU nucleotide sequence available in the GenBank, a family of divergent PGU genes is found in the species complex S. bayanus: S. bayanus var. uvarum, S. eubayanus, and hybrid taxon S. pastorianus. The PGU genes have different chromosome localization.     

First report of an atypical new <Emphasis Type="Italic">Aspergillus parasiticus</Emphasis> isolates with nucleotide insertion in <Emphasis Type="Italic">aflR</Emphasis> gene resembling to <Emphasis Type="Italic">A. sojae</Emphasis>

Aflatoxins are toxic and carcinogenic secondary metabolites produced primarily by the filamentous fungi Aspergillus flavus and Aspergillus parasiticus and cause toxin contamination in food chain worldwide. Aspergillus oryzae and Aspergillus sojae are highly valued as koji molds in the traditional preparation of fermented foods, such as miso, sake, and shoyu. Koji mold species are generally perceived of as being nontoxigenic and are generally recognized as safe (GRAS). Fungal isolates were collected from a California orchard and a few were initially identified to be A. sojae using β-tubulin gene sequences blasted against NCBI data base. These new isolates all produced aflatoxins B₁, B₂, G₁, and G₂ and were named as Pistachio Winter Experiment (PWE) strains. Thus, it is very important to further characterize these strains for food safety purposes. The full length of aflR gene of these new isolates was sequenced. Comparison of aflR DNA sequences of PWE, A. parasiticus and A. sojae, showed that the aflatoxigenic PWE strains had the six base insertion (CTCATG) similar to domesticated A. sojae, but a pre-termination codon TGA at nucleotide positions 1153–1155 was absent due to a nucleotide codon change from T to C. Colony morphology and scanning microscopic imaging of spore surfaces showed similarity of PWE strains to both A. parasiticus and A. sojae. Concordance analysis of multi locus DNA sequences indicated that PWE strains were closely linked between A. parasiticus and A. sojae. The finding documented the first report that such unique strains have been found in North America and in the world.     

Preferential retention,expression profile and potential functional diversity analysis of HD-Zip gene family in <Emphasis Type="Italic">Brassica rapa</Emphasis>

Homeodomain-Leu zipper (HD-Zip) gene family performs important biological functions related to organ development, photomorphogenesis and abiotic stress response in higher plants. However, systematic analysis of HD-Zip genes in Brassica rapa has not been performed. In the present study, a bioinformatics approach was used to identify and characterize the BraHD-Zip gene family in B. rapa. A total of 88 members were identified. All putative BraHD-Zip proteins contained a clear HD and LZ combined domain. Eighty-seven BraHD-Zips were non-randomly located on ten chromosomes. This gene family was mainly expanded following the whole genome triplication event and was preferentially over-retained relative to its neighboring genes in B. rapa. On phylogenetic analysis, the BraHD-Zips could be categorized into four distinct major groups (I–IV). Each group exhibited variant gene structures and motif distributions. Some syntenic orthologous gene pairs presented diverse expression profiles, which indicate that these gene pairs may be involved in the development of new functions during evolution. In summary, our analysis provided genome-wide insights into the expansion, preferential retention, expression profiles and functional diversity of BraHD-Zip genes following whole genome triplication in B. rapa.     

The high diversity of <Emphasis Type="Italic">Lotus corniculatus</Emphasis> endosymbionts in soils of northwest Spain

The diversity of rhizobia that establish symbiosis with Lotus corniculatus has scarcely been studied. Several species of Mesorhizobium are endosymbionts of this legume, including Mesorhizobium loti, the type species of this genus. We analysed the genetic diversity of strains nodulating Lotus corniculatus in Northwest Spain and ten different RAPD patterns were identified among 22 isolates. The phylogenetic analysis of the 16S rRNA gene showed that the isolated strains belong to four divergent phylogenetic groups within the genus Mesorhizobium. These phylogenetic groups are widely distributed worldwide and the strains nodulate L. corniculatus in several countries of Europe, America and Asia. Three of the groups include the currently described Mesorhizobium species M. loti, M. erdmanii and M. jarvisii which are L. corniculatus endosymbionts. An analysis of the recA and atpD genes showed that our strains belong to several clusters, one of them very closely related to M. jarvisii and the remanining ones phylogenetically divergent from all currently described Mesorhizobium species. Some of these clusters include L. corniculatus nodulating strains isolated in Europe, America and Asia, although the recA and atpD genes have been sequenced in only a few L. corniculatus endosymbionts. The results of this study revealed great phylogenetic diversity of strains nodulating L. corniculatus, allowing us to predict that even more diversity will be discovered as further ecosystems are investigated.     

Fine mapping of a <Emphasis Type="Italic">Phytophthora</Emphasis>-resistance gene <Emphasis Type="Italic">RpsWY</Emphasis> in soybean (<Emphasis Type="Italic">Glycine max</Emphasis> L.) by high-throughput genome-wide sequencing

Key message

Using a combination of phenotypic screening, genetic and statistical analyses, and high-throughput genome-wide sequencing, we have finely mapped a dominant Phytophthora resistance gene in soybean cultivar Wayao.

Abstract

Phytophthora root rot (PRR) caused by Phytophthora sojae is one of the most important soil-borne diseases in many soybean-production regions in the world. Identification of resistant gene(s) and incorporating them into elite varieties are an effective way for breeding to prevent soybean from being harmed by this disease. Two soybean populations of 191 F₂ individuals and 196 F_7:8 recombinant inbred lines (RILs) were developed to map Rps gene by crossing a susceptible cultivar Huachun 2 with the resistant cultivar Wayao. Genetic analysis of the F₂ population indicated that PRR resistance in Wayao was controlled by a single dominant gene, temporarily named RpsWY, which was mapped on chromosome 3. A high-density genetic linkage bin map was constructed using 3469 recombination bins of the RILs to explore the candidate genes by the high-throughput genome-wide sequencing. The results of genotypic analysis showed that the RpsWY gene was located in bin 401 between 4466230 and 4502773 bp on chromosome 3 through line 71 and 100 of the RILs. Four predicted genes (Glyma03g04350, Glyma03g04360, Glyma03g04370, and Glyma03g04380) were found at the narrowed region of 36.5 kb in bin 401. These results suggest that the high-throughput genome-wide resequencing is an effective method to fine map PRR candidate genes.