Polyploidy on islands - concerted evolution and gene loss amid chromosomal stasis.

BACKGROUND AND AIMS: Polyploidy is an important process that often generates genomic diversity within lineages, but it can also cause changes that result in loss of genomic material. Island lineages, while often polyploid, typically show chromosomal stasis but have not been investigated in detail regarding smaller-scale gene loss. Our aim was to investigate post-polyploidization genome dynamics in a chromosomally stable lineage of Malvaceae endemic to New Zealand. METHODS: We determined chromosome numbers and used fluorescence in situ hybridization to localize 18S and 5S rDNA. Gene sequencing of 18S rDNA, the internal transcribed spacers (ITS) with intervening 5.8S rDNA, and a low-copy nuclear gene, GBSSI-1, was undertaken to determine if gene loss occurred in the New Zealand lineage following polyploidy. KEY RESULTS: The chromosome number for all species investigated was 2n = 42, with the first published report for the monotypic Australian genus Asterotrichion. The five species investigated all had two 5S rDNA signals localized interstitially on the long arm of one of the largest chromosome pairs. All species, except Plagianthus regius, had two 18S rDNA signals localized proximally on the short arm of one of the smallest chromosome pairs. Plagianthus regius had two additional 18S rDNA signals on a separate chromosome, giving a total of four. Sequencing of nuclear ribosomal 18S rDNA and the ITS cistron indicated loss of historical ribosomal repeats. Phylogenetic analysis of a low-copy nuclear gene, GBSSI-1, indicated that some lineages maintained three copies of the locus, while others have lost one or two copies. CONCLUSIONS: Although island endemic lineages show chromosomal stasis, with no additional changes in chromosome number, they may undergo smaller-scale processes of gene loss and concerted evolution ultimately leading to further genome restructuring and downsizing.

A Malvaceae mystery: A mallow maelstrom of genome multiplications and maybe misleading methods?

Previous research suggests that Gossypium has undergone a 5- to 6-fold multiplication following its divergence from Theobroma. However, the number of events, or where they occurred in the Malvaceae phylogeny remains unknown. We analyzed transcriptomic and genomic data from representatives of eight of the nine Malvaceae subfamilies. Phylogenetic analysis of nuclear data placed Dombeya (Dombeyoideae) as sister to the rest of Malvadendrina clade, but the plastid DNA tree strongly supported Durio (Helicteroideae) in this position. Intraspecific Ks plots indicated that all sampled taxa, except Theobroma (Byttnerioideae), Corchorus (Grewioideae), and Dombeya (Dombeyoideae), have experienced whole genome multiplications (WGMs). Quartet analysis suggested WGMs were shared by Malvoideae-Bombacoideae and Sterculioideae-Tilioideae, but did not resolve whether these are shared with each other or Helicteroideae (Durio). Gene tree reconciliation and Bayesian concordance analysis suggested a complex history. Alternative hypotheses are suggested, each involving two independent autotetraploid and one allopolyploid event. They differ in that one entails an allopolyploid origin for the Durio lineage, whereas the other invokes an allopolyploid origin for Malvoideae-Bombacoideae. We highlight the need for more genomic information in the Malvaceae and improved methods to resolve complex evolutionary histories that may include allopolyploidy, incomplete lineage sorting, and variable rates of gene and genome evolution.

The compounding effects of high pollen limitation, selfing rates and inbreeding depression leave a New Zealand tree with few viable offspring.

BACKGROUND AND AIMS: Interactions between species are especially sensitive to environmental changes. The interaction between plants and pollinators is of particular interest given the potential current global decline in pollinators. Reduced pollinator services can be compensated for in some plant species by self-pollination. However, if inbreeding depression is high, selfed progeny could die prior to reaching adulthood, leading to cryptic recruitment failure. METHODS: To examine this scenario, pollinator abundance, pollen limitation, selfing rates and inbreeding depression were examined in 12 populations of varying disturbance levels in Sophora microphylla (Fabaceae), an endemic New Zealand tree species. KEY RESULTS: High pollen limitation was found in all populations (average of 58 % reduction in seed production, nine populations), together with high selfing rates (61 % of offspring selfed, six populations) and high inbreeding depression (selfed offspring 86 % less fit, six populations). Pollen limitation was associated with lower visitation rates by the two endemic bird pollinators. CONCLUSIONS: The results suggest that for these populations, over half of the seeds produced are genetically doomed. This reduction in the fitness of progeny due to reduced pollinator service is probably important to population dynamics of other New Zealand species. More broadly, the results suggest that measures of seed production or seedling densities may be a gross overestimate of the effective offspring production. This could lead to cryptic recruitment failure, i.e. a decline in successful reproduction despite high progeny production. Given the global extent of pollinator declines, cryptic recruitment failure may be widespread.

Gene silencing via DNA methylation in naturally occurring Tragopogon miscellus (Asteraceae) allopolyploids.

BACKGROUND: Hybridization coupled with whole-genome duplication (allopolyploidy) leads to a variety of genetic and epigenetic modifications in the resultant merged genomes. In particular, gene loss and gene silencing are commonly observed post-polyploidization. Here, we investigated DNA methylation as a potential mechanism for gene silencing in Tragopogon miscellus (Asteraceae), a recent and recurrently formed allopolyploid. This species, which also exhibits extensive gene loss, was formed from the diploids T. dubius and T. pratensis. RESULTS: Comparative bisulfite sequencing revealed CG methylation of parental homeologs for three loci (S2, S18 and TDF-44) that were previously identified as silenced in T. miscellus individuals relative to the diploid progenitors. One other locus (S3) examined did not show methylation, indicating that other transcriptional and post-transcriptional mechanisms are likely responsible for silencing that homeologous locus. CONCLUSIONS: These results indicate that Tragopogon miscellus allopolyploids employ diverse mechanisms, including DNA methylation, to respond to the potential shock of genome merger and doubling.

A roadmap of phylogenomic methods for studying polyploid plant genera.

Phylogenetic inference of polyploid species is the first step towards understanding their patterns of diversification. In this paper, we review the challenges and limitations of inferring species relationships of polyploid plants using traditional phylogenetic sequencing approaches, as well as the mischaracterization of the species tree from single or multiple gene trees. We provide a roadmap to infer interspecific relationships among polyploid lineages by comparing and evaluating the application of current phylogenetic, phylogenomic, transcriptomic, and whole-genome approaches using different sequencing platforms. For polyploid species tree reconstruction, we assess the following criteria: (1) the amount of prior information or tools required to capture the genetic region(s) of interest; (2) the probability of recovering homeologs for polyploid species; and (3) the time efficiency of downstream data analysis. Moreover, we discuss bioinformatic pipelines that can reconstruct networks of polyploid species relationships. In summary, although current phylogenomic approaches have improved our understanding of reticulate species relationships in polyploid-rich genera, the difficulties of recovering reliable orthologous genes and sorting all homeologous copies for allopolyploids remain a challenge. In the future, assembled long-read sequencing data will assist the recovery and identification of multiple gene copies, which can be particularly useful for reconstructing the multiple independent origins of polyploids.

Ploidy variation in Rhododendron subsection Maddenia and its implications for conservation.

Polyploidy, which is common in plants, can confound taxon recognition and hence conservation assessments. In the taxonomically complex genus Rhododendron, 25 % of the over 1,300 taxa are considered under threat and 27 % Near Threatened or Data Deficient, with their taxonomy needing to be resolved urgently. Although ploidy levels of Rhododendron taxa range from diploid (2x) to dodecaploid (12x) according to previous reports, the extent of polyploidy across the genus has not been examined. We first summarized the taxonomic distribution of polyploids in the genus based on the literature. Then as a case study, we estimated ploidy levels of 47 taxa in subsection Maddenia (subgenus Rhododendron, section Rhododendron) using flow cytometry, together with verification of meiotic chromosome counts for representative taxa. The summary of reported ploidy in Rhododendron indicates that polyploidy is most common in subgenera Pentanthera and Rhododendron. In subsection Maddenia, all examined taxa are diploids except for the R. maddenii complex that shows a high ploidy variation (2-8x, 12x). We investigated ploidy level of 12 taxa in subsection Maddenia for the first time, and estimated genome sizes of two Rhododendron species. Knowledge of ploidy levels will inform phylogenetic analysis of unresolved species complexes. Overall, our study of subsection Maddenia provides a model for examining multiple issues including taxonomic complexity, ploidy variation and geographic distribution in relation to biodiversity conservation.

Microsatellite markers for the endangered root holoparasite Dactylanthus taylorii (Balanophoraceae) from 454 pyrosequencing.

PREMISE OF THE STUDY: Microsatellite loci were isolated and developed as polymorphic markers for the New Zealand endemic root holoparasite Dactylanthus taylorii for use in population and conservation genetics studies. METHODS AND RESULTS: Shotgun 454 pyrosequencing was performed on genomic DNA pooled from three individuals of D. taylorii. From 61709 individual sequence reads, primers for 753 microsatellite loci were developed in silico and 72 of these were tested for consistent amplification and variability. Ten microsatellite loci were found to be polymorphic and consistently scorable when screened in 44 individuals from five geographically distant populations. The number of alleles per locus ranged from four to 16 with an average of 9.7, and average observed heterozygosity per locus was between 0.182 and 0.634. CONCLUSIONS: These polymorphic microsatellite markers establish an important resource for ongoing conservation initiatives and planned population genetic studies of D. taylorii.

Outlier analyses and genome-wide association study identify glgC and ERD6-like 4 as candidate genes for foliar water-soluble carbohydrate accumulation in Trifolium repens.

Increasing water-soluble carbohydrate (WSC) content in white clover is important for improving nutritional quality and reducing environmental impacts from pastoral agriculture. Elucidation of genes responsible for foliar WSC variation would enhance genetic improvement by enabling molecular breeding approaches. The aim of the present study was to identify single nucleotide polymorphisms (SNPs) associated with variation in foliar WSC in white clover. A set of 935 white clover individuals, randomly sampled from five breeding pools selectively bred for divergent (low or high) WSC content, were assessed with 14,743 genotyping-by-sequencing SNPs, using three outlier detection methods: PCAdapt, BayeScan and KGD-FST. These analyses identified 33 SNPs as discriminating between high and low WSC populations and putatively under selection. One SNP was located in the intron of ERD6-like 4, a gene coding for a sugar transporter located on the vacuole membrane. A genome-wide association study using a subset of 605 white clover individuals and 5,757 SNPs, identified a further 12 SNPs, one of which was associated with a starch biosynthesis gene, glucose-1-phosphate adenylyltransferase, glgC. Our results provide insight into genomic regions underlying WSC accumulation in white clover, identify candidate genomic regions for further functional validation studies, and reveal valuable information for marker-assisted or genomic selection in white clover.

Polyploidy on Islands: Its Emergence and Importance for Diversification.

Whole genome duplication or polyploidy is widespread among floras globally, but traditionally has been thought to have played a minor role in the evolution of island biodiversity, based on the low proportion of polyploid taxa present. We investigate five island systems (Juan Fernández, Galápagos, Canary Islands, Hawaiian Islands, and New Zealand) to test whether polyploidy (i) enhances or hinders diversification on islands and (ii) is an intrinsic feature of a lineage or an attribute that emerges in island environments. These island systems are diverse in their origins, geographic and latitudinal distributions, levels of plant species endemism (37% in the Galapagos to 88% in the Hawaiian Islands), and ploidy levels, and taken together are representative of islands more generally. We compiled data for vascular plants and summarized information for each genus on each island system, including the total number of species (native and endemic), generic endemicity, chromosome numbers, genome size, and ploidy levels. Dated phylogenies were used to infer lineage age, number of colonization events, and change in ploidy level relative to the non-island sister lineage. Using phylogenetic path analysis, we then tested how the diversification of endemic lineages varied with the direct and indirect effects of polyploidy (presence of polyploidy, time on island, polyploidization near colonization, colonizer pool size) and other lineage traits not associated with polyploidy (time on island, colonizer pool size, repeat colonization). Diploid and tetraploid were the most common ploidy levels across all islands, with the highest ploidy levels (>8x) recorded for the Canary Islands (12x) and New Zealand (20x). Overall, we found that endemic diversification of our focal island floras was shaped by polyploidy in many cases and certainly others still to be detected considering the lack of data in many lineages. Polyploid speciation on the islands was enhanced by a larger source of potential congeneric colonists and a change in ploidy level compared to overseas sister taxa.

New targets acquired: Improving locus recovery from the Angiosperms353 probe set.

PREMISE: Universal target enrichment kits maximize utility across wide evolutionary breadth while minimizing the number of baits required to create a cost-efficient kit. The Angiosperms353 kit has been successfully used to capture loci throughout the angiosperms, but the default target reference file includes sequence information from only 6-18 taxa per locus. Consequently, reads sequenced from on-target DNA molecules may fail to map to references, resulting in fewer on-target reads for assembly, and reducing locus recovery. METHODS: We expanded the Angiosperms353 target file, incorporating sequences from 566 transcriptomes to produce a 'mega353' target file, with each locus represented by 17-373 taxa. This mega353 file is a drop-in replacement for the original Angiosperms353 file in HybPiper analyses. We provide tools to subsample the file based on user-selected taxon groups, and to incorporate other transcriptome or protein-coding gene data sets. RESULTS: Compared to the default Angiosperms353 file, the mega353 file increased the percentage of on-target reads by an average of 32%, increased locus recovery at 75% length by 49%, and increased the total length of the concatenated loci by 29%. DISCUSSION: Increasing the phylogenetic density of the target reference file results in improved recovery of target capture loci. The mega353 file and associated scripts are available at: https://github.com/chrisjackson-pellicle/NewTargets.

Similar patterns of rDNA evolution in synthetic and recently formed natural populations of Tragopogon (Asteraceae) allotetraploids.

BACKGROUND: Tragopogon mirus and T. miscellus are allotetraploids (2n = 24) that formed repeatedly during the past 80 years in eastern Washington and adjacent Idaho (USA) following the introduction of the diploids T. dubius, T. porrifolius, and T. pratensis (2n = 12) from Europe. In most natural populations of T. mirus and T. miscellus, there are far fewer 35S rRNA genes (rDNA) of T. dubius than there are of the other diploid parent (T. porrifolius or T. pratensis). We studied the inheritance of parental rDNA loci in allotetraploids resynthesized from diploid accessions. We investigate the dynamics and directionality of these rDNA losses, as well as the contribution of gene copy number variation in the parental diploids to rDNA variation in the derived tetraploids. RESULTS: Using Southern blot hybridization and fluorescent in situ hybridization (FISH), we analyzed copy numbers and distribution of these highly reiterated genes in seven lines of synthetic T. mirus (110 individuals) and four lines of synthetic T. miscellus (71 individuals). Variation among diploid parents accounted for most of the observed gene imbalances detected in F1 hybrids but cannot explain frequent deviations from repeat additivity seen in the allotetraploid lines. Polyploid lineages involving the same diploid parents differed in rDNA genotype, indicating that conditions immediately following genome doubling are crucial for rDNA changes. About 19% of the resynthesized allotetraploid individuals had equal rDNA contributions from the diploid parents, 74% were skewed towards either T. porrifolius or T. pratensis-type units, and only 7% had more rDNA copies of T. dubius-origin compared to the other two parents. Similar genotype frequencies were observed among natural populations. Despite directional reduction of units, the additivity of 35S rDNA locus number is maintained in 82% of the synthetic lines and in all natural allotetraploids. CONCLUSIONS: Uniparental reductions of homeologous rRNA gene copies occurred in both synthetic and natural populations of Tragopogon allopolyploids. The extent of these rDNA changes was generally higher in natural populations than in the synthetic lines. We hypothesize that locus-specific and chromosomal changes in early generations of allopolyploids may influence patterns of rDNA evolution in later generations.

On the road to diploidization? Homoeolog loss in independently formed populations of the allopolyploid Tragopogon miscellus (Asteraceae).

BACKGROUND: Polyploidy (whole-genome duplication) is an important speciation mechanism, particularly in plants. Gene loss, silencing, and the formation of novel gene complexes are some of the consequences that the new polyploid genome may experience. Despite the recurrent nature of polyploidy, little is known about the genomic outcome of independent polyploidization events. Here, we analyze the fate of genes duplicated by polyploidy (homoeologs) in multiple individuals from ten natural populations of Tragopogon miscellus (Asteraceae), all of which formed independently from T. dubius and T. pratensis less than 80 years ago. RESULTS: Of the 13 loci analyzed in 84 T. miscellus individuals, 11 showed loss of at least one parental homoeolog in the young allopolyploids. Two loci were retained in duplicate for all polyploid individuals included in this study. Nearly half (48%) of the individuals examined lost a homoeolog of at least one locus, with several individuals showing loss at more than one locus. Patterns of loss were stochastic among individuals from the independently formed populations, except that the T. dubius copy was lost twice as often as T. pratensis. CONCLUSION: This study represents the most extensive survey of the fate of genes duplicated by allopolyploidy in individuals from natural populations. Our results indicate that the road to genome downsizing and ultimate genetic diploidization may occur quickly through homoeolog loss, but with some genes consistently maintained as duplicates. Other genes consistently show evidence of homoeolog loss, suggesting repetitive aspects to polyploid genome evolution.

Concerted evolution of rDNA in recently formed Tragopogon allotetraploids is typically associated with an inverse correlation between gene copy number and expression.

We analyzed nuclear ribosomal DNA (rDNA) transcription and chromatin condensation in individuals from several populations of Tragopogon mirus and T. miscellus, allotetraploids that have formed repeatedly within only the last 80 years from T. dubius and T. porrifolius and T. dubius and T. pratensis, respectively. We identified populations with no (2), partial (2), and complete (4) nucleolar dominance. It is probable that epigenetic regulation following allopolyploidization varies between populations, with a tendency toward nucleolar dominance by one parental homeologue. Dominant rDNA loci are largely decondensed at interphase while silent loci formed condensed heterochromatic regions excluded from nucleoli. Those populations where nucleolar dominance is fixed are epigenetically more stable than those with partial or incomplete dominance. Previous studies indicated that concerted evolution has partially homogenized thousands of parental rDNA units typically reducing the copy numbers of those derived from the T. dubius diploid parent. Paradoxically, despite their low copy number, repeats of T. dubius origin dominate rDNA transcription in most populations studied, i.e., rDNA units that are genetic losers (copy numbers) are epigenetic winners (high expression).

Microsatellite markers for Corybas (Orchidaceae) species in New Zealand.

PREMISE OF THE STUDY: Microsatellite markers were developed for New Zealand species of Corybas (Orchidaceae) to investigate population genetics and species delimitation. METHODS AND RESULTS: From sequencing a total genomic DNA library (using Illumina MiSeq), we developed 22 microsatellite markers for C. obscurus. The di- and trinucleotide repeat loci were initially trialed on individuals representing seven Corybas taxa (C. "rimutaka," C. confusus, C. hypogaeus, C. macranthus, C. obscurus, C. trilobus, and C. walliae) and had one to eight alleles per locus. Twelve polymorphic markers were further tested on six Corybas populations from three of the seven taxa (C. obscurus, C. "rimutaka," and C. trilobus). Observed and expected heterozygosities ranged from 0-1 and 0-0.859, respectively. The utility of these 12 loci was further validated in five related Corybas species (C. hypogaeus, C. obscurus, C. vitreus, C. walliae, and C. "rimutaka"; 38 individuals) representing populations from across the North and South Islands. The average value for genetic diversity among populations (F ST) of 0.439 shows differentiation among species. CONCLUSIONS: These markers will be useful for future studies aimed at delimiting species boundaries and examining the genetic diversity of the New Zealand Corybas species.

Evolution and expression of homeologous loci in Tragopogon miscellus (Asteraceae), a recent and reciprocally formed allopolyploid.

On both recent and ancient time scales, polyploidy (genome doubling) has been a significant evolutionary force in plants. Here, we examined multiple individuals from reciprocally formed populations of Tragopogon miscellus, an allotetraploid that formed repeatedly within the last 80 years from the diploids T. dubius and T. pratensis. Using cDNA-AFLPs followed by genomic and cDNA cleaved amplified polymorphic sequence (CAPS) analyses, we found differences in the evolution and expression of homeologous loci in T. miscellus. Fragment variation within T. miscellus, possibly attributable to reciprocal formation, comprised 0.6% of the cDNA-AFLP bands. Genomic and cDNA CAPS analyses of 10 candidate genes revealed that only one "transcript-derived fragment" (TDF44) showed differential expression of parental homeologs in T. miscellus; the T. pratensis homeolog was preferentially expressed by most polyploids in both populations. Most of the cDNA-AFLP polymorphisms apparently resulted from loss of parental fragments in the polyploids. Importantly, changes at the genomic level have occurred stochastically among individuals within the independently formed populations. Synthetic F(1) hybrids between putative diploid progenitors are additive of their parental genomes, suggesting that polyploidization rather than hybridization induces genomic changes in Tragopogon.

Cytonuclear Coordination Is Not Immediate upon Allopolyploid Formation in Tragopogon miscellus (Asteraceae) Allopolyploids.

Allopolyploids, formed by hybridization and chromosome doubling, face the immediate challenge of having duplicated nuclear genomes that interact with the haploid and maternally inherited cytoplasmic (plastid and mitochondrial) genomes. Most of our knowledge of the genomic consequences of allopolyploidy has focused on the fate of the duplicated nuclear genes without regard to their potential interactions with cytoplasmic genomes. As a step toward understanding the fates of nuclear-encoded subunits that are plastid-targeted, here we examine the retention and expression of the gene encoding the small subunit of Ribulose-1, 5-bisphosphate carboxylase/oxygenase (Rubisco; rbcS) in multiple populations of allotetraploid Tragopogon miscellus (Asteraceae). These polyploids formed recently (~80 years ago) and repeatedly from T. dubius and T. pratensis in the northwestern United States. Examination of 79 T. miscellus individuals from 10 natural populations, as well as 25 synthetic allotetraploids, including reciprocally formed plants, revealed a low percentage of naturally occurring individuals that show a bias in either gene (homeolog) loss (12%) or expression (16%), usually toward maintaining the maternal nuclear copy of rbcS. For individuals showing loss, seven retained the maternally derived rbcS homeolog only, while three had the paternally derived copy. All of the synthetic polyploid individuals examined (S0 and S1 generations) retained and expressed both parental homeologs. These results demonstrate that cytonuclear coordination does not happen immediately upon polyploid formation in Tragopogon miscellus.

Microsatellite markers for the New Zealand endemic Myosotis pygmaea species group (Boraginaceae) amplify across species.

PREMISE OF THE STUDY: Microsatellite loci were developed as polymorphic markers for the New Zealand endemic Myosotis pygmaea species group (Boraginaceae) for use in species delimitation and population and conservation genetic studies. METHODS AND RESULTS: Illumina MiSeq sequencing was performed on genomic DNA from seedlings of M. drucei. From trimmed paired-end sequences >400 bp, 484 microsatellite loci were identified. Twelve of 48 microsatellite loci tested were found to be polymorphic and consistently scorable when screened on 53 individuals from four populations representing the geographic range of M. drucei. They also amplify in all other species in the M. pygmaea species group, i.e., M. antarctica, M. brevis, M. glauca, and M. pygmaea, as well as 18 other Myosotis species. CONCLUSIONS: These 12 polymorphic microsatellite markers establish an important resource for research and conservation of the M. pygmaea species group and potentially other Southern Hemisphere Myosotis.

Tumor cell targeting by iron oxide nanoparticles is dominated by different factors in vitro versus in vivo.

Realizing the full potential of iron oxide nanoparticles (IONP) for cancer diagnosis and therapy requires selective tumor cell accumulation. Here, we report a systematic analysis of two key determinants for IONP homing to human breast cancers: (i) particle size and (ii) active vs passive targeting. In vitro, molecular targeting to the HER2 receptor was the dominant factor driving cancer cell association. In contrast, size was found to be the key determinant of tumor accumulation in vivo, where molecular targeting increased tumor tissue concentrations for 30 nm but not 100 nm IONP. Similar to the in vitro results, PEGylation did not influence in vivo IONP biodistribution. Thus, the results reported here indicate that the in vitro advantages of molecular targeting may not consistently extend to pre-clinical in vivo settings. These observations may have important implications for the design and clinical translation of advanced, multifunctional, IONP platforms.