Phylogeographic and demographic modeling analyses of the multiple origins of the rheophytic goldenrod Solidago yokusaiana Makino.

Understanding adaptation mechanisms is important in evolutionary biology. Parallel adaptation provides good opportunities to investigate adaptive evolution. To confirm parallel adaptation, it is effective to examine whether the phenotypic similarity has one or multiple origins and to use demographic modeling to consider the gene flow between ecotypes. Solidago yokusaiana is a rheophyte endemic to the Japanese Archipelago that diverged from Solidago virgaurea. This study examined the parallel origins of S. yokusaiana by distinguishing between multiple and single origins and subsequent gene flow. The haplotypes of noncoding chloroplast DNA and genotypes at 14 nuclear simple sequence repeat (nSSR) loci and single-nucleotide polymorphisms (SNPs) revealed by double-digest restriction-associated DNA sequencing (ddRADseq) were used for phylogeographic analysis; the SNPs were also used to model population demographics. Some chloroplast haplotypes were common to S. yokusaiana and its ancestor S. virgaurea. Also, the population genetic structures revealed by nSSR and SNPs did not correspond to the taxonomic species. The demographic modeling supported the multiple origins of S. yokusaiana in at least four districts and rejected a single origin with ongoing gene flow between the two species, implying that S. yokusaiana independently and repeatedly adapted to frequently flooding riversides.

Genetic consequences of being a dwarf: do evolutionary changes in life-history traits influence gene flow patterns in populations of the world's smallest goldenrod?

BACKGROUND AND AIMS: Contrasting life-history traits can evolve through generations of dwarf plant ecotypes, yet such phenotypic changes often involve decreased plant size and reproductive allocation, which can configure seed dispersal patterns and, subsequently, population demography. Therefore, evolutionary transitions to dwarfism can represent good study systems to test the roles of life-history traits in population demography by comparing genetic structure between related but phenotypically divergent ecotypes. METHODS: In this study, we examined an ecotypic taxon pair of the world's smallest goldenrod (stem height 2.6 cm) in alpine habitats and its closely related lowland taxon (30-40 cm) found on Yakushima Island, Japan. Genetic variation in chloroplast DNA sequences, nuclear microsatellites and genome-wide single-nucleotide polymorphisms were used to investigate 197 samples from 16 populations, to infer the population genetic demography and compare local genetic structure of the ecotypes. KEY RESULTS: We found a pronounced level of genetic differentiation among alpine dwarf populations, which were much less geographically isolated than their lowland counterparts. In particular, several neighbouring dwarf populations (located ~500 m apart) harboured completely different sets of chloroplast haplotypes and nuclear genetic clusters. Demographic modelling revealed that the dwarf populations have not exchanged genes at significant levels after population divergence. CONCLUSIONS: These lines of evidence suggest that substantial effects of genetic drift have operated on these dwarf populations. The low-growing stature and reduced fecundity (only 3.1 heads per plant) of the dwarf plants may have reduced gene flow and rare long-distance seed dispersal among habitat patches, although the effects of life-history traits require further evaluation using ecological approaches.

Temperature-related cline in the root mass fraction in East Asian wild radish along the Japanese archipelago.

Wild plants with a wide distribution, including those exposed to a wide variety of environmental conditions, may have variations in key functional traits relevant for agricultural applications. The East Asian wild radish (Raphanus sativus var. raphanistroides) is an appropriate model plant because it is widely distributed and has outstanding sink capacity as well as two cultivars within the species. Multiple common garden trials with 14 populations and three testing sites were conducted across the Japanese archipelago to quantify variations in yield and allocation. Significant inter-population variations and interaction effects with testing sites were detected for the root and shoot mass and the root mass fraction (RMF). While the rank order of the population changed drastically among sites and the variance components of genetic effects were small in yield traits (2.4%-4.7%), RMF displayed a large genetic variance (23.2%) and was consistently higher in the northern populations at all sites. Analyses revealed that the mean temperature of growing season of the seed origin was the most prominent factor explaining variation in RMF, irrespective of the sites. We concluded that the trait of resource allocation had a temperature-related cline and plants in cooler climates could invest more resources into their roots.

Phylogeographic analysis suggests two origins for the riparian azalea Rhododendron indicum (L.) Sweet.

Ecological speciation is an important factor in the diversification of plants. The distribution of the woody species Rhododendron indicum, which grows along rivers and is able to withstand water flow when rivers flood (i.e. it is a rheophyte), is disjunct, in contrast to the widespread distribution of its relative, Rhododendron kaempferi. This study aimed to elucidate the phylogenetic relationships between R. indicum and R. kaempferi and the evolutionary processes that gave rise to them. The sequences of three non-coding chloroplast DNA regions (total length 1977 bp) were obtained from 21 populations covering the ranges of the two species. In addition, genome-wide SNPs were genotyped from 20 populations using a genotyping by sequencing method. Leaf morphologies were measured for eight representative populations. Two chloroplast DNA haplotypes, which were detected in R. indicum, were shared between the two species. Genome-wide SNPs identified two lineages in R. indicum and these lineages did not constitute a monophyletic group. Each of these two lineages was related to geographically close populations of R. kaempferi. Leaf morphology, which is a characteristic feature in rheophytes, was not differentiated between the two lineages in R. indicum. The morphological similarity between the two heterogeneous lineages may be a result of parallel evolution from R. kaempferi or of introgressive hybridization between the species due to strong selective pressure imposed by flooding.

Phylogeographic analysis of the East Asian goldenrod (Solidago virgaurea complex, Asteraceae) reveals hidden ecological diversification with recurrent formation of ecotypes.

Background and Aims: The processes and mechanisms underlying lineage diversification are major topics in evolutionary biology. Eurasian goldenrod species of the Solidago virgaurea complex show remarkable morphological and ecological diversity in the Japanese Archipelago, with ecotypic taxa well adapted to specific environments (climate, edaphic conditions and disturbance regimes). The species complex is a suitable model to investigate the evolutionary processes of actively speciating plant groups, due to its ability to evolve in relation to environmental adaptation and its historical population dynamics. Methods: Two chloroplast markers, 18 nuclear microsatellite markers and ddRAD-sequencing were used to infer population genetic demography of S. virgaurea complex with its related species/genera. Key Results: Our analysis showed that populations in Japan form an evolutionary unit, which was genetically diverged from adjacent continental populations. The phylogenetic structure within the archipelago strongly corresponds to the geography, but interestingly there is no concordance between genetic structure and ecotypic boundaries; neighbouring populations of distinct ecotypes share a genetic background. Conclusions: We propose that the traits specific to the ecotypic entities are maintained by natural selection or are very recently generated and have little effect on the genomes, making genome-wide genetic markers unsuitable for detecting ecotypic differentiation. Furthermore, some sporadically distributed taxa (found as rheophytes and alpine plants) were repeatedly generated from a more widespread taxon in geographically distant areas by means of selection. Overall, this study showed that the goldenrod complex has a high ability to evolve, enabling rapid ecological diversification over a recent timeframe.

A Genome Scan for Genes Underlying Microgeographic-Scale Local Adaptation in a Wild Arabidopsis Species.

Adaptive divergence at the microgeographic scale has been generally disregarded because high gene flow is expected to disrupt local adaptation. Yet, growing number of studies reporting adaptive divergence at a small spatial scale highlight the importance of this process in evolutionary biology. To investigate the genetic basis of microgeographic local adaptation, we conducted a genome-wide scan among sets of continuously distributed populations of Arabidopsis halleri subsp. gemmifera that show altitudinal phenotypic divergence despite gene flow. Genomic comparisons were independently conducted in two distinct mountains where similar highland ecotypes are observed, presumably as a result of convergent evolution. Here, we established a de novo reference genome and employed an individual-based resequencing for a total of 56 individuals. Among 527,225 reliable SNP loci, we focused on those showing a unidirectional allele frequency shift across altitudes. Statistical tests on the screened genes showed that our microgeographic population genomic approach successfully retrieve genes with functional annotations that are in line with the known phenotypic and environmental differences between altitudes. Furthermore, comparison between the two distinct mountains enabled us to screen out those genes that are neutral or adaptive only in either mountain, and identify the genes involved in the convergent evolution. Our study demonstrates that the genomic comparison among a set of genetically connected populations, instead of the commonly-performed comparison between two isolated populations, can also offer an effective screening for the genetic basis of local adaptation.

Simultaneous evaluation of the effects of geographic, environmental and temporal isolation in ecotypic populations of Solidago virgaurea.

Early stages of ecological speciation can create populations with an ecology and reproduction timing distinct from those of related populations. Landscape genetic models incorporating environmental heterogeneity and population-specific reproductive traits enable the processes of population genetic differentiation to be inferred. We investigated genome-wide genetic variation in ecotypic populations of Solidago virgaurea sensu lato, a herbaceous plant inhabiting a wide range of habitats (woodlands, serpentine barrens and alpine grasslands) and displaying remarkable variation in flowering time. Simultaneous evaluation of environmental factors revealed an overwhelming effect of soil type differences on neutral genetic differentiation, compared with elevational differences. This result probably reflects the abrupt environmental changes generated by geological boundaries, whereas mountain slopes exhibit clinal changes, facilitating gene exchange between neighbouring populations. Temporal isolation was positively associated with genetic differentiation, with some early-flowering serpentine populations having allele frequencies distinct from adjacent nonserpentine populations. Overall, this study highlights the importance of ecological processes and of evolution of flowering time to promote genetic differentiation of S. virgaurea populations in a complex landscape.

Transient hybridization, not homoploid hybrid speciation, between ancient and deeply divergent conifers.

PREMISE OF THE STUDY: Homoploid hybrid speciation is receiving growing attention due the increasing recognition of its role in speciation. We investigate if individuals intermediate in morphology between the two species of the conifer genus Athrotaxis represent a homoploid hybrid species, A. laxifolia, or are spontaneous F1 hybrids. METHODS: A total of 1055 individuals of Athrotaxis cupressoides and A. selaginoides, morphologically intermediate individuals, and two putative hybrid swarms were sampled across the range of the genus and genotyped with 13 microsatellites. We used simulations to test the power of our data to identify the pure species, F1s, F2s, and backcross generations. KEY RESULTS: We found that Athrotaxis cupressoides and A. selaginoides are likely the most divergent congeneric conifers known, but the intermediates are F1 hybrids, sharing one allele each from A. cupressoides and A. selaginoides at six loci with completely species specific alleles. The hybrid swarms contain wide genetic variation with stronger affinities to the locally dominant species, A. selaginoides and A. selaginoides backcrosses outnumbering A. cupressoides backcrosses. In addition, we observed evidence for isolated advanced generation backcrosses within the range of the pure species. CONCLUSIONS: We conclude that, even though they can be large and long-lived, Athrotaxis hybrid swarms are on a trajectory of decline and will eventually be reabsorbed by the parental species. However, this process may take millennia and fossil evidence suggests that such events have occurred repeatedly since the early Quaternary. Given this timeline, our study highlights the many obstacles to homoploid hybrid speciation.

Wide distribution range of rhizobial symbionts associated with pantropical sea-dispersed legumes.

To understand the geographic distributions of rhizobia that associated with widely distributed wild legumes, 66 nodules obtained from 41 individuals including three sea-dispersed legumes (Vigna marina, Vigna luteola, and Canavalia rosea) distributed across the tropical and subtropical coastal regions of the world were studied. Partial sequences of 16S rRNA and nodC genes extracted from the nodules showed that only Bradyrhizobium and Sinorhizobium were associated with the pantropical legumes, and some of the symbiont strains were widely distributed over the Pacific. Horizontal gene transfer of nodulation genes were observed within the Bradyrhizobium and Sinorhizobium lineages. BLAST searches in GenBank also identified records of these strains from various legumes across the world, including crop species. However, one of the rhizobial strains was not found in GenBank, which implies the strain may have adapted to the littoral environment. Our results suggested that some rhizobia, which associate with the widespread sea-dispersed legume, distribute across a broad geographic range. By establishing symbiotic relationships with widely distributed rhizobia, the pantropical legumes may also be able to extend their range much further than other legume species.

A Receptor-Like Kinase, Related to Cell Wall Sensor of Higher Plants, is Required for Sexual Reproduction in the Unicellular Charophycean Alga, Closterium peracerosum-strigosum-littorale Complex.

Here, we cloned the CpRLK1 gene, which encodes a receptor-like protein kinase expressed during sexual reproduction, from the heterothallic Closterium peracerosum-strigosum-littorale complex, one of the closest unicellular alga to land plants. Mating-type plus (mt(+)) cells with knockdown of CpRLK1 showed reduced competence for sexual reproduction and formed an abnormally enlarged conjugation papilla after pairing with mt(-) cells. The knockdown cells were unable to release a naked gamete, which is indispensable for zygote formation. We suggest that the CpRLK1 protein is an ancient cell wall sensor that now functions to regulate osmotic pressure in the cell to allow proper gamete release.

AP2-type transcription factors determine stem cell identity in the moss Physcomitrella patens.

Stem cells are formed at particular times and positions during the development of multicellular organisms. Whereas flowering plants form stem cells only in the sporophyte generation, non-seed plants form stem cells in both the sporophyte and gametophyte generations. Although the molecular mechanisms underlying stem cell formation in the sporophyte generation have been extensively studied, only a few transcription factors involved in the regulation of gametophyte stem cell formation have been reported. The moss Physcomitrella patens forms a hypha-like body (protonema) and a shoot-like body (gametophore) from a protonema apical cell and a gametophore apical cell, respectively. These apical cells have stem cell characteristics and are formed as side branches of differentiated protonema cells. Here, we show that four AP2-type transcription factors orthologous to Arabidopsis thaliana AINTEGUMENTA, PLETHORA and BABY BOOM (APB) are indispensable for the formation of gametophore apical cells from protonema cells. Quadruple disruption of all APB genes blocked gametophore formation, even in the presence of cytokinin, which enhances gametophore apical cell formation in the wild type. All APB genes were expressed in emerging gametophore apical cells, but not in protonema apical cells. Heat-shock induction of an APB4 transgene driven by a heat-shock promoter increased the number of gametophores. Expression of all APB genes was induced by auxin but not by cytokinin. Thus, the APB genes function synergistically with cytokinin signaling to determine the identity of the two types of stem cells.

From ß- to α-proteobacteria: the origin and evolution of rhizobial nodulation genes nodIJ.

Although many α- and some ß-proteobacterial species are symbiotic with legumes, the evolutionary origin of nitrogen-fixing nodulation remains unclear. We examined α- and ß-proteobacteria whose genomes were sequenced using large-scale phylogenetic profiling and revealed the evolutionary origin of two nodulation genes. These genes, nodI and nodJ (nodIJ), play key roles in the secretion of Nod factors, which are recognized by legumes during nodulation. We found that only the nodulating ß-proteobacteria, including the novel strains isolated in this study, possess both nodIJ and their paralogous genes (DRA-ATPase/permease genes). Contrary to the widely accepted scenario of the α-proteobacterial origin of rhizobia, our exhaustive phylogenetic analysis showed that the entire nodIJ clade is included in the clade of Burkholderiaceae DRA-ATPase/permease genes, that is, the nodIJ genes originated from gene duplication in a lineage of the ß-proteobacterial family. After duplication, the evolutionary rates of nodIJ were significantly accelerated relative to those of homologous genes, which is consistent with their novel function in nodulation. The likelihood analyses suggest that this accelerated evolution is not associated with changes in either nonsynonymous/synonymous substitution rates or transition/transversion rates, but rather, in the GC content. Although the low GC content of the nodulation genes has been assumed to reflect past horizontal transfer events from donor rhizobial genomes with low GC content, no rhizobial genome with such low GC content has yet been found. Our results encourage a reconsideration of the origin of nodulation and suggest new perspectives on the role of the GC content of bacterial genes in functional adaptation.

Herbivore pressure by weevils associated with flower color polymorphism in Geranium thunbergii (Geraniaceae).

Although floral herbivory has recently received increased attention as an important factor influencing plant reproduction, relatively little is known about how its frequency and intensity vary depending on traits of host plants. Here we report that herbivore pressure by a weevil, Zacladus geranii, is associated with a flower color polymorphism of Geranium thunbergii (Geraniaceae). Pink and white flower color morphs have been reported in G. thunbergii, and we found in a three-year field survey in multiple populations that, generally, adult weevils more preferentially visited white flowers than pink flowers. Consistently, we found more severe damage by weevil larvae in white flowers. Overall herbivore pressure for G. thunbergii varied strongly between populations, and the difference seems to be partly explained by the co-occurrence of a related plant species, Geranium yezoense, in a population, as weevils preferred it to both color morphs of G. thunbergii, thereby relaxing overall herbivore pressure for G. thunbergii. Nonetheless, despite such high variability, the preference of weevils for white morphs over pink morphs of G. thunbergii was found across multiple populations. We discuss possible mechanisms causing the association between flower color and herbivore preference as well as its evolutionary consequences.

Development of mitochondrial DNA cytochrome c oxidase subunit I primer sets to construct DNA barcoding library using next-generation sequencing.

Insects are one of the most diverse eukaryotic groups on the planet, with one million or more species present, including those yet undescribed. The DNA barcoding system has been developed, which has aided in the identification of cryptic species and undescribed species. The mitochondrial cytochrome c oxidase I region (mtDNA COI) has been utilised for the barcoding analysis of insect taxa. Thereafter, next-generation sequencing (NGS) technology has been developed, allowing for rapid acquisition of massive amounts of sequence data for genetic analyses. Although NGS-based PCR primers designed to amplify the mtDNA COI region have been developed, their target regions were only a part of COI region and/or there were taxonomic bias for PCR amplification. As the mtDNA COI region is a traditional DNA marker for the DNA barcoding system, modified primers for this region would greatly contribute to taxonomic studies. In this study, we redesigned previously developed PCR primer sets that targetted the mtDNA COI barcoding region to improve amplification efficiency and to enable us to conduct sequencing analysis on NGS. As a result, the redesigned primer sets achieved a high success rate (> 85%) for species examined in this study, covering four insect orders (Coleoptera, Lepidoptera, Orthoptera and Odonata). Thus, by combining the primers with developed primer sets for 12S or 16S rRNA regions, we can conduct more detailed taxonomic, phylogeographic and conservation genetic studies using NGS.

Taxonomic review of the genus Rhinoncomimus (Coleoptera: Curculionidae: Ceutorhynchinae) with description of a new species from Yunnan, China.

The genus Rhinoncomimus Wagner, 1940 includes seven species from Eastern Asia. One new species, Rh. continuus sp. nov. from Yunnan, China, is described. Habitus photos, illustrations and descriptions of all species except Rh. rubripes Korotyaev, 2006 (a possible junior synonym of R. niger Chûjô and Morimoto, 1959) are provided in detail, as well as key to species and distribution maps. In addition, the host plant of the type species Rh. klapperichi Wagner, 1940, Polygonum hydropiper L. (Polygonaceae) is newly recorded.

Extensive Copy Number Variation Explains Genome Size Variation in the Unicellular Zygnematophycean Alga, Closterium peracerosum-strigosum-littorale Complex.

Genome sizes are known to vary within and among closely related species, but the knowledge about genomic factors contributing to the variation and their impacts on gene functions is limited to only a small number of species. This study identified a more than 2-fold heritable genome size variation among the unicellular Zygnematophycean alga, Closterium peracerosum-strigosum-littorale (C. psl.) complex, based on short-read sequencing analysis of 22 natural strains and F1 segregation analysis. Six de novo assembled genomes revealed that genome size variation is largely attributable to genome-wide copy number variation (CNV) among strains rather than mating type-linked genomic regions or specific repeat sequences such as rDNA. Notably, about 30% of genes showed CNV even between strains that can mate with each other. Transcriptome and gene ontology analysis demonstrated that CNV is distributed nonrandomly in terms of gene functions, such that CNV was more often observed in the gene set with stage-specific expression. Furthermore, in about 30% of these genes with CNV, the expression level does not increase proportionally with the gene copy number, suggesting presence of dosage compensation, which was overrepresented in genes involved in basic biological functions, such as translation. Nonrandom patterns in gene duplications and corresponding expression changes in terms of gene functions may contribute to maintaining the high level of CNV associated with extensive genome size variation in the C. psl. complex, despite its possible detrimental effects.

Evolution and biogeography of Pleurophyllum (Astereae, Asteraceae), a small genus of megaherbs endemic to the subantarctic islands.

PREMISE OF THE STUDY: The abundance of fossils in Antarctica suggests this continent was a center of diversification and a corridor for migration for many austral plant groups until the late Tertiary and may have played a pivotal role in shaping plant distributions in the southern hemisphere. Although the Antarctic flora was largely erased by glaciation during the Pleistocene, at least some Antarctic plant species found refuge on the subantarctic islands. METHODS: We used independent and combined analyses of ITS, ETS, trnK, and trnL DNA sequences to infer phylogenetic relations in Pleurophyllum, a small genus of three species that are endemic to the subantarctic islands of Australia and New Zealand. The inferred phylogeny provided a framework to reconstruct the origin and patterns of diversification in the genus. KEY RESULTS: We summarize support for the hypothesis that Pleurophyllum survived episodes of Pleistocene glaciation in the subantarctic islands and that its sisters dispersed northward in response to glacial advance. CONCLUSIONS: The distinctive flora of the subantarctic islands includes some of the last remnants of a once-diverse Antarctic flora. These plants may still retain distinctive features of their ancestors. Studies of endemic plants such as Pleurophyllum are the key to resolving this puzzle.

Reticulate evolution in the Crepidomanes minutum species complex (Hymenophyllaceae).

PREMISE OF THE STUDY: Hybridization is an important mechanism of speciation in plants and often results in complexes that comprise multiple sexual diploids and their hybrid offspring. However, the intricacy of these systems has prevented a thorough understanding of many groups. The Crepidomanes minutum species complex (Hymenophyllaceae) is a widely distributed, morphologically variable fern species complex. Although prior reports of apogamy and polyploidy suggest hybridization, it has never been the focus of a phylogenetic study. METHODS: Morphology, nuclear (gapCp), and chloroplast (rbcL) DNA sequences, cytology, field observation, and spore counts were used to infer phylogeny and trace hybrid origins. KEY RESULTS: The C. minutum species complex is composed of at least three major clades: the African clade, clade 1 (East Asia and the Pacific), and clade 2 (Southeast Asia and the South Pacific). Clades 1 and 2 differ strikingly in morphological variation (uniform in clade 1 vs. highly variable in clade 2) and occurrence of hybrids (rare in clade 1 vs. frequent in clade 2). Apogamy and polyploidy were confirmed as likely mechanisms of hybrid stabilization in clade 2. Despite the large genetic distance between clades 1 and 2, several specimens were observed with gapCp sequences from both; diploid genome size and sexual reproduction indicate maintenance of genetic diversity via introgression or incomplete lineage sorting, rather than ongoing hybridization, in these specimens. CONCLUSIONS: The C. minutum species complex is a reticulate network including multiple diploid lineages and their stabilized hybrid crosses. Additional sampling focused on reproductive mode and ploidy level is needed to delimit diploid species and hybrids.

Bacterial symbionts of the giant jewel stinkbug Eucorysses grandis (Hemiptera: Scutelleridae).

Microbiological characterization of gut symbiotic bacteria in a limited number of stinkbugs of the families Acanthosomatidae, Plataspidae, Pentatomidae, Scutelleridae, Parastrachiidae, Alydidae and Pyrrhocoridae has shown symbiotic association with midgut bacteria to be common in phytophagous taxa of these heteropteran insects. Here we investigated the midgut bacterial symbiont of Eucorysses grandis, a stinkbug of the family Scutelleridae. A specific gammaproteobacterium was consistently identified in insects from five different geographic origins. The bacterium was detected in 64 of 64 insects sampled from three host populations. Phylogenetic analyses revealed that the bacterium constitutes a distinct lineage in the Gammaproteobacteria, neither closely related to the gut symbiont of another scutellerid stinkbug, Cantao ocellatus, nor to gut symbionts of other stinkbugs. Diagnostic PCR, in situ hybridization and electron microscopy demonstrated that the bacterium is located extracelluarly, in the midgut fourth section, which possesses crypts. These results indicate that the primary gut symbionts have multiple evolutionary origins in the Scutelleridae. A Sodalis-allied facultative symbiont was also identified in some insects from natural populations. Biological aspects of the primary gut symbiont and the secondary Sodalis-allied symbiont are discussed.