Construction of a genetic linkage map for Camellia sinensis (tea)

High-density genetic map construction and QTL mapping of a zigzag-shaped stem trait in tea plant (Camellia sinensis)

The highly unique zigzag-shaped stem phenotype in tea plants boasts significant ornamental value and is exceptionally rare. To investigate the genetic mechanism behind this trait, we developed BC1 artificial hybrid populations. Our genetic analysis revealed the zigzag-shaped trait as a qualitative trait. Utilizing whole-genome resequencing, we constructed a high-density genetic map from the BC1 population, incorporating 5,250 SNP markers across 15 linkage groups, covering 3,328.51 cM with an average marker interval distance of 0.68 cM. A quantitative trait locus (QTL) for the zigzag-shaped trait was identified on chromosome 4, within a 61.2 to 97.2 Mb range, accounting for a phenotypic variation explained (PVE) value of 13.62%. Within this QTL, six candidate genes were pinpointed. To better understand their roles, we analyzed gene expression in various tissues and individuals with erect and zigzag-shaped stems. The results implicated CsXTH (CSS0035625) and CsCIPK14 (CSS0044366) as potential key contributors to the zigzag-shaped stem formation. These discoveries lay a robust foundation for future functional genetic mapping and tea plant genetic enhancement.

Molecular markers in tea plant (Camellia sinensis): Applications to evolution, genetic identification, and molecular breeding

Tea plants have a long cultivation history in the world, and the beverage (tea) made from its leaves is well known in the world. Due to the characteristics of self-incompatibility, long-term natural and artificial hybridization, tea plants have a very complex genetic background, which make the classification of tea plants unclear. Molecular marker, one type of genetic markers, has the advantages of stable inheritance, large amount of information, and high reliability. The development of molecular marker has facilitated the understanding of complex tea germplasm resources. So far, molecular markers had played important roles in the study of the origin and evolution, the preservation and identification of tea germplasms, and the excellent cultivars breeding of tea plants. However, the information is scattered, making it difficult to understand the advance of molecular markers in tea plants. In this paper, we summarized the development process and types of molecular markers in tea plants. In addition, the application advance of these molecular markers in tea plants was reviewed. Perspectives of molecular markers in tea plants were also systematically provided and discussed. The elaboration of molecular markers in this paper should help us to renew understanding of its application in tea plants.

Combating Climate Change in the Kenyan Tea Industry

Climate change triggered by global warming poses a major threat to agricultural systems globally. This phenomenon is characterized by emergence of pests and diseases, extreme weather events, such as prolonged drought, high intensity rains, hailstones and frosts, which are becoming more frequent ultimately impacting negatively to agricultural production including rain-fed tea cultivation. Kenya is predominantly an agricultural based economy, with the tea sector generating about 26% of the total export earnings and about 4% gross domestic product (GDP). In the recent years, however, the country has witnessed unstable trends in tea production associated with climate driven stresses. Toward mitigation and adaptation of climate change, multiple approaches for impact assessment, intensity prediction and adaptation have been advanced in the Kenyan tea sub-sector. Further, pressure on tea breeders to release improved climate-compatible cultivars for the rapidly deteriorating environment has resulted in the adoption of a multi-targeted approach seeking to understand the complex molecular regulatory networks associated with biotic and abiotic stresses adaptation and tolerance in tea. Genetic modeling, a powerful tool that assists in breeding process, has also been adopted for selection of tea cultivars for optimal performance under varying climatic conditions. A range of physiological and biochemical responses known to counteract the effects of environmental stresses in most plants that include lowering the rates of cellular growth and net photosynthesis, stomatal closure, and the accumulation of organic solutes such as sugar alcohols, or osmolytes have been used to support breeding programs through screening of new tea cultivars suitable for changing environment. This review describes simulation models combined with high resolution climate change scenarios required to quantify the relative importance of climate change on tea production. In addition, both biodiversity and ecosystem based approaches are described as a part of an overall adaptation strategy to mitigate adverse effects of climate change on tea in Kenya and gaps highlighted for urgent investigations.

Tea plant genomics: achievements, challenges and perspectives

Tea is among the world's most widely consumed non-alcoholic beverages and possesses enormous economic, health, and cultural values. It is produced from the cured leaves of tea plants, which are important evergreen crops globally cultivated in over 50 countries. Along with recent innovations and advances in biotechnologies, great progress in tea plant genomics and genetics has been achieved, which has facilitated our understanding of the molecular mechanisms of tea quality and the evolution of the tea plant genome. In this review, we briefly summarize the achievements of the past two decades, which primarily include diverse genome and transcriptome sequencing projects, gene discovery and regulation studies, investigation of the epigenetics and noncoding RNAs, origin and domestication, phylogenetics and germplasm utilization of tea plant as well as newly developed tools/platforms. We also present perspectives and possible challenges for future functional genomic studies that will contribute to the acceleration of breeding programs in tea plants.

Unraveling the Roles of Regulatory Genes during Domestication of Cultivated Camellia: Evidence and Insights from Comparative and Evolutionary Genomics

With the increasing power of DNA sequencing, the genomics-based approach is becoming a promising resolution to dissect the molecular mechanism of domestication of complex traits in trees. Genus Camellia possesses rich resources with a substantial value for producing beverage, ornaments, edible oil and more. Currently, a vast number of genetic and genomic research studies in Camellia plants have emerged and provided an unprecedented opportunity to expedite the molecular breeding program. In this paper, we summarize the recent advances of gene expression and genomic resources in Camellia species and focus on identifying genes related to key economic traits such as flower and fruit development and stress tolerances. We investigate the genetic alterations and genomic impacts under different selection programs in closely related species. We discuss future directions of integrating large-scale population and quantitative genetics and multiple omics to identify key candidates to accelerate the breeding process. We propose that future work of exploiting the genomic data can provide insights related to the targets of domestication during breeding and the evolution of natural trait adaptations in genus Camellia.

Biotechnological advances in tea (Camellia sinensis [L.] O. Kuntze): a review

This article presents a comprehensive review on the success and limitations of biotechnological approaches aimed at genetic improvement of tea with a purpose to explore possibilities to address challenging areas. Tea is a woody perennial tree with a life span of more than 100 years. Conventional breeding of tea is slow and limited primarily to selection which leads to narrowing down of its genetic base. Harnessing the benefits of wild relatives has been negligible due to low cross-compatibility, genetic drag and undesirable alleles for low yield. Additionally, being a recalcitrant species, in vitro propagation of tea is constrained too. Nevertheless, maneuvering with tissue/cell culture techniques, a considerable success has been achieved in the area of micropropagation, somatic embryogenesis as well as genetic transformation. Besides, use of molecular markers, "expressomics" (transcriptomics, proteomics, metabolomics), map-based cloning towards construction of physical maps, generation of expressed sequenced tags (ESTs) have facilitated the identification of QTLs and discovery of genes associated with abiotic or biotic stress tolerance and agronomic traits. Furthermore, the complete genome (or at least gene space) sequence of tea is expected to be accessible in the near future which will strengthen combinational approaches for improvement of tea. This review presents a comprehensive account of the success and limitations of the biotechnological tools and techniques hitherto applied to tea and its wild relatives. Expectedly, this will form a basis for making further advances aimed at genetic improvement of tea in particular and of economically important woody perennials in general.

Large-Scale SNP Discovery and Genotyping for Constructing a High-Density Genetic Map of Tea Plant Using Specific-Locus Amplified Fragment Sequencing (SLAF-seq)

Genetic maps are important tools in plant genomics and breeding. The present study reports the large-scale discovery of single nucleotide polymorphisms (SNPs) for genetic map construction in tea plant. We developed a total of 6,042 valid SNP markers using specific-locus amplified fragment sequencing (SLAF-seq), and subsequently mapped them into the previous framework map. The final map contained 6,448 molecular markers, distributing on fifteen linkage groups corresponding to the number of tea plant chromosomes. The total map length was 3,965 cM, with an average inter-locus distance of 1.0 cM. This map is the first SNP-based reference map of tea plant, as well as the most saturated one developed to date. The SNP markers and map resources generated in this study provide a wealth of genetic information that can serve as a foundation for downstream genetic analyses, such as the fine mapping of quantitative trait loci (QTL), map-based cloning, marker-assisted selection, and anchoring of scaffolds to facilitate the process of whole genome sequencing projects for tea plant.

Construction of a SSR-based genetic map and identification of QTLs for catechins content in tea plant (Camellia sinensis)

Catechins are the most important bioactive compounds in tea, and have been demonstrated to possess a wide variety of pharmacological activities. To characterize quantitative trait loci (QTLs) for catechins content in the tender shoots of tea plant, we constructed a moderately saturated genetic map using 406 simple sequence repeat (SSR) markers, based on a pseudo-testcross population of 183 individuals derived from an intraspecific cross of two Camellia sinensis varieties with diverse catechins composition. The map consisted of fifteen linkage groups (LGs), corresponding to the haploid chromosome number of tea plant (2n = 2x = 30). The total map length was 1,143.5 cM, with an average locus spacing of 2.9 cM. A total of 25 QTLs associated with catechins content were identified over two measurement years. Of these, nine stable QTLs were validated across years, and clustered into four main chromosome regions on LG03, LG11, LG12 and LG15. The population variability explained by each QTL was predominantly at moderate-to-high levels and ranged from 2.4% to 71.0%, with an average of 17.7%. The total number of QTL for each trait varied from four to eight, while the total population variability explained by all QTLs for a trait ranged between 38.4% and 79.7%. This is the first report on the identification of QTL for catechins content in tea plant. The results of this study provide a foundation for further cloning and functional characterization of catechin QTLs for utilization in improvement of tea plant.

Floral transcriptome sequencing for SSR marker development and linkage map construction in the tea plant (Camellia sinensis)

Despite the worldwide consumption and high economic importance of tea, the plant (Camellia sinensis) is not well studied in molecular biology. Under the few circumstances in which the plant is studied, C. sinensis flowers, which are important for reproduction and cross-breeding, receive less emphasis than investigation of its leaves or roots. Using high-throughput Illumina RNA sequencing, we analyzed a C. sinensis floral transcriptome, and 26.9 million clean reads were assembled into 75,531 unigenes averaging 402 bp. Among them, 50,792 (67.2%) unigenes were annotated with a BLAST search against the NCBI Non-Redundant (NR) database and 10,290 (16.67%) were detected that contained one or more simple sequence repeats (SSRs). From these SSR-containing sequences, 2,439 candidate SSR markers were developed and 720 were experimentally tested, validating 431 (59.9%) novel polymorphic SSR markers for C. sinensis. Then, a consensus SSR-based linkage map was constructed that covered 1,156.9 cM with 237 SSR markers distributed in 15 linkage groups. Both transcriptome information and the genetic map of C. sinensis presented here offer a valuable foundation for molecular biology investigations such as functional gene isolation, quantitative trait loci mapping, and marker-assisted selection breeding in this important species.

Genetic map of the primocane-fruiting and thornless traits of tetraploid blackberry

Blackberry primocane fruiting, fruiting on first-year canes, has the potential to expand blackberry production both seasonally and geographically. The incorporation of the primocane-fruiting trait into cultivars with desirable horticultural attributes is challenging due to its recessive nature and tetrasomic inheritance. Molecular marker-assisted selection has high potential to facilitate incorporation, because breeders already use morphological marker-assisted selection of seedlings without marginal cotyledonary hairs to identify progeny that will be thornless when mature. The development of a genetic linkage map with these two traits is the first step to utilizing molecular markers in breeding for thornless primocane-fruiting blackberry cultivars. A full-sib family segregating for thornlessness and primocane fruiting, from a cross between 'APF-12' and 'Arapaho', was used to construct the first genetic map of tetraploid blackberry. Segregation patterns of several dominant markers and the two phenotypic traits fit those expected uniquely with tetrasomic inheritance (e.g., 5:1, 11:1 and 35:1). Some loci showed significant double reduction frequencies, but genotypes that could have originated only from double reduction were not found. The map consists of seven linkage groups (LG) in each parent, consistent with the basic number of chromosomes (2n = 4x = 28). Naming of LG1-LG6 followed that of the recently revised system for raspberry using SSR markers in common between blackberry and raspberry, and LG7 was tentatively defined by default. The loci controlling primocane fruiting and thornlessness were not linked to each other; thornless/thorny, the S Locus, was mapped on LG4, and the primocane-/floricane-fruiting locus, named in this work the F Locus, on LG7.

Construction of a high-density reference linkage map of tea (Camellia sinensis)

A few linkage maps of tea have been constructed using pseudo-testcross theory based on dominant marker systems. However, dominant markers are not suitable as landmark markers across a wide range of materials. Therefore, we developed co-dominant SSR markers from genomic DNA and ESTs and constructed a reference map using these co-dominant markers as landmarks. A population of 54 F(1) clones derived from reciprocal crosses between 'Sayamakaori' and 'Kana-Ck17' was used for the linkage analysis. Maps of both parents were constructed from the F(1) population that was taken for BC(1) population. The order of most of the dominant markers in the parental maps was consistent. We constructed a core map by merging the linkage data for markers that detected polymorphisms in both parents. The core map contains 15 linkage groups, which corresponds to the basic chromosome number of tea. The total length of the core map is 1218 cM. Here, we present the reference map as a central core map sandwiched between the parental maps for each linkage group; the combined maps contain 441 SSRs, 7 CAPS, 2 STS and 674 RAPDs. This newly constructed linkage map can be used as a basic reference linkage map of tea.

An EST-derived SNP and SSR genetic linkage map of cassava (Manihot esculenta Crantz)

Cassava (Manihot esculenta Crantz) is one of the most important food security crops in the tropics and increasingly being adopted for agro-industrial processing. Genetic improvement of cassava can be enhanced through marker-assisted breeding. For this, appropriate genomic tools are required to dissect the genetic architecture of economically important traits. Here, a genome-wide SNP-based genetic map of cassava anchored in SSRs is presented. An outbreeder full-sib (F1) family was genotyped on two independent SNP assay platforms: an array of 1,536 SNPs on Illumina's GoldenGate platform was used to genotype a first batch of 60 F1. Of the 1,358 successfully converted SNPs, 600 which were polymorphic in at least one of the parents and was subsequently converted to KBiosciences' KASPar assay platform for genotyping 70 additional F1. High-precision genotyping of 163 informative SSRs using capillary electrophoresis was also carried out. Linkage analysis resulted in a final linkage map of 1,837 centi-Morgans (cM) containing 568 markers (434 SNPs and 134 SSRs) distributed across 19 linkage groups. The average distance between adjacent markers was 3.4 cM. About 94.2% of the mapped SNPs and SSRs have also been localized on scaffolds of version 4.1 assembly of the cassava draft genome sequence. This more saturated genetic linkage map of cassava that combines SSR and SNP markers should find several applications in the improvement of cassava including aligning scaffolds of the cassava genome sequence, genetic analyses of important agro-morphological traits, studying the linkage disequilibrium landscape and comparative genomics.

SSR and EST-SSR-based genetic linkage map of cassava (Manihot esculenta Crantz)

Simple sequence repeat (SSR) markers provide a powerful tool for genetic linkage map construction that can be applied for identification of quantitative trait loci (QTL). In this study, a total of 640 new SSR markers were developed from an enriched genomic DNA library of the cassava variety 'Huay Bong 60' and 1,500 novel expressed sequence tag-simple sequence repeat (EST-SSR) loci were developed from the Genbank database. To construct a genetic linkage map of cassava, a 100 F(1) line mapping population was developed from the cross Huay Bong 60 by 'Hanatee'. Polymorphism screening between the parental lines revealed that 199 SSRs and 168 EST-SSRs were identified as novel polymorphic markers. Combining with previously developed SSRs, we report a linkage map consisted of 510 markers encompassing 1,420.3 cM, distributed on 23 linkage groups with a mean distance between markers of 4.54 cM. Comparison analysis of the SSR order on the cassava linkage map and the cassava genome sequences allowed us to locate 284 scaffolds on the genetic map. Although the number of linkage groups reported here revealed that this F(1) genetic linkage map is not yet a saturated map, it encompassed around 88% of the cassava genome indicating that the map was almost complete. Therefore, sufficient markers now exist to encompass most of the genomes and efficiently map traits in cassava.

Construction, characterization, and preliminary BAC-end sequence analysis of a bacterial artificial chromosome library of the tea plant (Camellia sinensis)

We describe the construction and characterization of a publicly available BAC library for the tea plant, Camellia sinensis. Using modified methods, the library was constructed with the aim of developing public molecular resources to advance tea plant genomics research. The library consists of a total of 401,280 clones with an average insert size of 135 kb, providing an approximate coverage of 13.5 haploid genome equivalents. No empty vector clones were observed in a random sampling of 576 BAC clones. Further analysis of 182 BAC-end sequences from randomly selected clones revealed a GC content of 40.35% and low chloroplast and mitochondrial contamination. Repetitive sequence analyses indicated that LTR retrotransposons were the most predominant sequence class (86.93%-87.24%), followed by DNA retrotransposons (11.16%-11.69%). Additionally, we found 25 simple sequence repeats (SSRs) that could potentially be used as genetic markers.

A molecular linkage map of olive (Olea europaea L) based on RAPD, microsatellite, and SCAR markers

An integrated molecular linkage map of olive (Olea europaea L.) was constructed based on randomly amplified polymorphic DNA (RAPD), sequence characterized amplified region (SCAR), and microsatellite markers using the pseudo-testcross strategy. A mapping population of 104 individuals was generated from an F1 full-sib family of a cross between 'Frantoio' and 'Kalamata'. The hybridity of the mapping population was confirmed by genetic similarity and nonmetric multidimensional scaling. Twenty-three linkage groups were mapped for 'Kalamata', covering 759 cM of the genome with 89 loci and an average distance between loci of 11.5 cM. Twenty-seven linkage groups were mapped for 'Frantoio', covering 798 cM of the genome with 92 loci and an average distance between loci of 12.3 cM. For the integrated map, 15 linkage groups covered 879 cM of the genome with 101 loci and an average distance between loci of 10.2 cM. The size of the genomic DNA was estimated to be around 3000 cM. A sequence characterized amplified region marker linked to olive peacock disease resistance was mapped to linkage group 2 of the integrated map. These maps will be the starting point for studies on the structure, evolution, and function of the olive genome. When the mapping progeny pass through their juvenile phase and assume their adult characters, mapping morphological markers and identification of quantitative trait loci for adaptive traits will be the primary targets.

Construction of two genetic linkage maps in cultivated tetraploid alfalfa (Medicago sativa) using microsatellite and AFLP markers

BACKGROUND

Alfalfa (Medicago sativa) is a major forage crop. The genetic progress is slow in this legume species because of its autotetraploidy and allogamy. The genetic structure of this species makes the construction of genetic maps difficult. To reach this objective, and to be able to detect QTLs in segregating populations, we used the available codominant microsatellite markers (SSRs), most of them identified in the model legume Medicago truncatula from EST database. A genetic map was constructed with AFLP and SSR markers using specific mapping procedures for autotetraploids. The tetrasomic inheritance was analysed in an alfalfa mapping population.

RESULTS

We have demonstrated that 80% of primer pairs defined on each side of SSR motifs in M. truncatula EST database amplify with the alfalfa DNA. Using a F1 mapping population of 168 individuals produced from the cross of 2 heterozygous parental plants from Magali and Mercedes cultivars, we obtained 599 AFLP markers and 107 SSR loci. All but 3 SSR loci showed a clear tetrasomic inheritance. For most of the SSR loci, the double-reduction was not significant. For the other loci no specific genotypes were produced, so the significant double-reduction could arise from segregation distortion. For each parent, the genetic map contained 8 groups of four homologous chromosomes. The lengths of the maps were 2649 and 3045 cM, with an average distance of 7.6 and 9.0 cM between markers, for Magali and Mercedes parents, respectively. Using only the SSR markers, we built a composite map covering 709 cM.

CONCLUSIONS

Compared to diploid alfalfa genetic maps, our maps cover about 88-100% of the genome and are close to saturation. The inheritance of the codominant markers (SSR) and the pattern of linkage repulsions between markers within each homology group are consistent with the hypothesis of a tetrasomic meiosis in alfalfa. Except for 2 out of 107 SSR markers, we found a similar order of markers on the chromosomes between the tetraploid alfalfa and M. truncatula genomes indicating a high level of colinearity between these two species. These maps will be a valuable tool for alfalfa breeding and are being used to locate QTLs.

Heterologous nuclear and chloroplast microsatellite amplification and variation in tea, Camellia sinensis

The advantage of the cross transferability of heterologous chloroplast and nuclear microsatellite primers was taken to detect polymorphism among 24 tea (Camellia sinensis (L.) O. Kuntze) genotypes, including both the assamica and the sinensis varieties. Primer information was obtained from the closely related Camellia japonica species for four nuclear microsatellites, and from Nicotiana tabaccum for seven universal chloroplast microsatellites. All of the nuclear microsatellite loci tested generated an expected DNA fragment in tea, revealing between three and five alleles per locus. Four out of the seven chloroplast microsatellites primers amplified positively, and of these only one was polymorphic with three alleles, which is in agreement with the conserved nature of chloroplast microsatellites at the intraspecific level. A factorial correspondence analysis carried out on all genotypes and nuclear microsatellite alleles separated the assamica and sinensis genotypes into two groups, thus demonstrating the value of these markers in establishing the genetic relationship between tea varieties. Genetic diversity measured with nuclear microsatellites was higher than that measured with other types of molecular markers, offering prospects for their use in fingerprinting, mapping, and population genetic studies, whereas polymorphisms detected at a cpSSR locus will allow the determination of plastid inheritance in the species.

Excess heterozygosity contributes to genetic map expansion in pea recombinant inbred populations

Several plant genetic maps presented in the literature are longer than expected from cytogenetic data. Here we compare F(2) and RI maps derived from a cross between the same two parental lines and show that excess heterozygosity contributes to map inflation. These maps have been constructed using a common set of dominant markers. Although not generally regarded as informative for F(2) mapping, these allowed rapid map construction, and the resulting data analysis has provided information not otherwise obvious when examining a population from only one generation. Segregation distortion, a common feature of most populations and marker systems, found in the F(2) but not the RI, has identified excess heterozygosity. A few markers with a deficiency of heterozygotes were found to map to linkage group V (chromosome 3), which is known to form rod bivalents in this cross. Although the final map length was longer for the F(2) population, the mapped order of markers was generally the same in the F(2) and RI maps. The data presented in this analysis reconcile much of the inconsistency between map length estimates from chiasma counts and genetic data.

Inheritance of random amplified polymorphic DNA (RAPD-PCR) markers and their use in population genetic studies of Calomys musculinus (Rodentia, Muridae), the reservoir of Argentine hemorrhagic fever

In order to assess the reliability of RAPD markers in the estimation of the genetic structure of natural populations of the murid rodent Calomys musculinus (reservoir of Junin virus, ethiological agent of Argentine Hemorrhagic Fever), we have analyzed the heritability of RAPD bands in 10 parents and their offspring (33 individuals). Fourteen out of a total of 119 bands obtained were absent in the parental patterns, but consistently amplified in offspring from some families. These bands can be eliminated from analyses. Overall degree of band sharing between individuals, including non-parental bands, correctly grouped members of a family in the same cluster in a UPGMA tree, with a high bootstrap percentage. Results support the usefulness of RAPDs as hereditable markers. One hundred polymorphic RAPD loci were identified in three natural populations of C. musculinus. Mean expected heterozygosity in three natural populations ranged from 0.206 to 0.220. Allele frequency based and phenotype based measures of genetic differentiation among natural populations of C. musculinus gave similar results (Weir and Cockerham's theta = 0.133; Excoffier et al.'s phi = 0.127). These values were considerably higher than those found previously using allozymes as genetic markers, and are compatible with moderate to low levels of gene flow among populations.