SNPs in stress-responsive rice genes: validation, genotyping, functional relevance and population structure

ABSTRACT: BACKGROUND: Single nucleotide polymorphism (SNP) validation and large-scale genotyping are required to maximize the use of DNA sequence variation and determine the functional relevance of candidate genes for complex stress tolerance traits through genetic association in rice. We used the bead array platform-based Illumina GoldenGate assay to validate and genotype SNPs in a select set of stress-responsive genes to understand their functional relevance and study the population structure in rice. RESULTS: Of the 384 putative SNPs assayed, we successfully validated and genotyped 362 (94.3%). Of these 325 (84.6%) showed polymorphism among the 91 rice genotypes examined. Physical distribution, degree of allele sharing, admixtures and introgression, and amino acid replacement of SNPs in 263 abiotic and 62 biotic stress-responsive genes provided clues for identification and targeted mapping of trait-associated genomic regions. We assessed the functional and adaptive significance of validated SNPs in a set of contrasting drought tolerant upland and sensitive lowland rice genotypes by correlating their allelic variation with amino acid sequence alterations in catalytic domains and three-dimensional secondary protein structure encoded by stress-responsive genes. We found a strong genetic association among SNPs in the nine stress-responsive genes with upland and lowland ecological adaptation. Higher nucleotide diversity was observed in indica accessions compared with other rice sub-populations based on different population genetic parameters. The inferred ancestry of 16% among rice genotypes was derived from admixed populations with the maximum between upland aus and wild Oryza species. CONCLUSIONS: SNPs validated in biotic and abiotic stress-responsive rice genes can be used in association analyses to identify candidate genes and develop functional markers for stress tolerance in rice.     

Homoeology and phylogeny of the A,S, and D genomes of the triticinae

Summary Two-dimensional gel electrophoresis of seedling proteins has been undertaken on diploid representatives of the A, S, and DTriticum genomes. The 457 reproducible spots of the hexaploid Chinese Spring were taken into account in the comparisons. More than 82% of the hexaploid spots were found among the representatives of each of the three genomes. Thus, it was proposed that the homoeologous structural loci, on the three homoeologous genomes, can be present in the same allelic forms and that from 72% to 86% of the spots of the hexaploid pattern are in two or three doses. The dendrogram, deduced from the similarity indices computed between the 16 diploid accessions, was very consistent with our present knowledge of the genetic relationships in theTriticum genus. This result, in addition to previous ones, shows that the comparison of two-dimensional protein patterns represents a novel experimental approach for studying the phylogenetic relationships in theTriticinae.     

利用SSR标记分析云南、西藏和新疆小麦的遗传多样性

用185对SSR引物对52份中国西部特有小麦的遗传多样性进行了研究分析。在31份云南小麦材料中,共检测到488个等位变异,每一个SSR引物可检测到1至9个等位变异,平均为2.64个;平均PIC值为0.2764。在15份西藏小麦材料中,共检测到472个等位变异,每个引物可扩增出1到8个等位变异,平均为2.55个;平均PIC值为0.3082。在6份新疆小麦材料中,共检测到308个等位变异,每一个SSR引物可检测1到5个等位变异,平均为1.66个;平均PIC值为0.1944。185对SSR引物在云南、西藏和新疆小麦的21条染色体、7个部分同源群和3个染色体组上检测到的等位位点的多态性存在明显差异。云南、西藏和新疆小麦均以3B染色体较高,而1D染色体最低;在7个部分同源群中,均以第三部分同源群最高,第六部分同源群最低;在A、B和D染色体组上,均以B染色体组最高,D染色体组最低,A染色体组居中。利用185对SSR引物计算了云南、西藏和新疆小麦群体内及其群体间的遗传距离(GD)和平均遗传距离,结果显示,西藏小麦和云南小麦群体内的平均遗传距离要高于新疆小麦,而云南小麦和西藏小麦间的平均遗传距离低于两者与新疆小麦的平均遗传距离。聚类分析结果也表明,云南小麦和西藏小麦的亲缘关系较近,但两者与新疆小麦的亲缘关系相对较远。     

Identification of homologous, homoeologous and paralogous sequence variants in an outbreeding allopolyploid species based on comparison with progenitor taxa

The combination of homologous, homoeologous and paralogous classes of sequence variation presents major challenges for SNP discovery in outbreeding allopolyploid species. Previous in vitro gene-associated SNP discovery studies in the allotetraploid forage legume white clover (Trifolium repens L.) were vulnerable to such effects, leading to prohibitive levels of attrition during SNP validation. Identification of T. occidentale and T. pallescens as the putative diploid progenitors of white clover has permitted discrimination of the different sequence variant categories. Amplicons from selected abiotic stress tolerance-related genes were obtained using mapping family parents and individuals from each diploid species. Following cloning, progenitor comparison allowed tentative assignment of individual haplotypes to one or other sub-genome, as well as to gene copies within sub-genomes. A high degree of coincidence and identity between SNPs and HSVs was observed. Close similarity was observed between the genome of T. occidentale and one white clover sub-genome, but the affinity between T. pallescens and the other sub-genome was weaker, suggesting that a currently uncharacterised taxon may be the true second progenitor. Selected validated SNPs were attributed to individual sub-genomes by assignment to and naming of homoeologous linkage groups, providing the basis for improved genetic trait-dissection studies. The approach described in this study is broadly applicable to a range of allopolyploid taxa of equivocal ancestry.     

Mutation Scanning in Wheat by Exon Capture and Next-Generation Sequencing

Targeted Induced Local Lesions in Genomes (TILLING) is a reverse genetics approach to identify novel sequence variation in genomes, with the aims of investigating gene function and/or developing useful alleles for breeding. Despite recent advances in wheat genomics, most current TILLING methods are low to medium in throughput, being based on PCR amplification of the target genes. We performed a pilot-scale evaluation of TILLING in wheat by next-generation sequencing through exon capture. An oligonucleotide-based enrichment array covering ~2 Mbp of wheat coding sequence was used to carry out exon capture and sequencing on three mutagenised lines of wheat containing previously-identified mutations in the TaGA20ox1 homoeologous genes. After testing different mapping algorithms and settings, candidate SNPs were identified by mapping to the IWGSC wheat Chromosome Survey Sequences. Where sequence data for all three homoeologues were found in the reference, mutant calls were unambiguous; however, where the reference lacked one or two of the homoeologues, captured reads from these genes were mis-mapped to other homoeologues, resulting either in dilution of the variant allele frequency or assignment of mutations to the wrong homoeologue. Competitive PCR assays were used to validate the putative SNPs and estimate cut-off levels for SNP filtering. At least 464 high-confidence SNPs were detected across the three mutagenized lines, including the three known alleles in TaGA20ox1, indicating a mutation rate of ~35 SNPs per Mb, similar to that estimated by PCR-based TILLING. This demonstrates the feasibility of using exon capture for genome re-sequencing as a method of mutation detection in polyploid wheat, but accurate mutation calling will require an improved genomic reference with more comprehensive coverage of homoeologues.     

Sequence variation in two lignin biosynthesis genes,cinnamoyl CoA reductase (CCR) and cinnamyl alcohol dehydrogenase 2 (CAD2)

Cinnamoyl CoA reductase (CCR) and cinnamyl alcohol dehydrogenase 2 (CAD2) are genes which may influence variation in lignin content and composition within plants. Sequence variation within these genes may be responsible for changes in enzyme activity and/or specificity, which could cause variation in lignin content or composition. This study examines sequence variation within these two genes in Eucalyptus globulus, an important species used in pulp and paper-making. Twenty-one single nucleotide polymorphisms (SNPs) were identified in the exons of CCR, of which nine were neutral mutations and 12 were missense mutations. Six of the missense mutations affected highly conserved amino acids within the protein sequence of CCR. Eight SNPs were identified in the CAD2 exons, six of which were neutral mutations and two which were missense mutations. One of the missense mutations affected a highly conserved amino acid within the protein sequence. In addition, 32 SNPs were identified in the CCR introns along with four insertion/deletions and two polyA length variation regions. Polymorphism affecting highly conserved amino acids may alter enzyme function and this molecular variation may be linked to variation in lignin profiles. Selecting positive alleles which produce favourable lignin profiles would be advantageous in tree breeding programs.     

A SNARE-Like Protein and Biotin Are Implicated in Soybean Cyst Nematode Virulence

Phytoparasitic nematodes that are able to infect and reproduce on plants that are considered resistant are referred to as virulent. The mechanism(s) that virulent nematodes employ to evade or suppress host plant defenses are not well understood. Here we report the use of a genetic strategy (allelic imbalance analysis) to associate single nucleotide polymorphisms (SNPs) with nematode virulence genes in Heterodera glycines, the soybean cyst nematode (SCN). To accomplish this analysis, a custom SCN SNP array was developed and used to genotype SCN F3-derived populations grown on resistant and susceptible soybean plants. Three SNPs reproducibly showed allele imbalances between nematodes grown on resistant and susceptible plants. Two candidate SCN virulence genes that were tightly linked to the SNPs were identified. One SCN gene encoded biotin synthase (HgBioB), and the other encoded a bacterial-like protein containing a putative SNARE domain (HgSLP-1). The two genes mapped to two different linkage groups. HgBioB contained sequence polymorphisms between avirulent and virulent nematodes. However, the gene encoding HgSLP-1 had reduced copy number in virulent nematode populations and appears to produce multiple forms of the protein via intron retention and alternative splicing. We show that HgSLP-1 is an esophageal-gland protein that is secreted by the nematode during plant parasitism. Furthermore, in bacterial co-expression experiments, HgSLP-1 co-purified with the SCN resistance protein Rhg1 α-SNAP, suggesting that these two proteins physically interact. Collectively our data suggest that multiple SCN genes are involved in SCN virulence, and that HgSLP-1 may function as an avirulence protein and when absent it helps SCN evade host defenses.     

Genetic mapping of the wheat dimeric α-amylase inhibitor multi-gene family using allele-specific primers based on intergenomic SNPs

Single nucleotide polymorphisms (SNPs) identified in EST sequences can be used to map expressed genes. Though SNPs are useful markers for genetic mapping, SNP mapping of genes in common wheat is challenging because the genetic complement of wheat consists of three closely related genomes (designated A, B, and D), and most genes are present in triplicate sets. Mapping multi-gene family members is further complicated by the fact that it is difficult to distinguish SNP differences between the various paralogs from those between the different genomes. We have developed a PCR-based method for assigning wheat EST sequences to their proper genetic loci by first identifying and mapping SNPs that distinguish the three genomes. To develop this method, we focused on EST sequences encoding the dimeric α-amylase inhibitors (WDAI), The WDAI coding regions of hexaploid wheat were aligned. The sequences were classified into three groups based on nucleotide variations. Twenty-two SNPs were identified that distinguish the three groups. Group-specific primers based on these SNPs were designed to permit selective amplification of each group. The chromosomal location of each group was then determined using Group 3 ditelosomic lines of Chinese Spring. Groups 1 and 2 were assigned to chromosome locations 3DS and 3BS, respectively, whereas no sequence could be assigned to 3AS. A remarkable feature of this method is the ability to discriminate the location of homoeologous multigenes in the three genomes of wheat. This strategy can be useful for assigning unknown wheat EST sequences to specific chromosomes.     

Large-scale profiling and identification of potential regulatory mechanisms for allelic gene expression in colorectal cancer cells

Allelic variation in gene expression is common in humans and this variation is associated with phenotypic variation. In this study, we employed high-density single nucleotide polymorphism (SNP) chips containing 13,900 exonic SNPs to identify genes with allelic gene expression in cells from colorectal cancer cell lines. We found 2 monoallelically expressed genes (ERAP2 and MYLK4), 32 genes with an allelic imbalance in their expression, and 13 genes showing allele substitution by RNA editing. Among a total of 34 allelically expressed genes in colorectal cancer cells, 15 genes (44.1%) were associated with cis-acting eQTL, indicating that large portions of allelically expressed genes are regulated by cis-acting mechanisms of gene expression. In addition, potential regulatory variants present in the proximal promoter regions of genes showing either monoallelic expression or allelic imbalance were not tightly linked with coding SNPs, which were detected with allelic gene expression. These results suggest that multiple rare variants could be involved in the cis-acting regulatory mechanism of allelic gene expression. In the comparison with allelic gene expression data from Centre d'Etude du Polymorphisme Humain (CEPH) family B cells, 12 genes showed B-cell specific allelic imbalance and 1 noncoding SNP showed colorectal cancer cell-specific allelic imbalance. In addition, different patterns of allele substitution were observed between B cells and colorectal cancer cells. Overall, our study not only indicates that allelic gene expression is common in colorectal cancer cells, but our study also provides a better understanding of allele-specific gene expression in colorectal cancer cells.     

Identification of two new alleles, IGHV3-23*04 and IGHJ6*04, and the complete sequence of the IGHV3-h pseudogene in the human immunoglobulin locus and their prevalences in Danish Caucasians

More than 100 variable (V), 27 diversity (D), and six joining (J) genes are encoded in the human heavy chain locus, and many of these genes exists in different allelic forms. The number of genes and the allelic differences help to create diversity in the immunoglobulin receptors, a key feature of the adaptive immune system. We here report the identification of two novel and seemingly functional alleles of human heavy chain genes. The variable IGHV3-23*04 allele is found with an allele frequency of 0.21 amongst Danish Caucasians, whereas the novel joining IGHJ6*04 allele is rare (allele frequency 0.02). We also report the full sequence of IGHV3-h. The gene exists in two allelic forms but is only found in 58% of the Danish Caucasians studied. The methionine translation initiation codon is mutated, ATG→AAG, and we therefore propose that the gene is a pseudogene incapable of being translated.     

Isolation and characterization of SNP variation at 90 anonymous loci in the banded wren (<Emphasis Type="Italic">Thryothorus pleurostictus)</Emphasis>

Single nucleotide polymorphisms (SNPs) are becoming more commonly used as molecular markers in conservation studies. However, relatively few studies have employed SNPs for species with little or no existing sequence data, partly due to the practical challenge of locating appropriate SNP loci in these species. Here we describe an application of SNP discovery via shotgun cloning that requires no pre-existing sequence data and is readily applied to all taxa. Using this method, we isolated, cloned and screened for SNP variation at 90 anonymous sequence loci (51 kb total) from the banded wren (Thryothorus pleurostictus), a Central American species with minimal pre-existing sequence data and a documented paucity of microsatellite allelic variation. We identified 168 SNPs (a mean of one SNP/305 bp, with SNPs unevenly distributed across loci). Further characterization of variation at 41 of these SNP loci among 256 individuals including 37 parent–offspring families suggests that they provide substantial information for defining the genetic mating system of this species, and that SNPs may be generally useful for this purpose when other markers are problematic.     

Allelic variation in gene expression identified through computational analysis of the dbEST database

Differential expression between the two alleles of an individual and between people with different genotypes has been commonly observed. Quantitative differences in gene expression between people may provide the genetic basis for the phenotypic difference between individuals and may be the primary cause of complex diseases. In this paper, we developed a computational method to identify genes that displayed allelic variation in gene expression in human EST libraries. To model allele-specific gene expression, we first identified EST libraries in which both A and B alleles were expressed and then identified allelic variation in gene expression based on the EST counts for each allele using a binomial test. Among 1107 SNPs that had a sufficient number of ESTs for the analysis, 524 (47%) displayed allelic variation in at least one cDNA library. We verified experimentally the allelic variation in gene expression for 6 of these SNPs. The frequency of allelic variation observed in EST libraries was similar to the previous studies using the SNP chip and primer extension method. We found that genes that displayed allelic variation were distributed throughout the human genome and were enriched in certain chromosome regions. The SNPs and genes identified in this study will provide a rich source for evaluating the effects of those SNPs and associated haplotypes in human health and diseases.     

Identification of SNPs in Interferon Gamma,Interleukin-22, and Their Receptors and Associations with Health and Production-Related Traits in Canadian Holstein Bulls

Genetic variants in a number of immunoregulatory genes have been previously associated with health and production traits in dairy cattle. Therefore, in the following study, the genes coding interferon gamma (IFNG), IFNG receptor 1 and 2 domains, interleukin-22 (IL22), and IL22 receptor alpha 1, were investigated for single nucleotide polymorphisms (SNPs) in Holstein bulls. These SNPs, along with SNPs previously identified in IL10, IL10 receptor, and transforming growth factor beta 1 (TGFB1) genes, were evaluated for statistical associations to estimated breeding values for milk somatic cell score (SCS), a trait highly correlated to mastitis incidence, and various production-related traits, including milk yield, protein yield, fat yield, and lactation persistency. While no significant associations were found between these SNPs and SCS, SNPs in IL10 receptor beta subunit showed a significant effect on protein yield and lactation persistency. While there is evidence that IL10 plays an important role during lactation, it is also likely that the effects of SNPs in IL10 receptor beta subunit on protein yield and lactation persistency are due to linkage disequilibrium with a neighboring QTL.