Testing four candidate barcoding markers in temperate woody bamboos (Poaceae: Bambusoideae)

Abstract Bambusoideae is an important subfamily of the grass family Poaceae that has considerable economic, ecologic and cultural value. In addition, Bambusoideae species are important constituents of the forest vegetation in China. Because of the paucity of flower‐bearing specimens and homoplasies of morphological characters, it is difficult to identify species of Bambusoideae using morphology alone, especially in the case of temperate woody bamboos (i.e. Arundinarieae). To this end, DNA barcoding has shown great potential in identifying species. The present study is the first attempt to test the feasibility of four proposed DNA barcoding markers (matK, rbcL, trnH–psbA, and internal transcribed spacer [ITS]) in identifying 27 species of the temperate woody bamboos. Three plastid markers showed high levels of universality, whereas the universality of ITS was comparatively low. A single plastid marker provided low levels of discrimination success at both the genus and species levels (<12%). Among the combinations of plastid markers, the highest discriminatory power was obtained using the combination of rbcL+matK (14.8%). Using a combination of three markers did not increase species discrimination. The nuclear region ITS alone could identify 66.7% of species, although fewer taxa were included in the ITS analyses than in the plastid analyses. When ITS was integrated with a single or combination of plastid markers, the species discriminatory power was significantly improved. We suggest that a combination of rbcL+ ITS, which exhibited the highest species identification power of all combinations in the present study, could be used as a potential DNA barcode for temperate woody bamboos.     

Identification of putative orthologous genes for the phylogenetic reconstruction of temperate woody bamboos (Poaceae: Bambusoideae)

The temperate woody bamboos (Arundinarieae) are highly diverse in morphology but lack a substantial amount of genetic variation. The taxonomy of this lineage is intractable, and the relationships within the tribe have not been well resolved. Recent studies indicated that this tribe could have a complex evolutionary history. Although phylogenetic studies of the tribe have been carried out, most of these phylogenetic reconstructions were based on plastid data, which provide lower phylogenetic resolution compared with nuclear data. In this study, we intended to identify a set of desirable nuclear genes for resolving the phylogeny of the temperate woody bamboos. Using two different methodologies, we identified 209 and 916 genes, respectively, as putative single copy orthologous genes. A total of 112 genes was successfully amplified and sequenced by next‐generation sequencing technologies in five species sampled from the tribe. As most of the genes exhibited intra‐individual allele heterozygotes, we investigated phylogenetic utility by reconstructing the phylogeny based on individual genes. Discordance among gene trees was observed and, to resolve the conflict, we performed a range of analyses using BUCKy and HybTree. While caution should be taken when inferring a phylogeny from multiple conflicting genes, our analysis indicated that 74 of the 112 investigated genes are potential markers for resolving the phylogeny of the temperate woody bamboos.     

High-throughput sequencing of six bamboo chloroplast genomes: phylogenetic implications for temperate woody bamboos (Poaceae: Bambusoideae)

Background

Bambusoideae is the only subfamily that contains woody members in the grass family, Poaceae. In phylogenetic analyses, Bambusoideae, Pooideae and Ehrhartoideae formed the BEP clade, yet the internal relationships of this clade are controversial. The distinctive life history (infrequent flowering and predominance of asexual reproduction) of woody bamboos makes them an interesting but taxonomically difficult group. Phylogenetic analyses based on large DNA fragments could only provide a moderate resolution of woody bamboo relationships, although a robust phylogenetic tree is needed to elucidate their evolutionary history. Phylogenomics is an alternative choice for resolving difficult phylogenies.

Methodology/Principal Findings

Here we present the complete nucleotide sequences of six woody bamboo chloroplast (cp) genomes using Illumina sequencing. These genomes are similar to those of other grasses and rather conservative in evolution. We constructed a phylogeny of Poaceae from 24 complete cp genomes including 21 grass species. Within the BEP clade, we found strong support for a sister relationship between Bambusoideae and Pooideae. In a substantial improvement over prior studies, all six nodes within Bambusoideae were supported with ≥0.95 posterior probability from Bayesian inference and 5/6 nodes resolved with 100% bootstrap support in maximum parsimony and maximum likelihood analyses. We found that repeats in the cp genome could provide phylogenetic information, while caution is needed when using indels in phylogenetic analyses based on few selected genes. We also identified relatively rapidly evolving cp genome regions that have the potential to be used for further phylogenetic study in Bambusoideae.

Conclusions/Significance

The cp genome of Bambusoideae evolved slowly, and phylogenomics based on whole cp genome could be used to resolve major relationships within the subfamily. The difficulty in resolving the diversification among three clades of temperate woody bamboos, even with complete cp genome sequences, suggests that these lineages may have diverged very rapidly.     

Genetic Variation and Evolution of the Alpine Bamboos (Poaceae: Bambusoideae) using DNA Sequence Data

Thamnocalamus , Fargesia and Yushania, of the alpine bamboos and one species of Ampelocalamus as an out-group were studied. The results indicated that Thamnocalamus spathiflorus var. crassinodus and the Fargesia spathacea clade form the basal groups but bootstrap support was weak. Among the rest of the species, including species previously placed in Fargesia (plus Borinda) and Yushania, the F. yunnanensis subclade and the F. communis subclade were recognized but internal support for such groups was again low. The result indicated that, Fargesia and Yushania as delimited by morphological characters, are not monophyletic in the ITS phylogeny and require further resolution. We revealed relatively high levels of genetic variability in the alpine bamboos and showed that the ITS region could be used to improve generic delimitation of the woody bamboos in general. Received 18 September 2000/ Accepted in revised form 9 May 2001     

Identification of sex-linked SNP markers using RAD sequencing suggests ZW/ZZ sex determination in Pistacia vera L.

Background

Pistachio (Pistacia vera L.) is a dioecious species that has a long juvenility period. Therefore, development of marker-assisted selection (MAS) techniques would greatly facilitate pistachio cultivar-breeding programs. The sex determination mechanism is presently unknown in pistachio. The generation of sex-linked markers is likely to reduce time, labor, and costs associated with breeding programs, and will help to clarify the sex determination system in pistachio.

Results

Restriction site-associated DNA (RAD) markers were used to identify sex-linked markers and to elucidate the sex determination system in pistachio. Eight male and eight female F₁ progenies from a Pistacia vera L. Siirt × Bağyolu cross, along with the parents, were subjected to RAD sequencing in two lanes of a Hi-Seq 2000 sequencing platform. This generated 449 million reads, comprising approximately 37.7 Gb of sequences. There were 33,757 polymorphic single nucleotide polymorphism (SNP) loci between the parents. Thirty-eight of these, from 28 RAD reads, were detected as putative sex-associated loci in pistachio. Validation was performed by SNaPshot analysis in 42 mature F₁ progenies and in 124 cultivars and genotypes in a germplasm collection. Eight loci could distinguish sex with 100% accuracy in pistachio. To ascertain cost-effective application of markers in a breeding program, high-resolution melting (HRM) analysis was performed; four markers were found to perfectly separate sexes in pistachio. Because of the female heterogamety in all candidate SNP loci, we report for the first time that pistachio has a ZZ/ZW sex determination system. As the reported female-to-male segregation ratio is 1:1 in all known segregating populations and there is no previous report of super-female genotypes or female heteromorphic chromosomes in pistachio, it appears that the WW genotype is not viable.

Conclusion

Sex-linked SNP markers were identified and validated in a large germplasm and proved their suitability for MAS in pistachio. HRM analysis successfully validated the sex-linked markers for MAS. For the first time in dioecious pistachio, a female heterogamety ZW/ZZ sex determination system is suggested.     

Identification of a large SNP dataset in Larimichthys crocea using specific‐locus amplified fragment sequencing

The large yellow croaker, Larimichthys crocea, is a commercially important drum fish (Family: Sciaenidae) native to the East and South China Sea. Habitat deterioration and overfishing have led to significant population decline and the collapse of its fishery over the past decades. Today, the market supply of L. crocea depends solely on stocks produced in hatcheries and farms. Common issues that occur in the culture of L. crocea include germplasm degradation, precocious puberty, elevated disease susceptibility and growth retardation. In this study, we employed SLAF‐seq (specific‐locus amplified fragment sequencing) technology to identify single nucleotide polymorphism (SNP) loci across the L. crocea genome. Sixty samples were selected for SLAF analysis out of 1000 progeny in the same cohort of a cultured stock. Our analysis obtained a total of 151 253 SLAFs, of which 65.88% (99 652) were identified to be polymorphic, scoring a total of 710 567 putative SNPs. Further filtration resulted in a final panel of 1782 SNP loci. The data derived from this work could be beneficial for understanding the genetics of complex phenotypic traits as well as for developing marker‐selection‐assisted breeding programs in L. crocea.     

Non-monophyly of the woody bamboos (Bambuseae; Poaceae): a multi-gene region phylogenetic analysis of Bambusoideae s.s.

The taxonomy of Bambusoideae is in a state of flux and phylogenetic studies are required to help resolve systematic issues. Over 60 taxa, representing all subtribes of Bambuseae and related non-bambusoid grasses were sampled. A combined analysis of five plastid DNA regions, trnL intron, trnL-F intergenic spacer, atpB-rbcL intergenic spacer, rps16 intron, and matK, was used to study the phylogenetic relationships among the bamboos in general and the woody bamboos in particular. Within the BEP clade (Bambusoideae s.s., Ehrhartoideae, Pooideae), Pooideae were resolved as sister to Bambusoideae s.s. Tribe Bambuseae, the woody bamboos, as currently recognized were not monophyletic because Olyreae, the herbaceous bamboos, were sister to tropical Bambuseae. Temperate Bambuseae were sister to the group consisting of tropical Bambuseae and Olyreae. Thus, the temperate Bambuseae would be better treated as their own tribe Arundinarieae than as a subgroup of Bambuseae. Within the tropical Bambuseae, neotropical Bambuseae were sister to the palaeotropical and Austral Bambuseae. In addition, Melocanninae were found to be sister to the remaining palaeotropical and Austral Bambuseae. We discuss phylogenetic and morphological patterns of diversification and interpret them in a biogeographic context.     

RAD sequencing,morphometry and synecology clarify the taxonomy of the Melica ciliata (Poaceae) complex in France and Poland

Melica (Poaceae) consist of about 92 species distributed across temperate regions of the world. Within section Dalycum, Melica ciliata sensu lato forms a taxonomic complex of several species and subspecies with clinal morphological variation causing conflicting identifications. To resolve taxonomic confusion, we used three complementary approaches, through molecular, morphological, and phytoecological analyses. The double-digest restriction-associated DNA markers significantly support the monophyly of three taxa: (i) the Mediterranean Melica magnolii, (ii) the Eurasian Melica transsilvanica subsp. transsilvanica, and (iii) the west-European M. ciliata subsp. glauca. This differentiation is corroborated by the analysis of 22 morphometric variables. Furthermore, phytoecological analysis of 221 floristic inventories revealed habitat distinctions among these taxa. Our approach of integrative taxonomy argues for a specific distinction for these three taxa, and we include a key to separate these forms. These new molecular data on the section Dalycum, subsection Ciliatae, call for further phylogenetic analyses including samples of M. ciliata subsp. ciliata and other East-Mediterranean and South African taxa.     

Mass production of SNP markers in a nonmodel passerine bird through RAD sequencing and contig mapping to the zebra finch genome

Here, we present an adaptation of restriction‐site‐associated DNA sequencing (RAD‐seq) to the Illumina HiSeq2000 technology that we used to produce SNP markers in very large quantities at low cost per unit in the Réunion grey white‐eye (Zosterops borbonicus), a nonmodel passerine bird species with no reference genome. We sequenced a set of six pools of 18–25 individuals using a single sequencing lane. This allowed us to build around 600 000 contigs, among which at least 386 000 could be mapped to the zebra finch (Taeniopygia guttata) genome. This yielded more than 80 000 SNPs that could be mapped unambiguously and are evenly distributed across the genome. Thus, our approach provides a good illustration of the high potential of paired‐end RAD sequencing of pooled DNA samples combined with comparative assembly to the zebra finch genome to build large contigs and characterize vast numbers of informative SNPs in nonmodel passerine bird species in a very efficient and cost‐effective way.     

Genome‐wide sequence data suggest the possibility of pollinator sharing by host shift in dioecious figs (Moraceae,Ficus)

The obligate mutualism of figs and fig‐pollinating wasps has been one of the classic models used for testing theories of co‐evolution and cospeciation due to the high species‐specificity of these relationships. To investigate the species‐specificity between figs and fig pollinators and to further understand the speciation process in obligate mutualisms, we examined the genetic differentiation and phylogenetic relationships of four closely related fig‐pollinating wasp species (Blastophaga nipponica, Blastophaga taiwanensis, Blastophaga tannoensis and Blastophaga yeni) in Japan and Taiwan using genome‐wide sequence data, including mitochondrial DNA sequences. In addition, population structure was analysed for the fig wasps and their host species using microsatellite data. The results suggest that the three Taiwanese fig wasp species are a single panmictic population that pollinates three dioecious fig species, which are sympatrically distributed, have large differences in morphology and ecology and are also genetically differentiated. Our results illustrate the first case of pollinator sharing by host shift in the subgenus Ficus. On the other hand, there are strict genetic codivergences between allopatric populations of the two host–pollinator pairs. The possible processes that produce these pollinator‐sharing events are discussed based on the level and pattern of genetic differentiation in these figs and fig wasps.     

Genotyping by RAD sequencing enables mapping of fatty acid composition traits in perennial ryegrass (Lolium perenne (L.))

Perennial ryegrass (Lolium perenne L.) is the most important forage crop in temperate livestock agriculture. Its nutritional quality has significant impact on the quality of meat and milk for human consumption. Evidence suggests that higher energy content in forage can assist in reducing greenhouse gas emissions from ruminants. Increasing the fatty acid content (especially α‐linolenic acid, an omega‐3 fatty acid) may thus contribute to better forage, but little is known about the genetic basis of variation for this trait. To this end, quantitative trait loci (QTLs) were identified associated with major fatty acid content in perennial ryegrass using a population derived from a cross between the heterozygous and outbreeding high‐sugar grass variety AberMagic and an older variety, Aurora. A genetic map with 434 restriction‐associated DNA (RAD) and SSR markers was generated. Significant QTLs for the content of palmitic (C16:0) on linkage groups (LGs) 2 and 7; stearic (C18:0) on LGs 3, 4 and 7; linoleic (C18:2n‐6) on LGs 2 and 5; and α‐linolenic acids (C18:3n‐3) on LG 1 were identified. Two candidate genes (a lipase and a beta‐ketoacyl CoA synthase), both associated with C16:0, and separately with C18:2n‐6 and C18:0 contents, were identified. The physical positions of these genes in rice and their genetic positions in perennial ryegrass were consistent with established syntenic relationships between these two species. Validation of these associations is required, but the utility of RAD markers for rapid generation of genetic maps and QTL analysis has been demonstrated for fatty acid composition in a global forage crop.     

Differentiation in a geographical mosaic of plants coevolving with ants: phylogeny of the Leonardoxa africana complex (Fabaceae: Caesalpinioideae) using amplified fragment length polymorphism markers

Comprising four allopatric subspecies that exhibit various grades of ant-plant interactions, from diffuse to obligate and symbiotic associations, the Leonardoxa africana complex (Fabaceae, Caesalpinioideae) provides a good opportunity to investigate the evolutionary history of ant-plant mutualisms. A previous study of the L. africana complex based on chloroplast DNA noncoding sequences revealed a lack of congruence between clades suggested by morphological and plastid characters. In this study, we analysed phylogenetic relationships within the L. africana complex using a Bayesian probability approach on amplified fragment length polymorphism markers. The results reported permit partial validation of the four subspecies of L. africana previously defined by morphological and ecological markers. Incongruences between phylogenies based on chloroplast DNA and amplified fragment length polymorphism markers are discussed in the light of morphological and ecological data, and confronted with hypotheses of convergence, lineage sorting and introgression.     

On polyphyly of the former section Ochlopoa and the hybridogenic section Acroleucae (Poa,Poaceae): insights from molecular phylogenetic analyses

In this study, sequence data from the inert nuclear region ITS1‐5.8S rDNA‐ITS2 and the chloroplast region trnL–F, as well as a few morphological characters, are analysed to the relationships among known annual Poa (bluegrasses). It is shown that all taxa from the Poa annua aggregate distinguished by lemma characters and growth form have identical ITS and trnL–trnF sequences, all ITS sequences of this aggregate are the same as thethose of P. supina, and all trnL–trnF sequences are homologous with those of P. infirma. Furthermore, no differences were found between unusual morphotypes of P. supina with short spinules on their panicle branches and typical plants, but Siberian samples were found to have a slightly differentiated trnL–trnF region. These results suggest a hybrid origin of the Asian annual bluegrasses. Their maternal genome is close to that of P. sect. Homalopoa, but their ITS sequences are different. Some ITS sequences from annual Asian bluegrasses are resolved among representatives of P. sect. Stenopoa while for other (morphologically closely similar) species they fall in a clade with representatives of P. sect. Malacanthae. The latter group is distant from P. sect. Ochlopoa and is better treated as a separate section, viz P. sect. Acroleucae. The American annual bluegrasses are heterogeneous and also rather distant from P. sect. Ochlopoa. Poa chapmaniana, a species with cleistogamic flowers, is nested among the basal Subantarctic sections, far away from the taxa with which it has previously been considered related. It is indeed closer to P. sect. Ochlopoa than to other annual American bluegrasses. Thus, the studied annual species in fact belong to four independent evolutionary lines (or six including the separate genus Eremopoa and the Turkish Poa jubata), one of which, Acroleucae, has gone through three reticulation events. As in previous studies, our analysis did not support the generic status of P. sect. Ochlopoa.     

RAD sequencing reveals genomewide divergence between independent invasions of the European green crab (Carcinus maenas) in the Northwest Atlantic

Genomic studies of invasive species can reveal both invasive pathways and functional differences underpinning patterns of colonization success. The European green crab (Carcinus maenas) was initially introduced to eastern North America nearly 200 years ago where it expanded northwards to eastern Nova Scotia. A subsequent invasion to Nova Scotia from a northern European source allowed further range expansion, providing a unique opportunity to study the invasion genomics of a species with multiple invasions. Here, we use restriction‐site‐associated DNA sequencing‐derived SNPs to explore fine‐scale genomewide differentiation between these two invasions. We identified 9137 loci from green crab sampled from 11 locations along eastern North America and compared spatial variation to mitochondrial COI sequence variation used previously to characterize these invasions. Overall spatial divergence among invasions was high (pairwise F_ST ~0.001 to 0.15) and spread across many loci, with a mean F_ST ~0.052 and 52% of loci examined characterized by F_ST values >0.05. The majority of the most divergent loci (i.e., outliers, ~1.2%) displayed latitudinal clines in allele frequency highlighting extensive genomic divergence among the invasions. Discriminant analysis of principal components (both neutral and outlier loci) clearly resolved the two invasions spatially and was highly correlated with mitochondrial divergence. Our results reveal extensive cryptic intraspecific genomic diversity associated with differing patterns of colonization success and demonstrates clear utility for genomic approaches to delineating the distribution and colonization success of aquatic invasive species.     

Isolation and characterization of microsatellite markers for Achyranthes bidentata (Amaranthaceae) using next-generation sequencing platform

Achyranthes bidentata Blume (Amaranthaceae) is a perennial herb that is widely distributed in India, Java, China, and Japan. The natural resources of A. bidentata in its geo-authentic product area have rapidly declined in recent years because of the over-collection of its roots. To devise adequate conservation and management strategies for this species, its genetic diversity and population structure should be characterized. Roche 454 pyrosequencing combined with magnetic bead enrichment was used to develop microsatellite markers for A. bidentata. A total of 903 microsatellite loci were identified from 42,004 individual sequence reads. One hundred microsatellite loci were selected to test the primer amplification efficiency across 16 individuals from two A. bidentata populations. Of these tested markers, 8 yielded polymorphic amplification products, 29 yielded single alleles. For polymorphic primer pairs, the number of alleles per locus ranged from 2 to 4, with an average of 2.75. The observed and expected heterozygosities ranged from 0.353 to 0.671 and 0.250 to 0.938, respectively. The inbreeding coefficient varied from −0.692 to 0.627. This set of markers will provide useful tools for examining genetic diversity and population structure, and aid in better understanding of the conservation of A. bidentata.     

油棕种壳厚度控制基因SHELL的SNP分子标记开发

油棕属棕榈科多年生木本油料作物,果实含油量高达50%,且单位面积产油量高,享有"世界油王"美誉。油棕果实由外果皮、中果皮、内果皮(种壳)、种子四个部分组成,产油部分主要是中果皮和种子,其中种壳厚度是影响果实含油量的重要因素。SHELL基因控制种壳厚度,是一类MADS-box同源基因,SHELL基因在厚壳种和无壳种中的变异主要是第一个外显子上的两个SNP位点。该研究根据两个SNP位点进行特异标记开发,根据已知的油棕SHELL基因的序列,设计了4对SNP引物。4对SNP引物以2个SNP位点设计,每个SNP位点设计2对SNP标记,并均在引物3'端第二位引入强错配碱基。以2份薄壳种油棕材料和2份厚壳种油棕材料DNA为模板,扩增筛选油棕SHELL基因SNP引物。通过PCR扩增发现,设计的SHELL基因特异SNP标记EgSh(N)-f/EgSh(SNP)-2r能够鉴别油棕厚壳种和薄壳种。再用24株油棕树进行特异性验证,发现该标记能较准确地判断油棕的厚薄壳。该研究结果表明SNP标记EgSh(N)-f/EgSh(SNP)-2r可用来进行油棕种质资源早期分子鉴定,为高产油棕品种选育提供了技术支撑。     

Molecular phylogeny of the Canary Island lacertids (Gallotia): mitochondrial DNA restriction fragment divergence in relation to sequence divergence and geological time

Restriction fragment length polymorphisms of 6 base pair recognising endonucleases are used to reconstruct the phylogeny of the endemic Canary Island lacertid, Gallotia. The division into conventional species is upheld by this molecular analysis and the western Canary Island lizard (G. galloti) and eastern Canary Island lizard (G. atlantica) are hypothesized to be sister species. A more comprehensive study of the intraspecific relationships of G. galloti, based on nineteen restriction enzymes, indicates that there are distinct southern and northern lineages within this species. The phylogenetic analysis does not uphold the conventional subspecies, but suggests an alternative arrangement with one northern (La Palma, Tenerife) and one southern (Gomera, Hierro) subspecies. The inferred timing of molecular divergence of populations of G. galloti, based on RFLP analysis, is compatible with the geological timing for island origin and fossil data. Mantel tests show that mitochondrial RFLP divergence is correlated with mitrochondrial 12s rRNA and cytochrome oxidase I sequence divergence and highly correlated with mitochondrial cytochrome b sequence divergence.