The use of chloroplast genome sequences to solve phylogenetic incongruences in Polystachya Hook (Orchidaceae Juss).

BACKGROUND: Current evidence suggests that for more robust estimates of species tree and divergence times, several unlinked genes are required. However, most phylogenetic trees for non-model organisms are based on single sequences or just a few regions, using traditional sequencing methods. Techniques for massive parallel sequencing or next generation sequencing (NGS) are an alternative to traditional methods that allow access to hundreds of DNA regions. Here we use this approach to resolve the phylogenetic incongruence found in Polystachya Hook. (Orchidaceae), a genus that stands out due to several interesting aspects, including cytological (polyploid and diploid species), evolutionary (reticulate evolution) and biogeographical (species widely distributed in the tropics and high endemism in Brazil). The genus has a notoriously complicated taxonomy, with several sections that are widely used but probably not monophyletic. METHODS: We generated the complete plastid genome of 40 individuals from one clade within the genus. The method consisted in construction of genomic libraries, hybridization to RNA probes designed from available sequences of a related species, and subsequent sequencing of the product. We also tested how well a smaller sample of the plastid genome would perform in phylogenetic inference in two ways: by duplicating a fast region and analyzing multiple copies of this dataset, and by sampling without replacement from all non-coding regions in our alignment. We further examined the phylogenetic implications of non-coding sequences that appear to have undergone hairpin inversions (reverse complemented sequences associated with small loops). RESULTS: We retrieved 131,214 bp, including coding and non-coding regions of the plastid genome. The phylogeny was able to fully resolve the relationships among all species in the targeted clade with high support values. The first divergent species are represented by African accessions and the most recent ones are among Neotropical species. DISCUSSION: Our results indicate that using the entire plastid genome is a better option than screening highly variable markers, especially when the expected tree is likely to contain many short branches. The phylogeny inferred is consistent with the proposed origin of the genus, showing a probable origin in Africa, with later dispersal into the Neotropics, as evidenced by a clade containing all Neotropical individuals. The multiple positions of Polystachya concreta (Jacq.) Garay & Sweet in the phylogeny are explained by allotetraploidy. Polystachya estrellensis Rchb.f. can be considered a genetically distinct species from P. concreta and P. foliosa (Lindl.) Rchb.f., but the delimitation of P. concreta remains uncertain. Our study shows that NGS provides a powerful tool for inferring relationships at low taxonomic levels, even in taxonomically challenging groups with short branches and intricate morphology.

Allele phasing is critical to revealing a shared allopolyploid origin of Medicago arborea and M. strasseri (Fabaceae).

BACKGROUND: Whole genome duplication plays a central role in plant evolution. There are two main classes of polyploid formation: autopolyploids which arise within one species by doubling of similar homologous genomes; in contrast, allopolyploidy (hybrid polyploidy) arise via hybridization and subsequent doubling of nonhomologous (homoeologous) genomes. The distinction between polyploid origins can be made using gene phylogenies, if alleles from each genome can be correctly retrieved. We examined whether two closely related tetraploid Mediterranean shrubs (Medicago arborea and M. strasseri) have an allopolyploid origin - a question that has remained unsolved despite substantial previous research. We sequenced and analyzed ten low-copy nuclear genes from these and related species, phasing all alleles. To test the efficacy of allele phasing on the ability to recover the evolutionary origin of polyploids, we compared these results to analyses using unphased sequences. RESULTS: In eight of the gene trees the alleles inferred from the tetraploids formed two clades, in a non-sister relationship. Each of these clades was more closely related to alleles sampled from other species of Medicago, a pattern typical of allopolyploids. However, we also observed that alleles from one of the remaining genes formed two clades that were sister to one another, as is expected for autopolyploids. Trees inferred from unphased sequences were very different, with the tetraploids often placed in poorly supported and different positions compared to results obtained using phased alleles. CONCLUSIONS: The complex phylogenetic history of M. arborea and M. strasseri is explained predominantly by shared allotetraploidy. We also observed that an increase in woodiness is correlated with polyploidy in this group of species and present a new possibility that woodiness could be a transgressive phenotype. Correctly phased homoeologues are likely to be critical for inferring the hybrid origin of allopolyploid species, when most genes retain more than one homoeologue. Ignoring homoeologous variation by merging the homoeologues can obscure the signal of hybrid polyploid origins and produce inaccurate results.

Assignment of homoeologs to parental genomes in allopolyploids for species tree inference, with an example from Fumaria (papaveraceae).

There is a rising awareness that species trees are best inferred from multiple loci while taking into account processes affecting individual gene trees, such as substitution model error (failure of the model to account for the complexity of the data) and coalescent stochasticity (presence of incomplete lineage sorting [ILS]). Although most studies have been carried out in the context of dichotomous species trees, these processes operate also in more complex evolutionary histories involving multiple hybridizations and polyploidy. Recently, methods have been developed that accurately handle ILS in allopolyploids, but they are thus far restricted to networks of diploids and tetraploids. We propose a procedure that improves on this limitation by designing a workflow that assigns homoeologs to hypothetical diploid ancestral genomes prior to genome tree construction. Conflicting assignment hypotheses are evaluated against substitution model error and coalescent stochasticity. Incongruence that cannot be explained by stochastic mechanisms needs to be explained by other processes (e.g., homoploid hybridization or paralogy). The data can then be filtered to build multilabeled genome phylogenies using inference methods that can recover species trees, either in the face of substitution model error and coalescent stochasticity alone, or while simultaneously accounting for hybridization. Methods are already available for folding the resulting multilabeled genome phylogeny into a network. We apply the workflow to the reconstruction of the reticulate phylogeny of the plant genus Fumaria (Papaveraceae) with ploidal levels ranging from 2[Formula: see text] to 14[Formula: see text]. We describe the challenges in recovering nuclear NRPB2 homoeologs in high ploidy species while combining in vivo cloning and direct sequencing techniques. Using parametric bootstrapping simulations we assign nuclear homoeologs and chloroplast sequences (four concatenated loci) to their common hypothetical diploid ancestral genomes. As these assignments hinge on effective population size assumptions, we investigate how varying these assumptions impacts the recovered multilabeled genome phylogeny.

Stereopsidales--a new order of mushroom-forming fungi.

One new order, one new family, and one new combination are presented, as the result of molecular phylogenetic analyses. The new order Stereopsidales and the new family Stereopsidaceae are described incorporating Stereopsis radicans and S. globosa, formerly Clavulicium globosum. We show that not only do these species represent an old overlooked lineage, but both species harbor cryptic diversity. In addition, a third species, C. macounii, appears as a plausible sister to the new lineage, but there is conflict in the data. All specimens of S. radicans and S. globosa analysed here are from the South and Central Americas; several records of S. radicans have been made also from tropical Asia. We expect the true diversity in this group to be a lot higher than presented in this paper. Stereopsis radicans was formerly included in Polyporales, but a placement within that order is rejected by our data through SH tests. The dataset consisted of four nuclear markers: rpb2, tef1, LSU and SSU, each of which was analysed separately using maximum likelihood and Bayesian inference. Recombination detection tests indicate no plausible recombinations. The potential of S. radicans, S. globosa and C. macounii being amphitallic is briefly discussed.

Evolution of a complex disease resistance gene cluster in diploid Phaseolus and tetraploid Glycine.

We used a comparative genomics approach to investigate the evolution of a complex nucleotide-binding (NB)-leucine-rich repeat (LRR) gene cluster found in soybean (Glycine max) and common bean (Phaseolus vulgaris) that is associated with several disease resistance (R) genes of known function, including Rpg1b (for Resistance to Pseudomonas glycinea1b), an R gene effective against specific races of bacterial blight. Analysis of domains revealed that the amino-terminal coiled-coil (CC) domain, central nucleotide-binding domain (NB-ARC [for APAF1, Resistance genes, and CED4]), and carboxyl-terminal LRR domain have undergone distinct evolutionary paths. Sequence exchanges within the NB-ARC domain were rare. In contrast, interparalogue exchanges involving the CC and LRR domains were common, consistent with both of these regions coevolving with pathogens. Residues under positive selection were overrepresented within the predicted solvent-exposed face of the LRR domain, although several also were detected within the CC and NB-ARC domains. Superimposition of these latter residues onto predicted tertiary structures revealed that the majority are located on the surface, suggestive of a role in interactions with other domains or proteins. Following polyploidy in the Glycine lineage, NB-LRR genes have been preferentially lost from one of the duplicated chromosomes (homeologues found in soybean), and there has been partitioning of NB-LRR clades between the two homeologues. The single orthologous region in common bean contains approximately the same number of paralogues as found in the two soybean homeologues combined. We conclude that while polyploidization in Glycine has not driven a stable increase in family size for NB-LRR genes, it has generated two recombinationally isolated clusters, one of which appears to be in the process of decay.

Use of allele-specific sequencing primers is an efficient alternative to PCR subcloning of low-copy nuclear genes.

Direct Sanger sequencing of polymerase chain reaction (PCR)-amplified nuclear genes leads to polymorphic sequences when allelic variation is present. To overcome this problem, most researchers subclone the PCR products to separate alleles. An alternative is to directly sequence the separate alleles using allele-specific primers. We tested two methods to enhance the specificity of allele-specific primers for use in direct sequencing: using short primers and amplification refractory mutation system (ARMS) technique. By shortening the allele-specific primer to 15-13 nucleotides, the single mismatch in the ultimate base of the primer is enough to hinder the amplification of the nontarget allele in direct sequencing and recover only the targeted allele at high accuracy. The deliberate addition of a second mismatch, as implemented in the ARMS technique, was less successful and seems better suited for allele-specific amplification in regular PCR rather than in direct sequencing.

A molecular phylogeny of the orange subfamily(Rutaceae: Aurantioideae) using nine cpDNA sequences.

The breeding of new, high-quality citrus cultivars depends on dependable information about the relationships of taxa within the tribe Citreae; therefore, it is important to have a well-supported phylogeny of the relationships between species not only to advance breeding strategies, but also to advance conservation strategies for the wild taxa. The recent history of the systematics of Citrus (Rutaceae: Aurantioideae) and its allies, in the context of Rutaceae taxonomy as a whole, is reviewed. The most recent classification is tested using nine cpDNA sequence regions in representatives of all genera of the subfam. Aurantioideae (save Limnocitrus) and numerous species and hybrids referred to Citrus s.l. Aurantioideae are confirmed as monophyletic. Within Aurantioideae, tribe Clauseneae are not monophyletic unless Murraya s.s. and Merrillia are removed to Aurantieae. Within tribe Aurantieae, the three traditionally recognized subtribes are not monophyletic. Triphasiinae is not monophyletic unless Oxanthera is returned to Citrus (Citrinae). Balsamocitrinae is polyphyletic. Feroniella, traditionally considered allied closely to Limonia (=Feronia), is shown to be nested in Citrus. The proposed congenericity of Severinia and Atalantia is confirmed. The most recent circumscription of Citrus is strongly supported by this analysis, with hybrids appearing with their putative maternal parents. The genus was resolved into two clades, one comprising wild species from New Guinea, Australia, and New Caledonia (formerly Clymenia, Eremocitrus, Microcitrus, Oxanthera), but surprisingly also Citrus medica, traditionally believed to be native in India. The second clade is largely from the Asian mainland (including species formerly referred to Fortunella and Poncirus).

Differential accumulation of retroelements and diversification of NB-LRR disease resistance genes in duplicated regions following polyploidy in the ancestor of soybean.

The genomes of most, if not all, flowering plants have undergone whole genome duplication events during their evolution. The impact of such polyploidy events is poorly understood, as is the fate of most duplicated genes. We sequenced an approximately 1 million-bp region in soybean (Glycine max) centered on the Rpg1-b disease resistance gene and compared this region with a region duplicated 10 to 14 million years ago. These two regions were also compared with homologous regions in several related legume species (a second soybean genotype, Glycine tomentella, Phaseolus vulgaris, and Medicago truncatula), which enabled us to determine how each of the duplicated regions (homoeologues) in soybean has changed following polyploidy. The biggest change was in retroelement content, with homoeologue 2 having expanded to 3-fold the size of homoeologue 1. Despite this accumulation of retroelements, over 77% of the duplicated low-copy genes have been retained in the same order and appear to be functional. This finding contrasts with recent analyses of the maize (Zea mays) genome, in which only about one-third of duplicated genes appear to have been retained over a similar time period. Fluorescent in situ hybridization revealed that the homoeologue 2 region is located very near a centromere. Thus, pericentromeric localization, per se, does not result in a high rate of gene inactivation, despite greatly accelerated retrotransposon accumulation. In contrast to low-copy genes, nucleotide-binding-leucine-rich repeat disease resistance gene clusters have undergone dramatic species/homoeologue-specific duplications and losses, with some evidence for partitioning of subfamilies between homoeologues.

The reticulate history of Medicago (Fabaceae).

The phylogenetic history of Medicago was examined for 60 accessions from 56 species using two nuclear genes (CNGC5 and beta-cop) and one mitochondrial region (rpS14-cob). The results of several analyses revealed that extensive robustly supported incongruence exists among the nuclear genes, the cause of which we seek to explain. After rejecting several processes, hybridization and lineage sorting of ancestral polymorphisms remained as the most likely factors promoting incongruence. Using coalescence simulations, we rejected lineage sorting alone as an explanation of the differences among gene trees. The results indicate that hybridization has been common and ongoing among lineages since the origin of Medicago. Coalescence provides a good framework to test the causes of incongruence commonly seen among gene trees but requires knowledge of effective population sizes and generation times. We estimated the effective population size at 240,000 individuals and assumed a generation time of 1 year in Medicago (many are annual plants). A sensitivity analysis showed that our conclusions remain unchanged using a larger effective population size and/or longer generation time.