Development Of Single Nucleotide Polymorphisms Research Articles

PSXII-10 Development of a 70k Snp Genotyping Array for American Mink (Neogale vison)

Abstract Mink is an important fur animal species; however, little is known about the genomics of complex traits in mink. The development of single nucleotide polymorphism (SNP) array has accelerated genomics application in animal breeding by facilitating genotyping and associated analyses. The current project aimed to design a medium-density SNP array in American mink for genomic studies and selection in this species. For this purpose, whole genome sequence with 30x coverage was generated for 100 individuals from the Canadian Centre for Fur Animal Research (CCFAR) at the Dalhousie Faculty of Agriculture (Truro, NS, Canada) and Millbank Fur Farm (Rockwood, ON, Canada). Variant calling was done using Samtools and GATK pipelines. A total of 8,373,854 bi-allelic SNPs identified by both pipelines were used for an initial selection of 91K SNPs for submission to the Affymetrix scoring system. After filtering, 62,375 SNPs were selected for probe design and production of an Axiom Mink Genotyping Array 70K. This 70K SNP array was used to genotype 2,973 mink for validation. After quality control of genotypes, 47,159 SNPs were classified as high-resolution SNPs, and among them 24,237 were polymorphic. A preliminary genome-wide association study of 1,356 mink for body length at harvest identified a SNP (AX-647623051) significantly associated (P < 1.26 × 10-5) with the trait. In conclusion, this 70K mink SNP array could be used for genomic studies and potentially as a starting point for genomic selection programs in American mink.

Development of single nucleotide polymorphism (SNP) application for detection and genotyping of RSIV-type megalocytiviruses.

Red sea bream iridovirus (RSIV) belonging to the genus Megalocytivirus of the family Iridoviridae is the cause of serious mass mortality of cultured marine fishes. RSIV-type megalocytiviruses show extremely high nucleotide sequence identities. Thus, epidemiological studies on this virus are limited. This study developed two primer sets amplifying the regions possessing single nucleotide polymorphism (SNP) to determine the relationships and divergence of RSIV-type megalocytiviruses isolated from cultured marine fishes in Japan. The two regions were designed according to the genome sequences of the representative RSIV genotype II of megalocytivirus members in GenBank. The SNP 1 and 2 regions have sequences homologous to hypothetical protein ORF 24 and ORF 31, respectively, of RSIV (accession no. AP017456.1). By sequencing the regions, 53 polymorphic sites were identified. The phylogenetic analysis of 25 RSIV-type megalocytivirus isolates, classified into RSIV cluster, was clustered into eight haplotypes (seven haplotypes from Oita, two haplotypes from Ehime, and one haplotype shared between Oita and Ehime). These findings suggested that SNP in the RSIV genome is a powerful application for the detection and identification of RSIV-type megalocytiviruses.

384 Genetic and Genomic Analysis in Livestock with Increasing Datasets

Abstract In animal breeding and genetics, statistical methods have been used for several decades to identify animals with the best genetic potential. One of the foundations for computing accurate estimated breeding values (EBV) is the amount of data that is used in the evaluation system — as the more data points one animal has, the more accurate its EBV is going to be. However, the animal breeding and genetics field periodically faces a big data paradox, where efficient methods have to be developed to handle the amount of data collected over time, given the computing capacity becomes the limiting factor. For instance, running genetic evaluations based on phenotypes and pedigree for a million animals was impossible in 1970. Methods and algorithms evolved to a point where using data for millions of animals was not a problem, until genomic information became available. After the development of single nucleotide polymorphism (SNP) chips for livestock in 2008, genomic information started being used in addition to phenotypes and pedigree to further improve accuracy of EBV. However, each animal is genotyped for around 50,000 SNP, which makes this data dense and difficult to work with. Over 790,000 Angus and 3.4 million Holstein cattle have been genotyped in the US as of April 2020. As the amount of new data considerably increases every week, most of the genomic evaluations are done on a weekly basis. Given that data have to be processed, the computation of EBV cannot take more than four days, which can be challenging depending on the model. In this talk we will discuss the challenges and solutions for successful genomic evaluations in large livestock populations. Finally, perspectives on the use of whole-genome sequence data and high-throughput phenotypes in genomic analysis will be summarized.

Development of Single Nucleotide Polymorphism (SNP) Markers in Tropical Crops

Understanding genetic diversity, association studies, evolution analysis, quantitative trait loci, marker-assisted selection and genome-wide association in tropical crops are important for improving plant characteristics in order to increase food sustainability in tropical countries. Single nucleotide polymorphism (SNP) marker is becoming the most popular molecular marker for those studies. By using SNP marker, genes associated with important traits can be identified efficiently compared to the other molecular markers. This review describes about how SNP can be discovered in the plant genomes and the application of SNP in plant breeding, especially in tropical crops such as rice, maize, peas, potato, tomato, cassava, taro, etc. Keywords: food sustainability, plant breeding, SNP marker, tropical crops

Development of Single Nucleotide Polymorphism (SNP) Markers for Analysis of Population Structure and Invasion Pathway in the Coconut Leaf Beetle Brontispa longissima (Gestro) Using Restriction Site-Associated DNA (RAD) Genotyping in Southern China

To determine population genomic structure through high-throughput sequencing techniques has revolutionized research on non-model organisms. The coconut leaf beetle, Brontispa longissima (Gestro), is a widely distributed pest in Southern China. Here, we used restriction site-associated DNA (RAD) genotyping to assess the invasion pathway by detecting and estimating the degree of genetic differentiation among 51 B. longissima accessions collected from Southern China. A total of 10,127 SNPs were obtained, the screened single nucleotide polymorphism (SNP) information was used to construct the phylogenetic tree, FST analysis, principal component analysis, and population structure analysis. Genetic structure analysis was used to infer the population structure; the result showed that all accessions were divided into Hainan population and non-Hainan population. The Hainan population remained stable, only the Sansha population differentiated, and the non-Hainan populations have gradually differentiated into smaller sub-populations. We concluded that there are two sources of invasion of B. longissima into mainland China: Taiwan and Hainan. With the increase of the invasion time, the Hainan population was relatively stable, and the Taiwan population was differentiated into three sub-populations. Based on the unrooted phylogenetic tree, we infer that Taiwan and Hainan are the two invasive base points. The Taiwan population invaded Fujian, Guangdong, and Guangxi, while the Hainan population invaded Yunnan and Sansha. Our results provide strong evidence for the utility of RAD sequencing (RAD-seq) in population genetics studies, and our generated SNP resource could provide a valuable tool for population genomics studies of B. longissima in the future.

Development of a set of SNP markers for population genetics of the red gorgonian (Paramuricea clavata), an emblematic species of the Mediterranean coralligenous

Transcriptome sequencing was used for the development of single nucleotide polymorphisms (SNP) for the red gorgonian (Paramuricea clavata). A total of 20,736 SNPs were identified, and 1718 had a coverage of over 100 reads. Of the 480 SNPs tested, 347 SNPs were successfully genotyped at 95 samples from the NW Mediterranean using a MassARRAY System. This set of markers will be of great value for population genetics and phylogeography.

Development of a set of SNP markers for population genetics studies of Ipe (Handroanthus sp.), a valuable tree genus from Latin America

A combination of restriction associated DNA sequencing (RADSeq) and low coverage MiSeq genome sequencing was used for the development of single nucleotide polymorphisms (SNP) and INDEL (insertion/deletions) genetic markers for Ipe (Handroanthus sp.). Of the 402 putative loci identified, 389 SNPs and INDELs (315 nuclear SPNs, six chloroplast INDELs, 15 chloroplast SNPs, 12 mitochondrial INDELs and 41 mitochondrial SNPs) were successfully genotyped at 93 individuals from Brazil, Bolivia and French Guiana using a MassARRAY® iPLEX™ platform. This set of markers will be invaluable for population genetics, phylogeography and DNA fingerprinting studies

Dealing with paralogy in RADseq data: in silico detection and single nucleotide polymorphism validation in Robinia pseudoacacia L.

The RADseq technology allows researchers to efficiently develop thousands of polymorphic loci across multiple individuals with little or no prior information on the genome. However, many questions remain about the biases inherent to this technology. Notably, sequence misalignments arising from paralogy may affect the development of single nucleotide polymorphism (SNP) markers and the estimation of genetic diversity. We evaluated the impact of putative paralog loci on genetic diversity estimation during the development of SNPs from a RADseq dataset for the nonmodel tree species Robinia pseudoacacia L. We sequenced nine genotypes and analyzed the frequency of putative paralogous RAD loci as a function of both the depth of coverage and the mismatch threshold allowed between loci. Putative paralogy was detected in a very variable number of loci, from 1% to more than 20%, with the depth of coverage having a major influence on the result. Putative paralogy artificially increased the observed degree of polymorphism and resulting estimates of diversity. The choice of the depth of coverage also affected diversity estimation and SNP validation: A low threshold decreased the chances of detecting minor alleles while a high threshold increased allelic dropout. SNP validation was better for the low threshold (4×) than for the high threshold (18×) we tested. Using the strategy developed here, we were able to validate more than 80% of the SNPs tested by means of individual genotyping, resulting in a readily usable set of 330 SNPs, suitable for use in population genetics applications.

Development of single nucleotide polymorphism markers in the large and complex rubber tree genome using next-generation sequence data

The development of single nucleotide polymorphism (SNP) markers provides the opportunity to improve many areas of plant breeding and population genetics. Unfortunately, for species such as the rubber tree (Hevea brasiliensis), the use of next-generation sequencing for genomic SNP discovery is very difficult because of the large genome size and the abundance of repeated sequences. Access to a set of validated SNP markers is a significant advantage for rubber researchers who wish to apply SNPs in scientific research. Here, we performed genomic sequencing of H. brasiliensis and generated 10,993,648 short reads, which were assembled into 10,071 contigs (N50 = 3078) by a de novo assembly strategy. A total of 2446 contigs presented no hits in the current H. brasiliensis genome assembly and may therefore be considered novel genomic sequences of rubber tree. A total of 143 putative polymorphic positions were selected, gene annotations were available for 58.7 % of the markers, and all of the sequences could be anchored to the released H. brasiliensis genome. These SNPs were validated in eight genotypes of H. brasiliensis and 15 F1 plants from a mapping population, resulting in 30 (20.9 %) positions correctly classified. The analysis revealed key candidate genes responsible for defence mechanisms and provided markers for further genetic improvement of Hevea in breeding programmes.

Scrimer: designing primers from transcriptome data.

With the rise of next-generation sequencing methods, it has become increasingly possible to obtain genomewide sequence data even for nonmodel species. Such data are often used for the development of single nucleotide polymorphism (SNP) markers, which can subsequently be screened in a larger population sample using a variety of genotyping techniques. Many of these techniques require appropriate locus-specific PCR and genotyping primers. Currently, there is no publicly available software for the automated design of suitable PCR and genotyping primers from next-generation sequence data. Here we present a pipeline called Scrimer that automates multiple steps, including adaptor removal, read mapping, selection of SNPs and multiple primer design from transcriptome data. The designed primers can be used in conjunction with several widely used genotyping methods such as SNaPshot or MALDI-TOF genotyping. Scrimer is composed of several reusable modules and an interactive bash workflow that connects these modules. Even the basic steps are presented, so the workflow can be executed in a step-by-step manner. The use of standard formats throughout the pipeline allows data from various sources to be plugged in, as well as easy inspection of intermediate results with visualization tools of the user's choice.

Genotyping and development of single-nucleotide polymorphism (SNP) markers associated with blast resistance genes in rice using GoldenGate assay

Development and large-scale genotyping of single-nucleotide polymorphism (SNP) is required to use identified sequence variation in the alleles of different genes to determine their functional relevance to the candidate gene(s). In the present study, Illumina GoldenGate assay was used to validate and genotype SNPs in a set of six major rice blast resistance genes, viz. Pi-ta, Piz(t), Pi54, Pi9, Pi5(1) and Pib, distributed over five chromosomes, to understand their functional relevance and study the population structure in rice. All the selected SNPs loci (96) of six blast (Magnaporthe oryzae) resistance genes were genotyped successfully in 92 rice lines with an overall genotype call rate of 92.0 % and minimum GenTrain cutoff score of ≥0.448. The highest genotyped SNPs were found in japonica type (97.1 %) rice lines, followed by indica (92.12 %), indica basmati (91.84 %) and minimum in case of wild species (82.0 %). Among the genotyped loci, the highest score (98.68 %) was observed in case of Piz(t), followed by Pi-ta, Pi5(1), Pib, Pi54 and Pi9. Polymorphism was obtained in 87.5 % SNPs loci producing 7,728 genotype calls. Minor allele frequency ranged from 0.01 to 0.49 and has good differentiating power for distinguishing different rice accessions. Population structure analysis revealed that a set of genotypes from four rice subpopulations had “admix” ancestry (>26 %) with more than one genetic background of indica, japonica and wild types. SNPs markers were validated in a set of 92 rice lines and converted into CAPS markers which can be used in blast resistance breeding programme.

English

Accurate identification of individual genotypes is important for cacao ( Theobroma cacao L.) breeding, germplasm conservation and seed propagation. The development of single nucleotide polymorphism (SNP) markers in cacao offers an effective way to use a high-throughput genotyping system for cacao genotype verification. In the present study, high-throughput genotyping with SNP markers was used to fingerprint 160 cacao trees in the germplasm collection at the Cocoa Research Institute of Ghana (CRIG). These accessions had been originally introduced from international germplasm collections. The multilocus SNP profiles, generated by the Sequenom Mass Spectrometry platform, were compared with the SNP profiles of reference trees maintained in the international cacao collections. The comparison unambiguously identified mislabeled trees. For materials introduced as hybrid seeds without an available reference genotype, parentage analysis and model-based assignment were applied to verify their recorded parentage and genetic background. Our study shows that a small set of polymorphic SNP markers can provide a robust and accurate result for cacao genotype identification. This protocol can be applied for large-scale genotyping of cacao as well as for many other crops. Keywords: Cacao, conservation, chocolate, DNA fingerprint, molecular marker, tropical plant, off-type, true-to-type, West Africa. African Journal of Biotechnology , Vol 13(21), 2127-2136

Transcriptome-based single nucleotide polymorphism markers for genome mapping in Japanese pear (Pyrus pyrifolia Nakai)

The development of single nucleotide polymorphism (SNP) markers in Japanese pear (Pyrus pyrifolia Nakai) offers the opportunity to use DNA markers for marker-assisted selection in breeding programs because of their high abundance, codominant inheritance, and potential for automated high-throughput analysis. We developed a 1,536-SNP bead array without a reference genome sequence from more than 44,000 base changes on the basis of a large-scale expressed sequence tag (EST) analysis combined with 454 genome sequencing data of Japanese pear ‘Housui’. Among the 1,536 SNPs on the array, 756 SNPs were genotyped, and 609 SNP loci were mapped to linkage groups on a genetic linkage map of ‘Housui’, based on progeny of an interspecific cross between European pear (Pyrus communis L.) ‘Bartlett’ and ‘Housui’. The newly constructed genetic linkage map consists of 951 loci, comprising 609 new SNPs, 110 pear genomic simple sequence repeats (SSRs), 25 pear EST–SSRs, 127 apple SSRs, 61 pear SNPs identified by the “potential intron polymorphism” method, and 19 other loci. The map covers 22 linkage groups spanning 1341.9 cM with an average distance of 1.41 cM between markers and is anchored to reference genetic linkage maps of European pears and apples. A total of 514 contigs containing mapped SNP loci showed significant similarity to known proteins by functional annotation analysis.

Genome-wide transcriptome study in wheat identified candidate genes related to processing quality, majority of them showing interaction (quality x development) and having temporal and spatial distributions

BackgroundThe cultivated bread wheat (Triticum aestivum L.) possesses unique flour quality, which can be processed into many end-use food products such as bread, pasta, chapatti (unleavened flat bread), biscuit, etc. The present wheat varieties require improvement in processing quality to meet the increasing demand of better quality food products. However, processing quality is very complex and controlled by many genes, which have not been completely explored. To identify the candidate genes whose expressions changed due to variation in processing quality and interaction (quality x development), genome-wide transcriptome studies were performed in two sets of diverse Indian wheat varieties differing for chapatti quality. It is also important to understand the temporal and spatial distributions of their expressions for designing tissue and growth specific functional genomics experiments.ResultsGene-specific two-way ANOVA analysis of expression of about 55 K transcripts in two diverse sets of Indian wheat varieties for chapatti quality at three seed developmental stages identified 236 differentially expressed probe sets (10-fold). Out of 236, 110 probe sets were identified for chapatti quality. Many processing quality related key genes such as glutenin and gliadins, puroindolines, grain softness protein, alpha and beta amylases, proteases, were identified, and many other candidate genes related to cellular and molecular functions were also identified. The ANOVA analysis revealed that the expression of 56 of 110 probe sets was involved in interaction (quality x development). Majority of the probe sets showed differential expression at early stage of seed development i.e. temporal expression. Meta-analysis revealed that the majority of the genes expressed in one or a few growth stages indicating spatial distribution of their expressions. The differential expressions of a few candidate genes such as pre-alpha/beta-gliadin and gamma gliadin were validated by RT-PCR. Therefore, this study identified several quality related key genes including many other genes, their interactions (quality x development) and temporal and spatial distributions.ConclusionsThe candidate genes identified for processing quality and information on temporal and spatial distributions of their expressions would be useful for designing wheat improvement programs for processing quality either by changing their expression or development of single nucleotide polymorphisms (SNPs) markers.

A simple strategy for development of single nucleotide polymorphisms from non-model species and its application in Panax.

Single nucleotide polymorphisms (SNPs) are widely employed in the studies of population genetics, molecular breeding and conservation genetics. In this study, we explored a simple route to develop SNPs from non-model species based on screening the library of single copy nuclear genes (SCNGs). Through application of this strategy in Panax, we identified 160 and 171 SNPs from P. quinquefolium and P. ginseng, respectively. Our results demonstrated that both P. ginseng and P. quinquefolium possessed a high level of nucleotide diversity. The number of haplotype per locus ranged from 1 to 12 for P. ginseng and from 1 to 9 for P. quinquefolium, respectively. The nucleotide diversity of total sites (πT) varied between 0.000 and 0.023 for P. ginseng and 0.000 and 0.035 for P. quinquefolium, respectively. These findings suggested that this approach is well suited for SNP discovery in non-model organisms and is easily employed in standard genetics laboratory studies.

Construction of high-density genetic linkage maps for orange-spotted grouper Epinephelus coioidesusing multiplexed shotgun genotyping

BackgroundOrange-spotted grouper, Epinephelus coioides, is one of the most valuable fish species in China. Commercial production of orange-spotted grouper could be increased by developing higher growth rates and improving commercially important traits. Information on genetic markers associated with quantitative trait loci (QTL) can be used in breeding programs to identify and select individuals carrying desired traits. A high-density genetic linkage map is the basis for QTL study, and multiplexed shotgun genotyping (MSG) facilitates the development of single nucleotide polymorphisms (SNPs) and genotyping. In this study, the first high-density genetic linkage maps for groupers were generated on the basis of the MSG method.ResultsThe sex-averaged map contained a total of 4,608 SNPs, which spanned 1581.7 cM, with a mean distance between SNPs of 0.34 cM. The 4,608 SNPs were located in 2,849 unique locations on the linkage map, with an average inter-location space at 0.56 cM. There were 2,516 SNPs on the female map, and the number of unique locus was 1,902. However, the male map contained more numbers of SNP (2,939) and unique locations (2,005). The total length of the female and male maps was 1,370.9 and 1,335.5 cM, respectively.ConclusionsThe high-resolution genetic linkage maps will be very useful for QTL analyses and marker-assisted selection (MAS) for economically important traits in molecular breeding of the orange-spotted grouper.

Development of single nucleotide polymorphisms in Phaseolus vulgaris and related Phaseolus spp

In this study, new single nucleotide polymorphism (SNP) markers were developed for common bean (Phaseolus vulgaris L.) and related Phaseolus species. The applied strategy presents new and interesting aspects, such as the choice of accessions used, which was aimed at capturing a large portion of the genetic diversity present in the common bean, with particular focus on wild and domesticated materials from Mesoamerica and the identification of loci for sequencing. Indeed, the primer pairs for 34 loci were designed with the main strategy being to search for single-copy orthologous genes among the legumes (for use in other legume species and comparative analyses). The 10 remaining loci were selected as being near to domestication quantitative trait loci or detected as putatively under selection during domestication in previous studies. To provide an efficient and inexpensive genotyping platform for geneticists and breeders, we used sequence data to develop 60 new SNP markers for KASPar assay genotyping. The same sample was also genotyped with SNP markers developed for common bean in other studies for the same assay. This allowed testing for systematic bias according to the criteria chosen to select the genotypes in which the genetic diversity is surveyed during SNP discovery. Finally, we show that most of the SNP markers worked well in a set of accessions of other species belonging to the Phaseolus genus. The genetic resources developed will be very useful not only for breeding, but also for biodiversity conservation management and evolutionary studies on legumes.

Génétique des schizophrénies : mise en perspective des schizophrénies à début précoce et autres pathologies du développement

Schizophrenia (SCZ) is a severe brain disorder characterized by hallucinations, delusions, flat and/or inappropriate affect and cognitive impairment. The lifetime risk is about 0.5% with heritability of 65–85%. The prevalence of early-onset schizophrenia (defined here as before 15 years of age) has not been well studied, but is likely to be 5–10% of all cases. The rarity of early-onset SCZ has made it difficult to study. We focus on genetic studies of adults with schizophrenia, highlighting results for early-onset schizophrenia where available. Prior to the past 5years, studies failed to find replicable association or linkage between SCZ and specific genes when appropriate statistical corrections for multiple testing were used. Many false positive results were probably reported using the candidate gene approach. Recently, the development of single nucleotide polymorphism (SNP) “chips” has permitted large genome-wide association study (GWAS) analyses that suggest that across all age groups, a proportion of genetic risk can be attributed to a large number of common SNP, each with a very small effect on risk (odds ratios of 1.1 or less). The greatest known genetic effect is conferred by the 1.5–3 Mb 22q.11.2 deletions, which occurs in ∼ 1/4000–1/6000 births with SCZ developing in 20–30% of carriers. Large SNP and aCGH microarray studies have now identified associations between SCZ and other rare, large copy number variations (CNV, insertions and deletions) with high odds ratios (5–10), including deletions of 1q21, 2p16.3 (neurexin-1 gene), 3q29 and 15q13.3, and duplications of 16p11.2. Some of these CNV are also associated with autism or other developmental disorders as well as epilepsy or intellectual deficiency, suggesting some overlap in the mechanisms that contribute to risks of these disorders. Based on preliminary data from larger-scale analyses in progress, approximately 1–2% of cases carry a CNV that has been clearly associated with SCZ (ORs 4–12). Whole exome and genome sequencing studies of large adult samples will be the next steps to identify rarer SCZ-associated mutations, including point mutations and smaller as well as rarer CNV. Genetic findings are beginning to contribute to an understanding of biological mechanisms of SCZ risk and may lead to new approaches to treatment.

Single nucleotide polymorphism discovery in albacore and Atlantic bluefin tuna provides insights into worldwide population structure

The optimal management of the commercially important, but mostly over-exploited, pelagic tunas, albacore (Thunnus alalunga Bonn., 1788) and Atlantic bluefin tuna (BFT; Thunnus thynnus L., 1758), requires a better understanding of population structure than has been provided by previous molecular methods. Despite numerous studies of both species, their population structures remain controversial. This study reports the development of single nucleotide polymorphisms (SNPs) in albacore and BFT and the application of these SNPs to survey genetic variability across the geographic ranges of these tunas. A total of 616 SNPs were discovered in 35 albacore tuna by comparing sequences of 54 nuclear DNA fragments. A panel of 53 SNPs yielded FST values ranging from 0.0 to 0.050 between samples after genotyping 460 albacore collected throughout the distribution of this species. No significant heterogeneity was detected within oceans, but between-ocean comparisons (Atlantic, Pacific and Indian oceans along with Mediterranean Sea) were significant. Additionally, a 17-SNP panel was developed in Atlantic BFT by cross-species amplification in 107 fish. This limited number of SNPs discriminated between samples from the two major spawning areas of Atlantic BFT (FST =0.116). The SNP markers developed in this study can be used to genotype large numbers of fish without the need for standardizing alleles among laboratories.

Regional scale distribution of imidazolinone herbicide-resistant alleles in red rice ( Oryza sativa L.) determined through SNP markers

Red rice is the main weed in rice paddy fields. Imidazolinone herbicides in resistant rice cultivars currently provide a unique opportunity to control red rice in large-scale rice fields. However, the continuous use of this technology has resulted in imidazolinone-resistant red rice biotypes. This study aimed to identify the mechanism of herbicide resistance and the frequency and spatial distribution of the known imidazolinone herbicide-resistant alleles in red rice. The nucleotide sequence of the ALS gene indicated that the G 654E, S 653D and A 122T mutations are present in the imidazolinone herbicide-resistant rice cultivars IRGA 422 CL, SATOR CL and PUITÁ INTA CL, respectively. This information and the nucleotide sequence surrounding these mutations were used for the development of single nucleotide polymorphism (SNP) molecular markers to identify the possible mutations that confer herbicide resistance in red rice. This analysis was carried out in a total of 481 plants from 38 populations collected as individuals that escaped control with the herbicides imazethapyr and imazapic in rice paddy fields in Southern Brazil. The G 654E mutation was the most frequent, being found in 100% and 90.9% of the populations in the 2006/2007 and 20007/2008 seasons, respectively. In addition, the S 653D and A 122T mutations were also present either alone or as double or triple mutations in some plants. Target site insensitivity is the predominant mechanism of resistance in red rice resistant to imidazolinone herbicides in Southern Brazil. The high frequency of the S 653D mutation, the same mutation responsible for the resistance in the rice cultivar largely used in Southern Brazil, indicates that gene flow is occurring from the rice cultivar to red rice. Management practices related to increasing crop sanitation and decreasing of herbicide selection pressure through crop rotation should be enforced to prevent the evolution of herbicide resistance in red rice.