SNP frequency, haplotype structure and linkage disequilibrium in elite maize inbred lines

Maize mutant screens: from classical methods to new CRISPR-based approaches

Mutations play a pivotal role in shaping the trajectory and outcomes of a species evolution and domestication. Maize (Zea mays) has been a major staple crop and model for genetic research for more than 100 yr. With the arrival of site-directed mutagenesis and genome editing (GE) driven by the Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR), maize mutational research is once again in the spotlight. If we combine the powerful physiological and genetic characteristics of maize with the already available and ever increasing toolbox of CRISPR-Cas, prospects for its future trait engineering are very promising. This review aimed to give an overview of the progression and learnings of maize screening studies analyzing forward genetics, natural variation and reverse genetics to focus on recent GE approaches. We will highlight how each strategy and resource has contributed to our understanding of maize natural and induced trait variability and how this information could be used to design the next generation of mutational screenings.

Genome-wide association study reveals SNP markers controlling drought tolerance and related agronomic traits in chickpea across multiple environments

Chickpea, renowned for its exceptional nutritional value, stands as a crucial crop, serving as a dietary staple in various parts of the world. However, its productivity faces a significant challenge in the form of drought stress. This challenge highlights the urgent need to find genetic markers linked to drought tolerance for effective breeding programs. The primary objective of this study is to identify genetic markers associated with drought tolerance to facilitate effective breeding programs. To address this, we cultivated 185 chickpea accessions in two distinct locations in Lebanon over a two-year period, subjecting them to both irrigated and rain-fed environments. We assessed 11 drought-linked traits, including morphology, growth, yield, and tolerance score. SNP genotyping revealed 1344 variable SNP markers distributed across the chickpea genome. Genetic diversity across populations originating from diverse geographic locations was unveiled by the PCA, clustering, and structure analysis indicating that these genotypes have descend from five or four distinct ancestors. A genome-wide association study (GWAS) revealed several marker trait associations (MTAs) associated with the traits evaluated. Within the rainfed conditions, 11 significant markers were identified, each associated with distinct chickpea traits. Another set of 11 markers exhibited associations in both rainfed and irrigated environments, reflecting shared genetic determinants across these conditions for the same trait. The analysis of linkage disequilibrium (LD) highlighted two genomic regions with notably strong LD, suggesting significant interconnections among several investigated traits. This was further investigated by the correlation between major markers associated with these traits. Gene annotation of the identified markers has unveiled insights into 28 potential genes that play a role in influencing various chickpea drought-linked traits. These traits encompass crucial aspects such as blooming organ development, plant growth, seed weight, starch metabolism, drought regulation, and height index. Among the identified genes are CPN60-2, hsp70, GDSL(GELP), AHL16, NAT3, FAB1B, bZIP, and GL21. These genes collectively contribute to the multifaceted response of chickpea plants to drought stress. Our identified genetic factors exert their influence in both irrigated and rainfed environments, emphasizing their importance in shaping chickpea characteristics.

ZmARF1 positively regulates low phosphorus stress tolerance via modulating lateral root development in maize

Phosphorus (P) deficiency is one of the most critical factors for plant growth and productivity, including its inhibition of lateral root initiation. Auxin response factors (ARFs) play crucial roles in root development via auxin signaling mediated by genetic pathways. In this study, we found that the transcription factor ZmARF1 was associated with low inorganic phosphate (Pi) stress-related traits in maize. This superior root morphology and greater phosphate stress tolerance could be ascribed to the overexpression of ZmARF1. The knock out mutant zmarf1 had shorter primary roots, fewer root tip number, and lower root volume and surface area. Transcriptomic data indicate that ZmLBD1, a direct downstream target gene, is involved in lateral root development, which enhances phosphate starvation tolerance. A transcriptional activation assay revealed that ZmARF1 specifically binds to the GC-box motif in the promoter of ZmLBD1 and activates its expression. Moreover, ZmARF1 positively regulates the expression of ZmPHR1, ZmPHT1;2, and ZmPHO2, which are key transporters of Pi in maize. We propose that ZmARF1 promotes the transcription of ZmLBD1 to modulate lateral root development and Pi-starvation induced (PSI) genes to regulate phosphate mobilization and homeostasis under phosphorus starvation. In addition, ZmERF2 specifically binds to the ABRE motif of the promoter of ZmARF1 and represses its expression. Collectively, the findings of this study revealed that ZmARF1 is a pivotal factor that modulates root development and confers low-Pi stress tolerance through the transcriptional regulation of the biological function of ZmLBD1 and the expression of key Pi transport proteins.

Genome-wide association study of maize resistance to Pythium aristosporum stalk rot

Stalk rot, a severe and widespread soil-borne disease in maize, globally reduces yield and quality. Recent documentation reveals that Pythium aristosporum has emerged as one of the dominant causal agents of maize stalk rot. However, a previous study of maize stalk rot disease resistance mechanisms and breeding had mainly focused on other pathogens, neglecting P. aristosporum. To mitigate crop loss, resistance breeding is the most economical and effective strategy against this disease. This study involved characterizing resistance in 295 inbred lines using the drilling inoculation method and genotyping them via sequencing. By combining with population structure, disease resistance phenotype, and genome-wide association study (GWAS), we identified 39 significant single-nucleotide polymorphisms (SNPs) associated with P. aristosporum stalk rot resistance by utilizing six statistical methods. Bioinformatics analysis of these SNPs revealed 69 potential resistance genes, among which Zm00001d051313 was finally evaluated for its roles in host defense response to P. aristosporum infection. Through virus-induced gene silencing (VIGS) verification and physiological index determination, we found that transient silencing of Zm00001d051313 promoted P. aristosporum infection, indicating a positive regulatory role of this gene in maize's antifungal defense mechanism. Therefore, these findings will help advance our current understanding of the underlying mechanisms of maize defense to Pythium stalk rot.

Whole-genome sequencing and variant discovery of Citrus reticulata "Kinnow" from Pakistan

Citrus is a source of nutritional and medicinal advantages, cultivated worldwide with major groups of sweet oranges, mandarins, grapefruits, kumquats, lemons and limes. Pakistan produces all major citrus groups with mandarin (Citrus reticulata) being the prominent group that includes local commercial cultivars Feutral's Early, Dancy, Honey, and Kinnow. The present study designed to understand the genetic architecture of this unique variety of Citrus reticulata 'Kinnow.' The whole-genome resequencing and variant calling was performed to map the genomic variability that might be responsible for its particular characteristics like taste, seedlessness, juice content, thickness of peel, and shelf-life. A total of 139,436,350 raw sequence reads were generated with 20.9 Gb data in Fastq format having 98% effectiveness and 0.2% base call error rate. Overall, 3,503,033 SNPs, 176,949 MNPs, 323,287 INS, and 333,083 DEL were identified using the GATK4 variant calling pipeline against Citrus clementina. Furthermore, g:Profiler was applied for annotating the newly found variants, harbor genes/transcripts and their involved pathways. A total of 73,864 transcripts harbors 4,336,352 variants, most of the observed variants were predicted in non-coding regions and 1009 transcripts were found well annotated by different databases. Out of total aforementioned transcripts, 588 involved in biological processes, 234 in molecular functions and 167 transcripts in cellular components. In a nutshell, 18,153 high impact variants and 216 genic variants found in the current study, which may be used after its functional validation for marker-assisted breeding programs of "Kinnow" to propagate its valued traits for the improvement of contemporary citrus varieties in the region.

Genotyping by Sequencing (GBS) for Genome-Wide SNP Identification in Plants

Marker-assisted selection has played a pivotal role in developing several elite varieties in the past two decades. Molecular markers employed in plant breeding programs have recently shifted from microsatellites or simple sequence repeats (SSRs) to single nucleotide polymorphisms (SNPs) due to the ubiquity of SNP markers in the genome and the availability of various high-throughput SNP genotyping platforms. Rapid advances in sequencing technologies and the reduction in sequencing cost have facilitated SNP discovery in several plant species including non-model organisms with little or no genomic resources. Despite the lower cost of sequencing, genome complexity reduction approaches are still useful for SNP identification because many applications do not require every base of the genome to be sequenced. Genotyping-by-sequencing (GBS) is a quick and affordable reduced representation method that can simultaneously identify and genotype a large number of SNPs that has been successfully applied to a wide range of plant species. This chapter describes a robust two-enzyme GBS method for SNP discovery and genotyping that has been verified in non-model plant species.

cgSSR marker-based genome-wide association study identified genomic regions for panicle characters and yield in rice (Oryza sativa L.)

BACKGROUND

To improve production efficiency, positive alleles corresponding to yield-related attributes must be accumulated in a single elite background. We designed and used cgSSR markers, which are superior to random SSR markers in genome-wide association study, to identify genomic regions that contribute to panicle characters and grain yield in this study.

RESULTS

As evidenced by the high polymorphic information content value and gene diversity coefficient, the new cgSSR markers were determined to be highly informative. These cgSSR markers were employed to generate genotype data for an association panel evaluated for four panicle characters and grain yield over three seasons. For five traits, 17 significant marker-trait associations on six chromosomes were discovered. The percentage of phenotypic variance that could be explained ranged from 4% to 13%. Unrelated gene-derived markers had a strong association with target traits as well.

CONCLUSION

Trait-associated cgSSR markers derived from corresponding or related genes ensure their utility in direct allele selection, while other linked markers aid in allele selection indirectly by altering the phenotype of interest. Through a marker-assisted breeding approach, these marker-trait associations can be leveraged to accumulate favourable alleles for yield enhancement in rice. © 2022 Society of Chemical Industry.

Identification of Genomic Regions Associated with Agronomic and Disease Resistance Traits in a Large Set of Multiple DH Populations

Breeding maize lines with the improved level of desired agronomic traits under optimum and drought conditions as well as increased levels of resistance to several diseases such as maize lethal necrosis (MLN) is one of the most sustainable approaches for the sub-Saharan African region. In this study, 879 doubled haploid (DH) lines derived from 26 biparental populations were evaluated under artificial inoculation of MLN, as well as under well-watered (WW) and water-stressed (WS) conditions for grain yield and other agronomic traits. All DH lines were used for analyses of genotypic variability, association studies, and genomic predictions for the grain yield and other yield-related traits. Genome-wide association study (GWAS) using a mixed linear FarmCPU model identified SNPs associated with the studied traits i.e., about seven and eight SNPs for the grain yield; 16 and 12 for anthesis date; seven and eight for anthesis silking interval; 14 and 5 for both ear and plant height; and 15 and 5 for moisture under both WW and WS environments, respectively. Similarly, about 13 and 11 SNPs associated with gray leaf spot and turcicum leaf blight were identified. Eleven SNPs associated with senescence under WS management that had depicted drought-stress-tolerant QTLs were identified. Under MLN artificial inoculation, a total of 12 and 10 SNPs associated with MLN disease severity and AUDPC traits, respectively, were identified. Genomic prediction under WW, WS, and MLN disease artificial inoculation revealed moderate-to-high prediction accuracy. The findings of this study provide useful information on understanding the genetic basis for the MLN resistance, grain yield, and other agronomic traits under MLN artificial inoculation, WW, and WS conditions. Therefore, the obtained information can be used for further validation and developing functional molecular markers for marker-assisted selection and for implementing genomic prediction to develop superior elite lines.

Multiomics approach reveals a role of translational machinery in shaping maize kernel amino acid composition

Maize (Zea mays) seeds are a good source of protein, despite being deficient in several essential amino acids. However, eliminating the highly abundant but poorly balanced seed storage proteins has revealed that the regulation of seed amino acids is complex and does not rely on only a handful of proteins. In this study, we used two complementary omics-based approaches to shed light on the genes and biological processes that underlie the regulation of seed amino acid composition. We first conducted a genome-wide association study to identify candidate genes involved in the natural variation of seed protein-bound amino acids. We then used weighted gene correlation network analysis to associate protein expression with seed amino acid composition dynamics during kernel development and maturation. We found that almost half of the proteome was significantly reduced during kernel development and maturation, including several translational machinery components such as ribosomal proteins, which strongly suggests translational reprogramming. The reduction was significantly associated with a decrease in several amino acids, including lysine and methionine, pointing to their role in shaping the seed amino acid composition. When we compared the candidate gene lists generated from both approaches, we found a nonrandom overlap of 80 genes. A functional analysis of these genes showed a tight interconnected cluster dominated by translational machinery genes, especially ribosomal proteins, further supporting the role of translation dynamics in shaping seed amino acid composition. These findings strongly suggest that seed biofortification strategies that target the translation machinery dynamics should be considered and explored further.

High-Throughput and Accurate Determination of Transgene Copy Number and Zygosity in Transgenic Maize: From DNA Extraction to Data Analysis

It is vital to develop high-throughput methods to determine transgene copy numbers initially and zygosity during subsequent breeding. In this study, the target sequence of the previously reported endogenous reference gene hmg was analyzed using 633 maize inbred lines, and two SNPs were observed. These SNPs significantly increased the PCR efficiency, while the newly developed hmg gene assay (hmg-taq-F2/R2) excluding these SNPs reduced the efficiency into normal ranges. The TaqMan amplification efficiency of bar and hmg with newly developed primers was calculated as 0.993 and 1.000, respectively. The inter-assay coefficient of variation (CV) values for the bar and hmg genes varied from 1.18 to 2.94%. The copy numbers of the transgene bar using new TaqMan assays were identical to those using dPCR. Significantly, the precision of one repetition reached 96.7% of that of three repetitions of single-copy plants analyzed by simple random sampling, and the actual accuracy reached 95.8%, confirmed by T₁ and T₂ progeny. With the high-throughput DNA extraction and automated data analysis procedures developed in this study, nearly 2700 samples could be analyzed within eight hours by two persons. The combined results suggested that the new hmg gene assay developed here could be a universal maize reference gene system, and the new assay has high throughput and high accuracy for large-scale screening of maize varieties around the world.

Development of high-resolution multiple-SNP arrays for genetic analyses and molecular breeding through genotyping by target sequencing and liquid chip

Genotyping platforms, as critical supports for genomics, genetics, and molecular breeding, have been well implemented at national institutions/universities in developed countries and multinational seed companies that possess high-throughput, automatic, large-scale, and shared facilities. In this study, we integrated an improved genotyping by target sequencing (GBTS) system with capture-in-solution (liquid chip) technology to develop a multiple single-nucleotide polymorphism (mSNP) approach in which mSNPs can be captured from a single amplicon. From one 40K maize mSNP panel, we developed three types of markers (40K mSNPs, 251K SNPs, and 690K haplotypes), and generated multiple panels with various marker densities (1K-40K mSNPs) by sequencing at different depths. Comparative genetic diversity analysis was performed with genic versus intergenic markers and di-allelic SNPs versus non-typical SNPs. Compared with the one-amplicon-one-SNP system, mSNPs and within-mSNP haplotypes are more powerful for genetic diversity detection, linkage disequilibrium decay analysis, and genome-wide association studies. The technologies, protocols, and application scenarios developed for maize in this study will serve as a model for the development of mSNP arrays and highly efficient GBTS systems in animals, plants, and microorganisms.

Gene-based mapping of trehalose biosynthetic pathway genes reveals association with source- and sink-related yield traits in a spring wheat panel

Trehalose 6‐phosphate (T6P) signalling regulates carbon use and allocation and is a target to improve crop yields. However, the specific contributions of trehalose phosphate synthase (TPS) and trehalose phosphate phosphatase (TPP) genes to source‐ and sink‐related traits remain largely unknown. We used enrichment capture sequencing on TPS and TPP genes to estimate and partition the genetic variation of yield‐related traits in a spring wheat (Triticum aestivum) breeding panel specifically built to capture the diversity across the 75,000 CIMMYT wheat cultivar collection. Twelve phenotypes were correlated to variation in TPS and TPP genes including plant height and biomass (source), spikelets per spike, spike growth and grain filling traits (sink) which showed indications of both positive and negative gene selection. Individual genes explained proportions of heritability for biomass and grain‐related traits. Three TPS1 homologues were particularly significant for trait variation. Epistatic interactions were found within and between the TPS and TPP gene families for both plant height and grain‐related traits. Gene‐based prediction improved predictive ability for grain weight when gene effects were combined with the whole‐genome markers. Our study has generated a wealth of information on natural variation of TPS and TPP genes related to yield potential which confirms the role for T6P in resource allocation and in affecting traits such as grain number and size confirming other studies which now opens up the possibility of harnessing natural genetic variation more widely to better understand the contribution of native genes to yield traits for incorporation into breeding programmes.

Genome-Wide Scanning Enabled SNP Discovery, Linkage Disequilibrium Patterns and Population Structure in a Panel of Fonio (Digitaria exilis [Kippist] Stapf) Germplasm

White fonio (Digitaria exilis) is a staple food for millions of people in arid and semi-arid areas of West Africa. Knowledge about nutritional and health benefits, insights into morphological diversity, and the recent development of genomic resources call for a better understanding of the genetic structure of the extant germplasm gathered throughout the region in order to set up a robust breeding program. We assessed the genetic diversity and population structure of 259 fonio individuals collected from six countries from West Africa (Nigeria, Benin, Guinea, Mali, Burkina Faso and Niger) in this study using 688 putative out of 21,324 DArTseq-derived SNP markers. Due to the inbreeding and small population size, the results revealed a substantial level of genetic variability. Furthermore, two clusters were found irrespective of the geographic origins of accessions. Moreover, the high level of linkage disequilibrium (LD) between loci observed resulted from the mating system of the crop, which is often associated with a low recombination rate. These findings fill the gaps about the molecular diversity and genetic structure of the white fonio germplasm in West Africa. This was required for the application of genomic tools that can potentially speed up the genetic gain in fonio millet breeding for complex traits such as yield, and other nutrient contents.

Genetic diversity and selection signatures in maize landraces compared across 50 years of in situ and ex situ conservation

Genomics-based, longitudinal comparisons between ex situ and in situ agrobiodiversity conservation strategies can contribute to a better understanding of their underlying effects. However, landrace designations, ambiguous common names, and gaps in sampling information complicate the identification of matching ex situ and in situ seed lots. Here we report a 50-year longitudinal comparison of the genetic diversity of a set of 13 accessions from the state of Morelos, Mexico, conserved ex situ since 1967 and retrieved in situ from the same donor families in 2017. We interviewed farmer families who donated in situ landraces to understand their germplasm selection criteria. Samples were genotyped by sequencing, producing 74,739 SNPs. Comparing the two sample groups, we show that ex situ and in situ genome-wide diversity was similar. In situ samples had 3.1% fewer SNPs and lower pairwise genetic distances (F_st 0.008-0.113) than ex situ samples (F_st 0.031-0.128), but displayed the same heterozygosity. Despite genome-wide similarities across samples, we could identify several loci under selection when comparing in situ and ex situ seed lots, suggesting ongoing evolution in farmer fields. Eight loci in chromosomes 3, 5, 6, and 10 showed evidence of selection in situ that could be related with farmers' selection criteria surveyed with focus groups and interviews at the sampling site in 2017, including wider kernels and larger ear size. Our results have implications for ex situ collection resampling strategies and the in situ conservation of threatened landraces.

Zea mays RNA-seq estimated transcript abundances are strongly affected by read mapping bias

Background

Genetic variation for gene expression is a source of phenotypic variation for natural and agricultural species. The common approach to map and to quantify gene expression from genetically distinct individuals is to assign their RNA-seq reads to a single reference genome. However, RNA-seq reads from alleles dissimilar to this reference genome may fail to map correctly, causing transcript levels to be underestimated. Presently, the extent of this mapping problem is not clear, particularly in highly diverse species. We investigated if mapping bias occurred and if chromosomal features associated with mapping bias. Zea mays presents a model species to assess these questions, given it has genotypically distinct and well-studied genetic lines.

Results

In Zea mays, the inbred B73 genome is the standard reference genome and template for RNA-seq read assignments. In the absence of mapping bias, B73 and a second inbred line, Mo17, would each have an approximately equal number of regulatory alleles that increase gene expression. Remarkably, Mo17 had 2–4 times fewer such positively acting alleles than did B73 when RNA-seq reads were aligned to the B73 reference genome. Reciprocally, over one-half of the B73 alleles that increased gene expression were not detected when reads were aligned to the Mo17 genome template. Genes at dissimilar chromosomal ends were strongly affected by mapping bias, and genes at more similar pericentromeric regions were less affected. Biased transcript estimates were higher in untranslated regions and lower in splice junctions. Bias occurred across software and alignment parameters.

Conclusions

Mapping bias very strongly affects gene transcript abundance estimates in maize, and bias varies across chromosomal features. Individual genome or transcriptome templates are likely necessary for accurate transcript estimation across genetically variable individuals in maize and other species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07577-3.

Genetic diversity and population structure of advanced clones selected over forty years by a potato breeding program in the USA

Knowledge regarding genetic diversity and population structure of breeding materials is essential for crop improvement. The Texas A&M University Potato Breeding Program has a collection of advanced clones selected and maintained in-vitro over a 40-year period. Little is known about its genetic makeup and usefulness for the current breeding program. In this study, 214 potato clones were genotyped with the Infinium Illumina 22 K V3 Potato Array. After filtering, a total of 10,106 single nucleotide polymorphic (SNP) markers were used for analysis. Heterozygosity varied by SNP, with an overall average of 0.59. Three groups of tetraploid clones primarily based on potato market classes, were detected using STRUCTURE software and confirmed by discriminant analysis of principal components. The highest coefficient of differentiation observed between the groups was 0.14. Signatures of selection were uncovered in genes controlling potato flesh and skin color, length of plant cycle and tuberization, and carbohydrate metabolism. A core set of 43 clones was obtained using Core Hunter 3 to develop a sub-collection that retains similar genetic diversity as the whole population, minimize redundancies, and facilitates long-term conservation of genetic resources. The comprehensive molecular characterization of our breeding clone bank collection contributes to understanding the genetic diversity of existing potato resources. This analysis could be applied to other breeding programs and assist in the selection of parents, fingerprinting, protection, and management of the breeding collections.

Genetic variation and population structure in China summer maize germplasm

Maize (Zea mays L.) germplasm in China Summer maize ecological region (CSM) or central corn-belt of China is diverse but has not been systematically characterized at molecular level. In this study, genetic variation, genome diversity, linkage disequilibrium patterns, population structure, and characteristics of different heterotic groups were studied using 525,141 SNPs obtained by Genotyping-By-Sequencing (GBS) for 490 inbred lines collected from researchers at CSM region. The SNP density is lower near centromere, but higher near telomere region of maize chromosome, the degree of linkage disequilibrium (r2) vary at different chromosome regions. Majority of the inbred lines (66.05%) show pairwise relative kinship near zero, indicating a large genetic diversity in the CSM breeding germplasm. Using 4849 tagSNPs derived from 3618 haplotype blocks, the 490 inbred lines were delineated into 3 supergroups, 6 groups, and 10 subgroups using ADMIXTURE software. A procedure of assigning inbred lines into heterotic groups using genomic data and tag-SNPs was developed and validated. Genome differentiation among different subgroups measured by Fst, and the genetic diversity within each subgroup measured by GD are both large. The share of heterotic groups that have significant North American germplasm contribution: P, SS, IDT, and X, accounts about 54% of the CSM breeding germplasm collection and has increased significantly in the last two decades. Two predominant types of heterotic pattern in CSM region are: M-Reid group × TSPT group, and X subgroup × Local subgroups.

Founder effects shape linkage disequilibrium and genomic diversity of a partially clonal invader

The genomic variation of an invasive species may be affected by complex demographic histories and evolutionary changes during the invasion. Here, we describe the relative influence of bottlenecks, clonality, and population expansion in determining genomic variability of the widespread red macroalga Agarophyton vermiculophyllum. Its introduction from mainland Japan to the estuaries of North America and Europe coincided with shifts from predominantly sexual to partially clonal reproduction and rapid adaptive evolution. A survey of 62,285 SNPs for 351 individuals from 35 populations, aligned to 24 chromosome-length scaffolds indicate that linkage disequilibrium (LD), observed heterozygosity (H_o ), Tajima's D, and nucleotide diversity (Pi) were greater among non-native than native populations. Evolutionary simulations indicate LD and Tajima's D were consistent with a severe population bottleneck. Also, the increased rate of clonal reproduction in the non-native range could not have produced the observed patterns by itself but may have magnified the bottleneck effect on LD. Elevated marker diversity in the genetic source populations could have contributed to the increased H_o and Pi observed in the non-native range. We refined the previous invasion source region to a ~50 km section of northeastern Honshu Island. Outlier detection methods failed to reveal any consistently differentiated loci shared among invaded regions, probably because of the complex A. vermiculophyllum demographic history. Our results reinforce the importance of demographic history, specifically founder effects, in driving genomic variation of invasive populations, even when localized adaptive evolution and reproductive system shifts are observed.

Mining Grapevine Downy Mildew Susceptibility Genes: A Resource for Genomics-Based Breeding and Tailored Gene Editing

Several pathogens continuously threaten viticulture worldwide. Until now, the investigation on resistance loci has been the main trend to understand the interaction between grapevine and the mildew causal agents. Dominantly inherited gene-based resistance has shown to be race-specific in some cases, to confer partial immunity, and to be potentially overcome within a few years since its introgression. Recently, on the footprint of research conducted in Arabidopsis, putative genes associated with downy mildew susceptibility have been discovered also in the grapevine genome. In this work, we deep-sequenced four putative susceptibility genes-namely VvDMR6.1, VvDMR6.2, VvDLO1, VvDLO2-in 190 genetically diverse grapevine genotypes to discover new sources of broad-spectrum and recessively inherited resistance. Identified Single Nucleotide Polymorphisms were screened in a bottleneck analysis from the genetic sequence to their impact on protein structure. Fifty-five genotypes showed at least one impacting mutation in one or more of the scouted genes. Haplotypes were inferred for each gene and two of them at the VvDMR6.2 gene were found significantly more represented in downy mildew resistant genotypes. The current results provide a resource for grapevine and plant genetics and could corroborate genomic-assisted breeding programs as well as tailored gene editing approaches for resistance to biotic stresses.

Identification, Association of Natural Variation and Expression Analysis of ZmNAC9 Gene Response to Low Phosphorus in Maize Seedling Stage

Phosphorus (P) is an essential macroelement supporting maize productivity and low-P stress is a limiting factor of maize growth and yield. Improving maize plant tolerance to low P through molecular breeding is an effective alternative to increase crop productivity. In this study, a total of 111 diverse maize inbred lines were used to identify the favorable alleles and nucleotide diversity of candidate ZmNAC9, which plays an important role in response to low P and regulation in root architecture. A significant difference was found under low- and sufficient-P conditions for each of the 22 seedling traits, and a total of 41 polymorphic sites including 32 single nucleotide polymorphisms (SNPs) and 9 insertion and deletions (InDels) were detected in ZmNAC9 among 111 inbred lines. Among the 41 polymorphic studied sites, a total of 39 polymorphic sites were associated with 20 traits except for the dry weight of shoots and forks, of which six sites were highly significantly associated with a diverse number of low-P tolerant root trait index values by using a mixed linear model (MLM) at −log10 P = 3.61. In addition, 29 polymorphic sites under P-sufficient and 32 polymorphic sites under P-deficient conditions were significantly associated with a diverse number of seedling traits, of which five polymorphic sites (position S327, S513, S514, S520, and S827) were strongly significantly associated with multiple seedling traits under low-P and normal-P conditions. Among highly significant sites, most of the sites were associated with root traits under low-P, normal-P, and low-P trait index values. Linkage disequilibrium (LD) was strong at (r2 > 1.0) in 111 inbred lines. Furthermore, the effect of five significant sites was verified for haplotypes in 111 lines and the favorable allele S520 showed a positive effect on the dry weight of roots under the low-P condition. Furthermore, the expression pattern confirmed that ZmNAC9 was highly induced by low P in the roots of the P-tolerant 178 inbred line. Moreover, the subcellular localization of ZmNAC9 encoded by protein was located in the cytoplasm and nucleus. Haplotypes carrying more favorable alleles exhibited superior effects on phenotypic variation and could be helpful in developing molecular markers in maize molecular breeding programs. Taken together, the finding of this study might lead to further functions of ZmNAC9 and genes that might be involved in responses to low-P stress in maize.

Development and validation of high-throughput and low-cost STARP assays for genes underpinning economically important traits in wheat

We developed and validated 56 gene-specific semi-thermal asymmetric reverse PCR (STARP) markers for 46 genes of important wheat quality, biotic and abiotic stress resistance, grain yield, and adaptation-related traits for marker-assisted selection in wheat breeding. Development of high-throughput, low-cost, gene-specific molecular markers is important for marker-assisted selection in wheat breeding. In this study, we developed 56 gene-specific semi-thermal asymmetric reverse PCR (STARP) markers for wheat quality, tolerance to biotic and abiotic stresses, grain yield, and adaptation-related traits. The STARP assays were validated by (1) comparison of the assays with corresponding diagnostic STS/CAPS markers on 40 diverse wheat cultivars and (2) characterization of allelic effects based on the phenotypic and genotypic data of three segregating populations and 305 diverse wheat accessions from China and 13 other countries. The STARP assays showed the advantages of high-throughput, accuracy, flexibility, simple assay design, low operational costs, and platform compatibility. The state-of-the-art assays of this study provide a robust and reliable molecular marker toolkit for wheat breeding programs.

Transcriptional dynamics of Zn-accumulation in developing kernels of maize reveals important Zn-uptake mechanisms

In the present study, transcriptomic analysis of 10-days old baby kernels of two contrasting maize genotypes, namely VQL-2 (high kernel Zn accumulator) and CM-145 (low kernel Zn accumulator), under low- and optimum- soil Zn conditions generated 1948 differentially expressed transcripts. Among these, 666 and 437 transcripts were up-regulated and down-regulated respectively in VQL-2; whereas, 437 and 408 transcripts were up-regulated and down-regulated respectively in CM-145. Remarkably, 135 transcription factors and 77 known Zn transporters expressed differentially. By comparing the transcripts differentially expressed between the optimum-Zn and low-Zn libraries of the contrasting genotypes, we identified 21,986 and 26,871 SNPs, respectively. Similarly, 6810 and 8192 InDels were found between optimum- and low-Zn growing conditions, respectively. Further, 21 differentially expressed genes were co-localized with already known QTLs associated with Zn uptake, such as qZn10, CQZnK9-1 and YNZnK6. These findings will be useful to develop high Zn-accumulator maize through marker-assisted breeding in future.

Association mapping for yield traits in Paeonia rockii based on SSR markers within transcription factors of comparative transcriptome

BACKGROUND

Allelic variation underlying the quantitative traits in plants is caused by the extremely complex regulation process. Tree peony originated in China is a peculiar ornamental, medicinal and oil woody plant. Paeonia rockii, one of tree peony species, is a precious emerging woody oil crop. However, in this valuable plant, the study of functional loci associated with yield traits has rarely been identified. Therefore, to explore the genetic architecture of 24 yield quantitative traits, the association mapping was first reported in 420 unrelated cultivated P. rockii individuals based on the next-generation sequencing (NGS) and single-molecule long-read sequencing (SMLRS).

RESULTS

The developed 58 pairs of polymorphic expressed sequence tag-simple sequence repeat (EST-SSR) markers from 959 candidate transcription factors (TFs) associated with yield were used for genotyping the 420 P. rockii accessions. We observed a high level of genetic diversity (polymorphic information content, PIC = 0.514) and low linkage disequilibrium (LD) between EST-SSRs. Moreover, four subpopulations in the association population were revealed by STRUCTURE analyses. Further, single-marker association analysis identified 141 significant associations, involving 17 quantitative traits and 41 EST-SSRs. These loci were mainly from AP2, TCP, MYB, HSF, bHLH, GATA, and B3 gene families and showed a small proportion of the phenotypic variance (3.79 to 37.45%).

CONCLUSIONS

Our results summarize a valuable collection of functional loci associated with yield traits in P. rockii, and provide a precious resource that reveals allelic variation underlying quantitative traits in Paeonia and other woody oil crops.

Genome-wide association mapping of QTLs implied in potato virus Y population sizes in pepper: evidence for widespread resistance QTL pyramiding

In this study, we looked for genetic factors in the pepper (Capsicum annuum) germplasm that control the number of potato virus Y (PVY) particles entering the plant (i.e. effective population size at inoculation) and the PVY accumulation at the systemic level (i.e. census population size). Using genotyping-by-sequencing (GBS) in a core collection of 256 pepper accessions, we obtained 10 307 single nucleotide polymorphisms (SNPs) covering the whole genome. Genome-wide association studies (GWAS) detected seven SNPs significantly associated with the virus population size at inoculation and/or systemic level on chromosomes 4, 6, 9 and 12. Two SNPs on chromosome 4 associated with both PVY population sizes map closely to the major resistance gene pvr2 encoding the eukaryotic initiation factor 4E. No obvious candidates for resistance were identified in the confidence intervals for the other chromosomes. SNPs detected on chromosomes 6 and 12 colocalized with resistance quantitative trait loci (QTLs) previously identified with a biparental population. These results show the efficiency of GBS and GWAS in C. annuum, indicate highly consistent results between GWAS and classical QTL mapping, and suggest that resistance QTLs identified with a biparental population are representative of a much larger collection of pepper accessions. Moreover, the resistance alleles at these different loci were more frequently combined than expected by chance in the core collection, indicating widespread pyramiding of resistance QTLs and widespread combination of resistance QTLs and major effect genes. Such pyramiding may increase resistance efficiency and/or durability.

Genome-Wide SNP discovery and genomic characterization in avocado (Persea americana Mill.)

Modern crop breeding is based on the use of genetically and phenotypically diverse plant material and, consequently, a proper understanding of population structure and genetic diversity is essential for the effective development of breeding programs. An example is avocado, a woody perennial fruit crop native to Mesoamerica with an increasing popularity worldwide. Despite its commercial success, there are important gaps in the molecular tools available to support on-going avocado breeding programs. In order to fill this gap, in this study, an avocado 'Hass' draft assembly was developed and used as reference to study 71 avocado accessions which represent the three traditionally recognized avocado horticultural races or subspecies (Mexican, Guatemalan and West Indian). An average of 5.72 M reads per individual and a total of 7,108 single nucleotide polymorphism (SNP) markers were produced for the 71 accessions analyzed. These molecular markers were used in a study of genetic diversity and population structure. The results broadly separate the accessions studied according to their botanical race in four main groups: Mexican, Guatemalan, West Indian and an additional group of Guatemalan × Mexican hybrids. The high number of SNP markers developed in this study will be a useful genomic resource for the avocado community.

Transcriptome sequencing, de novo assembly, characterisation of wild accession of blackgram (Vigna mungo var. silvestris) as a rich resource for development of molecular markers and validation of SNPs by high resolution melting (HRM) analysis

BACKGROUND

Blackgram [Vigna mungo (L.) Hepper], is an important legume crop of Asia with limited genomic resources. We report a comprehensive set of genic simple sequence repeat (SSR) and single nucleotide polymorphism (SNPs) markers using Illumina MiSeq sequencing of transcriptome and its application in genetic variation analysis and mapping.

RESULTS

Transcriptome sequencing of immature seeds of wild blackgram, V. mungo var. silvestris by Illumina MiSeq technology generated 1.9 × 10⁷ reads, which were assembled into 40,178 transcripts (TCS) with an average length of 446 bp covering 2.97 GB of the genome. A total of 38,753 CDS (Coding sequences) were predicted from 40,178 TCS and 28,984 CDS were annotated through BLASTX and mapped to GO and KEGG database resulting in 140 unique pathways. The tri-nucleotides were most abundant (39.9%) followed by di-nucleotide (30.2%). About 60.3 and 37.6% of SSR motifs were present in the coding sequences (CDS) and untranslated regions (UTRs) respectively. Among SNPs, the most abundant substitution type were transitions (Ts) (61%) followed by transversions (Tv) type (39%), with a Ts/Tv ratio of 1.58. A total of 2306 DEGs were identified by RNA Seq between wild and cultivar and validation was done by quantitative reverse transcription polymerase chain reaction. In this study, we genotyped SNPs with a validation rate of 78.87% by High Resolution Melting (HRM) Assay.

CONCLUSION

In the present study, 1621genic-SSR and 1844 SNP markers were developed from immature seed transcriptome sequence of blackgram and 31 genic-SSR markers were used to study genetic variations among different blackgram accessions. Above developed markers contribute towards enriching available genomic resources for blackgram and aid in breeding programmes.

Quantitative Trait Nucleotides Impacting the Technological Performances of Industrial Saccharomyces cerevisiae Strains

The budding yeast Saccharomyces cerevisiae is certainly the prime industrial microorganism and is related to many biotechnological applications including food fermentations, biofuel production, green chemistry, and drug production. A noteworthy characteristic of this species is the existence of subgroups well adapted to specific processes with some individuals showing optimal technological traits. In the last 20 years, many studies have established a link between quantitative traits and single-nucleotide polymorphisms found in hundreds of genes. These natural variations constitute a pool of QTNs (quantitative trait nucleotides) that modulate yeast traits of economic interest for industry. By selecting a subset of genes functionally validated, a total of 284 QTNs were inventoried. Their distribution across pan and core genome and their frequency within the 1,011 Saccharomyces cerevisiae genomes were analyzed. We found that 150 of the 284 QTNs have a frequency lower than 5%, meaning that these variants would be undetectable by genome-wide association studies (GWAS). This analysis also suggests that most of the functional variants are private to a subpopulation, possibly due to their adaptive role to specific industrial environment. In this review, we provide a literature survey of their phenotypic impact and discuss the opportunities and the limits of their use for industrial strain selection.

Identification and validation of SNP markers linked to seed toxicity in Jatropha curcas L

Edible/non-toxic varieties of Jatropha curcas L. are gaining increasing attention, providing both oil as biofuel feedstock or even as edible oil and the seed kernel meal as animal feed ingredient. They are a viable alternative to the limitation posed by the presence of phorbol esters in toxic varieties. Accurate genotyping of toxic/non-toxic accessions is critical to breeding management. The aim of this study was to identify SNP markers linked to seed toxicity in J. curcas. For SNP discovery, NGS technology was used to sequence the whole genomes of a toxic and non-toxic parent along with a bulk of 51 toxic and 30 non-toxic F₂ plants. To ascertain the association between SNP markers and seed toxicity trait, candidate SNPs were genotyped on 672 individuals segregating for seed toxicity and two collections of J. curcas composed of 96 individuals each. In silico SNP discovery approaches led to the identification of 64 candidate SNPs discriminating non-toxic and toxic samples. These SNPs were mapped on Chromosome 8 within the Linkage Group 8 previously identified as a genomic region important for phorbol ester biosynthesis. The association study identified two new SNPs, SNP_J22 and SNP_J24 significantly linked to low toxicity with R² values of 0.75 and 0.54, respectively. Our study released two valuable SNP markers for high-throughput, marker-assisted breeding of seed toxicity in J. curcas.

Genome-wide association study (GWAS) of leaf cuticular wax components in Camelina sativa identifies genetic loci related to intracellular wax transport

BACKGROUND

It is important to explore renewable alternatives (e.g. biofuels) that can produce energy sources to help reduce reliance on fossil oils, and reduce greenhouse gases and waste solids resulted from fossil oils consumption. Camelina sativa is an oilseed crop which has received increasing attention due to its short life cycle, broader adaptation regions, high oil content, high level of omega-3 unsaturated fatty acids, and low-input requirements in agriculture practices. To expand its Camelina production areas into arid regions, there is a need to breed for new drought-tolerant cultivars. Leaf cuticular wax is known to facilitate plant development and growth under water-limited conditions. Dissecting the genetic loci underlying leaf cuticular waxes is important to breed for cultivars with improved drought tolerance.

RESULTS

Here we combined phenotypic data and single nucleotide polymorphism (SNP) data from a spring C. sativa diversity panel using genotyping-by-sequencing (GBS) technology, to perform a large-scale genome-wide association study (GWAS) on leaf wax compositions. A total of 42 SNP markers were significantly associated with 15 leaf wax traits including major wax components such as total primary alcohols, total alkanes, and total wax esters as well as their constituents. The vast majority of significant SNPs were associated with long-chain carbon monomers (carbon chain length longer than C28), indicating the important effects of long-chain carbon monomers on leaf total wax biosynthesis. These SNP markers are located on genes directly or indirectly related to wax biosynthesis such as maintaining endoplasmic reticulum (ER) morphology and enabling normal wax secretion from ER to plasma membrane or Golgi network-mediated transport.

CONCLUSIONS

These loci could potentially serve as candidates for the genetic control involved in intracellular wax transport that might directly or indirectly facilitate leaf wax accumulation in C. sativa and can be used in future marker-assisted selection (MAS) to breed for the cultivars with high wax content to improve drought tolerance.

Natural variation of ZmHKT1 affects root morphology in maize at the seedling stage

Eight variants in ZmHKT1 promoter were significantly associated with root diameter, four haplotypes based on these significant variants were found, and Hap2 has the largest root diameter. Roots play an important role in uptake of water, nutrients and plant anchorage. Identification of gene and corresponding SNPs associated with root traits would enable develop maize lines with better root traits that might help to improve capacity for absorbing nutrients and water acquisition. The genomic sequences of a salt tolerance gene ZmHKT1 was resequenced in 349 maize inbred lines, and the association between nucleotide polymorphisms and seedling root traits was detected. A total of 269 variants in ZmHKT1 were identified, including 226 single nucleotide polymorphisms and 43 insertions and deletions. The gene displayed high level of nucleotide diversity, especially in non-genic regions. A total of 19 variations in untranslated region of ZmHKT1 were found to be associated with six seedling traits. Eight variants in promoter region were significantly associated with average root diameter (ARD), four haplotypes were found based on these significant variants, and Hap2 has the largest ARD. Two SNPs in high-linkage disequilibrium (SNP-415 and SNP 2169) with pleiotropic effects were significantly associated with plant height, root surface area, root volume, and shoot dry weight. This result revealed that ZmHKT1 was an important contributor to the phenotypic variations of seedling root traits in maize, these significant variants could use to develop functional markers to improve root traits.

Novel Features and Considerations for ERA and Regulation of Crops Produced by Genome Editing

Genome editing describes a variety of molecular biology applications enabling targeted and precise alterations of the genomes of plants, animals and microorganisms. These rapidly developing techniques are likely to revolutionize the breeding of new crop varieties. Since genome editing can lead to the development of plants that could also have come into existence naturally or by conventional breeding techniques, there are strong arguments that these cases should not be classified as genetically modified organisms (GMOs) and be regulated no differently from conventionally bred crops. If a specific regulation would be regarded necessary, the application of genome editing for crop development may challenge risk assessment and post-market monitoring. In the session “Plant genome editing—any novel features to consider for ERA and regulation?” held at the 14th ISBGMO, scientists from various disciplines as well as regulators, risk assessors and potential users of the new technologies were brought together for a knowledge-based discussion to identify knowledge gaps and analyze scenarios for the introduction of genome-edited crops into the environment. It was aimed to enable an open exchange forum on the regulatory approaches, ethical aspects and decision-making considerations.

Enhancing (crop) plant photosynthesis by introducing novel genetic diversity

Although some elements of the photosynthetic light reactions might appear to be ideal, the overall efficiency of light conversion to biomass has not been optimized during evolution. Because crop plants are depleted of genetic diversity for photosynthesis, efforts to enhance its efficiency with respect to light conversion to yield must generate new variation. In principle, three sources of natural variation are available: (i) rare diversity within extant higher plant species, (ii) photosynthetic variants from algae, and (iii) reconstruction of no longer extant types of plant photosynthesis. Here, we argue for a novel approach that outsources crop photosynthesis to a cyanobacterium that is amenable to adaptive evolution. This system offers numerous advantages, including a short generation time, virtually unlimited population sizes and high mutation rates, together with a versatile toolbox for genetic manipulation. On such a synthetic bacterial platform, 10 000 years of (crop) plant evolution can be recapitulated within weeks. Limitations of this system arise from its unicellular nature, which cannot reproduce all aspects of crop photosynthesis. But successful establishment of such a bacterial host for crop photosynthesis promises not only to enhance the performance of eukaryotic photosynthesis but will also reveal novel facets of the molecular basis of photosynthetic flexibility.This article is part of the themed issue 'Enhancing photosynthesis in crop plants: targets for improvement'.

Genetic relatedness of previously Plant-Variety-Protected commercial maize inbreds

The emergence of high-throughput, high-density genotyping methods combined with increasingly powerful computing systems has created opportunities to further discover and exploit the genes controlling agronomic performance in elite maize breeding populations. Understanding the genetic basis of population structure in an elite set of materials is an essential step in this genetic discovery process. This paper presents a genotype-based population analysis of all maize inbreds whose Plant Variety Protection certificates had expired as of the end of 2013 (283 inbreds) as well as 66 public founder inbreds. The results provide accurate population structure information and allow for important inferences in context of the historical development of North American elite commercial maize germplasm. Genotypic data was obtained via genotyping-by-sequencing on 349 inbreds. After filtering for missing data, 77,314 high-quality markers remained. The remaining missing data (average per individual was 6.22 percent) was fully imputed at an accuracy of 83 percent. Calculation of linkage disequilibrium revealed that the average r2 of 0.20 occurs at approximately 1.1 Kb. Results of population genetics analyses agree with previously published studies that divide North American maize germplasm into three heterotic groups: Stiff Stalk, Non-Stiff Stalk, and Iodent. Principal component analysis shows that population differentiation is indeed very complex and present at many levels, yet confirms that division into three main sub-groups is optimal for population description. Clustering based on Nei's genetic distance provides an additional empirical representation of the three main heterotic groups. Overall fixation index (FST), indicating the degree of genetic divergence between the three main heterotic groups, was 0.1361. Understanding the genetic relationships and population differentiation of elite germplasm may help breeders to maintain and potentially increase the rate of genetic gain, resulting in higher overall agronomic performance.

Cloning, molecular evolution and functional characterization of ZmbHLH16, the maize ortholog of OsTIP2 (OsbHLH142)

The transcription factor ZmbHLH16, the maize ortholog of OsTIP2 (OsbHLH142), was isolated in the present study. Tissue expression analysis showed that ZmbHLH16 is preferentially expressed in male reproductive organs. Subcellular location analysis of ZmbHLH16 via rice protoplast indicated that it is located in the nucleus. Through nucleotide variation analysis, 36 polymorphic sites in ZmbHLH16, including 23 single nucleotide polymorphisms and 13 InDels, were detected among 78 maize inbred lines. Neutrality tests and linkage disequilibrium analysis showed that ZmbHLH16 experienced no significant evolutionary pressure. Yeast one-hybrid experiment showed that the first 80 residues in the N-terminus of ZmbHLH16 had transactivation activity, whereas the full length did not. Genome-wide coexpression analysis showed that 395 genes were coexpressed with ZmbHLH16. Among these genes, the transcription factor ZmbHLH51 had similar expression pattern and identical subcellular localization to those of ZmbHLH16. Subsequently, the interaction between ZmbHLH51 and ZmbHLH16 was verified by yeast two-hybrid experiment. Through yeast two-hybrid analysis of series truncated ZmbHLH16 fragments, we found not only the typical bHLH domain [175-221 amino acids (a.a.)], but also that the 81-160 a.a. and 241-365 a.a. of ZmbHLH16 could interact with ZmbHLH51. All these results lay the foundation for further understanding the functions of ZmbHLH16.

Dominant Allele Phylogeny and Constitutive Subgenome Haplotype Inference in Bananas Using Mitochondrial and Nuclear Markers

Cultivated bananas (Musa spp.) have undergone domestication patterns involving crosses of wild progenitors followed by long periods of clonal propagation. Majority of cultivated bananas are polyploids with different constitutive subgenomes and knowledge on phylogenies to their progenitors at the species and subspecies levels is essential. Here, the mitochondrial (NAD1) and nuclear (CENH3) markers were used to phylogenetically position cultivated banana genotypes to diploid progenitors. The CENH3 nuclear marker was used to identify a minimum representative haplotype number in polyploids and diploid bananas based on single nucleotide polymorphisms. The mitochondrial marker NAD1 was observed to be ideal in differentiating bananas of different genomic constitutions based on size of amplicons as well as sequence. The genotypes phylogenetically segregated based on the dominant genome; AAB genotypes grouped with AA and AAA, and the ABB together with BB. Both markers differentiated banana sections, but could not differentiate subspecies within the A genomic group. On the basis of CENH3 marker, a total of 13 haplotypes (five in both diploid and triploid, three in diploids, and rest unique to triploids) were identified from the genotypes tested. The presence of haplotypes, which were common in diploids and triploids, stipulate possibility of a shared ancestry in the genotypes involved in this study. Furthermore, the presence of multiple haplotypes in some diploid bananas indicates their being heterozygous. The haplotypes identified in this study are of importance because they can be used to check the level of homozygozity in breeding lines as well as to track segregation in progenies.

Report on the development of putative functional SSR and SNP markers in passion fruits

BACKGROUND

Passionflowers Passiflora edulis and Passiflora alata are diploid, outcrossing and understudied fruit bearing species. In Brazil, passion fruit cultivation began relatively recently and has earned the country an outstanding position as the world's top producer of passion fruit. The fruit's main economic value lies in the production of juice, an essential exotic ingredient in juice blends. Currently, crop improvement strategies, including those for underexploited tropical species, tend to incorporate molecular genetic approaches. In this study, we examined a set of P. edulis transcripts expressed in response to infection by Xanthomonas axonopodis, (the passion fruit's main bacterial pathogen that attacks the vines), aiming at the development of putative functional markers, i.e. SSRs (simple sequence repeats) and SNPs (single nucleotide polymorphisms).

RESULTS

A total of 210 microsatellites were found in 998 sequences, and trinucleotide repeats were found to be the most frequent (31.4%). Of the sequences selected for designing primers, 80.9% could be used to develop SSR markers, and 60.6% SNP markers for P. alata. SNPs were all biallelic and found within 15 gene fragments of P. alata. Overall, gene fragments generated 10,003 bp. SNP frequency was estimated as one SNP every 294 bp. Polymorphism rates revealed by SSR and SNP loci were 29.4 and 53.6%, respectively.

CONCLUSIONS

Passiflora edulis transcripts were useful for the development of putative functional markers for P. alata, suggesting a certain level of sequence conservation between these cultivated species. The markers developed herein could be used for genetic mapping purposes and also in diversity studies.

The impacts of genetic polymorphisms in genes of base excision repair pathway on the efficacy and acute toxicities of (chemo)radiotherapy in patients with nasopharyngeal carcinoma

Purpose

To explore whether polymorphisms in base excision repair (BER) pathway genes are predictors of (chemo)radiotherapy outcome in patients with nasopharyngeal carcinoma (NPC).

Methods

We genotyped five potentially functional single nucleotide polymorphisms (SNPs) of three genes in the BER pathway in 174 NPC patients who were treated with (chemo)radiotherapy. Sequenom MassArray was used for SNPs analysis. The efficacy at the end of radiotherapy and at 3 months after radiotherapy was evaluated by Response Evaluation Criteria in Solid Tumors (RECIST). Acute radiation toxicity was scored using Radiation Therapy Oncology Group and the European Organization for Research and Treatment of Cancer (RTOG/EORTC) acute radiation morbidity scoring criteria. Logistic regression was employed to assess the multivariate analyses.

Results

We found that the wide genotype GG of X-ray repair cross-complementing 1 (XRCC1) rs25489 (GG vs GA: OR=3.833, 95%CI=1.512-9.714, P=0.005; GG vs GA+AA: OR=3.610, 95%CI=1.496-8.713, P=0.004) and the wide genotype CC of 8-oxoguanine DNA glycosylase (OGG1) rs1052133 (CC vs GG: OR=0.263, 95%CI=0.073-0.951, P=0.042; CC vs CG+GG: OR=0.454, 95%CI=0.195-1.053, P=0.066) were positively and negatively associated with primary tumor efficacy at the end of radiotherapy, respectively. By contrast, no association was found between BER gene polymorphisms and the treatment outcomes at 3 months post-treatment or the treatment-related acute toxicities.

Conclusions

The SNPs of the BER genes may act as biomarkers for the curative effect of (chemo)radiotherapy. Further study with long-time follow-up and large population is needed for accurate assessment.

Comparison of traditional and new generation DNA markers declares high genetic diversity and differentiated population structure of wild almond species

Wild almond species as sources of genetic variation may have crucial importance in breeding. A total of 389 accessions of 18 species have been analysed using inter-retrotransposon amplified polymorphism (IRAP), retrotransposon-microsatellite amplified polymorphism (REMAP), sequence-specific amplification polymorphism (S-SAP), amplified fragment length polymorphism (AFLP), inter simple sequence repeat (ISSR) and simple sequence repeats (SSR). Retrotransposon markers indicated the presence and movement of some Ty3-gypsy and Ty1-copia-elements in almond genome. Since transposable elements are associated with large-scale genome alterations, REMAP produced more reliable phylogenetic inferences than AFLP where homoplasy may affect clustering. In addition, high resolution melting (HRM) analysis was developed to detect SNPs. HRM analysis revealed 1:189 bp frequency of SNPs in exon positions, and the transition-to-transversion proportion was 1.84:1. The low transition bias suggests low methylation levels in almond genome. The polymorphic information content (PIC) was the highest for SSR markers, while SNPs had an average PIC of 0.59, which is close to the values of the rest of the markers. Huge genetic diversity, fragmented population structure and footprints of human selection was confirmed by merging information from all marker strategies. Considering time, cost and performance HRM can be a marker of choice in future studies of Prunus diversity.

De Novo Assembly, Annotation, and Characterization of Root Transcriptomes of Three Caladium Cultivars with a Focus on Necrotrophic Pathogen Resistance/Defense-Related Genes

Roots are vital to plant survival and crop yield, yet few efforts have been made to characterize the expressed genes in the roots of non-model plants (root transcriptomes). This study was conducted to sequence, assemble, annotate, and characterize the root transcriptomes of three caladium cultivars (Caladium × hortulanum) using RNA-Seq. The caladium cultivars used in this study have different levels of resistance to Pythiummyriotylum, the most damaging necrotrophic pathogen to caladium roots. Forty-six to 61 million clean reads were obtained for each caladium root transcriptome. De novo assembly of the reads resulted in approximately 130,000 unigenes. Based on bioinformatic analysis, 71,825 (52.3%) caladium unigenes were annotated for putative functions, 48,417 (67.4%) and 31,417 (72.7%) were assigned to Gene Ontology (GO) and Clusters of Orthologous Groups (COG), respectively, and 46,406 (64.6%) unigenes were assigned to 128 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. A total of 4518 distinct unigenes were observed only in Pythium-resistant "Candidum" roots, of which 98 seemed to be involved in disease resistance and defense responses. In addition, 28,837 simple sequence repeat sites and 44,628 single nucleotide polymorphism sites were identified among the three caladium cultivars. These root transcriptome data will be valuable for further genetic improvement of caladium and related aroids.

Development of intron polymorphism markers in major latex-like protein gene for locality-level and cultivar identification of Salvia miltiorrhiza

BACKGROUND

Salvia miltiorrhiza (Danshen) is one of the most widely used medicinal herbs in traditional Chinese medicine. Locality-level and cultivar identification is of great importance not only for protecting highest therapeutic effectiveness of Daodi Danshen, but also for the genetic conservation and utilization of existing S. miltiorrhiza populations.

RESULTS

Intron polymorphisms including SNPs (single nucleotide polymorphisms) and indels were exploited in major latex-like protein (MLP) gene. Based on these markers, genetic relationships among S. miltiorrhiza cultivar and populations in different locations were evaluated by constructing a dendrogram. Moreover, S. miltiorrhiza specimens from Laiwu region were geographically distinguishable by the developed SNP marker. A 204 bp-indel marker was exploited for the first space breeding cultivar Luyuan Danshen-1 (LD-1), and an effective real-time PCR assay was successfully developed for fast screening of LD-1 among local landraces.

CONCLUSIONS

MLP intron is a valuable DNA barcode for intra-specific study of S. miltiorrhiza populations, and the developed markers can serve as a useful tool for molecular identification of LD-1 cultivar and geographically distinct populations of S. miltiorrhiza.

Genome-Wide Association Study Dissecting the Genetic Architecture Underlying the Branch Angle Trait in Rapeseed (Brassica napus L.)

The rapeseed branch angle is an important morphological trait because an adequate branch angle enables more efficient light capture under high planting densities. Here, we report that the average angle of the five top branches provides a reliable representation of the average angle of all branches. Statistical analyses revealed a significantly positive correlation between the branch angle and multiple plant-type and yield-related traits. The 60 K Brassica Infinium^® single nucleotide polymorphism (SNP) array was utilized to genotype an association panel with 520 diverse accessions. A genome-wide association study was performed to determine the genetic architecture of branch angle, and 56 loci were identified as being significantly associated with the branch angle trait via three models, including a robust, novel, nonparametric Anderson-Darling (A-D) test. Moreover, these loci explained 51.1% of the phenotypic variation when a simple additive model was applied. Within the linkage disequilibrium (LD) decay ranges of 53 loci, we observed plausible candidates orthologous to documented Arabidopsis genes, such as LAZY1, SGR2, SGR4, SGR8, SGR9, PIN3, PIN7, CRK5, TIR1, and APD7. These results provide insight into the genetic basis of the branch angle trait in rapeseed and might facilitate marker-based breeding for improvements in plant architecture.

Single nucleotide polymorphism in sugar pathway and disease resistance genes in sugarcane

Single nucleotide polymorphism in sugar pathway and disease resistance genes showing genetic association with sugar content and red rot resistance would be useful in marker-assisted genetic improvement of sugarcane. Validation and genotyping of potential sequence variants in candidate genes are necessary to understand their functional significance and trait association potential. We discovered, characterized, validated and genotyped SNPs and InDels in sugar pathway and disease resistance genes of Saccharum complex and sugarcane varieties using amplicon sequencing and CAPS assays. The SNPs were abundant in the non-coding 3'UTRs than 5'UTRs and coding sequences depicting a strong bias toward C to T transition substitutions than transversions. Sequencing of cloned amplicons validated 61.6 and 45.2 % SNPs detected in silico in 21 sugar pathway and 16 disease resistance genes, respectively. Sixteen SNPs in four sugar pathway genes and 10 SNPs in nine disease resistance genes were validated through cost-effective CAPS assay. Functional and adaptive significance of SNP and protein haplotypes identified in sugar pathway and disease resistance genes was assessed by correlating their allelic variation with missense amino acid substitutions in the functional domains, alteration in protein structure models and possible modulation of catalytic enzyme activity in contrasting high and low sugar and moderately red rot resistant and highly susceptible sugarcane genotypes. A strong genetic association of five SNPs in the sugar pathway and disease resistance genes, and an InDel marker in the promoter sequence of sucrose synthase-2 gene, with sugar content and red rot resistance, was evident. The functionally relevant SNPs and InDels, detected and validated in sugar pathway and disease resistance genes, and genic CAPS markers designed, would be of immense use in marker-assisted genetic improvement of sugarcane for sugar content and disease resistance.

Molecular Evolution and Association of Natural Variation in ZmARF31 with Low Phosphorus Tolerance in Maize

Low-phosphorus (P) stress is one of the major factors constraining plant growth and yield. Improving plant tolerance to P starvation through molecular breeding is an efficient alternative to increase grain production. In the study, 331 diverse maize inbreds were used to detect nucleotide diversity and favorable alleles of ZmARF31, which plays a key role in low P responses and root architecture regulation. Significant phenotypic variation was found in each of 11 tested traits under both P and no-P treatments, and 30 single nucleotide polymorphisms (SNPs) and 14 insertion-deletions (InDels) were detected in ZmARF31 among the 331 maize inbreds. The 5'-untranslated region (UTR) of ZmARF31 showed a small linkage disequilibrium (LD) block under significant purifying selection, whereas the 3'-UTR showed the most abundant diversity and a larger LD block. Thirty, fourteen, and nine natural variations were identified in ZmARF31 that were associated with P-deficiency-tolerance traits (P ≤ 0.01) by using the general linear model (GLM), GLM incorporated with population structure, and mixed linear model, respectively. Four SNPs were significantly associated with the total dry weight (TDW) in the three models, of which SNPs S410 and S462 were located in a complete LD block. A further verification conducted in a recombinant inbred line population revealed that favorable allele G/G of non-synonymous mutation S410 and favorable allele with a 38 bp insertion of InDel S1442 exhibited positive genetic effects on the TDW and total root tips, respectively. Expression analysis further confirmed that ZmARF31 was highly expressed in the roots of low-P-tolerant inbred 178. The protein encoded by ZmARF31 was located both in the nucleus and cytoplasm. Haplotypes carrying more favorable alleles showed a greater effect on phenotypic variation than single loci. Such haplotypes should be helpful to develop valuable genetic markers and perform maize molecular breeding.

Exploiting genotyping by sequencing to characterize the genomic structure of the American cranberry through high-density linkage mapping

Background

The application of genotyping by sequencing (GBS) approaches, combined with data imputation methodologies, is narrowing the genetic knowledge gap between major and understudied, minor crops. GBS is an excellent tool to characterize the genomic structure of recently domesticated (~200 years) and understudied species, such as cranberry (Vaccinium macrocarpon Ait.), by generating large numbers of markers for genomic studies such as genetic mapping.

Results

We identified 10842 potentially mappable single nucleotide polymorphisms (SNPs) in a cranberry pseudo-testcross population wherein 5477 SNPs and 211 short sequence repeats (SSRs) were used to construct a high density linkage map in cranberry of which a total of 4849 markers were mapped. Recombination frequency, linkage disequilibrium (LD), and segregation distortion at the genomic level in the parental and integrated linkage maps were characterized for first time in cranberry. SSR markers, used as the backbone in the map, revealed high collinearity with previously published linkage maps. The 4849 point map consisted of twelve linkage groups spanning 1112 cM, which anchored 2381 nuclear scaffolds accounting for ~13 Mb of the estimated 470 Mb cranberry genome. Bin mapping identified 592 and 672 unique bins in the parentals and a total of 1676 unique marker positions in the integrated map. Synteny analyses comparing the order of anchored cranberry scaffolds to their homologous positions in kiwifruit, grape, and coffee genomes provided initial evidence of homology between cranberry and closely related species.

Conclusions

GBS data was used to rapidly saturate the cranberry genome with markers in a pseudo-testcross population. Collinearity between the present saturated genetic map and previous cranberry SSR maps suggests that the SNP locations represent accurate marker order and chromosome structure of the cranberry genome. SNPs greatly improved current marker genome coverage, which allowed for genome-wide structure investigations such as segregation distortion, recombination, linkage disequilibrium, and synteny analyses. In the future, GBS can be used to accelerate cranberry molecular breeding through QTL mapping and genome-wide association studies (GWAS).

Electronic supplementary material

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

A multiparental cross population for mapping QTL for agronomic traits in durum wheat (Triticum turgidum ssp. durum)

Multiparental cross designs for mapping quantitative trait loci (QTL) provide an efficient alternative to biparental populations because of their broader genetic basis and potentially higher mapping resolution. We describe the development and deployment of a recombinant inbred line (RIL) population in durum wheat (Triticum turgidum ssp. durum) obtained by crossing four elite cultivars. A linkage map spanning 2664 cM and including 7594 single nucleotide polymorphisms (SNPs) was produced by genotyping 338 RILs. QTL analysis was carried out by both interval mapping on founder haplotype probabilities and SNP bi-allelic tests for heading date and maturity date, plant height and grain yield from four field experiments. Sixteen QTL were identified across environments and detection methods, including two yield QTL on chromosomes 2BL and 7AS, with the former mapped independently from the photoperiod response gene Ppd-B1, while the latter overlapped with the vernalization locus VRN-A3. Additionally, 21 QTL with environment-specific effects were found. Our results indicated a prevalence of environment-specific QTL with relatively small effect on the control of grain yield. For all traits, functionally different QTL alleles in terms of direction and size of genetic effect were distributed among parents. We showed that QTL results based on founder haplotypes closely matched functional alleles at known heading date loci. Despite the four founders, only 2.1 different functional haplotypes were estimated per QTL, on average. This durum wheat population provides a mapping resource for detailed genetic dissection of agronomic traits in an elite background typical of breeding programmes.

Genetic Analysis in Maize Foundation Parents with Mapping Population and Testcross Population: Ye478 Carried More Favorable Alleles and Using QTL Information Could Improve Foundation Parents

The development of maize foundation parents is an important part of genetics and breeding research, and applying new genetic information to produce foundation parents has been challenging. In this study, we focused on quantitative trait loci (QTLs) and general combining ability (GCA) of Ye478, a widely used foundation parent in China. We developed three sets of populations for QTL mapping and to analyze the GCA for some agronomic traits. The assessment of 15 traits resulted in the detection of 251 QTLs in six tested environments, with 119 QTLs identified through a joint analysis across all environments. Further, analyses revealed that most favorable alleles for plant type-related traits were from Ye478, and more than half of the favorable alleles for yield-related traits were from R08, another foundation parent used in southwestern China, suggesting that different types of foundation parents carried different favorable alleles. We observed that the GCA for most traits (e.g., plant height and 100-kernel weight) was maintained in the inbred lines descended from the foundation parents. Additionally, the continuous improvement in the GCA of the descendants of the foundation parents was consistent with the main trend in maize breeding programs. We identified three significant genomic regions that were highly conserved in three Ye478 descendants, including the stable QTL for plant height. The GCA for the traits in the F₇ generation revealed that the QTLs for the given traits per se were affected by additive effects in the same way in different populations.