The impact of intragenic recombination on phylogenetic reconstruction at the sectional level in Eucalyptus when using a single copy nuclear gene (cinnamoyl CoA reductase)

Choosing and Using Introns in Molecular Phylogenetics

Introns are now commonly used in molecular phylogenetics in an attempt to recover gene trees that are concordant with species trees, but there are a range of genomic, logistical and analytical considerations that are infrequently discussed in empirical studies that utilize intron data. This review outlines expedient approaches for locus selection, overcoming paralogy problems, recombination detection methods and the identification and incorporation of LVHs in molecular systematics. A range of parsimony and Bayesian analytical approaches are also described in order to highlight the methods that can currently be employed to align sequences and treat indels in subsequent analyses. By covering the main points associated with the generation and analysis of intron data, this review aims to provide a comprehensive introduction to using introns (or any non-coding nuclear data partition) in contemporary phylogenetics.

Rice Chloroplast Genome Variation Architecture and Phylogenetic Dissection in Diverse Oryza Species Assessed by Whole-Genome Resequencing

BACKGROUND

Chloroplast genome variations have been detected, despite its overall conserved structure, which has been valuable for plant population genetics and evolutionary studies. Here, we described chloroplast variation architecture of 383 rice accessions from diverse regions and different ecotypes, in order to mine the rice chloroplast genome variation architecture and phylogenetic.

RESULTS

A total of 3677 variations across the chloroplast genome were identified with an average density of 27.33 per kb, in which wild rice showing a higher variation density than cultivated groups. Chloroplast genome nucleotide diversity investigation indicated a high degree of diversity in wild rice than in cultivated rice. Genetic distance estimation revealed that African rice showed a low level of breeding and connectivity with the Asian rice, suggesting the big distinction of them. Population structure and principal component analysis revealed the existence of clear clustering of African and Asian rice, as well as the indica and japonica in Asian cultivated rice. Phylogenetic analysis based on maximum likelihood and Bayesian inference methods and the population splits test suggested and supported the independent origins of indica and japonica within Asian cultivated rice. In addition, the African cultivated rice was thought to be domesticated differently from Asian cultivated rice.

CONCLUSIONS

The chloroplast genome variation architecture in Asian and African rice are different, as well as within Asian or African rice. Wild rice and cultivated rice also have distinct nucleotide diversity or genetic distance. In chloroplast level, the independent origins of indica and japonica within Asian cultivated rice were suggested and the African cultivated rice was thought to be domesticated differently from Asian cultivated rice. These results will provide more candidate evidence for the further rice chloroplast genomic and evolution studies.

Transcriptionally active LTR retrotransposons in Eucalyptus genus are differentially expressed and insertionally polymorphic

BACKGROUND

In Eucalyptus genus, studies on genome composition and transposable elements (TEs) are particularly scarce. Nearly half of the recently released Eucalyptus grandis genome is composed by retrotransposons and this data provides an important opportunity to understand TE dynamics in Eucalyptus genome and transcriptome.

RESULTS

We characterized nine families of transcriptionally active LTR retrotransposons from Copia and Gypsy superfamilies in Eucalyptus grandis genome and we depicted genomic distribution and copy number in two Eucalyptus species. We also evaluated genomic polymorphism and transcriptional profile in three organs of five Eucalyptus species. We observed contrasting genomic and transcriptional behavior in the same family among different species. RLC_egMax_1 was the most prevalent family and RLC_egAngela_1 was the family with the lowest copy number. Most families of both superfamilies have their insertions occurring <3 million years, except one Copia family, RLC_egBianca_1. Protein theoretical models suggest different properties between Copia and Gypsy domains. IRAP and REMAP markers suggested genomic polymorphisms among Eucalyptus species. Using EST analysis and qRT-PCRs, we observed transcriptional activity in several tissues and in all evaluated species. In some families, osmotic stress increases transcript values.

CONCLUSION

Our strategy was successful in isolating transcriptionally active retrotransposons in Eucalyptus, and each family has a particular genomic and transcriptional pattern. Overall, our results show that retrotransposon activity have differentially affected genome and transcriptome among Eucalyptus species.

A six nuclear gene phylogeny of Citrus (Rutaceae) taking into account hybridization and lineage sorting

Background

Genus Citrus (Rutaceae) comprises many important cultivated species that generally hybridize easily. Phylogenetic study of a group showing extensive hybridization is challenging. Since the genus Citrus has diverged recently (4–12 Ma), incomplete lineage sorting of ancestral polymorphisms is also likely to cause discrepancies among genes in phylogenetic inferences. Incongruence of gene trees is observed and it is essential to unravel the processes that cause inconsistencies in order to understand the phylogenetic relationships among the species.

Methodology and Principal Findings

(1) We generated phylogenetic trees using haplotype sequences of six low copy nuclear genes. (2) Published simple sequence repeat data were re-analyzed to study population structure and the results were compared with the phylogenetic trees constructed using sequence data and coalescence simulations. (3) To distinguish between hybridization and incomplete lineage sorting, we developed and utilized a coalescence simulation approach. In other studies, species trees have been inferred despite the possibility of hybridization having occurred and used to generate null distributions of the effect of lineage sorting alone (by coalescent simulation). Since this is problematic, we instead generate these distributions directly from observed gene trees. Of the six trees generated, we used the most resolved three to detect hybrids. We found that 11 of 33 samples appear to be affected by historical hybridization. Analysis of the remaining three genes supported the conclusions from the hybrid detection test.

Conclusions

We have identified or confirmed probable hybrid origins for several Citrus cultivars using three different approaches–gene phylogenies, population structure analysis and coalescence simulation. Hybridization and incomplete lineage sorting were identified primarily based on differences among gene phylogenies with reference to null expectations via coalescence simulations. We conclude that identifying hybridization as a frequent cause of incongruence among gene trees is critical to correctly infer the phylogeny among species of Citrus.

Chloroplast genome sequence confirms distinctness of Australian and Asian wild rice

Cultivated rice (Oryza sativa) is an AA genome Oryza species that was most likely domesticated from wild populations of O. rufipogon in Asia. O. rufipogon and O. meridionalis are the only AA genome species found within Australia and occur as widespread populations across northern Australia. The chloroplast genome sequence of O. rufipogon from Asia and Australia and O. meridionalis and O. australiensis (an Australian member of the genus very distant from O. sativa) was obtained by massively parallel sequencing and compared with the chloroplast genome sequence of domesticated O. sativa. Oryza australiensis differed in more than 850 sites single nucleotide polymorphism or indel from each of the other samples. The other wild rice species had only around 100 differences relative to cultivated rice. The chloroplast genomes of Australian O. rufipogon and O. meridionalis were closely related with only 32 differences. The Asian O. rufipogon chloroplast genome (with only 68 differences) was closer to O. sativa than the Australian taxa (both with more than 100 differences). The chloroplast sequences emphasize the genetic distinctness of the Australian populations and their potential as a source of novel rice germplasm. The Australian O. rufipogon may be a perennial form of O. meridionalis.

Coalescent simulations reveal hybridization and incomplete lineage sorting in Mediterranean Linaria

We examined the phylogenetic history of Linaria with special emphasis on the Mediterranean sect. Supinae (44 species). We revealed extensive highly supported incongruence among two nuclear (ITS, AGT1) and two plastid regions (rpl32-trnL^UAG, trnS-trnG). Coalescent simulations, a hybrid detection test and species tree inference in *BEAST revealed that incomplete lineage sorting and hybridization may both be responsible for the incongruent pattern observed. Additionally, we present a multilabelled *BEAST species tree as an alternative approach that allows the possibility of observing multiple placements in the species tree for the same taxa. That permitted the incorporation of processes such as hybridization within the tree while not violating the assumptions of the *BEAST model. This methodology is presented as a functional tool to disclose the evolutionary history of species complexes that have experienced both hybridization and incomplete lineage sorting. The drastic climatic events that have occurred in the Mediterranean since the late Miocene, including the Quaternary-type climatic oscillations, may have made both processes highly recurrent in the Mediterranean flora.

The influence of pre- and post-zygotic barriers on interspecific Corymbia hybridization

BACKGROUND AND AIMS

Corymbia species from different sections hybridize readily, with some of increasing economic importance to plantation forestry. This study explores the locations of reproductive barriers between interspecific Corymbia hybrids and investigates the reproductive success of a wide taxonomic range of C. torelliana hybrid crosses.

METHODS

Pollen, pistil and embryo development were investigated for four C. torelliana crosses (×C. torelliana, ×C. citriodora subsp. citriodora, ×C. tessellaris and ×C. intermedia) using fluorescent and standard microscopy to identify the locations of interspecific reproductive isolating barriers. Corymbia torelliana was also crossed with 16 taxa, representing six of the seven Corymbia sections, both Corymbia subgenera and one species each from the related genera, Angophora and Eucalyptus. All crosses were assessed for capsule and seed yields.

KEY RESULTS

Interspecific C. torelliana hybridization was controlled by pre-zygotic reproductive isolating barriers inhibiting pollen adhesion to the stigma, pollen germination, pollen tube growth in the style and pollen tube penetration of the micropyle. Corymbia torelliana (subgenus Blakella, sect. Torellianae) was successfully hybridized with Corymbia species from subgenus Blakella, particularly C. citriodora subsp. citriodora, C. citriodora subsp. variegata, C. henryi (sect. Maculatae) and C. tessellaris (sect. Abbreviatae), and subgenus Corymbia, particularly C. clarksoniana and C. erythrophloia (sect. Septentrionales). Attempted intergeneric hybrids between C. torelliana and either Angophora floribunda or Eucalyptus pellita were unsuccessful.

CONCLUSIONS

Corymbia hybrids were formed between species from different sections and subgenera, but not with species from the related genera Angophora or Eucalyptus. Reproductive isolation between the interspecific Corymbia hybrid crosses was controlled by early- and late-acting pre-zygotic isolating barriers, with reproductive success generally decreasing with increasing taxonomic distance between parent species. These findings support the monophyly of Corymbia and the close relationships of infrageneric clades. The hybridizing propensity of Corymbia species provides opportunities for breeding but suggests risks of environmental gene flow.

Molecular evolution of Adh and LEAFY and the phylogenetic utility of their introns in Pyrus (Rosaceae)

Background

The genus Pyrus belongs to the tribe Pyreae (the former subfamily Maloideae) of the family Rosaceae, and includes one of the most important commercial fruit crops, pear. The phylogeny of Pyrus has not been definitively reconstructed. In our previous efforts, the internal transcribed spacer region (ITS) revealed a poorly resolved phylogeny due to non-concerted evolution of nrDNA arrays. Therefore, introns of low copy nuclear genes (LCNG) are explored here for improved resolution. However, paralogs and lineage sorting are still two challenges for applying LCNGs in phylogenetic studies, and at least two independent nuclear loci should be compared. In this work the second intron of LEAFY and the alcohol dehydrogenase gene (Adh) were selected to investigate their molecular evolution and phylogenetic utility.

Results

DNA sequence analyses revealed a complex ortholog and paralog structure of Adh genes in Pyrus and Malus, the pears and apples. Comparisons between sequences from RT-PCR and genomic PCR indicate that some Adh homologs are putatively nonfunctional. A partial region of Adh1 was sequenced for 18 Pyrus species and three subparalogs representing Adh1-1 were identified. These led to poorly resolved phylogenies due to low sequence divergence and the inclusion of putative recombinants. For the second intron of LEAFY, multiple inparalogs were discovered for both LFY1int2 and LFY2int2. LFY1int2 is inadequate for phylogenetic analysis due to lineage sorting of two inparalogs. LFY2int2-N, however, showed a relatively high sequence divergence and led to the best-resolved phylogeny. This study documents the coexistence of outparalogs and inparalogs, and lineage sorting of these paralogs and orthologous copies. It reveals putative recombinants that can lead to incorrect phylogenetic inferences, and presents an improved phylogenetic resolution of Pyrus using LFY2int2-N.

Conclusions

Our study represents the first phylogenetic analyses based on LCNGs in Pyrus. Ancient and recent duplications lead to a complex structure of Adh outparalogs and inparalogs in Pyrus and Malus, resulting in neofunctionalization, nonfunctionalization and possible subfunctionalization. Among all investigated orthologs, LFY2int2-N is the best nuclear marker for phylogenetic reconstruction of Pyrus due to suitable sequence divergence and the absence of lineage sorting.

Comparative and phylogenomic analyses of cinnamoyl-CoA reductase and cinnamoyl-CoA-reductase-like gene family in land plants

The biosynthesis of monolignols, the main components of lignin, involves many intermediates and enzymes. The cinnamoyl-CoA reductase (CCR) enzyme catalyzes the conversion of cinnamoyl-CoAs to cinnamaldehydes, i.e. the first specific step in lignin synthesis. The CCR and CCR-like gene family was studied partially in several plant species. This is a comprehensive study of the CCR and CCR-like gene family including genome organization, gene structure, phylogeny across land plant species, and, expression profiling in Populus. Analysis of amino acid motifs enabled the identification of sequence variations in the CCR catalytic site and annotates CCR and CCR-like genes. CCR and CCR-like genes were distributed in three major phylogenetic classes of which one includes the bona fide CCR genes. The other two classes include CCR and CCR-like, of which several genes present a high similarity to cinnamyl alcohol dehydrogenase, or dihydroflavonol reductase (DFR) genes. All CCR, CCR-like, and DFR classes were deeply rooted in the phylogeny of land plants suggesting that their evolution preceded the evolution of lycophytes. Over two thirds of CCR and CCR-like Populus genes were physically distributed on duplicated regions. This suggests that these duplication/retention processes contributed significantly to the size of the CCR and CCR-like gene family. The Populus CCR and CCR-like genes showed six expression patterns in the tissues studied with a preferential expression of PoptrCCR12 in xylem. The other genes present divergent expression profiles with some preferentially expressed in leaves, bark, or both. Several CCR and CCR-like genes were induced or repressed under various abiotic stresses suggesting that their duplication was followed by the evolution of divergent expression profiles and divergence of functions.

Population genetic analysis and phylogeny reconstruction in Eucalyptus (Myrtaceae) using high-throughput, genome-wide genotyping

A set of over 8000 Diversity Arrays Technology (DArT) markers was tested for its utility in high-resolution population and phylogenetic studies across a range of Eucalyptus taxa. Small-scale population studies of Eucalyptus camaldulensis, Eucalyptus cladocalyx, Eucalyptus globulus, Eucalyptus grandis, Eucalyptus nitens, Eucalyptus pilularis and Eucalyptus urophylla demonstrated the potential of genome-wide genotyping with DArT markers to differentiate species, to identify interspecific hybrids and to resolve biogeographic disjunctions within species. The population genetic studies resolved geographically partitioned clusters in E. camaldulensis, E. cladocalyx, E. globulus and E. urophylla that were congruent with previous molecular studies. A phylogenetic study of 94 eucalypt species provided results that were largely congruent with traditional taxonomy and ITS-based phylogenies, but provided more resolution within major clades than had been obtained previously. Ascertainment bias (the bias introduced in a phylogeny from using markers developed in a small sample of the taxa that are being studied) was not detected. DArT offers an unprecedented level of resolution for population genetic, phylogenetic and evolutionary studies across the full range of Eucalyptus species.

Origins and introgression of polyploid species in Mentzelia section Trachyphytum (Loasaceae)

PREMISE OF THE STUDY

Polyploid speciation has been important in plant evolution. However, the conditions that favor the origination and persistence of polyploids are still not well understood. Here, we examine origins of 16 polyploid species in Mentzelia section Trachyphytum. •

METHODS

We used phylogeny reconstructions based on DNA sequences from plastid regions and the nuclear gene isocitrate dehydrogenase (idh) to construct hypotheses of introgression and polyploidization. •

KEY RESULTS

Molecular data suggest that homoploid hybridization has been surprisingly common in Trachyphytum. Diploid species had unequal involvement in polyploid origins, but most polyploid taxa had allopolyploid origins from extant progenitors. A few polyploids with extreme phenotypes did not appear to have extant progenitors. We infer that the progenitors of these species were derived from extinct diploid lineages or ancestral lineages of multiple extant diploids. In agreement with other recent studies, we recovered molecular evidence of multiple phylogenetically distinct origins for several polyploid taxa, including the widespread octoploid M. albicaulis. •

CONCLUSIONS

Evidence of high levels of introgression and allopolyploidy suggests that hybridization has played an important role in the evolution of Trachyphytum. Although idh sequences exhibited complicated evolution, including gene duplication, deletion, and recombination, they provided a higher percentage of informative characters for phylogeny reconstruction than the most variable plastid regions, allowing tests of hypotheses regarding polyploid origins. Given the necessity for rapidly evolving low-copy nuclear genes, researchers studying hybridization and polyploidy may increasingly turn to complex sequence data.

Reticulate or tree-like chloroplast DNA evolution in Sileneae (Caryophyllaceae)?

Despite sampling of up to 25kb of chloroplast DNA sequence from 24 species in Sileneae a number of nodes in the phylogeny remain poorly supported and it is not expected that additional sequence sampling will converge to a reliable phylogenetic hypothesis in these parts of the tree. The main reason for this is probably a combination of rapid radiation and substitution rate heterogeneity. Poor resolution among closely related species are often explained by low levels of variation in chloroplast data, but the problem with our data appear to be high levels of homoplasy. Tree-like cpDNA evolution cannot be rejected, but apparent incongruent patterns between different regions are evaluated with the possibility of ancient interspecific chloroplast recombination as explanatory model. However, several major phylogenetic relationships, previously not recognized, are confidently resolved, e.g. the grouping of the two SW Anatolian taxa S. cryptoneura and S. sordida strongly disagrees with previous studies on nuclear DNA sequence data, and indicate a possible case of homoploid hybrid origin. The closely related S. atocioides and S. aegyptiaca form a sister group to Lychnis and the rest of Silene, thus suggesting that Silene may be paraphyletic, despite recent revisions based on molecular data.

An AFLP marker approach to lower-level systematics in Eucalyptus (Myrtaceae)

Genus Eucalyptus, with over 700 species, presents a number of systematic difficulties including taxa that hybridize or intergrade across environmental gradients. To date, no DNA marker has been found capable of resolving phylogeny below the sectional level in the major subgenera. Molecular markers are needed to support taxonomic revision, assess the extent of genetic divergence at lower taxonomic levels, and inform conservation efforts. We examined the utility of 930 amplified fragment length polymorphisms (AFLPs) for analyzing relationships among Tasmanian taxa of subgenus Symphyomyrtus section Maidenaria. Phenetic and cladistic analyses resolved species into clusters demonstrating significant genetic partitioning, largely concordant with series defined in the most recent taxonomic revision of Eucalyptus. Some departures from current taxonomy were noted, indicating possible cases of morphological convergence and character reversion. Although the resolution obtained using AFLP was greatly superior to that of single sequence markers, the data demonstrated high homoplasy and incomplete resolution of closely related species. The results of this study and others are consistent with recent speciation and reticulate evolution in Maidenaria. We conclude that a combination of phylogenetic and population genetic approaches using multiple molecular markers offers the best prospects for understanding taxonomic relationships below the sectional level in Eucalyptus.

Robust inference of positive selection from recombining coding sequences

MOTIVATION

Accurate detection of positive Darwinian selection can provide important insights to researchers investigating the evolution of pathogens. However, many pathogens (particularly viruses) undergo frequent recombination and the phylogenetic methods commonly applied to detect positive selection have been shown to give misleading results when applied to recombining sequences. We propose a method that makes maximum likelihood inference of positive selection robust to the presence of recombination. This is achieved by allowing tree topologies and branch lengths to change across detected recombination breakpoints. Further improvements are obtained by allowing synonymous substitution rates to vary across sites.

RESULTS

Using simulation we show that, even for extreme cases where recombination causes standard methods to reach false positive rates >90%, the proposed method decreases the false positive rate to acceptable levels while retaining high power. We applied the method to two HIV-1 datasets for which we have previously found that inference of positive selection is invalid owing to high rates of recombination. In one of these (env gene) we still detected positive selection using the proposed method, while in the other (gag gene) we found no significant evidence of positive selection.

AVAILABILITY

A HyPhy batch language implementation of the proposed methods and the HIV-1 datasets analysed are available at http://www.cbio.uct.ac.za/pub_support/bioinf06. The HyPhy package is available at http://www.hyphy.org, and it is planned that the proposed methods will be included in the next distribution. RDP2 is available at http://darwin.uvigo.es/rdp/rdp.html