High-resolution organellar genome analysis of Triticum and Aegilops sheds new light on cytoplasm evolution in wheat

Characterizing chloroplast genomes and inferring maternal divergence of the Triticum-Aegilops complex

The Triticum (wheat)-Aegilops (goatgrass) complex has been extensively studied, but the evolutionary history of polyploid wheats has not been fully elucidated. The chloroplast (cp) with maternal inheritance and homoplasy can simplify the sequence-based evolutionary inferences, but informative inferences would require a complete and accurate cp genome sequence. In this study, 16 cp genomes representing five Aegilops and 11 Triticum species and subspecies were sequenced, assembled and annotated, yielding five novel circular cp genome sequences. Analyzing the assembled cp genomes revealed no marked differences in genome structure and gene arrangement across the assayed species. A polymorphism analysis of 72 published cp genome sequences representing 10 Aegilops and 15 Triticum species and subspecies detected 1183 SNPs and 1881 SSRs. More than 80% SNPs detected resided on the downstream and upstream gene regions and only 2.78% or less SNPs were predicted to be deleterious. The largest nucleotide diversity was observed in the short single-copy genomic region. Relatively weak selection pressure on cp coding genes was detected. Different phylogenetic analyses confirmed that the maternal divergence of the Triticum-Aegilops complex had three deep lineages each representing a diploid species with nuclear A, B, or D genome. Dating the maternal divergence yielded age estimates of divergence that matched well with those reported previously. The divergence between emmer and bread wheats occurred at 8200-11,200 years ago. These findings are useful for further genomic studies, provide insight into cp genome evolvability and allow for better understanding of the maternal divergence of the Triticum-Aegilops complex.

Cultivated Tomato (Solanum lycopersicum L.) Suffered a Severe Cytoplasmic Bottleneck during Domestication: Implications from Chloroplast Genomes

In various crops, genetic bottlenecks occurring through domestication can limit crop resilience to biotic and abiotic stresses. In the present study, we investigated nucleotide diversity in tomato chloroplast genome through sequencing seven plastomes of cultivated accessions from the Campania region (Southern Italy) and two wild species among the closest (Solanum pimpinellifolium) and most distantly related (S. neorickii) species to cultivated tomatoes. Comparative analyses among the chloroplast genomes sequenced in this work and those available in GenBank allowed evaluating the variability of plastomes and defining phylogenetic relationships. A dramatic reduction in genetic diversity was detected in cultivated tomatoes, nonetheless, a few de novo mutations, which still differentiated the cultivated tomatoes from the closest wild relative S. pimpinellifolium, were detected and are potentially utilizable as diagnostic markers. Phylogenetic analyses confirmed that S. pimpinellifolium is the closest ancestor of all cultivated tomatoes. Local accessions all clustered together and were strictly related with other cultivated tomatoes (S. lycopersicum group). Noteworthy, S. lycopersicum var. cerasiforme resulted in a mixture of both cultivated and wild tomato genotypes since one of the two analyzed accessions clustered with cultivated tomato, whereas the other with S. pimpinellifolium. Overall, our results revealed a very reduced cytoplasmic variability in cultivated tomatoes and suggest the occurrence of a cytoplasmic bottleneck during their domestication.

The chloroplast view of the evolution of polyploid wheat

Polyploid wheats comprise four species: Triticum turgidum (AABB genomes) and T. aestivum (AABBDD) in the Emmer lineage, and T. timopheevii (AAGG) and T. zhukovskyi (AAGGA(m) A(m) ) in the Timopheevi lineage. Genetic relationships between chloroplast genomes were studied to trace the evolutionary history of the species. Twenty-five chloroplast genomes were sequenced, and 1127 plant accessions were genotyped, representing 13 Triticum and Aegilops species. The A. speltoides (SS genome) diverged before the divergence of T. urartu (AA), A. tauschii (DD) and the Aegilops species of the Sitopsis section. Aegilops speltoides forms a monophyletic clade with the polyploid Emmer and Timopheevi wheats, which originated within the last 0.7 and 0.4 Myr, respectively. The geographic distribution of chloroplast haplotypes of the wild tetraploid wheats and A. speltoides illustrates the possible geographic origin of the Emmer lineage in the southern Levant and the Timopheevi lineage in northern Iraq. Aegilops speltoides is the closest relative of the diploid donor of the chloroplast (cytoplasm), as well as the B and G genomes to Timopheevi and Emmer lineages. Chloroplast haplotypes were often shared by species or subspecies within major lineages and between the lineages, indicating the contribution of introgression to the evolution and domestication of polyploid wheats.

Sequencing of chloroplast genomes from wheat, barley, rye and their relatives provides a detailed insight into the evolution of the Triticeae tribe

Using Roche/454 technology, we sequenced the chloroplast genomes of 12 Triticeae species, including bread wheat, barley and rye, as well as the diploid progenitors and relatives of bread wheat Triticum urartu, Aegilops speltoides and Ae. tauschii. Two wild tetraploid taxa, Ae. cylindrica and Ae. geniculata, were also included. Additionally, we incorporated wild Einkorn wheat Triticum boeoticum and its domesticated form T. monococcum and two Hordeum spontaneum (wild barley) genotypes. Chloroplast genomes were used for overall sequence comparison, phylogenetic analysis and dating of divergence times. We estimate that barley diverged from rye and wheat approximately 8–9 million years ago (MYA). The genome donors of hexaploid wheat diverged between 2.1–2.9 MYA, while rye diverged from Triticum aestivum approximately 3–4 MYA, more recently than previously estimated. Interestingly, the A genome taxa T. boeoticum and T. urartu were estimated to have diverged approximately 570,000 years ago. As these two have a reproductive barrier, the divergence time estimate also provides an upper limit for the time required for the formation of a species boundary between the two. Furthermore, we conclusively show that the chloroplast genome of hexaploid wheat was contributed by the B genome donor and that this unknown species diverged from Ae. speltoides about 980,000 years ago. Additionally, sequence alignments identified a translocation of a chloroplast segment to the nuclear genome which is specific to the rye/wheat lineage. We propose the presented phylogeny and divergence time estimates as a reference framework for future studies on Triticeae.

Genetic structure of Thai rice and rice accessions obtained from the International Rice Research Institute

BACKGROUND

Although the genetic structure of rice germplasm has been characterized worldwide, few studies investigated germplasm from Thailand, the world's largest exporter of rice. Thailand and the International Rice Research Institute (IRRI) have diverse collections of rice germplasm, which could be used to develop breeding lines with desirable traits. This study aimed to investigate the level of genetic diversity and structures of Thai and selected IRRI germplasm. Understanding the genetic structure and relationships among these germplasm will be useful for parent selection used in rice breeding programs.

RESULTS

From the 98 InDel markers tested for single copy and polymorphism, 19 markers were used to evaluate 43 Thai and 57 IRRI germplasm, including improved cultivars, breeding lines, landraces, and 5 other Oryza species. The Thai accessions were selected from all rice ecologies such as irrigated, deep water, upland, and rainfed lowland ecosystems. The IRRI accessions were groups of germplasm having agronomic desirable traits, including temperature-sensitive genetic male sterility (TGMS), new plant type, early flowering, and biotic and abiotic stress resistances. Most of the InDel markers were genes with diverse functions. These markers produced the total of 127 alleles for all loci, with a mean of 6.68 alleles per locus, and a mean Polymorphic Information Content (PIC) of 0.440. Genetic diversity of Thai rice were 0.3665, 0.4479 and 0.3972 for improved cultivars, breeding lines, and landraces, respectively, while genetic diversity of IRRI improved and breeding lines were 0.3272 and 0.2970, respectively. Cluster, structure, and differentiation analyses showed six distinct groups: japonica, TGMS, deep-water, IRRI germplasm, Thai landraces and breeding lines, and other Oryza species.

CONCLUSIONS

Thai and IRRI germplasm were significantly different. Thus, they can be used to broaden the genetic base and trait improvements. Cluster, structure, and differentiation analyses showed concordant results having six distinct groups, in agreement with their development, and ecologies.

Storage products and transcriptional analysis of the endosperm of cultivated wheat and two wild wheat species

The starch and protein in wheat (Triticum aestivum L.) endosperm provide 20 percent of the calories eaten by humans and were heavily selected for during domestication. We examined the main storage products and gene expression patterns that may embody compositional differences between two wild species Aegilops crassa and Aegilops tauschii and cultivated bread wheat. The storage product profiles differed significantly with T. aestivum accumulating twice as much carbon as the wild species, while the latter had 1.5 to 2-fold more total nitrogen per seed. Transcriptional analyses of endosperms of similar fresh weight were compared using a cDNA macroarray. Aegilops tauschii, and especially Ae. crassa had stronger hybridizations with storage protein sequences, but while there were differences in transcripts for starch biosynthetic genes, they were less dramatic. Of these, we cloned the Starch Branching Enzymes (SBE) IIa promoter region and the genomic clone of the Brittle-1 (Bt1) ADPglucose transporter. While Ae. crassa SBEIIa sequence was more divergent than that of Ae. tauschii’s compared to bread wheat, there were no sequence polymorphisms that would explain the observed expression differences in Bt1 between these species. Furthermore, while there were nucleotide differences between Bt1 in Ae. crassa and bread wheat, they were synonymous at the amino acid level. Some of transcriptional differences identified here, however, deserve further examination as part of a strategy to manipulate wheat starch and protein composition.

Characterization of Coffea chloroplast microsatellites and evidence for the recent divergence of C. arabica and C. eugenioides chloroplast genomes

Comparative sequencing of >7 kb of highly variable chloroplast genome regions (atpB–rbcL, trnS–trnG, rpl22–rps19, and rps19–rpl2 spacers; introns in atpF, trnG, trnK, and rpl16) with microsatellites known from other angiosperms was carried out in Coffea . Samples comprised 8 diploid species of Coffea, 5 individuals of tetraploid C. arabica representing geographically distant wild populations from Ethiopia, 2 commercial cultivars of C. arabica, and Psilanthus leroyi and Ixora coccinea as outgroups. Phylogeny reconstruction using maximum parsimony and Bayesian inference resulted in congruent topologies with high support for C. arabica and C. eugenioides being sisters. Partitioned analyses showed that all regions except the atpB–rbcL spacer resolved this sister-group, although this was often unsupported. The large sequence data set further shows that chloroplast genomes of C. arabica and C. eugenioides each possess apomorphies, indicating that not C. eugenioides but an ancestor or close relative of C. eugenioides is the maternal parent of C. arabica. Seven variable chloroplast microsatellites were characterized in Coffea. Most microsatellites are poly(A/T) stretches, whereas one in the trnS–trnG spacer has an (AT)n motif. Most strikingly, all individuals of C. arabica possess identical sequences, suggesting a single chloroplast haplotype. This can be explained by a recent origin of C. arabica in a unique allopolyploidization event, or by severe bottleneck effects in the evolutionary history of the species. Reconstruction of the evolution of microstructural mutations shows much higher levels of homoplasy in microsatellite loci than in other parts of spacers and introns. Microsatellites are inferred to evolve by insertion and deletion of 1 to 3 motif copies in one step.

Independent Wheat B and G Genome Origins in Outcrossing Aegilops Progenitor Haplotypes

The origin of modern wheats involved alloploidization among related genomes. To determine if Aegilops speltoides was the donor of the B and G genomes in AABB and AAGG tetraploids, we used a 3-tiered approach. Using 70 amplified fragment length polymorphism (AFLP) loci, we sampled molecular diversity among 480 wheat lines from their natural habitats encompassing all S genome Aegilops, the putative progenitors of wheat B and G genomes. Fifty-nine Aegilops representatives for S genome diversity were compared at 375 AFLP loci with diploid, tetraploid, and 11 nulli-tetrasomic Triticum aestivum wheat lines. B genome-specific markers allowed pinning the origin of the B genome to S chromosomes of A. speltoides, while excluding other lineages. The outbreeding nature of A. speltoides influences its molecular diversity and bears upon inferences of B and G genome origins. Haplotypes at nuclear and chloroplast loci ACC1, G6PDH, GPT, PGK1, Q, VRN1, and ndhF for approximately 70 Aegilops and Triticum lines (0.73 Mb sequenced) reveal both B and G genomes of polyploid wheats as unique samples of A. speltoides haplotype diversity. These have been sequestered by the AABB Triticum dicoccoides and AAGG Triticum araraticum lineages during their independent origins.

Quantification of genetic relationships among A genomes of wheats

The genetic relationships of A genomes of Triticum urartu (Au) and Triticum monococcum (Am) in polyploid wheats are explored and quantified by AFLP fingerprinting. Forty-one accessions of A-genome diploid wheats, 3 of AG-genome wheats, 19 of AB-genome wheats, 15 of ABD-genome wheats, and 1 of the D-genome donor Ae. tauschii have been analysed. Based on 7 AFLP primer combinations, 423 bands were identified as potentially A genome specific. The bands were reduced to 239 by eliminating those present in autoradiograms of Ae. tauschii, bands interpreted as common to all wheat genomes. Neighbour-joining analysis separates T. urartu from T. monococcum. Triticum urartu has the closest relationship to polyploid wheats. Triticum turgidum subsp. dicoccum and T. turgidum subsp. durum lines are included in tightly linked clusters. The hexaploid spelts occupy positions in the phylogenetic tree intermediate between bread wheats and T. turgidum. The AG-genome accessions cluster in a position quite distant from both diploid and other polyploid wheats. The estimates of similarity between A genomes of diploid and polyploid wheats indicate that, compared with Am, Au has around 20% higher similarity to the genomes of polyploid wheats. Triticum timo pheevii AG genome is molecularly equidistant from those of Au and Am wheats.

Phylogenetic relationships of Triticum and Aegilops and evidence for the origin of the A, B, and D genomes of common wheat (Triticum aestivum)

Common wheat (Triticum aestivum) has for decades been a textbook example of the evolution of a major crop species by allopolyploidization. Using a sophisticated extension of the PCR technique, we have successfully isolated two single-copy nuclear genes, DMC1 and EF-G, from each of the three genomes found in hexaploid wheat (BA(u)D) and from the two genomes of the tetraploid progenitor Triticum turgidum (BA(u)). By subjecting these sequences to phylogenetic analysis together with sequences from representatives of all the diploid Triticeae genera we are able for the first time to provide simultaneous and strongly supported evidence for the D genome being derived from Aegilops tauschii, the A(u) genome being derived from Triticum urartu, and the hitherto enigmatic B genome being derived from Aegilops speltoides. Previous problems of identifying the B genome donor may be associated with a higher diversification rate of the B genome compared to the A(u) genome in the polyploid wheats. The phylogenetic hypothesis further suggests that neither Triticum, Aegilops, nor Triticum plus Aegilops are monophyletic.

Molecular evidence for multiple polyploidization and lineage recombination in the Chrysanthemum indicum polyploid complex (Asteraceae)

The Chrysanthemum indicum polyploid complex comprises morphologically differentiated diploids, tetraploids and hybrids between C. indicum and C. lavandulifolium. The relationships between species and cytotypes within this complex remain poorly understood. Random amplified polymorphic DNAs (RAPDs), intersimple sequence repeats (ISSRs) and chloroplast SSR markers were used to elucidate the genetic diversity and relationships of the C. indicum polyploid complex. Molecular analysis of three diploid and nine tetraploid populations provided strong evidence for recurrent origins and lineage recombination in the C. indicum polyploid complex. The high similarity in molecular marker profiles and cpDNA haplotypes between the diploids and tetraploids distributed in the Shen-Nong-Jia Mountain area of China suggested an autopolyploid origin of the tetraploids, while the tetraploids from other populations may have originated via allopolyploidization. Lineage recombination was revealed by the extensive sharing of chloroplast haplotypes and genetic markers among the tetraploid populations with different origins. Multiple differentiation and hybridization/polyploidization cycles have led to an evolutionary reticulation in the C. indicum polyploid complex, and resulted in the difficulties in systematic classification.

Chloroplast and nuclear microsatellite analysis of Aegilops cylindrica

Aegilops cylindrica Host (2n = 4x = 28, genome CCDD) is an allotetraploid formed by hybridization between the diploid species Ae. tauschii Coss. (2n = 2x = 14, genome DD) and Ae. markgrafii (Greuter) Hammer (2n = 2x = 14, genome CC). Previous research has shown that Ae. tauschii contributed its cytoplasm to Ae. cylindrica. However, our analysis with chloroplast microsatellite markers showed that 1 of the 36 Ae. cylindrica accessions studied, TK 116 (PI 486249), had a plastome derived from Ae. markgrafii rather than Ae. tauschii. Thus, Ae. markgrafii has also contributed its cytoplasm to Ae. cylindrica. Our analysis of chloroplast and nuclear microsatellite markers also suggests that D-type plastome and the D genome in Ae. cylindrica were closely related to, and were probably derived from, the tauschii gene pool of Ae. tauschii. A determination of the likely source of the C genome and the C-type plastome in Ae. cylindrica was not possible.

Genealogical use of chloroplast DNA variation for intraspecific studies of Aegilops tauschii Coss

Intraspecific patterns of chloroplast DNA variation was studied in Aegilops tauschii Coss., the D-genome progenitor of bread wheat. Nucleotide sequences of ten chloroplast microsatellite loci were analyzed for 63 accessions that cover the central part of the species distribution. As is often the case with nuclear microsatellites, those of chloroplasts of Ae. tauschii bear complex mutations. Several types of mutations other than change in the microsatellite repeat number were found, including base substitutions and length mutations in flanking regions. In total, eight mutations were present in the flanking regions of four loci. Most mutations in the flanking regions of microsatellite repeats are associated with biallelic polymorphisms. Phylogeographic analyses showed that such biallelic polymorphisms are useful to investigate intraspecific patterns of monophyletic lineage divergence. In contrast, most microsatellite repeat sites are multiallelic, variable within intraspecific lineages, and useful to compare degrees of genetic diversity between lineages. These findings show that the chloroplast genome harbors evolutionary variations informative for intraspecific studies of Ae. tauschii and can be analyzed by genealogical approaches.