SNP discovery and development of a high-density genotyping array for sunflower

Development and characterization of a new sunflower source of resistance to race G of Orobanche cumana Wallr. derived from Helianthus anomalus

KEY MESSAGE

A new OrAnom1 gene introgressed in cultivated sunflower from wild Helianthus anomalus confers late post-attachment resistance to Orobanche cumana race G and maps to a target interval in Chromosome 4 where two receptor-like kinases (RLKs) have been identified in the H. anomalus genome as putative candidates. Sunflower broomrape is a parasitic weed that infects sunflower (Helianthus annuus L.) roots causing severe yield losses. Breeding for resistance is the most effective and sustainable control method. In this study, we report the identification, introgression, and genetic and physiological characterization of a new sunflower source of resistance to race G of broomrape developed from the wild annual sunflower H. anomalus (accession PI 468642). Crosses between PI 468642 and the susceptible line P21 were carried out, and the genetic study was conducted in BC1F1, BC1F2, and its derived BC1F3 populations. A BC1F5 germplasm named ANOM1 was developed through selection for race G resistance and resemblance to cultivated sunflower. The resistant trait showed monogenic and dominant inheritance. The gene, named OrAnom1, was mapped to Chromosome 4 within a 1.2 cM interval and co-segregated with 7 SNP markers. This interval corresponds to a 1.32 Mb region in the sunflower reference genome, housing a cluster of receptor-like kinase and receptor-like protein (RLK-RLP) genes. Notably, the analysis of the H. anomalus genome revealed the absence of RLPs in the OrAnom1 target region but featured two RLKs as possible OrAnom1 candidates. Rhizotron and histological studies showed that OrAnom1 determines a late post-attachment resistance mechanism. Broomrape can establish a vascular connection with the host, but parasite growth is stopped before tubercle development, showing phenolic compounds accumulation and tubercle necrosis. ANOM1 will contribute to broadening the genetic basis of broomrape resistance in the cultivated sunflower pool and to a better understanding of the molecular basis of the sunflower-broomrape interaction.

Event-level risk for negative alcohol consequences in emerging adults: The role of affect, motivation, and context

OBJECTIVE

Decades of research has found support for the motivational model of alcohol use, such that positive/negative affect are indirectly associated with drinking behavior through drinking motives. However, research on event-level drinking motives is in its nascent stage, and studies have yet to consider how drinking context plays a role in the motivational pathway to both event- and person-level drinking behavior. Therefore, the present study seeks to test whether drinking context mediates the effect of affect and motivation on drinking outcomes at both the event- and person-level.

METHOD

Data for this Stage 1 Registered Report will come from a recently completed ecological momentary assessment (EMA) study in emerging adults. The study collected data on 131 emerging adults, of whom 107 reported event-level social and solitary drinking during the EMA period. Multilevel structural equation modeling will be used to test whether predrinking affect is associated with predrinking motives, and whether drinking context (social vs. solitary drinking) mediates the effect of drinking motives on drinking outcomes. Models will parse within-/between-person variance, allowing the present study to test whether drinking context serves as a mechanism of risk in the motivational model at the event-level, or solely at the between-person level. Findings will inform personalized interventions and motivational models of drinking behavior. (PsycInfo Database Record (c) 2024 APA, all rights reserved).

On the usefulness of mock genomes to define heterotic pools, testers, and hybrid predictions in orphan crops

The advances in genomics in recent years have increased the accuracy and efficiency of breeding programs for many crops. Nevertheless, the adoption of genomic enhancement for several other crops essential in developing countries is still limited, especially for those that do not have a reference genome. These crops are more often called orphans. This is the first report to show how the results provided by different platforms, including the use of a simulated genome, called the mock genome, can generate in population structure and genetic diversity studies, especially when the intention is to use this information to support the formation of heterotic groups, choice of testers, and genomic prediction of single crosses. For that, we used a method to assemble a reference genome to perform the single-nucleotide polymorphism (SNP) calling without needing an external genome. Thus, we compared the analysis results using the mock genome with the standard approaches (array and genotyping-by-sequencing (GBS)). The results showed that the GBS-Mock presented similar results to the standard methods of genetic diversity studies, division of heterotic groups, the definition of testers, and genomic prediction. These results showed that a mock genome constructed from the population's intrinsic polymorphisms to perform the SNP calling is an effective alternative for conducting genomic studies of this nature in orphan crops, especially those that do not have a reference genome.

Transcriptomics of developing wild sunflower seeds from the extreme ends of a latitudinal gradient differing in seed oil composition

Seed oil composition, an important agronomic trait in cultivated sunflower, varies latitudinally across the native range of its wild progenitor. This pattern is thought to be driven by selection for a higher proportion of saturated fatty acids in southern populations compared with northern populations, likely due to the different temperatures experienced during seed germination. To investigate whether these differences in fatty acid composition between northern and southern populations correspond to transcriptional variation in the expression of genes involved in fatty acid metabolism, we sequenced RNA from developing seeds of sunflowers from Texas, USA, and Saskatchewan, Canada (the extreme ends of sunflower's latitudinal range) grown in a common garden. We found 4,741 genes to be differentially expressed between Texas and Canada, including several genes involved in lipid metabolism. Several differentially expressed lipid metabolism genes also colocalized with known oil quantitative trait loci (QTL). The genes producing stearoyl-ACP-desaturases (SAD) were of particular interest because of their known role in the conversion of fully saturated into unsaturated fatty acids. Two SAD genes were more highly expressed in seeds from Canadian populations, consistent with the observation of increased levels of unsaturated fatty acids in seeds from that region. We also constructed a gene co-expression network to investigate regional variation in network modules. The results of this analysis revealed regional differentiation for eight of 12 modules but no clear relationship with oil biosynthesis. Overall, the differential expression of SAD genes offers a partial explanation for the observed differences in seed oil composition between Texas and Canada, while the expression patterns of other metabolic genes suggest complex regulation of fatty acid production and usage across latitudes.

Applications of Genomic Tools in Plant Breeding: Crop Biofortification

Crop breeding has mainly been focused on increasing productivity, either directly or by decreasing the losses caused by biotic and abiotic stresses (that is, incorporating resistance to diseases and enhancing tolerance to adverse conditions, respectively). Quite the opposite, little attention has been paid to improve the nutritional value of crops. It has not been until recently that crop biofortification has become an objective within breeding programs, through either conventional methods or genetic engineering. There are many steps along this long path, from the initial evaluation of germplasm for the content of nutrients and health-promoting compounds to the development of biofortified varieties, with the available and future genomic tools assisting scientists and breeders in reaching their objectives as well as speeding up the process. This review offers a compendium of the genomic technologies used to explore and create biodiversity, to associate the traits of interest to the genome, and to transfer the genomic regions responsible for the desirable characteristics into potential new varieties. Finally, a glimpse of future perspectives and challenges in this emerging area is offered by taking the present scenario and the slow progress of the regulatory framework as the starting point.

Genetic and physiological characterization of sunflower resistance provided by the wild-derived OrDeb2 gene against highly virulent races of Orobanche cumana Wallr

Or_Deb2 confers post-attachment resistance to Orobanche cumana and is located in a 1.38 Mbp genomic interval containing a cluster of receptor-like kinase and receptor-like protein genes with nine high-confidence candidates. Sunflower broomrape is a holoparasitic angiosperm that parasitizes on sunflower roots, severely constraining crop yield. Breeding for resistance is the most effective method of control. Or_Deb2 is a dominant resistance gene introgressed into cultivated sunflower from a wild-related species that confers resistance to highly virulent broomrape races. The objectives of this study were as follows: (i) locate Or_Deb2 into the sunflower genome and determine putative candidate genes and (ii) characterize its underlying resistance mechanism. A segregating population from a cross between the sunflower resistant line DEB2, carrying Or_Deb2, and a susceptible line was phenotyped for broomrape resistance in four experiments, including different environments and two broomrape races (F_GV and G_TK). This population was also densely genotyped with microsatellite and SNP markers, which allowed locating Or_Deb2 within a 0.9 cM interval in the upper half of Chromosome 4. This interval corresponded to a 1.38 Mbp genomic region of the sunflower reference genome that contained a cluster of genes encoding LRR (leucine-rich repeat) receptor-like proteins lacking a cytoplasmic kinase domain and receptor-like kinases with one or two kinase domains and lacking an extracellular LRR region, which were valuable candidates for Or_Deb2. Rhizotron and histological studies showed that Or_Deb2 determines a post-attachment resistance response that blocks O. cumana development mainly at the cortex before the establishment of host-parasite vascular connections. This study will contribute to understand the interaction between crops and parasitic weeds, to establish durable breeding strategies based on genetic resistance and provide useful tools for marker-assisted selection and Or_Deb2 map-based cloning.

Construction of transgenic detection system of Brassica napus L. based on single nucleotide polymorphism chip

Brassica napus L. is a vital oil crop in China. As auxiliary tools for rapeseed breeding, transgenic technologies play a considerable role in heterosis, variety improvement, and pest resistance. Research on transgenic detection technologies is of great significance for the introduction, supervision, and development of transgenic rapeseed in China. However, the transgenic detection methods currently in use are complex and time-consuming, with low output. A single nucleotide polymorphism (SNP) chip can effectively overcome such limitations. In the present study, we collected 40 transgenic elements and designed 291 probes. The probe sequences were submitted to Illumina Company, and the Infinium chip technology was used to prepare SNP chips. In the present Brassica napus transgenic detection experiment, 84 high-quality probes of 17 transgenic elements were preliminarily screened, and genotyping effect was optimised for the probe signal value. Ultimately, a transgenic detection system for B. napus was developed. The developed system has the advantages of simple operation, minimal technical errors, and stable detection outcomes. A transgenic detection sensitivity test revealed that the probe designed could accurately detect 1% of transgenic samples and had high detection sensitivity. In addition, in repeatability tests, the CaMV35S promoter coefficient of variation was approximately 3.58%. Therefore, the SNP chip had suitable repeatability in transgene detection. The SNP chip developed could be used to construct transgenic detection systems for B. napus.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-03062-6.

A TILLING by sequencing approach to identify induced mutations in sunflower genes

The Targeting Induced Local Lesions in Genomes (TILLING) technology is a reverse genetic strategy broadly applicable to every kind of genome and represents an attractive tool for functional genomic and agronomic applications. It consists of chemical random mutagenesis followed by high-throughput screening of point mutations in targeted genomic regions. Although multiple methods for mutation discovery in amplicons have been described, next-generation sequencing (NGS) is the tool of choice for mutation detection because it quickly allows for the analysis of a large number of amplicons. The aim of the present work was to screen a previously generated sunflower TILLING population and identify alterations in genes involved in several important and complex physiological processes. Twenty-one candidate sunflower genes were chosen as targets for the screening. The TILLING by sequencing strategy allowed us to identify multiple mutations in selected genes and we subsequently validated 16 mutations in 11 different genes through Sanger sequencing. In addition to addressing challenges posed by outcrossing, our detection and validation of mutations in multiple regulatory loci highlights the importance of this sunflower population as a genetic resource.

Role of NGS and SNP genotyping methods in sugarcane improvement programs

Sugarcane (Saccharum spp.) is one of the most economically significant crops because of its high sucrose content and it is a promising biomass feedstock for biofuel production. Sugarcane genome sequencing and analysis is a difficult task due to its heterozygosity and polyploidy. Long sequence read technologies, PacBio Single-Molecule Real-Time (SMRT) sequencing, the Illumina TruSeq, and the Oxford Nanopore sequencing could solve the problem of genome assembly. On the applications side, next generation sequencing (NGS) technologies played a major role in the discovery of single nucleotide polymorphism (SNP) and the development of low to high throughput genotyping platforms. The two mainstream high throughput genotyping platforms are the SNP microarray and genotyping by sequencing (GBS). This paper reviews the NGS in sugarcane genomics, genotyping methodologies, and the choice of these methods. Array-based SNP genotyping is robust, provides consistent SNPs, and relatively easier downstream data analysis. The GBS method identifies large scale SNPs across the germplasm. A combination of targeted GBS and array-based genotyping methods should be used to increase the accuracy of genomic selection and marker-assisted breeding.

A 62K genic-SNP chip array for genetic studies and breeding applications in pigeonpea (Cajanus cajan L. Millsp.)

Pigeonpea is the second most important pulse legume crop for food and nutritional security of South Asia that requires accelerated breeding using high throughput genomic tools. Single nucleotide polymorphisms (SNPs) are highly suitable markers for this purpose because of their bi-allelic nature, reproducibility and high abundance in the genome. Here we report on development and use of a pigeonpea 62 K SNP chip array ‘CcSNPnks’ for Affymetrix GeneTitan^® platform. The array was designed after filtering 645,662 genic-SNPs identified by re-sequencing of 45 diverse genotypes and has 62,053 SNPs from 9629 genes belonging to five different categories, including 4314 single-copy genes unique to pigeonpea, 4328 single-copy genes conserved between soybean and pigeonpea, 156 homologs of agronomically important cloned genes, 746 disease resistance and defense response genes and 85 multi-copy genes of pigeonpea. This fully genic chip has 28.94% exonic, 33.04% intronic, 27.56% 5′UTR and 10.46% 3′UTR SNPs and incorporates multiple SNPs per gene allowing gene haplotype network analysis. It was used successfully for the analysis of genetic diversity and population structure of 95 pigeonpea varieties and high resolution mapping of 11 yield related QTLs for number of branches, pod bearing length and number of seeds per pod in a biparental RIL population. As an accurate high-density genotyping tool, ‘CcSNPnks’ chip array will be useful for high resolution fingerprinting, QTL mapping and genome wide as well as gene-based association studies in pigeonpea.

Korean soybean core collection: Genotypic and phenotypic diversity population structure and genome-wide association study

A core collection is a subset that represents genetic diversity of the total collection. Soybean (Glycine max (L.) Merr.) is one of major food and feed crops. It is the world's most cultivated annual herbaceous legume. Constructing a core collection for soybean could play a pivotal role in conserving and utilizing its genetic variability for research and breeding programs. To construct and evaluate a Korean soybean core collection, genotypic and phenotypic data as well as population structure, were analyzed. The Korean soybean core collection consisted of 430 accessions selected from 2,872 collections based on Affymetrix Axiom® 180k SoyaSNP array data. The core collection represented 99% of genotypic diversity of the total collection. Analysis of population structure clustered the core collection into five subpopulations. Accessions from South Korea and North Korea were distributed across five subpopulations. Analysis of molecular variance indicated that only 2.01% of genetic variation could be explained by geographic origins while 16.18% of genetic variation was accounted for by subpopulations. Genome-wide association study (GWAS) for days to flowering, flower color, pubescent color, and growth habit confirmed that the core collection had the same genetic diversity for tested traits as the total collection. The Korean soybean core collection was constructed based on genotypic information of the 180k SNP data. Size and phenotypic diversity of the core collection accounted for approximately 14.9% and 18.1% of the total collection, respectively. GWAS of core and total collections successfully confirmed loci associated with tested traits. Consequently, the present study showed that the Korean soybean core collection could provide fundamental and practical material and information for both soybean genetic research and breeding programs.

Development and validation of an SNP genotyping array and construction of a high-density linkage map in castor

Castor is a commercially important oilseed crop that provides raw materials for several industries. Currently, the availability of genomic resources for castor is very limited. In this study, genome-wide SNPs were discovered in castor via whole-genome sequencing of 14 diverse lines to an average of 34X coverage. A total of 2,179,759 putative SNPs were detected, and a genotyping array was designed with 6,000 high-quality SNPs representing 2,492 scaffolds of the draft castor genome (87.5% genome coverage). The array was validated by genotyping a panel of 314 inbred castor lines, which resulted in 5,025 scorable SNPs with a high call rate (98%) and reproducibility (100%). Using this array, a consensus linkage map consisting of 1,978 SNP loci was constructed with an average inter-marker distance of 0.55 cM. The genome-wide SNP data, the genotyping array and the dense linkage map are valuable genomic tools for promoting high-throughput genomic research and molecular breeding in castor.

Development and Application of High-Density Axiom Cajanus SNP Array with 56K SNPs to Understand the Genome Architecture of Released Cultivars and Founder Genotypes

As one of the major outputs of next-generation sequencing (NGS), a large number of genome-wide single-nucleotide polymorphisms (SNPs) have been developed in pigeonpea [ (L.) Huth.]. However, SNPs require a genotyping platform or assay to be used in different evolutionary studies or in crop improvement programs. Therefore, we developed an Axiom SNP array with 56K SNPs uniformly distributed across the genome and assessed its utility in a genetic diversity study. From the whole-genome resequencing (WGRS) data on 104 pigeonpea lines, ∼2 million sequence variations (SNPs and insertion-deletions [InDels]) were identified, from which a subset of 56,512 unique and informative sequence variations were selected to develop the array. The Axiom SNP array developed was used for genotyping 103 pigeonpea lines encompassing 63 cultivars released between 1960 and 2014 and 40 breeding, germplasm, and founder lines. Genotyping data thus generated on 103 pigeonpea lines provided 51,201 polymorphic SNPs and InDels. Genetic diversity analysis provided in-depth insights into the genetic architecture and trends in temporal diversity in pigeonpea cultivars. Therefore, the continuous use of the high-density Axiom SNP array developed will accelerate high-resolution trait mapping, marker-assisted breeding, and genomic selection efforts in pigeonpea.

Sunflower Hybrid Breeding: From Markers to Genomic Selection

In sunflower, molecular markers for simple traits as, e.g., fertility restoration, high oleic acid content, herbicide tolerance or resistances to Plasmopara halstedii, Puccinia helianthi, or Orobanche cumana have been successfully used in marker-assisted breeding programs for years. However, agronomically important complex quantitative traits like yield, heterosis, drought tolerance, oil content or selection for disease resistance, e.g., against Sclerotinia sclerotiorum have been challenging and will require genome-wide approaches. Plant genetic resources for sunflower are being collected and conserved worldwide that represent valuable resources to study complex traits. Sunflower association panels provide the basis for genome-wide association studies, overcoming disadvantages of biparental populations. Advances in technologies and the availability of the sunflower genome sequence made novel approaches on the whole genome level possible. Genotype-by-sequencing, and whole genome sequencing based on next generation sequencing technologies facilitated the production of large amounts of SNP markers for high density maps as well as SNP arrays and allowed genome-wide association studies and genomic selection in sunflower. Genome wide or candidate gene based association studies have been performed for traits like branching, flowering time, resistance to Sclerotinia head and stalk rot. First steps in genomic selection with regard to hybrid performance and hybrid oil content have shown that genomic selection can successfully address complex quantitative traits in sunflower and will help to speed up sunflower breeding programs in the future. To make sunflower more competitive toward other oil crops higher levels of resistance against pathogens and better yield performance are required. In addition, optimizing plant architecture toward a more complex growth type for higher plant densities has the potential to considerably increase yields per hectare. Integrative approaches combining omic technologies (genomics, transcriptomics, proteomics, metabolomics and phenomics) using bioinformatic tools will facilitate the identification of target genes and markers for complex traits and will give a better insight into the mechanisms behind the traits.

Revolution in Genotyping Platforms for Crop Improvement

In the past decade, the application of high-throughput sequencing to crop genotyping has given rise to novel platforms capable of genotyping tens of thousands of genome-wide DNA markers. Coupled with the decreasing costs of sequencing, this rapid increase in markers allows accelerated and highly accurate genotyping of entire crop populations and diversity sets using single nucleotide polymorphisms (SNPs). These revolutionary advances accelerate crop improvement by facilitating a more precise connection of phenotype to genotype through association studies, linkage mapping and diversity analysis. The platforms driving the advances in genotyping are array technologies and genotyping by sequencing (GBS) methods, which include both low-coverage whole genome resequencing (skim sequencing) and reduced representation sequencing (RRS) approaches. Here, we outline and compare these genotyping platforms and provide a perspective on the promising future of crop genotyping. While SNP arrays provide high quality, simple handling, and unchallenging analysis, the lower cost of RRS and the greater data volume produced by skim sequencing suggest that use of GBS will become more prevalent in crop genomics as sequencing costs decrease and data analysis becomes more streamlined. Graphical Abstract.

Main and epistatic QTL analyses for Sclerotinia Head Rot resistance in sunflower

Sclerotinia Head Rot (SHR), a disease caused by Sclerotinia sclerotiorum, is one of the most limiting factors in sunflower production. In this study, we identified genomic loci associated with resistance to SHR to support the development of assisted breeding strategies. We genotyped 114 Recombinant Inbred Lines (RILs) along with their parental lines (PAC2 -partially resistant-and RHA266 -susceptible-) by using a 384 single nucleotide polymorphism (SNP) Illumina Oligo Pool Assay to saturate a sunflower genetic map. Subsequently, we tested these lines for SHR resistance using assisted inoculations with S. sclerotiorum ascospores. We also conducted a randomized complete-block assays with three replicates to visually score disease incidence (DI), disease severity (DS), disease intensity (DInt) and incubation period (IP) through four field trials (2010-2014). We finally assessed main effect quantitative trait loci (M-QTLs) and epistatic QTLs (E-QTLs) by composite interval mapping (CIM) and mixed-model-based composite interval mapping (MCIM), respectively. As a result of this study, the improved map incorporates 61 new SNPs over candidate genes. We detected a broad range of narrow sense heritability (h2) values (1.86-59.9%) as well as 36 M-QTLs and 13 E-QTLs along 14 linkage groups (LGs). On LG1, LG10, and LG15, we repeatedly detected QTLs across field trials; which emphasizes their putative effectiveness against SHR. In all selected variables, most of the identified QTLs showed high determination coefficients, associated with moderate to high heritability values. Using markers shared with previous Sclerotinia resistance studies, we compared the QTL locations in LG1, LG2, LG8, LG10, LG11, LG15 and LG16. This study constitutes the largest report of QTLs for SHR resistance in sunflower. Further studies focusing on the regions in LG1, LG10, and LG15 harboring the detected QTLs are necessary to identify causal alleles and contribute to unraveling the complex genetic basis governing the resistance.

Using semantics for representing experimental protocols

Background

An experimental protocol is a sequence of tasks and operations executed to perform experimental research in biological and biomedical areas, e.g. biology, genetics, immunology, neurosciences, virology. Protocols often include references to equipment, reagents, descriptions of critical steps, troubleshooting and tips, as well as any other information that researchers deem important for facilitating the reusability of the protocol. Although experimental protocols are central to reproducibility, the descriptions are often cursory. There is the need for a unified framework with respect to the syntactic structure and the semantics for representing experimental protocols.

Results

In this paper we present “SMART Protocols ontology”, an ontology for representing experimental protocols. Our ontology represents the protocol as a workflow with domain specific knowledge embedded within a document. We also present the Sample Instrument Reagent Objective (SIRO) model, which represents the minimal common information shared across experimental protocols. SIRO was conceived in the same realm as the Patient Intervention Comparison Outcome (PICO) model that supports search, retrieval and classification purposes in evidence based medicine. We evaluate our approach against a set of competency questions modeled as SPARQL queries and processed against a set of published and unpublished protocols modeled with the SP Ontology and the SIRO model. Our approach makes it possible to answer queries such as Which protocols use tumor tissue as a sample.

Conclusion

Improving reporting structures for experimental protocols requires collective efforts from authors, peer reviewers, editors and funding bodies. The SP Ontology is a contribution towards this goal. We build upon previous experiences and bringing together the view of researchers managing protocols in their laboratory work. Website: https://smartprotocols.github.io/.

Crop Breeding Chips and Genotyping Platforms: Progress, Challenges, and Perspectives

There is a rapidly rising trend in the development and application of molecular marker assays for gene mapping and discovery in field crops and trees. Thus far, more than 50 SNP arrays and 15 different types of genotyping-by-sequencing (GBS) platforms have been developed in over 25 crop species and perennial trees. However, much less effort has been made on developing ultra-high-throughput and cost-effective genotyping platforms for applied breeding programs. In this review, we discuss the scientific bottlenecks in existing SNP arrays and GBS technologies and the strategies to develop targeted platforms for crop molecular breeding. We propose that future practical breeding platforms should adopt automated genotyping technologies, either array or sequencing based, target functional polymorphisms underpinning economic traits, and provide desirable prediction accuracy for quantitative traits, with universal applications under wide genetic backgrounds in crops. The development of such platforms faces serious challenges at both the technological level due to cost ineffectiveness, and the knowledge level due to large genotype-phenotype gaps in crop plants. It is expected that such genotyping platforms will be achieved in the next ten years in major crops in consideration of (a) rapid development in gene discovery of important traits, (b) deepened understanding of quantitative traits through new analytical models and population designs, (c) integration of multi-layer -omics data leading to identification of genes and pathways responsible for important breeding traits, and (d) improvement in cost effectiveness of large-scale genotyping. Crop breeding chips and genotyping platforms will provide unprecedented opportunities to accelerate the development of cultivars with desired yield potential, quality, and enhanced adaptation to mitigate the effects of climate change.

Genome-wide analysis of allele frequency change in sunflower crop-wild hybrid populations evolving under natural conditions

Crop-wild hybridization occurs in numerous plant species and could alter the genetic structure and evolutionary dynamics of wild populations. Studying crop-derived alleles in wild populations is also relevant to assessing/mitigating the risks associated with transgene escape. To date, crop-wild hybridization has generally been examined via short-term studies, typically within a single generation, focusing on few traits or genetic markers. Little is known about patterns of selection on crop-derived alleles over multiple generations, particularly at a genome-wide scale. Here, we documented patterns of natural selection in an experimental crop × wild sunflower population that was allowed to evolve under natural conditions for two generations at two locations. Allele frequencies at a genome-wide collection of SNPs were tracked across generations, and a common garden experiment was conducted to compare trait means between generations. These data allowed us to identify instances of selection on crop-derived alleles/traits and, in concert with QTL mapping results, test for congruence between our genotypic and phenotypic results. We found that natural selection overwhelmingly favours wild alleles and phenotypes. However, crop alleles in certain genomic regions can be favoured, and these changes often occurred in parallel across locations. We did not, however, consistently observe close agreement between our genotypic and phenotypic results. For example, when a trait evolved towards the wild phenotype, wild QTL alleles associated with that trait did not consistently increase in frequency. We discuss these results in the context of crop allele introgression into wild populations and implications for the management of GM crops.

An innovative SNP genotyping method adapting to multiple platforms and throughputs

An innovative genotyping method designated as semi-thermal asymmetric reverse PCR (STARP) was developed for genotyping individual SNPs with improved accuracy, flexible throughputs, low operational costs, and high platform compatibility. Multiplex chip-based technology for genome-scale genotyping of single nucleotide polymorphisms (SNPs) has made great progress in the past two decades. However, PCR-based genotyping of individual SNPs still remains problematic in accuracy, throughput, simplicity, and/or operational costs as well as the compatibility with multiple platforms. Here, we report a novel SNP genotyping method designated semi-thermal asymmetric reverse PCR (STARP). In this method, genotyping assay was performed under unique PCR conditions using two universal priming element-adjustable primers (PEA-primers) and one group of three locus-specific primers: two asymmetrically modified allele-specific primers (AMAS-primers) and their common reverse primer. The two AMAS-primers each were substituted one base in different positions at their 3' regions to significantly increase the amplification specificity of the two alleles and tailed at 5' ends to provide priming sites for PEA-primers. The two PEA-primers were developed for common use in all genotyping assays to stringently target the PCR fragments generated by the two AMAS-primers with similar PCR efficiencies and for flexible detection using either gel-free fluorescence signals or gel-based size separation. The state-of-the-art primer design and unique PCR conditions endowed STARP with all the major advantages of high accuracy, flexible throughputs, simple assay design, low operational costs, and platform compatibility. In addition to SNPs, STARP can also be employed in genotyping of indels (insertion-deletion polymorphisms). As vast variations in DNA sequences are being unearthed by many genome sequencing projects and genotyping by sequencing, STARP will have wide applications across all biological organisms in agriculture, medicine, and forensics.

Development and Evaluation of a High Density Genotyping 'Axiom_Arachis' Array with 58 K SNPs for Accelerating Genetics and Breeding in Groundnut

Single nucleotide polymorphisms (SNPs) are the most abundant DNA sequence variation in the genomes which can be used to associate genotypic variation to the phenotype. Therefore, availability of a high-density SNP array with uniform genome coverage can advance genetic studies and breeding applications. Here we report the development of a high-density SNP array 'Axiom_Arachis' with 58 K SNPs and its utility in groundnut genetic diversity study. In this context, from a total of 163,782 SNPs derived from DNA resequencing and RNA-sequencing of 41 groundnut accessions and wild diploid ancestors, a total of 58,233 unique and informative SNPs were selected for developing the array. In addition to cultivated groundnuts (Arachis hypogaea), fair representation was kept for other diploids (A. duranensis, A. stenosperma, A. cardenasii, A. magna and A. batizocoi). Genotyping of the groundnut 'Reference Set' containing 300 genotypes identified 44,424 polymorphic SNPs and genetic diversity analysis provided in-depth insights into the genetic architecture of this material. The availability of the high-density SNP array 'Axiom_Arachis' with 58 K SNPs will accelerate the process of high resolution trait genetics and molecular breeding in cultivated groundnut.

Range-wide phenotypic and genetic differentiation in wild sunflower

BACKGROUND

Divergent phenotypes and genotypes are key signals for identifying the targets of natural selection in locally adapted populations. Here, we used a combination of common garden phenotyping for a variety of growth, plant architecture, and seed traits, along with single-nucleotide polymorphism (SNP) genotyping to characterize range-wide patterns of diversity in 15 populations of wild sunflower (Helianthus annuus L.) sampled along a latitudinal gradient in central North America. We analyzed geographic patterns of phenotypic diversity, quantified levels of within-population SNP diversity, and also determined the extent of population structure across the range of this species. We then used these data to identify significantly over-differentiated loci as indicators of genomic regions that likely contribute to local adaptation.

RESULTS

Traits including flowering time, plant height, and seed oil composition (i.e., percentage of saturated fatty acids) were significantly correlated with latitude, and thus differentiated northern vs. southern populations. Average pairwise FST was found to be 0.21, and a STRUCTURE analysis identified two significant clusters that largely separated northern and southern individuals. The significant FST outliers included a SNP in HaFT2, a flowering time gene that has been previously shown to co-localize with flowering time QTL, and which exhibits a known cline in gene expression.

CONCLUSIONS

Latitudinal differentiation in both phenotypic traits and SNP allele frequencies is observed across wild sunflower populations in central North America. Such differentiation may play an important adaptive role across the range of this species, and could facilitate adaptation to a changing climate.

Identification and characterization of contrasting sunflower genotypes to early leaf senescence process combining molecular and physiological studies (Helianthus annuus L.)

Leaf senescence is a complex mechanism ruled by multiple genetic and environmental variables that affect crop yields. It is the last stage in leaf development, is characterized by an active decline in photosynthetic rate, nutrients recycling and cell death. The aim of this work was to identify contrasting sunflower inbred lines differing in leaf senescence and to deepen the study of this process in sunflower. Ten sunflower genotypes, previously selected by physiological analysis from 150 inbred genotypes, were evaluated under field conditions through physiological, cytological and molecular analysis. The physiological measurement allowed the identification of two contrasting senescence inbred lines, R453 and B481-6, with an increase in yield in the senescence delayed genotype. These findings were confirmed by cytological and molecular analysis using TUNEL, genomic DNA gel electrophoresis, flow sorting and gene expression analysis by qPCR. These results allowed the selection of the two most promising contrasting genotypes, which enables future studies and the identification of new biomarkers associated to early senescence in sunflower. In addition, they allowed the tuning of cytological techniques for a non-model species and its integration with molecular variables.

New resources for genetic studies in Populus nigra: genome-wide SNP discovery and development of a 12k Infinium array

Whole genome resequencing of 51 Populus nigra (L.) individuals from across Western Europe was performed using Illumina platforms. A total number of 1 878 727 SNPs distributed along the P. nigra reference sequence were identified. The SNP calling accuracy was validated with Sanger sequencing. SNPs were selected within 14 previously identified QTL regions, 2916 expressional candidate genes related to rust resistance, wood properties, water-use efficiency and bud phenology and 1732 genes randomly spread across the genome. Over 10 000 SNPs were selected for the construction of a 12k Infinium Bead-Chip array dedicated to association mapping. The SNP genotyping assay was performed with 888 P. nigra individuals. The genotyping success rate was 91%. Our high success rate was due to the discovery panel design and the stringent parameters applied for SNP calling and selection. In the same set of P. nigra genotypes, linkage disequilibrium throughout the genome decayed on average within 5-7 kb to half of its maximum value. As an application test, ADMIXTURE analysis was performed with a selection of 600 SNPs spread throughout the genome and 706 individuals collected along 12 river basins. The admixture pattern was consistent with genetic diversity revealed by neutral markers and the geographical distribution of the populations. These newly developed SNP resources and genotyping array provide a valuable tool for population genetic studies and identification of QTLs through natural-population based genetic association studies in P. nigra.

Development of cleaved amplified polymorphic sequence markers and a CAPS-based genetic linkage map in watermelon (Citrullus lanatus [Thunb.] Matsum. and Nakai) constructed using whole-genome re-sequencing data

Cleaved amplified polymorphic sequence (CAPS) markers are useful tools for detecting single nucleotide polymorphisms (SNPs). This study detected and converted SNP sites into CAPS markers based on high-throughput re-sequencing data in watermelon, for linkage map construction and quantitative trait locus (QTL) analysis. Two inbred lines, Cream of Saskatchewan (COS) and LSW-177 had been re-sequenced and analyzed by Perl self-compiled script for CAPS marker development. 88.7% and 78.5% of the assembled sequences of the two parental materials could map to the reference watermelon genome, respectively. Comparative assembled genome data analysis provided 225,693 and 19,268 SNPs and indels between the two materials. 532 pairs of CAPS markers were designed with 16 restriction enzymes, among which 271 pairs of primers gave distinct bands of the expected length and polymorphic bands, via PCR and enzyme digestion, with a polymorphic rate of 50.94%. Using the new CAPS markers, an initial CAPS-based genetic linkage map was constructed with the F2 population, spanning 1836.51 cM with 11 linkage groups and 301 markers. 12 QTLs were detected related to fruit flesh color, length, width, shape index, and brix content. These newly CAPS markers will be a valuable resource for breeding programs and genetic studies of watermelon.

Diversity analysis and genomic prediction of Sclerotinia resistance in sunflower using a new 25 K SNP genotyping array

KEY MESSAGE

We have developed a SNP array for sunflower containing more than 25 K markers, representing single loci mostly in or near transcribed regions of the genome. The array was successfully applied to genotype a diversity panel of lines, hybrids, and mapping populations and represented well the genetic diversity of cultivated sunflower. Results of PCoA and population substructure analysis underlined the complexity of the genetic composition of current elite breeding material. The performance of this genotyping platform for genome-based prediction of phenotypes and detection of QTL with improved resolution could be demonstrated based on the re-evaluation of a population segregating for resistance to Sclerotinia midstalk rot. Given our results, the newly developed 25 K SNP array is expected to be of great utility for the most important applications in genome-based sunflower breeding and research.

ABSTRACT

Genotyping with a large number of molecular markers is a prerequisite to conduct genome-based genetic analyses with high precision. Here, we report the design and performance of a 25 K SNP genotyping array for sunflower (Helianthus annuus L.). SNPs were discovered based on variant calling in de novo assembled, UniGene-based contigs of sunflower derived from whole genome sequencing and amplicon sequences originating from four and 48 inbred lines, respectively. After inclusion of publically available transcriptome-derived SNPs, in silico design of the Illumina(®) Infinium iSelect HD BeadChip yielded successful assays for 22,299 predominantly haplotype-specific SNPs. The array was validated in a sunflower diversity panel including inbred lines, open-pollinated varieties, introgression lines, landraces, recombinant inbred lines, and F2 populations. Validation provided 20,502 high-quality bi-allelic SNPs with stable cluster performance whereby each SNP marker represents a single locus mostly in or near transcribed regions of the sunflower genome. Analyses of population structure and quantitative resistance to Sclerotinia midstalk rot demonstrate that this array represents a significant improvement over currently available genomic tools for genetic diversity analyses, genome-wide marker-trait association studies, and genetic mapping in sunflower.

Transcriptome Analysis of Gerbera hybrida Including in silico Confirmation of Defense Genes Found

For the ornamental crop Gerbera hybrida, breeding at the moment is done using conventional methods. As this has drawbacks in breeding speed and efficiency, especially for complex traits like disease resistance, we set out to develop genomic resources. The leaf and flower bud transcriptomes of four parents, used to generate two gerbera populations, were sequenced using Illumina paired-end sequencing. In total, 36,770 contigs with an average length of 1397 bp were generated and these have been the starting point for SNP identification and annotation. The consensus contig sequences were used to map reads of individual parents, to identify genotype specific SNPs, and to assess the presence of common SNPs between genotypes. Comparison with the non-redundant protein database (nr) showed that 29,146 contigs gave BLAST hits. Of sequences with blast results, 73.3% obtained a clear gene ontology (GO) annotation. EST contigs coding for enzymes were found in Kyoto Encyclopedia of Genes and Genomes maps (KEGG). Through, these annotated data and KEGG molecular interaction network, transcripts associated with the phenylpropanoid metabolism, other secondary metabolite biosynthesis pathways, phytohormone biosynthesis and signal transduction were analyzed in more detail. Identifying genes involved in these processes could provide genetic and genomic resources for studying the mechanism of disease resistance in gerbera.

Development of two major resources for pea genomics: the GenoPea 13.2K SNP Array and a high-density, high-resolution consensus genetic map

Single nucleotide polymorphism (SNP) arrays represent important genotyping tools for innovative strategies in both basic research and applied breeding. Pea is an important food, feed and sustainable crop with a large (about 4.45 Gbp) but not yet available genome sequence. In the present study, 12 pea recombinant inbred line populations were genotyped using the newly developed GenoPea 13.2K SNP Array. Individual and consensus genetic maps were built providing insights into the structure and organization of the pea genome. Largely collinear genetic maps of 3918-8503 SNPs were obtained from all mapping populations, and only two of these exhibited putative chromosomal rearrangement signatures. Similar distortion patterns in different populations were noted. A total of 12 802 transcript-derived SNP markers placed on a 15 079-marker high-density, high-resolution consensus map allowed the identification of ohnologue-rich regions within the pea genome and the localization of local duplicates. Dense syntenic networks with sequenced legume genomes were further established, paving the way for the identification of the molecular bases of important agronomic traits segregating in the mapping populations. The information gained on the structure and organization of the genome from this research will undoubtedly contribute to the understanding of the evolution of the pea genome and to its assembly. The GenoPea 13.2K SNP Array and individual and consensus genetic maps are valuable genomic tools for plant scientists to strengthen pea as a model for genetics and physiology and enhance breeding.

Single-copy gene based 50 K SNP chip for genetic studies and molecular breeding in rice

Single nucleotide polymorphism (SNP) is the most abundant DNA sequence variation present in plant genomes. Here, we report the design and validation of a unique genic-SNP genotyping chip for genetic and evolutionary studies as well as molecular breeding applications in rice. The chip incorporates 50,051 SNPs from 18,980 different genes spanning 12 rice chromosomes, including 3,710 single-copy (SC) genes conserved between wheat and rice, 14,959 SC genes unique to rice, 194 agronomically important cloned rice genes and 117 multi-copy rice genes. Assays with this chip showed high success rate and reproducibility because of the SC gene based array with no sequence redundancy and cross-hybridisation problems. The usefulness of the chip in genetic diversity and phylogenetic studies of cultivated and wild rice germplasm was demonstrated. Furthermore, its efficacy was validated for analysing background recovery in improved mega rice varieties with submergence tolerance developed through marker-assisted backcross breeding.

A bi-filtering method for processing single nucleotide polymorphism array data improves the quality of genetic map and accuracy of quantitative trait locus mapping in doubled haploid populations of polyploid Brassica napus

Background

Single nucleotide polymorphism (SNP) markers have a wide range of applications in crop genetics and genomics. Due to their polyploidy nature, many important crops, such as wheat, cotton and rapeseed contain a large amount of repeat and homoeologous sequences in their genomes, which imposes a huge challenge in high-throughput genotyping with sequencing and/or array technologies. Allotetraploid Brassica napus (AACC, 2n = 4x = 38) comprises of two highly homoeologous sub-genomes derived from its progenitor species B. rapa (AA, 2n = 2x = 20) and B. oleracea (CC, 2n = 2x = 18), and is an ideal species to exploit methods for reducing the interference of extensive inter-homoeologue polymorphisms (mHemi-SNPs and Pseudo-simple SNPs) between closely related sub-genomes.

Results

Based on a recent B. napus 6K SNP array, we developed a bi-filtering procedure to identify unauthentic lines in a DH population, and mHemi-SNPs and Pseudo-simple SNPs in an array data matrix. The procedure utilized both monomorphic and polymorphic SNPs in the DH population and could effectively distinguish the mHemi-SNPs and Pseudo-simple SNPs that resulted from superposition of the signals from multiple SNPs. Compared with conventional procedure for array data processing, the bi-filtering method could minimize the pseudo linkage relationship caused by the mHemi-SNPs and Pseudo-simple SNPs, thus improving the quality of SNP genetic map. Furthermore, the improved genetic map could increase the accuracies of mapping of QTLs as demonstrated by the ability to eliminate non-real QTLs in the mapping population.

Conclusions

The bi-filtering analysis of the SNP array data represents a novel approach to effectively assigning the multi-loci SNP genotypes in polyploid B. napus and may find wide applications to SNP analyses in polyploid crops.

Electronic supplementary material

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

RNA-Seq Uncovers SNPs and Alternative Splicing Events in Asian Lotus (Nelumbo nucifera)

RNA-Seq is an efficient way to comprehensively identify single nucleotide polymorphisms (SNPs) and alternative splicing (AS) events from the expressed genes. In this study, we conducted transcriptome sequencing of four Asian lotus (Nelumbo nucifera) cultivars using Illumina HiSeq2000 platform to identify SNPs and AS events in lotus. A total of 505 million pair-end RNA-Seq reads were generated from four cultivars, of which 86% were mapped to the lotus reference genome. Using the four sets of data together, a total of 357,689 putative SNPs were identified with an average density of one SNP per 2.2 kb. These SNPs were located in 1,253 scaffolds and 15,016 expressed genes. A/G and C/T were the two major types of SNPs in the Asian lotus transcriptome. In parallel, a total of 177,540 AS events were detected in the four cultivars and were distributed in 64% of the expressed genes of lotus. The predominant type of AS events was alternative 5’ first exon, which accounted for 41.2% of all the observed AS events, and exon skipping only accounted for 4.3% of all AS. Gene Ontology analysis was conducted to analyze the function of the genes containing SNPs and AS events. Validation of selected SNPs and AS events revealed that 74% of SNPs and 80% of AS events were reliable, which indicates that RNA-Seq is an efficient approach to uncover gene-associated SNPs and AS events. A large number of SNPs and AS events identified in our study will facilitate further genetic and functional genomics research in lotus.

Association mapping in sunflower (Helianthus annuus L.) reveals independent control of apical vs. basal branching

BACKGROUND

Shoot branching is an important determinant of plant architecture and influences various aspects of growth and development. Selection on branching has also played an important role in the domestication of crop plants, including sunflower (Helianthus annuus L.). Here, we describe an investigation of the genetic basis of variation in branching in sunflower via association mapping in a diverse collection of cultivated sunflower lines.

RESULTS

Detailed phenotypic analyses revealed extensive variation in the extent and type of branching within the focal population. After correcting for population structure and kinship, association analyses were performed using a genome-wide collection of SNPs to identify genomic regions that influence a variety of branching-related traits. This work resulted in the identification of multiple previously unidentified genomic regions that contribute to variation in branching. Genomic regions that were associated with apical and mid-apical branching were generally distinct from those associated with basal and mid-basal branching. Homologs of known branching genes from other study systems (i.e., Arabidopsis, rice, pea, and petunia) were also identified from the draft assembly of the sunflower genome and their map positions were compared to those of associations identified herein. Numerous candidate branching genes were found to map in close proximity to significant branching associations.

CONCLUSIONS

In sunflower, variation in branching is genetically complex and overall branching patterns (i.e., apical vs. basal) were found to be influenced by distinct genomic regions. Moreover, numerous candidate branching genes mapped in close proximity to significant branching associations. Although the sunflower genome exhibits localized islands of elevated linkage disequilibrium (LD), these non-random associations are known to decay rapidly elsewhere. The subset of candidate genes that co-localized with significant associations in regions of low LD represents the most promising target for future functional analyses.

Relocation of a rust resistance gene R 2 and its marker-assisted gene pyramiding in confection sunflower (Helianthus annuus L.)

The rust resistance gene R 2 was reassigned to linkage group 14 of the sunflower genome. DNA markers linked to R 2 were identified and used for marker-assisted gene pyramiding in a confection type genetic background. Due to the frequent evolution of new pathogen races, sunflower rust is a recurring threat to sunflower production worldwide. The inbred line Morden Cross 29 (MC29) carries the rust resistance gene, R 2 , conferring resistance to numerous races of rust fungus in the US, Canada, and Australia, and can be used as a broad-spectrum resistance resource. Based on phenotypic assessments and SSR marker analyses on the 117 F2 individuals derived from a cross of HA 89 with MC29 (USDA), R 2 was mapped to linkage group (LG) 14 of the sunflower, and not to the previously reported location on LG9. The closest SSR marker HT567 was located at 4.3 cM distal to R 2 . Furthermore, 36 selected SNP markers from LG14 were used to saturate the R 2 region. Two SNP markers, NSA_002316 and SFW01272, flanked R 2 at a genetic distance of 2.8 and 1.8 cM, respectively. Of the three closely linked markers, SFW00211 amplified an allele specific for the presence of R 2 in a marker validation set of 46 breeding lines, and SFW01272 was also shown to be diagnostic for R 2 . These newly developed markers, together with the previously identified markers linked to the gene R 13a , were used to screen 524 F2 individuals from a cross of a confection R 2 line and HA-R6 carrying R 13a . Eleven homozygous double-resistant F2 plants with the gene combination of R 2 and R 13a were obtained. This double-resistant line will be extremely useful in confection sunflower, where few rust R genes are available, risking evolution of new virulence phenotypes and further disease epidemics.

Population structure and genetic diversity characterization of a sunflower association mapping population using SSR and SNP markers

BACKGROUND

Argentina has a long tradition of sunflower breeding, and its germplasm is a valuable genetic resource worldwide. However, knowledge of the genetic constitution and variability levels of the Argentinean germplasm is still scarce, rendering the global map of cultivated sunflower diversity incomplete. In this study, 42 microsatellite loci and 384 single nucleotide polymorphisms (SNPs) were used to characterize the first association mapping population used for quantitative trait loci mapping in sunflower, along with a selection of allied open-pollinated and composite populations from the germplasm bank of the National Institute of Agricultural Technology of Argentina. The ability of different kinds of markers to assess genetic diversity and population structure was also evaluated.

RESULTS

The analysis of polymorphism in the set of sunflower accessions studied here showed that both the microsatellites and SNP markers were informative for germplasm characterization, although to different extents. In general, the estimates of genetic variability were moderate. The average genetic diversity, as quantified by the expected heterozygosity, was 0.52 for SSR loci and 0.29 for SNPs. Within SSR markers, those derived from non-coding regions were able to capture higher levels of diversity than EST-SSR. A significant correlation was found between SSR and SNP- based genetic distances among accessions. Bayesian and multivariate methods were used to infer population structure. Evidence for the existence of three different genetic groups was found consistently across data sets (i.e., SSR, SNP and SSR + SNP), with the maintainer/restorer status being the most prevalent characteristic associated with group delimitation.

CONCLUSION

The present study constitutes the first report comparing the performance of SSR and SNP markers for population genetics analysis in cultivated sunflower. We show that the SSR and SNP panels examined here, either used separately or in conjunction, allowed consistent estimations of genetic diversity and population structure in sunflower breeding materials. The generated knowledge about the levels of diversity and population structure of sunflower germplasm is an important contribution to this crop breeding and conservation.

Methods for the design, implementation, and analysis of illumina infinium™ SNP assays in plants

The advent of Next-Generation sequencing-by-synthesis technologies has fuelled SNP discovery, genotyping, and screening of populations in myriad ways for many species, including various plant species. One technique widely applied to screening a large number of SNP markers over a large number of samples is the Illumina Infinium™ assay.

De novo transcriptome sequence assembly and identification of AP2/ERF transcription factor related to abiotic stress in parsley (Petroselinum crispum)

Parsley is an important biennial Apiaceae species that is widely cultivated as herb, spice, and vegetable. Previous studies on parsley principally focused on its physiological and biochemical properties, including phenolic compound and volatile oil contents. However, little is known about the molecular and genetic properties of parsley. In this study, 23,686,707 high-quality reads were obtained and assembled into 81,852 transcripts and 50,161 unigenes for the first time. Functional annotation showed that 30,516 unigenes had sequence similarity to known genes. In addition, 3,244 putative simple sequence repeats were detected in curly parsley. Finally, 1,569 of the identified unigenes belonged to 58 transcription factor families. Various abiotic stresses have a strong detrimental effect on the yield and quality of parsley. AP2/ERF transcription factors have important functions in plant development, hormonal regulation, and abiotic response. A total of 88 putative AP2/ERF factors were identified from the transcriptome sequence of parsley. Seven AP2/ERF transcription factors were selected in this study to analyze the expression profiles of parsley under different abiotic stresses. Our data provide a potentially valuable resource that can be used for intensive parsley research.

Genetic mapping, marker assisted selection and allelic relationships for the Pu 6 gene conferring rust resistance in sunflower

Rust resistance in the sunflower line P386 is controlled by Pu 6 , a gene which was reported to segregate independently from other rust resistant genes, such as R 4 . The objectives of this work were to map Pu 6 , to provide and validate molecular tools for its identification, and to determine the linkage relationship of Pu 6 and R 4 . Genetic mapping of Pu 6 with six markers covered 24.8 cM of genetic distance on the lower end of linkage Group 13 of the sunflower consensus map. The marker most closely linked to Pu 6 was ORS316 at 2.5 cM in the distal position. ORS316 presented five alleles when was assayed with a representative set of resistant and susceptible lines. Allelism test between Pu 6 and R 4 indicated that both genes are linked at a genetic distance of 6.25 cM. This is the first confirmation based on an allelism test that at least two members of the R adv /R 4 /R 11 / R 13a /R 13b /Pu 6 cluster of genes are at different loci. A fine elucidation of the architecture of this complex locus will allow designing and constructing completely new genomic regions combining genes from different resistant sources and the elimination of the linkage drag around each resistant gene.

Genome wide SNP identification in chickpea for use in development of a high density genetic map and improvement of chickpea reference genome assembly

BACKGROUND

In the whole genome sequencing, genetic map provides an essential framework for accurate and efficient genome assembly and validation. The main objectives of this study were to develop a high-density genetic map using RAD-Seq (Restriction-site Associated DNA Sequencing) genotyping-by-sequencing (RAD-Seq GBS) and Illumina GoldenGate assays, and to examine the alignment of the current map with the kabuli chickpea genome assembly.

RESULTS

Genic single nucleotide polymorphisms (SNPs) totaling 51,632 SNPs were identified by 454 transcriptome sequencing of Cicer arietinum and Cicer reticulatum genotypes. Subsequently, an Illumina GoldenGate assay for 1,536 SNPs was developed. A total of 1,519 SNPs were successfully assayed across 92 recombinant inbred lines (RILs), of which 761 SNPs were polymorphic between the two parents. In addition, the next generation sequencing (NGS)-based GBS was applied to the same population generating 29,464 high quality SNPs. These SNPs were clustered into 626 recombination bins based on common segregation patterns. Data from the two approaches were used for the construction of a genetic map using a population derived from an intraspecific cross. The map consisted of 1,336 SNPs including 604 RAD recombination bins and 732 SNPs from Illumina GoldenGate assay. The map covered 653 cM of the chickpea genome with an average distance between adjacent markers of 0.5 cM. To date, this is the most extensive genetic map of chickpea using an intraspecific population. The alignment of the map with the CDC Frontier genome assembly revealed an overall conserved marker order; however, a few local inconsistencies within the Cicer arietinum pseudochromosome 1 (Ca1), Ca5 and Ca8 were detected. The map enabled the alignment of 215 unplaced scaffolds from the CDC Frontier draft genome assembly. The alignment also revealed varying degrees of recombination rates and hotspots across the chickpea genome.

CONCLUSIONS

A high-density genetic map using RAD-Seq GBS and Illumina GoldenGate assay was developed and aligned with the existing kabuli chickpea draft genome sequence. The analysis revealed an overall conserved marker order, although some localized inversions between draft genome assembly and the genetic map were detected. The current analysis provides an insight of the recombination rates and hotspots across the chickpea genome.

A high-throughput SNP array in the amphidiploid species Brassica napus shows diversity in resistance genes

Single-nucleotide polymorphisms (SNPs)are molecular markers based on nucleotide variation and can be used for genotyping assays across populations and to track genomic inheritance. SNPs offer a comprehensive genotyping alternative to whole-genome sequencing for both agricultural and research purposes including molecular breeding and diagnostics, genome evolution and genetic diversity analyses, genetic mapping, and trait association studies. Here genomic SNPs were discovered between four cultivars of the important amphidiploid oilseed species Brassica napus and used to develop a B. napus Infinium™ array containing 5,306 SNPs randomly dispersed across the genome. Assay success was high, with >94 % of these producing a reproducible, polymorphic genotype in the 1,070 samples screened. Although the assay was designed to B. napus, successful SNP amplification was achieved in the B. napus progenitor species, Brassica rapa and Brassica oleracea, and to a lesser extent in the related species Brassica nigra. Phylogenetic analysis was consistent with the expected relationships between B. napus individuals. This study presents an efficient custom SNP assay development pipeline in the complex polyploid Brassica genome and demonstrates the utility of the array for high-throughput genotyping in a number of related Brassica species. It also demonstrates the utility of this assay in genotyping resistance genes on chromosome A7, which segregate amongst the 1,070 samples.

Identification of candidate genes associated with leaf senescence in cultivated sunflower (Helianthus annuus L.)

Cultivated sunflower (Helianthus annuus L.), an important source of edible vegetable oil, shows rapid onset of senescence, which limits production by reducing photosynthetic capacity under specific growing conditions. Carbon for grain filling depends strongly on light interception by green leaf area, which diminishes during grain filling due to leaf senescence. Transcription factors (TFs) regulate the progression of leaf senescence in plants and have been well explored in model systems, but information for many agronomic crops remains limited. Here, we characterize the expression profiles of a set of putative senescence associated genes (SAGs) identified by a candidate gene approach and sunflower microarray expression studies. We examined a time course of sunflower leaves undergoing natural senescence and used quantitative PCR (qPCR) to measure the expression of 11 candidate genes representing the NAC, WRKY, MYB and NF-Y TF families. In addition, we measured physiological parameters such as chlorophyll, total soluble sugars and nitrogen content. The expression of Ha-NAC01, Ha-NAC03, Ha-NAC04, Ha-NAC05 and Ha-MYB01 TFs increased before the remobilization rate increased and therefore, before the appearance of the first physiological symptoms of senescence, whereas Ha-NAC02 expression decreased. In addition, we also examined the trifurcate feed-forward pathway (involving ORE1, miR164, and ethylene insensitive 2) previously reported for Arabidopsis. We measured transcription of Ha-NAC01 (the sunflower homolog of ORE1) and Ha-EIN2, along with the levels of miR164, in two leaves from different stem positions, and identified differences in transcription between basal and upper leaves. Interestingly, Ha-NAC01 and Ha-EIN2 transcription profiles showed an earlier up-regulation in upper leaves of plants close to maturity, compared with basal leaves of plants at pre-anthesis stages. These results suggest that the H. annuus TFs characterized in this work could play important roles as potential triggers of leaf senescence and thus can be considered putative candidate genes for senescence in sunflower.

Selection on crop-derived traits and QTL in sunflower (Helianthus annuus) crop-wild hybrids under water stress

Locally relevant conditions, such as water stress in irrigated agricultural regions, should be considered when assessing the risk of crop allele introgression into wild populations following hybridization. Although research in cultivars has suggested that domestication traits may reduce fecundity under water stress as compared to wild-like phenotypes, this has not been investigated in crop-wild hybrids. In this study, we examine phenotypic selection acting on, as well as the genetic architecture of vegetative, reproductive, and physiological characteristics in an experimental population of sunflower crop-wild hybrids grown under wild-like low water conditions. Crop-derived petiole length and head diameter were favored in low and control water environments. The direction of selection differed between environments for leaf size and leaf pressure potential. Interestingly, the additive effect of the crop-derived allele was in the direction favored by selection for approximately half the QTL detected in the low water environment. Selection favoring crop-derived traits and alleles in the low water environment suggests that a subset of these alleles would be likely to spread into wild populations under water stress. Furthermore, differences in selection between environments support the view that risk assessments should be conducted under multiple locally relevant conditions.

A high-density SNP Map of sunflower derived from RAD-sequencing facilitating fine-mapping of the rust resistance gene R12

A high-resolution genetic map of sunflower was constructed by integrating SNP data from three F₂ mapping populations (HA 89/RHA 464, B-line/RHA 464, and CR 29/RHA 468). The consensus map spanned a total length of 1443.84 cM, and consisted of 5,019 SNP markers derived from RAD tag sequencing and 118 publicly available SSR markers distributed in 17 linkage groups, corresponding to the haploid chromosome number of sunflower. The maximum interval between markers in the consensus map is 12.37 cM and the average distance is 0.28 cM between adjacent markers. Despite a few short-distance inversions in marker order, the consensus map showed high levels of collinearity among individual maps with an average Spearman's rank correlation coefficient of 0.972 across the genome. The order of the SSR markers on the consensus map was also in agreement with the order of the individual map and with previously published sunflower maps. Three individual and one consensus maps revealed the uneven distribution of markers across the genome. Additionally, we performed fine mapping and marker validation of the rust resistance gene R₁₂, providing closely linked SNP markers for marker-assisted selection of this gene in sunflower breeding programs. This high resolution consensus map will serve as a valuable tool to the sunflower community for studying marker-trait association of important agronomic traits, marker assisted breeding, map-based gene cloning, and comparative mapping.

Molecular evolution of candidate genes for crop-related traits in sunflower (Helianthus annuus L.)

Evolutionary analyses aimed at detecting the molecular signature of selection during crop domestication and/or improvement can be used to identify genes or genomic regions of likely agronomic importance. Here, we describe the DNA sequence-based characterization of a pool of candidate genes for crop-related traits in sunflower. These genes, which were identified based on homology to genes of known effect in other study systems, were initially sequenced from a panel of improved lines. All genes that exhibited a paucity of sequence diversity, consistent with the possible effects of selection during the evolution of cultivated sunflower, were then sequenced from a panel of wild sunflower accessions an outgroup. These data enabled formal tests for the effects of selection in shaping sequence diversity at these loci. When selection was detected, we further sequenced these genes from a panel of primitive landraces, thereby allowing us to investigate the likely timing of selection (i.e., domestication vs. improvement). We ultimately identified seven genes that exhibited the signature of positive selection during either domestication or improvement. Genetic mapping of a subset of these genes revealed co-localization between candidates for genes involved in the determination of flowering time, seed germination, plant growth/development, and branching and QTL that were previously identified for these traits in cultivated × wild sunflower mapping populations.

Chromosomal evolution and patterns of introgression in helianthus

Knowledge of the nature and extent of karyotypic differences between species provides insight into the evolutionary history of the genomes in question and, in the case of closely related species, the potential for genetic exchange between taxa. We constructed high-density genetic maps of the silverleaf sunflower (Helianthus argophyllus) and Algodones Dune sunflower (H. niveus ssp. tephrodes) genomes and compared them to a consensus map of cultivated sunflower (H. annuus) to identify chromosomal rearrangements between species. The genetic maps of H. argophyllus and H. niveus ssp. tephrodes included 17 linkage groups each and spanned 1337 and 1478 cM, respectively. Comparative analyses revealed greater divergence between H. annuus and H. niveus ssp. tephrodes (13 inverted segments, 18 translocated segments) than between H. annuus and H. argophyllus (10 inverted segments, 8 translocated segments), consistent with their known phylogenetic relationships. Marker order was conserved across much of the genome, with 83 and 64% of the H. argophyllus and H. niveus ssp. tephrodes genomes, respectively, being syntenic with H. annuus. Population genomic analyses between H. annuus and H. argophyllus, which are sympatric across a portion of the natural range of H. annuus, revealed significantly elevated genetic structure in rearranged portions of the genome, indicating that such rearrangements are associated with restricted gene flow between these two species.

A high-density SNP genotyping array for rice biology and molecular breeding

A high-density single nucleotide polymorphism (SNP) array is critically important for geneticists and molecular breeders. With the accumulation of huge amounts of genomic re-sequencing data and available technologies for accurate SNP detection, it is possible to design high-density and high-quality rice SNP arrays. Here we report the development of a high-density rice SNP array and its utility. SNP probes were designed by screening more than 10 000 000 SNP loci extracted from the re-sequencing data of 801 rice varieties and an array named RiceSNP50 was produced on the Illumina Infinium platform. The array contained 51 478 evenly distributed markers, 68% of which were within genic regions. Several hundred rice plants with parent/F1 relationships were used to generate a high-quality cluster file for accurate SNP calling. Application tests showed that this array had high genotyping accuracy, and could be used for different objectives. For example, a core collection of elite rice varieties was clustered with fine resolution. Genome-wide association studies (GWAS) analysis correctly identified a characterized QTL. Further, this array was successfully used for variety verification and trait introgression. As an accurate high-throughput genotyping tool, RiceSNP50 will play an important role in both functional genomics studies and molecular breeding.

Single-nucleotide polymorphism markers from de-novo assembly of the pomegranate transcriptome reveal germplasm genetic diversity

Pomegranate is a valuable crop that is grown commercially in many parts of the world. Wild species have been reported from India, Turkmenistan and Socotra. Pomegranate fruit has a variety of health-beneficial qualities. However, despite this crop's importance, only moderate effort has been invested in studying its biochemical or physiological properties or in establishing genomic and genetic infrastructures. In this study, we reconstructed a transcriptome from two phenotypically different accessions using 454-GS-FLX Titanium technology. These data were used to explore the functional annotation of 45,187 fully annotated contigs. We further compiled a genetic-variation resource of 7,155 simple-sequence repeats (SSRs) and 6,500 single-nucleotide polymorphisms (SNPs). A subset of 480 SNPs was sampled to investigate the genetic structure of the broad pomegranate germplasm collection at the Agricultural Research Organization (ARO), which includes accessions from different geographical areas worldwide. This subset of SNPs was found to be polymorphic, with 10.7% loci with minor allele frequencies of (MAF<0.05). These SNPs were successfully used to classify the ARO pomegranate collection into two major groups of accessions: one from India, China and Iran, composed of mainly unknown country origin and which was more of an admixture than the other major group, composed of accessions mainly from the Mediterranean basin, Central Asia and California. This study establishes a high-throughput transcriptome and genetic-marker infrastructure. Moreover, it sheds new light on the genetic interrelations between pomegranate species worldwide and more accurately defines their genetic nature.

Genome-wide DNA polymorphisms in seven rice cultivars of temperate and tropical japonica groups

Elucidation of the rice genome is expected to broaden our understanding of genes related to the agronomic characteristics and the genetic relationship among cultivars. In this study, we conducted whole-genome sequencings of 6 cultivars, including 5 temperate japonica cultivars and 1 tropical japonica cultivar (Moroberekan), by using next-generation sequencing (NGS) with Nipponbare genome as a reference. The temperate japonica cultivars contained 2 sake brewing (Yamadanishiki and Gohyakumangoku), 1 landrace (Kameji), and 2 modern cultivars (Koshihikari and Norin 8). Almost >83% of the whole genome sequences of the Nipponbare genome could be covered by sequenced short-reads of each cultivar, including Omachi, which has previously been reported to be a temperate japonica cultivar. Numerous single nucleotide polymorphisms (SNPs), insertions, and deletions were detected among the various cultivars and the Nipponbare genomes. Comparison of SNPs detected in each cultivar suggested that Moroberekan had 5-fold more SNPs than the temperate japonica cultivars. Success of the 2 approaches to improve the efficacy of sequence data by using NGS revealed that sequencing depth was directly related to sequencing coverage of coding DNA sequences: in excess of 30× genome sequencing was required to cover approximately 80% of the genes in the rice genome. Further, the contigs prepared using the assembly of unmapped reads could increase the value of NGS short-reads and, consequently, cover previously unavailable sequences. These approaches facilitated the identification of new genes in coding DNA sequences and the increase of mapping efficiency in different regions. The DNA polymorphism information between the 7 cultivars and Nipponbare are available at NGRC_Rices_Build1.0 (http://www.nodai-genome.org/oryza_sativa_en.html).