Distribution of <Emphasis Type="Italic">Divo</Emphasis> in <Emphasis Type="Italic">Coffea</Emphasis> genomes,a poorly described family of angiosperm LTR-Retrotransposons

Coffea arabica (the Arabica coffee) is an allotetraploid species originating from a recent hybridization between two diploid species: C. canephora and C. eugenioides. Transposable elements can drive structural and functional variation during the process of hybridization and allopolyploid formation in plants. To learn more about the evolution of the C. arabica genome, we characterized and studied a new Copia LTR-Retrotransposon (LTR-RT) family in diploid and allotetraploid Coffea genomes called Divo. It is a complete and relatively compact LTR-RT element (~5 kb), carrying typical Gag and Pol Copia type domains. Reverse Trancriptase (RT) domain-based phylogeny demonstrated that Divo is a new and well-supported family in the Bianca lineage, but strictly restricted to dicotyledonous species. In C. canephora, Divo is expressed and showed a genomic distribution along gene rich and gene poor regions. The copy number, the molecular estimation of insertion time and the analysis at orthologous locations of insertions in diploid and allotetraploid coffee genomes suggest that Divo underwent a different and recent transposition activity in C. arabica and C. canephora when compared to C. eugenioides. The analysis of this novel LTR-RT family represents an important step toward uncovering the genome structure and evolution of C. arabica allotetraploid genome.     

A panel of the most suitable reference genes for RT-qPCR expression studies of coffee: screening their stability under different conditions

The reliability of analyses using real-time quantitative polymerase chain reaction (RT-qPCR) depends on the selection of appropriate reference genes to correct for sample-to-sample and run-to-run variations. The aim of the present study was to select the most suitable reference genes for gene expression analyses in tissue samples from coffee, Coffea arabica L. (Arabica) grown under well-watered (WW) and water-deficit (WD) conditions and C. canephora Pierre ex A. Froehner (Robusta) grown under WW conditions. Expression profiles and stabilities were evaluated for 12 reference genes in different tissues from C. arabica and for 8 genes in tissues from C. canephora. The web-based RefFinder tool, which combines the geNorm, NormFinder, Bestkeeper, and Delta-Ct algorithms, was employed to assess the stability of the tested genes. The most stable reference genes identified for all tissues grouped (WW/WD) of C. arabica were clathrin adaptor protein medium subunit (AP47), ubiquitin (UBQ), 60S ribosomal protein L39 (RPL39), and elongation factor 1α (EF1α), while class III alcohol dehydrogenase (ADH2), β-actin (ACT), glyceraldehyde 3-phosphate dehydrogenase (GAPDH), and ubiquitin (UBQ) genes were the most stable for all tissues grouped (WW) of C. canephora tissues. Validation by the expression level analysis of CaACO-like demonstrated that the use of the best and the worst set of reference genes produced different expression results. The results reinforce the general assumption that there is no universal reference gene and that it is essential to select the most appropriate gene for each individual experiment to apply adequate normalization procedures of RT-qPCR data.     

Nucleotide Diversity of the Coding and Promoter Regions of <Emphasis Type="Italic">DREB1D</Emphasis>, a Candidate Gene for Drought Tolerance in <Emphasis Type="Italic">Coffea</Emphasis> Species

Climate change is posing a major challenge to coffee production worldwide leading to a need for the development of coffee cultivars with increased drought tolerance. In several plant species, the use of DREB genes in crop improvement has achieved promising results to desiccation tolerance engineering. Recent studies reported CcDREB1D specific patterns of expression in Coffea canephora and functional evidence of this gene involvement in drought stress responses. However, knowledge on natural diversity of this gene is largely unknown. In this context, this study aimed at evaluating the sequence variability of the DREB1D gene in several Coffea genotypes. Nucleotide variation in promoters and coding regions of this gene were evaluated in a population consisting of 38 genotypes of C. canephora, C. arabica and C. eugenioides, most of them characterized by different phenotypes (tolerance vs. susceptibility) in relation to drought. The genetic diversity of the loci revealed different haplotypes for the promoter and coding regions. In particular, our findings suggest association between drought tolerance and the genetic variations on DREB1D promoter regions, but not with those from its corresponding coding regions. Gene expression studies revealed up-regulated expression of DREB1D gene upon drought mainly in leaves of drought-tolerant clones of C. canephora, and in response to drought, high, and low temperatures in leaves of C. arabica, suggesting a key role of this gene in coffee responses to abiotic stress.     

Population structure and genetic diversity of coffee progenies derived from Catuaí and Híbrido de Timor revealed by genome-wide SNP marker

The use of single nucleotide polymorphism (SNP) molecular markers has provided advances in selection methodologies used in breeding programs of different crops, reducing cost and time of cultivar release. Despite the great economic and social importance of Coffea arabica, studies with SNP markers are scarce and a small number of SNP are available for this species, when compared with other crops of agronomic importance. Thus, the objective of this study was to identify and validate SNP molecular markers for the species Coffea arabica and to introduce these markers to genetic breeding by means of an accurate analysis of the diversity and genetic structure of breeding populations of this species. After quality filtering, 11,187 SNP markers were selected from the coffee population obtained from crosses between the genotypes Catuaí and Híbrido de Timor. A great number of markers were distributed in the 11 chromosomes, within transcribed regions, and were used to estimate the genetic dissimilarity among the individuals of the breeding population. Dendrogram analysis and a Bayesian approach demonstrated the formation of two groups and the discrimination of all genotypes evaluated. The expressive number of SNP molecular markers distributed throughout C. arabica genome was efficient to discriminate all the accessions evaluated in the experiment, clustering them according to their genealogies. This work identified mixtures within the progenies. The genotyping data also provided detailed information about the parental genotypes and led to the identification of new candidate parents to be introduced to the breeding program. The study discussed population structure and its consequence in obtaining improved varieties of C. arabica.     

Development and evaluation of a genome‐wide Coffee 8.5K SNP array and its application for high‐density genetic mapping and for investigating the origin of Coffea arabica L.

Coffee species such as Coffea canephora P. (Robusta) and C. arabica L. (Arabica) are important cash crops in tropical regions around the world. C. arabica is an allotetraploid (2n = 4x = 44) originating from a hybridization event of the two diploid species C. canephora and C. eugenioides (2n = 2x = 22). Interestingly, these progenitor species harbour a greater level of genetic variability and are an important source of genes to broaden the narrow Arabica genetic base. Here, we describe the development, evaluation and use of a single‐nucleotide polymorphism (SNP) array for coffee trees. A total of 8580 unique and informative SNPs were selected from C. canephora and C. arabica sequencing data, with 40% of the SNP located in annotated genes. In particular, this array contains 227 markers associated to 149 genes and traits of agronomic importance. Among these, 7065 SNPs (~82.3%) were scorable and evenly distributed over the genome with a mean distance of 54.4 Kb between markers. With this array, we improved the Robusta high‐density genetic map by adding 1307 SNP markers, whereas 945 SNPs were found segregating in the Arabica mapping progeny. A panel of C. canephora accessions was successfully discriminated and over 70% of the SNP markers were transferable across the three species. Furthermore, the canephora‐derived subgenome of C. arabica was shown to be more closely related to C. canephora accessions from northern Uganda than to other current populations. These validated SNP markers and high‐density genetic maps will be useful to molecular genetics and for innovative approaches in coffee breeding.     

Marker-assisted selection provides arabica coffee with genes from other <Emphasis Type="Italic">Coffea</Emphasis> species targeting on multiple resistance to rust and coffee berry disease

Selecting superior genotypes is facilitated by marker-assisted selection (MAS), which is particularly suitable for transferring disease resistance alleles because it nullifies environmental effects and allows selection of resistant individuals in the absence of the pathogen or race, enabling preventive breeding. Molecular markers linked to two major genes (S_H3 and S_H?), conferring resistance to coffee rust, and those linked to the Ck-1 gene, conferring resistance to coffee berry disease (CBD), have previously been identified. These markers were validated and used in a progeny of crosses between Indian selections with Coffea arabica cultivars. Eleven resistant individuals homozygous for S_H3 were identified by MAS. Of these, seven carry S_H? from Híbrido de Timor and the gene introduced from Coffea liberica (S_H3). S_H? was characterized as derived from Coffea canephora. Thus, it was possible to identify C. arabica genotypes carrying important genes for rust resistance introgressed from other coffee species. MAS also allowed identification of sources of CBD resistance for use in preventive breeding for resistance to this serious disease. Using two validated molecular markers, two coffee plants carrying Ck-1 were identified: the UFV 328-60 genotype (F₂) was resistant and homozygous based on both molecular markers but exhibited no markers related to S_H3 and S_H?, and the UFV 317-12 genotype (F₁) was resistant and homozygous but resistant and heterozygous based on CBD-Sat207 and CBD-Sat235, respectively. Along with possessing Ck-1, the latter carries S_H?. Overall, plants carrying different genes for resistance to rust and CBD were identified. These plants are important sources for gene pyramiding in breeding programs aimed at multiple and durable resistance.     

Specific features in using SNP markers developed for allopolyploid wheat

In this work, we analyzed 54 domestic cultivars of hexaploid (common) wheat Triticum aestivum L. (AABBDD genome) and accessions of tetraploid wheats of the Timopheevi group (AAGG) and rye Secale cereale (RR) using 21 SNP markers for common wheat. It was demonstrated that application of the SNP markers developed and verified for particular common wheat cultivars in allele-specific PCR analysis of other cultivars with different geographic origins could lead to an incorrect estimation of the similarity between the genotypes tested. The studied SNP markers of common wheat are inappropriate for analyzing genomes of other cereal species, in particular, T. timopheevii wheats and rye S. cereale.     

Coffee volatiles induced after mechanical injury and beetle herbivory attract the coffee berry borer and two of its parasitoids

The coffee berry borer (CBB), Hypothenemus hampei (Ferrari), is the most important insect pest of coffee worldwide. In this study, we used headspace solid-phase microextraction coupled with gas chromatography–mass spectrometry to sample and identify volatile compounds from Robusta coffee berries, Coffea canephora Pierre ex Froehner, infested with CBB and with mechanical damage. Furthermore, we evaluated the behavioral responses of the CBB and two of its parasitoids, Prorops nasuta Waterstone and Phymastichus coffea LaSalle, to three selected coffee volatile compounds in a Y-tube olfactometer. We found in the effluvia of red coffee berry compounds not previously reported for this coffee species. Our results show that Robusta coffee berries release induced volatiles either by insect herbivory or by mechanical damage. Small amount of butyl acetate, unknown compound 2, α-longipinene, longiborneol and longiborneol acetate are produced only in infested coffee berries fruits. Quantitatively, nine compounds account for the difference between healthy berries, infested, or mechanically damaged berries. Trans-ocimene, 4,8-dimethyl-3,7-nonadien-2-ol, α-copaene and kaurene increased amount levels in infested berries, while amount of methyl salicylate and linalool increased in mechanically damaged coffee berries. The olfactometric bioassays showed that CBB females and its two parasitoids were attracted to methyl salicylate. In addition, H. hampei and P. nasuta were attracted to linalool, and P. nasuta and P. coffea were attracted to trans-ocimene.     

Survey of Philippine coffee beans for the presence of ochratoxigenic fungi

In 2012 to 2014, Philippine green coffee beans from Coffea arabica in Benguet and Ifugao; Coffea canephora var. Robusta in Abra, Cavite, and Ifugao; and Coffea liberica and Coffea excelsea from Cavite were collected and assessed for the distribution of fungi with the potential to produce ochratoxin A (OTA). The presence of fungal species was evaluated both before and after surface sterilization. There were remarkable ecological and varietal differences in the population of OTA-producing species from the five provinces. Aspergillus ochraceus, A. westerdijkiae, and Penicillium verruculosum were detected from Arabica in Benguet and Ifugao while Aspergillus carbonarius, Aspergillus niger, and Aspergillus japonicus were isolated in Excelsa, Liberica, and Robusta varieties from Abra, Cavite, and Davao. Contamination by Aspergillus and Penicillium species was found on 59 and 19 %, respectively, of the 57 samples from five provinces. After disinfection with 1 % sodium hypochlorite, the levels of infection by Aspergillus and Penicillium fell to 40 and 17 %, respectively. A total of 1184 fungal isolates were identified to species level comprising Aspergillus sections Circumdati (four species), Clavati (one), Flavi (one), Fumigati (one), Nigri (three), and Terrie (one). Within section Circumdati, 70 % of A. ochraceus produced OTA as high as 16238 ng g^?1 while 40 % of A. westerdijkiae produced maximum OTA of 36561 ng g^?1 in solid agar. Within section Nigri, 16.76 % of A. niger produced OTA at the highest 18439 ng g^?1, 10 % of A. japonicus at maximum level of 174 ng g^?1, and 21.21 % of A. carbonarius yielded maximum OTA of 1900 ng g^?1. Of the 12 species of Penicillium isolated, P. verruculosum was ochratoxigenic, with a maximum OTA production of 12 ng g^?1.     

Proteinase inhibitors from <Emphasis Type="Italic">Cajanus platycarpus</Emphasis> accessions active against pod borer <Emphasis Type="Italic">Helicoverpa armigera</Emphasis>

Cajanus platycarpus, a wild relative of Cajanus cajan, is an important source for various agronomically desirable traits, including resistance towards pod borer, Helicoverpa armigera. In the present study, the inhibitory activity of proteinase inhibitors (PIs) present in crude protein extracted from different accessions of C. platycarpus and cultivars of C. cajan was evaluated against H. armigera under in vitro and in vivo conditions. The PIs active against H. armigera gut trypsin-like proteinases (HGPs), referred to as ‘HGPIs’, were more pronounced in mature dry seeds of C. platycarpus accessions when compared with cultivars, which is also evident through gelatin activity staining studies. Therefore, the inhibitory activity of HGPIs was further evaluated in various plant organs of C. platycarpus accessions, such as leaves, flowers, pods, developing seeds at 8–10 days (DAP-I), 18–20 days (DAP-II), and 28–32 days after pollination (DAP-III). However, the HGPI activity was more pronounced in mature dry seeds > DAP-III > DAP-II > DAP-I > flowers > pods > leaves. The observed quantitative allocation of HGPIs closely resembled “Optimal Defense Theory”. Further, bioassays demonstrated that there was a significant reduction in the body weight of the larvae fed upon crude PI extracts of C. platycarpus accessions with concomitant increase in mortality rate and the formation of larval–pupal intermediates. Nevertheless, such changes were not observed when the larvae were fed on crude PI extracts of C. cajan cultivars. These results suggest that the PI gene(s) from C. platycarpus accessions could be exploited in the management of H. armigera by introgression into C. cajan cultivars.     

Genome-wide association study reveals genetic architecture of coleoptile length in wheat

Key message

Eight QTL for coleoptile length were identified in a genome-wide association study on a set of 893 wheat accessions, four of which are novel loci.

Abstract

Wheat cultivars with long coleoptiles are preferred in wheat-growing regions where deep planting is practiced. However, the wide use of gibberellic acid (GA)-insensitive dwarfing genes, Rht-B1b and Rht-D1b, makes it challenging to breed dwarf wheat cultivars with long coleoptiles. To understand the genetic basis of coleoptile length, we performed a genome-wide association study on a set of 893 landraces and historical cultivars using 5011 single nucleotide polymorphism (SNP) markers. Structure analysis revealed four subgroups in the association panel. Association analysis results suggested that Rht-B1b and Rht-D1b genes significantly reduced coleoptile length, and eight additional quantitative trait loci (QTL) for coleoptile length were also identified. These QTL explained 1.45–3.18 and 1.36–3.11% of the phenotypic variation in 2015 and 2016, respectively, and their allelic substitution effects ranged from 0.31 to 1.75 cm in 2015, and 0.63–1.55 cm in 2016. Of the eight QTL, QCL.stars-1BS1, QCL.stars-2DS1, QCL.stars-4BS2, and QCL.stars-5BL1 are likely novel loci for coleoptile length. The favorable alleles in each accession ranged from two to eight with an average of 5.8 at eight loci in the panel, and more favorable alleles were significantly associated with longer coleoptile, suggesting that QTL pyramiding is an effective approach to increase wheat coleoptile length.

     

Combination of newly developed SNP and InDel markers for genotyping the <Emphasis Type="Italic">Cf-9</Emphasis> locus conferring disease resistance to leaf mold disease in the tomato

The Cf-9 gene in the tomato is known to confer resistance against leaf mold disease caused by Cladosporium fulvum, and a gene-based marker targeted to the Cf-9 allele has been widely used as a crop protection approach. However, we found this marker to be misleading in genotyping. Therefore, we developed new single-nucleotide polymorphism (SNP) and insertion and deletion (InDel) markers targeted to the Cf-9 allele in order to increase genotyping accuracy and facilitate high-throughput screening. The DNA sequences of reported Cf-9, cf-9, Cf-0, and closely related Cf-4 alleles were compared, and two functional and non-synonymous SNPs were found to distinguish the Cf-9 resistance allele from the cf-9, Cf-0, and Cf-4 alleles. An SNP marker including these two SNPs was developed and applied to the genotyping of 33 tomato cultivars by high-resolution melting analysis. Our SNP marker was able to select all three Cf-9 genotypes (resistant, heterozygous, and susceptible alleles). Interestingly, two cultivars were grouped separately from these three genotypes. To further examine this outgroup, we preformed polymerase chain reaction (PCR) on two InDel regions identified by sequence comparison of the Cf-9 and Cf-4 genes. The band patterns revealed that these two cultivars carried Cf-4 rather than Cf-9 alleles and that three cultivars classified in the Cf-9 resistance group actually carried both Cf-9 and Cf-4 genes. To determine whether these genotyping results were consistent with disease resistance phenotypes, we examined the induction of a hypersensitive response by transiently expressing the corresponding effector genes, and found that the results matched perfectly with the genotyping results. These findings indicate that the combination of our SNP and InDel markers allows resistant Cf-9 alleles to be distinguished from cf-9 and Cf-4 alleles, which will be useful for marker-assisted selection of tomato cultivars resistant to C. fulvum.     

Study on chloroplast DNA diversity of cultivated and wild pears (<Emphasis Type="Italic">Pyrus</Emphasis> L.) in Northern China

Eight pairs of chloroplast DNA (cpDNA) universal primers selected from 34 pairs were used to assess the genetic diversity of 132 pear accessions in Northern China. Among them, six amplified cpDNA fragments showed genetic diversity. A total of 24 variable sites, including 1 singleton variable site and 23 parsimony informative sites, as well as 21 insertion-deletion fragments, were obtained from the combined cpDNA sequences (5309–5535 bp). Two trnL-trnF-487 haplotypes, five trnL-trnF-413 haplotypes, five rbcL haplotypes, six trnS-psbC haplotypes, eight accD-psaI haplotypes and 12 rps16-trnQ haplotypes were identified among the individuals. Twenty-one haplotypes were identified based on the combined fragments. The values of nucleotide diversity (P_i), average number of nucleotide differences (k) and haplotype diversity (H_d) were 0.00070, 3.56408 and 0.7960, respectively. No statistical significance was detected in Tajima’s D test. Remarkably, the important cpDNA haplotypes and their representing accessions were identified clearly in this study. H_19 was considered as one of the ancient haplotypes and was a divergent centre. H_16 was the most common haplotype of the wild accessions. H_2 was the haplotype representing the most pear germplasm resources (46 cultivars and two wild Ussurian Pear accessions), followed by haplotype H_5 (30 cultivars, two wild Ussurian Pear accessions and four sand pears in outgroups) representing the cultivars ‘Dangshan Suli’ and ‘Yali’, which harbour the largest and the second largest cultivation areas in China. More importantly, this study demonstrated, for the first time, the supposed evolution routes of Pyrus based on cpDNA divergence in the background of pear phylogeny in Northern China.     

Single nucleotide polymorphism tightly linked to a major QTL on chromosome 7A for both kernel length and kernel weight in wheat

Thousand-kernel weight (TKW) is one of the major components of grain yield in wheat (Triticum aestivum). Identifying major quantitative trait loci (QTLs) for TKW and developing effective markers are prerequisite for success in marker-assisted selection (MAS) to improve wheat yield through breeding. This study mapped a major QTL, designated as TaTKW-7AL, for increasing TKW on the long arm of chromosome 7A of ‘Clark’ to a 1.3-cM interval between single nucleotide polymorphism (SNP) markers IWB13913 and IWA5913. This QTL explained 19.7 % of the phenotypic variation for TKW. A QTL for increasing kernel length (KL), one of the major components of TKW, was mapped in the same interval as TaTKW-7AL, suggesting that increased TKW by the QTL in ‘Clark’ is most likely due to the increased KL. Association analysis on a diversity panel of 200 US winter wheat accessions also identified a haplotype of three SNP markers (IWB13913, IWB6693 and IWA5913) that were tightly associated with the both KL and TKW. The analysis of allele frequencies of the haplotype in the diversity panel suggested that the favorable allele of TaTKW-7AL has not been strongly selected for in practice and has potential to be used to improve grain yield in US hard winter wheat breeding. Two user-friendly flanking KASPar markers, IWB13913 and IWA5913, were developed for MAS of TaTKW-7AL.     

Molecular mapping of powdery mildew resistance gene <Emphasis Type="Italic">PmSGD</Emphasis> in Chinese wheat landrace Shangeda using RNA-seq with bulk segregant analysis

Wheat powdery mildew, caused by the fungal pathogen Blumeria graminis f. sp. tritici (Bgt), is one of the most devastating diseases of wheat in China and causes serious yield losses. Resistance genes are urgently needed by wheat breeding programs to combat this disease. In the present study, genetic analysis of powdery mildew resistance was conducted on segregated F₂ and F_2:3 populations derived from the cross of Shangeda (providing good resistance to powdery mildew) and Chancellor (susceptible to powdery mildew). The results showed that the resistance of Shangeda to E09 was controlled by a single recessive gene, tentatively designated as PmSGD. In addition, RNA sequencing of the parental lines Shangeda and Chancellor and the corresponding bulked pools derived from homozygous resistant or susceptible F_2:3 lines was implemented to identify single-nucleotide polymorphisms (SNPs). The PmSGD gene was estimated to be located in the 240–250-Mb region of chromosome 7B based on the characteristics of putative SNP loci distributed on 21 wheat chromosomes. Among the developed SNP markers, 17 (57%) markers were linked to PmSGD flanked by SNP2-57 and SNP2-46, with genetic distances of 0.4 and 0.8 cM, respectively. The reaction patterns of Shangeda and cultivars (lines) carrying the Pm5e, Pmhym, mlxbd, and PmTm4 genes to 22 Bgt isolates indicated that PmSGD may be allelic or very closely linked to those genes. All of the SNP loci linked to PmSGD were used to test 38 cultivars with known Pm gene(s), and the results suggested that these SNP loci are useful for pyramiding PmSGD by marker-assisted selection.