Double migration of the endangered Tricyrtis formosana (Liliaceae) in Japan

The Ryukyu Islands of Japan are a biodiversity hotspot due to geographical and historical factors. Tricyrtis formosana is a perennial herbaceous plant that commonly found in Taiwan. But only a few populations have been identified in a limited habitat on Iriomote Island, while populations of unknown origin occur near human settlements in an area on the main island of Okinawa. To better understand these populations of the phylogenetic uniqueness and intrinsic vulnerability, we conducted comparative analyses including (1) phylogeny and population structure with MIG-seq data, (2) photosynthesis-related traits of plants grown under common conditions and (3) transcriptome analysis to detect deleterious variations. Results revealed that T. formosana was split into two clades by the congeners and that Iriomote and Okinawa populations independently derived from ancestral Taiwanese populations in each clade. Photosynthetic efficiency was lowest in the Iriomote population, followed by Okinawa and Taiwan. Transcriptome analysis showed that the Iriomote population accumulated more deleterious variations, suggesting intrinsic vulnerability. These results indicate that each T. formosana population in Japan is phylogenetically unique and has been independently dispersed from Taiwan, and that the Iriomote population presents a high conservation difficulty with a unique photosynthesis-related characteristic and a larger amount of deleterious variations.

Temporal and interspecific dietary variation in wintering ducks in agricultural landscapes

Farmlands are becoming more important as waterfowl foraging habitats, while natural wetlands are being lost globally. However, it is unclear how waterfowl coexist in agricultural landscapes by resource partitioning. We evaluated the diets of seven sympatric dabbling ducks foraging in rice paddy and lotus fields around Lake Kasumigaura, the second largest lake in Japan, during two wintering seasons (from November to February) by faecal DNA metabarcoding using chloroplast trnL and mitochondrial CO1 region sequences. We examined 420 faecal samples and found different patterns of dietary diversity and composition among the duck species. The pattern also differed between plant and invertebrate food. Dietary niche partitioning was clear in plant food. Large-bodied ducks intensively use crop plants, and other ducks might mediate competition by using terrestrial and aquatic plants that are suitable for their foraging behaviours or microhabitats. Dietary segregation among species was the most apparent in February, when the abundance of foraging ducks was the largest. This study illustrated the complex pattern of dietary niche partitioning of dabbling ducks in agricultural landscapes, which might be difficult to evaluate by conventional approaches. The availability of crop plants, as well as other plant food resources in flooted areas and farmland dikes, may enable ducks to coexist by spatial or behavioural resource partitioning.

Phylogeographical evidence for historical long-distance dispersal in the flightless stick insect Ramulus mikado

Exploring how organisms overcome geographical barriers to dispersal is a fundamental question in biology. Passive long-distance dispersal events, although infrequent and unpredictable, have a considerable impact on species range expansions. Despite limited active dispersal capabilities, many stick insect species have vast geographical ranges, indicating that passive long-distance dispersal is vital for their distribution. A potential mode of passive dispersal in stick insects is via the egg stage within avian digestive tracts, as suggested by experimental evidence. However, detecting such events under natural conditions is challenging due to their rarity. Therefore, to indirectly assess the potential of historical avian-mediated dispersal, we examined the population genetic structure of the flightless stick insect Ramulus mikado across Japan, based on a multifaceted molecular approach [cytochrome oxidase subunit I (COI) haplotypes, nuclear simple sequence repeat markers and genome-wide single nucleotide polymorphisms]. Subsequently, we identified unique phylogeographic patterns, including the discovery of identical COI genotypes spanning considerable distances, which substantiates the notion of passive long-distance genotypic dispersal. Overall, all the molecular data revealed the low and mostly non-significant genetic differentiation among populations, with identical or very similar genotypes across distant populations. We propose that long-distance dispersal facilitated by birds is the plausible explanation for the unique phylogeographic pattern observed in this flightless stick insect.

Genetic and demographic signatures accompanying the evolution of the selfing syndrome in Daphne kiusiana, an evergreen shrub

BACKGROUND AND AIMS

The evolution of mating systems from outcrossing to self-fertilization is a common transition in flowering plants. This shift is often associated with the 'selfing syndrome', which is characterized by less visible flowers with functional changes to control outcrossing. In most cases, the evolutionary history and demographic dynamics underlying the evolution of the selfing syndrome remain poorly understood.

METHODS

Here, we characterize differences in the demographic genetic consequences and associated floral-specific traits between two distinct geographical groups of a wild shrub, Daphne kiusiana, endemic to East Asia; plants in the eastern region (southeastern Korea and Kyushu, Japan) exhibit smaller and fewer flowers compared to those of plants in the western region (southwestern Korea). Genetic analyses were conducted using nuclear microsatellites and chloroplast DNA (multiplexed phylogenetic marker sequencing) datasets.

KEY RESULTS

A high selfing rate with significantly increased homozygosity characterized the eastern lineage, associated with lower levels of visibility and herkogamy in the floral traits. The two lineages harboured independent phylogeographical histories. In contrast to the western lineage, the eastern lineage showed a gradual reduction in the effective population size with no signs of a severe bottleneck despite its extreme range contraction during the last glacial period.

CONCLUSIONS

Our results suggest that the selfing-associated morphological changes in D. kiusiana are of relatively old origin (at least 100 000 years ago) and were driven by directional selection for efficient self-pollination. We provide evidence that the evolution of the selfing syndrome in D. kiusiana is not strongly associated with a severe population bottleneck.

The contrary conservation situations of two local critically endangered species, Vaccinium emarginatum (Ericaceae) and Elatostema platyphyllum (Urticaceae), growing on the eastern edge of the distribution

As biodiversity loss continues, there is an urgent need to develop efficient conservation measures to protect diversity with limited conservation resources. Conservation targets have generally been selected based on their population size, but more detailed assessments clarifying the phylogenetic genetic status, history, and phylogenetic uniqueness of rare species is crucial to set more appropriate and effective conservation measures. In Japan, the Act on Conservation of Endangered Species of Wild Fauna and Flora designated endangered plants with high conservation priority, but >40% of these species also grow overseas. We conducted comparative analyses based on ddRADseq and MIG-seq to evaluate the population conservation status and value of Vaccinium emarginatum and Elatostema platyphyllum which are growing across national borders at the eastern edge of their species distribution range. The analyses revealed contrasting conservation status between the two species; the Japanese population of V. emarginatum had lower genetic diversity at the individual level and phylogenetically differentiated from Taiwanese populations, while that of E. platyphyllum had higher diversity at the individual level and is a relatively recent migrant with little phylogenetical differentiation from Taiwanese populations. The two species, which share the common feature of being critically rare in Japan, showed contrasting genetic/phylogenetic characteristics. This study provided useful information for appropriate conservation measures based on species’ phylogenetic traits and genetic diversity.

Effects of fruit dimorphism on genetic structure and gene flow in the coastal shrub Scaevola taccada

BACKGROUND AND AIMS

Plant propagules often possess specialized morphologies that facilitate dispersal across specific landscapes. In the fruit dimorphism of a coastal shrub, Scaevola taccada, individual plants produce either cork-morph or pulp-morph fruits. The former is buoyant and common on sandy beaches, whereas the latter does not float, is bird-dispersed, and is common on elevated sites such as slopes on sea cliffs and behind rocky shores. We hypothesized that beach populations bridge the heterogeneous landscapes by serving as a source of both fruit types, while dispersal is biased for the pulp morph on elevated sites within the islands and for the cork morph between beaches of different islands. Based on this hypothesis, we predicted that populations in elevated sites would diverge genetically over time due to isolation by distance, whereas beach populations would maintain high genetic similarity via current gene flow.

METHODS

The genetic structure and gene flow in S. taccada were evaluated by investigating genome-wide single nucleotide polymorphisms in plants from 17 sampling sites on six islands (belonging to the Ryukyu, Daito and Ogasawara Islands) in Japan.

KEY RESULTS

Geographical isolation was detected among the three distant island groups. Analyses within the Ryukyu Islands suggested that sandy beach populations were characterized by genetic admixture, whereas populations in elevated sites were relatively isolated between the islands. Pairwise FST values between islands were lowest between sandy beaches, intermediate between sandy beaches and elevated sites, and highest between elevated sites.

CONCLUSIONS

Dispersal across the ocean by cork morphs is sufficiently frequent to prevent genetic divergence between beaches of different islands. Stronger genetic isolation of elevated sites between islands suggests that bird dispersal by pulp morphs is restricted mainly within islands. These contrasting patterns of gene flow realized by fruit dimorphism provide evidence that fruit characteristics can strongly mediate genetic structure.

Demography and selection analysis of the incipient adaptive radiation of a Hawaiian woody species

Ecological divergence in a species provides a valuable opportunity to study the early stages of speciation. We focused on Metrosideros polymorpha, a unique example of the incipient radiation of woody species, to examine how an ecological divergence continues in the face of gene flow. We analyzed the whole genomes of 70 plants collected throughout the island of Hawaii, which is the youngest island with the highest altitude in the archipelago and encompasses a wide range of environments. The continuous M. polymorpha forest stands on the island of Hawaii were differentiated into three genetic clusters, each of which grows in a distinctive environment and includes substantial genetic and phenotypic diversity. The three genetic clusters showed signatures of selection in genomic regions encompassing genes relevant to environmental adaptations, including genes associated with light utilization, oxidative stress, and leaf senescence, which are likely associated with the ecological differentiation of the species. Our demographic modeling suggested that the glaberrima cluster in wet environments maintained a relatively large population size and two clusters split: polymorpha in the subalpine zone and incana in dry and hot conditions. This ecological divergence possibly began before the species colonized the island of Hawaii. Interestingly, the three clusters recovered genetic connectivity coincidentally with a recent population bottleneck, in line with the weak reproductive isolation observed in the species. This study highlights that the degree of genetic differentiation between ecologically-diverged populations can vary depending on the strength of natural selection in the very early phases of speciation.

Knowledge about how genetic barriers are formed between populations in distinct environments is valuable to understand the processes of speciation and conserve biodiversity. Metrosideros polymorpha, an endemic woody species in the Hawaiian Islands, is a good system to study developing genetic barriers in a species, because it colonized the diverse environments and diversified the morphology for a relatively short period of time. We analyzed the genomes of 70 M. polymorpha plants from a broad range of environments on the island of Hawaii to infer the current and past genetic barriers among them. Currently, M. polymorpha plants growing in different environments have substantially different genomes, especially at the genomic regions with genes putatively controlling physiology to fit in distinct environment. However, in its history, they had hybridized with one another, possibly because plants formerly growing in different environments came into close contact due to the climate changes. It is suggested that genetic barriers can easily strengthen or weaken depending on environments splitting the ecology of a species before reproductive isolation becomes complete.

The genome of Shorea leprosula (Dipterocarpaceae) highlights the ecological relevance of drought in aseasonal tropical rainforests

Hyperdiverse tropical rainforests, such as the aseasonal forests in Southeast Asia, are supported by high annual rainfall. Its canopy is dominated by the species-rich tree family of Dipterocarpaceae (Asian dipterocarps), which has both ecological (e.g., supports flora and fauna) and economical (e.g., timber production) importance. Recent ecological studies suggested that rare irregular drought events may be an environmental stress and signal for the tropical trees. We assembled the genome of a widespread but near threatened dipterocarp, Shorea leprosula, and analyzed the transcriptome sequences of ten dipterocarp species representing seven genera. Comparative genomic and molecular dating analyses suggested a whole-genome duplication close to the Cretaceous-Paleogene extinction event followed by the diversification of major dipterocarp lineages (i.e. Dipterocarpoideae). Interestingly, the retained duplicated genes were enriched for genes upregulated by no-irrigation treatment. These findings provide molecular support for the relevance of drought for tropical trees despite the lack of an annual dry season.

Development of microsatellite markers for the geographically parthenogenetic stick insect Phraortes elongatus (Insecta: Phasmatodea)

Many plant and animal species exhibit geographic parthenogenesis, wherein unisexual (= parthenogenetic) lineages are more common in their marginal habitats such as high latitude or altitudes than their closely related bisexual counterparts. The Japanese stick insect, Phraortes elongatus (Thunberg) (Insecta: Phasmatodea), is known as a geographically parthenogenetic species due to the existence of both bisexual and unisexual populations. Here, we developed microsatellite markers to infer the genetic variation among populations of P. elongatus. Totally, 13 primer pairs were developed for the species, and they were tested on 47 samples collected from both a bisexual population and a unisexual population. All 13 loci were polymorphic in the bisexual population, whereas no loci were polymorphic in the unisexual population. The loss of variation in the unisexual population implies automixis with terminal fusion or gamete duplication as the mode of parthenogenesis. The markers developed in this study will be helpful for further comprehensive analysis of the genetic diversity and gene flow between bisexual and parthenogenetic lineages of P. elongatus.

Directional seed and pollen dispersal and their separate effects on anisotropy of fine-scale spatial genetic structure among seedlings in a dioecious, wind-pollinated, and wind-dispersed tree species, Cercidiphyllum japonicum

Prevailing directions of seed and pollen dispersal may induce anisotropy of the fine-scale spatial genetic structure (FSGS), particularly in wind-dispersed and wind-pollinated species. To examine the separate effects of directional seed and pollen dispersal on FSGS, we conducted a population genetics study for a dioecious, wind-pollinated, and wind-dispersed tree species, Cercidiphyllum japonicum Sieb. et Zucc, based on genotypes at five microsatellite loci of 281 adults of a population distributed over a ca. 80 ha along a stream and 755 current-year seedlings. A neighborhood model approach with exponential-power-von Mises functions indicated shorter seed dispersal (mean = 69.1 m) and much longer pollen dispersal (mean = 870.6 m), effects of dispersal directions on the frequencies of seed and pollen dispersal, and the directions with most frequent seed and pollen dispersal (prevailing directions). Furthermore, the distance of effective seed dispersal within the population was estimated to depend on the dispersal direction and be longest at the direction near the prevailing direction. Therefore, patterns of seed and pollen dispersal may be affected by effective wind directions during the period of respective dispersals. Isotropic FSGS and spatial sibling structure analyses indicated a significant FSGS among the seedlings generated by the limited seed dispersal, but anisotropic analysis for the seedlings indicated that the strength of the FSGS varied with directions between individuals and was weakest at a direction near the directions of the most frequent and longest seed dispersal but far from the prevailing direction of pollen dispersal. These results suggest that frequent and long-distance seed dispersal around the prevailing direction weakens the FSGS around the prevailing direction. Therefore, spatially limited but directional seed dispersal would determine the existence and direction of FSGS among the seedlings.

Inferring historical survivals of climate relicts: the effects of climate changes, geography, and population-specific factors on herbaceous hydrangeas

Climate relicts hold considerable importance because they have resulted from numerous historical changes. However, there are major interspecific variations among the ways by which they survived climate changes. Therefore, investigating the factors and timing that affected population demographics can expand our understanding of how climate relicts responded to historical environmental changes. Here, we examined herbaceous hydrangeas of genus Deinanthe in East Asia, which show limited distributions and a remarkable disjunction between Japan and central China. Chloroplast genome and restriction site-associated DNA sequencing revealed that speciation event occurred in the late Miocene (ca. 7-9 Mya) in response to global climate change. Two lineages apparently remained not branched until the middle Quaternary, and afterwards started to diverge to regional population groups. The narrow endemic species in central China showed lower genetic diversity (He = 0.082), as its population size rapidly decreased during the Holocene due to isolation in montane refugia. Insular populations in the three Japanese islands (He = 0.137-0.160) showed a genetic structure that was inconsistent with sea barriers, indicating that it was shaped in the glacial period when its range retreated to coastal refugia on the exposed sea floor. Demographic modelling by stairway-plot analysis reconstructed variable responses of Japanese populations: some experienced glacial bottlenecks in refugial isolation, while post-glacial range expansion seemingly exerted founder effects on other populations. Overall, this study demonstrated the involvement of not just one, but multiple factors, such as the interplay between climate changes, geography, and other population-specific factors, that determine the demographics of climate relicts.

Detection of Endangered Aquatic Plants in Rapid Streams Using Environmental DNA

Podostemaceae are a eudicot family of plants that grow on rapid streams and waterfalls. Two genera and six species of this family are distributed in Japan, all of which are threatened with extinction. It is difficult to find these species from the river side and it takes much effort to investigate their distribution. In this study, we attempted to determine the presence and absence of the Podostemaceae species by environmental DNA (eDNA) metabarcoding. Four species of Podostemaceae were detected near four known habitats, and the detected species were in perfect agreement with the results of a past survey that was based on visual observation. The marker used in this study had sufficient resolution to distinguish all six Podostemaceae species distributed in Japan and detected multiple species growing in a site. These results show that eDNA metabarcoding can quickly detect rare aquatic plants that are difficult to find by visual observation and can provide important information regarding their conservation.

Genomic insights on the contribution of balancing selection and local adaptation to the long-term survival of a widespread living fossil tree, Cercidiphyllum japonicum

'Living fossils' are testimonies of long-term sustained ecological success, but how demographic history and natural selection contributed to their survival, resilience, and persistence in the face of Quaternary climate fluctuations remains unclear. To better understand the interplay between demographic history and selection in shaping genomic diversity and evolution of such organisms, we assembled the whole genome of Cercidiphyllum japonicum, a widespread East Asian Tertiary relict tree, and resequenced 99 individuals of C. japonicum and its sister species, Cercidiphyllum magnificum (Central Japan). We dated this speciation event to the mid-Miocene, and the intraspecific lineage divergence of C. japonicum (China vs Japan) to the Early Pliocene. Throughout climatic upheavals of the late Tertiary/Quaternary, population bottlenecks greatly reduced the genetic diversity of C. japonicum. However, this polymorphism loss was likely counteracted by, first, long-term balancing selection at multiple chromosomal and heterozygous gene regions, potentially reflecting overdominance, and, second, selective sweeps at stress response and growth-related genes likely involved in local adaptation. Our findings contribute to a better understanding of how living fossils have survived climatic upheaval and maintained an extensive geographic range; that is, both types of selection could be major factors contributing to the species' survival, resilience, and persistence.

Corrigendum: Development of microsatellite markers for the annual andromonoecious herb Commelina communis f. ciliata (Commelinaceae) [Genes Genet. Syst. (2019) 94, p. 133-138]

Table 1 on p. 135 should be replaced with the corrected Table 1 shown bellow.

Genetic consequences of plant edaphic specialization to solfatara fields: Phylogenetic and population genetic analysis of Carex angustisquama (Cyperaceae)

Edaphic specialization is one of the main drivers of plant diversification and has multifaceted effects on population dynamics. Carex angustisquama is a sedge plant growing only on heavily acidified soil in solfatara fields, where only extremophytes can survive. Because of the lack of closely related species in similar habitats and its disjunct distribution, the species offers ideal settings to investigate the effects of adaptation to solfatara fields and of historical biogeography on the genetic consequences of plant edaphic specialization to solfatara fields. Here, genome-wide single nucleotide polymorphisms were used to reveal the phylogenetic origin of C. angustisquama, and 16 expressed sequence tag-simple sequence repeat markers were employed to infer population demography of C angustisquama. Molecular phylogenetic analysis strongly indicated that C. angustisquama formed a monophyletic clade with Carex doenitzii, a species growing on nonacidified soil in the sympatric subalpine zone. The result of population genetic analysis showed that C. angustisquama has much lower genetic diversity than the sister species, and notably, all 16 loci were completely homozygous in most individuals of C. angustisquama. Approximate Bayesian computation analysis supported the model that assumed hierarchical declines of population size through its evolutionary sequence. We propose that the edaphic specialist in solfatara fields has newly attained the adaptation to solfatara fields in the process of speciation. Furthermore, we found evidence of a drastic reduction in genetic diversity in C. angustisquama, suggesting that the repeated founder effects associated with edaphic specialization and subsequent population demography lead to the loss of genetic diversity of this extremophyte in solfatara fields.

Development of microsatellite markers for the endangered sleeper Eleotris oxycephala (Perciformes: Eleotridae)

The amphidromous sleeper Eleotris oxycephala (Perciformes: Eleotridae) is mainly distributed along the Kuroshio Current in East Asia, and this current is thought to be the main driver of the species' dispersal. Due to anthropogenic environmental changes in rivers, E. oxycephala is ranked as a threatened or near-threatened species in the red lists of 12 prefectures in Japan. Moreover, there is concern that the species' dispersal pattern could be changed due to fluctuations in the Kuroshio Current caused by global warming. In this study, 40 microsatellite markers were developed for E. oxycephala, and their suitability was tested on 43 individuals from two populations of E. oxycephala from Kanagawa and Miyazaki Prefectures. The number of alleles, expected heterozygosity and fixation index at each locus were 2-10 (mean = 5.350), 0.034-0.860 (mean = 0.650) and -0.261-0.448 (mean = 0.065), respectively. Furthermore, there was a lack of genetic difference between the two populations (F_ST = 0.008, F'_ST = 0.024), indicating widespread gene flow via the Kuroshio Current. These markers will be useful to evaluate the genetic structure and infer population demographic history of E. oxycephala populations, which may assist in the conservation of this species.

New microsatellite markers recognize differences in tandem repeats among four related Gastrodia species (Orchidaceae)

Gastrodia is the most species-rich genus among mycoheterotrophic plants, and is thus an essential taxon to understand the mechanism of species diversification in mycoheterotrophs. In this study, we developed microsatellite markers with high transferability for four Gastrodia species to examine genetic differentiation and similarity among species, populations and individuals. The 12 microsatellite markers developed from a G. fontinalis library showed high transferability for the ramets that identified G. nipponica, G. kuroshimensis and G. takeshimensis. In addition to the high transferability of these markers, we observed low allele variation within a sampled population of each species and allele differences among the four species. The 12 markers described here will be useful for investigating the genetic differences among and within the Gastrodia species, which evolved by a limitation of gene flow.

Development of microsatellite markers for the annual andromonoecious herb Commelina communis f. ciliata (Commelinaceae)

Commelina communis f. ciliata (Commelinaceae), a newly distinguished taxon, is an annual andromonoecious herb exhibiting a mixed mating system, the details of which remain unclear. We developed microsatellite markers for use in exploring the evolution of andromonoecy and mixed mating in the species. Fifteen microsatellite loci were developed using next-generation sequencing. The primer sets were used to evaluate 65 C. communis f. ciliata individuals from three populations in Japan; we found 1-13 alleles per locus and the expected heterozygosity ranged from 0.00 to 0.76. The markers are potentially useful to examine intra- and interspecies genetic structure and the mixed mating strategy of Commelina species via paternity analysis.

Endangered island endemic plants have vulnerable genomes

Loss of genetic diversity is known to decrease the fitness of species and is a critical factor that increases extinction risk. However, there is little evidence for higher vulnerability and extinction risk in endangered species based on genomic differences between endangered and non-endangered species. This is true even in the case of functional loci, which are more likely to relate to the fitness of species than neutral loci. Here, we compared the genome-wide genetic diversity, proportion of duplicated genes (PD), and accumulation of deleterious variations of endangered island endemic (EIE) plants from four genera with those of their non-endangered (NE) widespread congeners. We focused on exhaustive sequences of expressed genes obtained by RNA sequencing. Most EIE species exhibited significantly lower genetic diversity and PD than NE species. Additionally, all endangered species accumulated deleterious variations. Our findings provide new insights into the genomic traits of EIE species.

Genomic reconstruction of 100 000-year grassland history in a forested country: population dynamics of specialist forbs

Grassland ecosystems worldwide have been extensively converted to other land uses and are globally imperiled. Because many grasslands have been maintained by human activities, understanding their origin and history is fundamentally important to better contemporary management. However, existing methods to reconstruct past vegetation can produce contrasting views on grassland history. Here, we inferred demographic histories of 40 populations of four grassland forb species throughout Japan using high-resolution genome sequences and model-flexible demographic simulation based on the site frequency spectrum. Although two species showed a slight decline in population size between 100 000–10 000 years ago, our results suggest that population sizes of studied species have been maintained within the range of 0.5–2.0 times the most recent estimates for at least 100 000 years across Japan. Our results suggest that greater than 90% declines in Japanese grasslands and subsequent losses of grassland species in the last 100 years are geologically and biologically important and will have substantial consequences for Japanese biota and culture. People have had critical roles in maintaining disturbance-dependent grassland ecosystems and biota in this warm and wet forested country. In these contexts, disturbances associated with forest harvesting and traditional extensive farming have the potential to maintain grassland ecosystems and can provide important opportunities to reconcile resource production and conservation of grassland biodiversity.

Development of microsatellite markers for the completely cleistogamous species Gastrodia takeshimensis (Orchidaceae) that are transferable to its chasmogamous sister G. nipponica

We developed microsatellite markers to compare the genetic variation between the putatively cleistogamous Gastrodia takeshimensis (Orchidaceae) and its chasmogamous sister species G. nipponica. We expected low genetic variation in G. takeshimensis in view of its hypothesized cleistogamy. Eighteen primer pairs were developed from a G. takeshimensis genomic DNA library, and their characteristics were tested for G. takeshimensis and G. nipponica. Seven loci were polymorphic in G. nipponica, whereas all loci showed no polymorphism in G. takeshimensis. Genetic diversity was thus not detected in G. takeshimensis, and it seems to have been lost by repeated selfing in the completely closed flower. The 18 markers described here will be useful for investigating the genetic variation between a cleistogamous species and its chasmogamous sister species.

Development of microsatellite markers for an endangered fern in the Ryukyus, Plagiogyria koidzumii (Plagiogyriaceae)

We developed 10 microsatellite markers for Plagiogyria koidzumii, a critically endangered fern species found on Iriomotejima Island in the Ryukyu Archipelago, Japan and in Taiwan. These markers showed polymorphism among 65 wild individuals from Iriomotejima Island; the number of alleles per locus was 2-14, and mean observed and expected heterozygosity in the largest population were 0.276 and 0.277, respectively. A genetic structure analysis using these markers indicated clear genetic differentiation even within the narrow geographic range (ca. 10 × 8 km) on Iriomotejima Island. These microsatellite markers should be valuable for measuring genetic diversity and comparing genetic structure within and between populations.

Turnover of southern cypresses in the post-Gondwanan world: extinction, transoceanic dispersal, adaptation and rediversification

Cupressaceae subfamily Callitroideae has been an important exemplar for vicariance biogeography, but its history is more than just disjunctions resulting from continental drift. We combine fossil and molecular data to better assess its extinction and, sometimes, rediversification after past global change. Key fossils were reassessed and their phylogenetic placement for calibration was determined using trait mapping and Bayes Factors. Five vicariance hypotheses were tested by comparing molecular divergence times with the timing of tectonic rifting. The role of adaptation to fire (serotiny) in its spread across a drying Australia was tested for Callitris. Our findings suggest that three transoceanic disjunctions within the Callitroideae probably arose from long-distance dispersal. A signature of extinction, centred on the end-Eocene global climatic chilling and drying, is evident in lineages-through-time plots and in the fossil record. Callitris, the most diverse extant callitroid genus, suffered extinctions but surviving lineages adapted and re-radiated into dry, fire-prone biomes that expanded in the Neogene. Serotiny, a key adaptation to fire, likely evolved in Callitris coincident with the biome shift. Both extinction and adaptive shifts have probably played major roles in this chronicle of turnover and renewal, but better understanding of biogeographical history requires improved taxonomy of fossils.

Development and characterization of EST-SSR markers for Carex angustisquama (Cyperaceae), an extremophyte in solfatara fields

PREMISE OF THE STUDY

Expressed sequence tag-simple sequence repeat (EST-SSR) markers were developed for Carex angustisquama (Cyperaceae) to investigate the evolutionary history of this plant that is endemic to solfatara fields in northern Japan.

METHODS AND RESULTS

Using RNA-Seq data generated by the Illumina HiSeq 2000, 20 EST-SSR markers were developed. Polymorphisms were assessed in C. angustisquama and the closely related species C. doenitzii and C. podogyna. In C. angustisquama, many loci were monomorphic within populations; the average number of alleles ranged from one to five, and levels of expected heterozygosity ranged from 0.000 to 0.580, while all markers were polymorphic in a population of C. doenitzii. This indicates that low genetic polymorphism of C. angustisquama is likely due to the species' population dynamics, rather than to null alleles at the developed markers.

CONCLUSIONS

These markers will be used to assess genetic diversity and structure and to investigate evolutionary history in future studies of C. angustisquama and related species.

Historical changes in grassland area determined the demography of semi-natural grassland butterflies in Japan

Semi-natural grassland areas expanded worldwide several thousand years ago following an increase in anthropogenic activities. However, semi-natural grassland habitat areas have been declining in recent decades due to changes in landuse, which have caused a loss of grassland biodiversity. Reconstructing historical and recent demographic changes in semi-natural grassland species will help clarify the factors affecting their population decline. Here we quantified past and recent demographic histories of Melitaea ambigua (Lepidoptera; Nymphalidae), an endangered grassland butterfly species in Japan. We examined changes in demography over the past 10,000 years based on 1378 bp of mitochondrial COI gene. We then examined changes in its genetic diversity and structure during the last 30 years using nine microsatellite DNA markers. The effective population size of M. ambigua increased about 3000-6000 years ago. In contrast, the genetic diversity and effective population sizes of many populations significantly declined from the 1980s to 2010s, which is consistent with a recent decline in the species population size. Our data suggest that the M. ambigua demography can be traced to changes in area covered by semi-natural grasslands throughout the Holocene.

Development of microsatellite markers for the endangered orchid Calanthe izu-insularis (Orchidaceae)

Microsatellite markers were developed for the endangered orchid Calanthe izu-insularis (Orchidaceae). This species is unique to the Izu Islands in Japan. Unfortunately, its population size has decreased because of excessive collection for horticultural purposes. In addition, although natural hybridization between C. izu-insularis and C. discolor var. discolor has been reported, morphological differences between C. izu-insularis and the hybridized individuals remain unclear. Using next-generation sequencing, 11 polymorphic microsatellite markers were developed. All developed markers could amplify C. aristulifera and nine markers could amplify C. d. var. discolor, two other orchid species that are also endangered in Japan. The number of alleles and expected heterozygosity at each locus were 1-6 (mean, 2.35) and 0.00-0.79 (mean, 0.30), respectively. These microsatellite markers will help conservation geneticists in their investigation of the proportion of pure C. izu-insularis individuals in the Izu Islands.

DNA barcoding reveals seasonal shifts in diet and consumption of deep-sea fishes in wedge-tailed shearwaters

The foraging ecology of pelagic seabirds is difficult to characterize because of their large foraging areas. In the face of this difficulty, DNA metabarcoding may be a useful approach to analyze diet compositions and foraging behaviors. Using this approach, we investigated the diet composition and its seasonal variation of a common seabird species on the Ogasawara Islands, Japan: the wedge-tailed shearwater Ardenna pacifica. We collected fecal samples during the prebreeding (N = 73) and rearing (N = 96) periods. The diet composition of wedge-tailed shearwater was analyzed by Ion Torrent sequencing using two universal polymerase chain reaction primers for the 12S and 16S mitochondrial DNA regions that targeted vertebrates and mollusks, respectively. The results of a BLAST search of obtained sequences detected 31 and 1 vertebrate and mollusk taxa, respectively. The results of the diet composition analysis showed that wedge-tailed shearwaters frequently consumed deep-sea fishes throughout the sampling season, indicating the importance of these fishes as a stable food resource. However, there was a marked seasonal shift in diet, which may reflect seasonal changes in food resource availability and wedge-tailed shearwater foraging behavior. The collected data regarding the shearwater diet may be useful for in situ conservation efforts. Future research that combines DNA metabarcoding with other tools, such as data logging, may provide further insight into the foraging ecology of pelagic seabirds.

Development of microsatellite markers for partially and putative fully mycoheterotrophic varieties of Pyrola japonica sensu lato (Ericaceae)

We developed microsatellite markers to compare the genetic variation and reproductive biology between the partially mycoheterotrophic Pyrola japonica var. japonica and the putative fully mycoheterotrophic P. japonica var. subaphylla. Fifteen primer pairs were developed for P. japonica sensu lato and they were tested on 77 ramets from three populations of the two varieties. Thirteen loci were polymorphic in at least one of the two var. japonica populations, whereas only four loci were polymorphic in the var. subaphylla population. The considerably lower genetic variation of the var. subaphylla population may be attributed to frequent selfing and/or inbreeding. The markers developed in this study will be useful for comparing the genetic diversity of P. japonica s. l. populations and measuring gene flow within and between populations and varieties.

Effects of pex1 disruption on wood lignin biodegradation, fruiting development and the utilization of carbon sources in the white-rot Agaricomycete Pleurotus ostreatus and non-wood decaying Coprinopsis cinerea

Peroxisomes are well-known organelles that are present in most eukaryotic organisms. Mutant phenotypes caused by the malfunction of peroxisomes have been shown in many fungi. However, these have never been investigated in Agaricomycetes, which include white-rot fungi that degrade wood lignin in nature almost exclusively and play an important role in the global carbon cycle. Based on the results of a forward genetics study to identify mutations causing defects in the ligninolytic activity of the white-rot Agaricomycete Pleurotus ostreatus, we report phenotypes of pex1 disruptants in P. ostreatus, which are defective in two major features of white-rot Agaricomycetes: lignin biodegradation and mushroom formation. Pex1 disruption was also shown to cause defects in the hyphal growth of P. ostreatus on certain sawdust and minimum media. We also demonstrated that pex1 is essential for fruiting initiation in the non-wood decaying Agaricomycete Coprinopsis cinerea. However, unlike P. ostreatus, significant defects in hyphal growth on the aforementioned agar medium were not observed in C. cinerea. This result, together with previous C. cinerea genetic studies, suggests that the regulation mechanisms for the utilization of carbon sources are altered during the evolution of Agaricomycetes or Agaricales.

Application of a simplified method of chloroplast enrichment to small amounts of tissue for chloroplast genome sequencing

PREMISE OF THE STUDY

High-throughput sequencing of genomic DNA can recover complete chloroplast genome sequences, but the sequence data are usually dominated by sequences from nuclear/mitochondrial genomes. To overcome this deficiency, a simple enrichment method for chloroplast DNA from small amounts of plant tissue was tested for eight plant species including a gymnosperm and various angiosperms.

METHODS

Chloroplasts were enriched using a high-salt isolation buffer without any step gradient procedures, and enriched chloroplast DNA was sequenced by multiplexed high-throughput sequencing.

RESULTS

Using this simple method, significant enrichment of chloroplast DNA-derived reads was attained, allowing deep sequencing of chloroplast genomes. As an example, the chloroplast genome of the conifer Callitris sulcata was assembled, from which polymorphic microsatellite loci were isolated successfully.

DISCUSSION

This chloroplast enrichment method from small amounts of plant tissue will be particularly useful for studies that use sequencers with relatively small throughput and that cannot use large amounts of tissue (e.g., for endangered species).

Identification of two mutations that cause defects in the ligninolytic system through an efficient forward genetics in the white-rot agaricomycete Pleurotus ostreatus

White-rot fungi play an important role in the global carbon cycle because they are the species that almost exclusively biodegrade wood lignin in nature. Lignin peroxidases (LiPs), manganese peroxidases (MnPs) and versatile peroxidases (VPs) are considered key players in the ligninolytic system. Apart from LiPs, MnPs and VPs, however, only few other factors involved in the ligninolytic system have been investigated using molecular genetics, implying the existence of unidentified elements. By combining classical genetic techniques with next-generation sequencing technology, they successfully showed an efficient forward genetics approach to identify mutations causing defects in the ligninolytic system of the white-rot fungus Pleurotus ostreatus. In this study, they identified two genes - chd1 and wtr1 - mutations in which cause an almost complete loss of Mn²⁺ -dependent peroxidase activity. The chd1 gene encodes a putative chromatin modifier, and wtr1 encodes an agaricomycete-specific protein with a putative DNA-binding domain. The chd1-1 mutation and targeted disruption of wtr1 hamper the ability of P. ostreatus to biodegrade wood lignin. Examination of the effects of the aforementioned mutation and disruption on the expression of certain MnP/VP genes suggests that a complex mechanism underlies the ligninolytic system in P. ostreatus.

Chloroplast DNA sequencing and detailed microsatellite genotyping of all remnant populations suggests that only a single genet survives of the critically endangered plant Rehmannia japonica

Rehmannia japonica (Thunb.) Makino ex T. Yamaz. is an endangered perennial herb species in Japan. Although earlier the Japanese considered it a variety of R. glutinosa, recent Japanese taxonomists have consistently regarded it as an independent species. According to the historical literature, Rehmannia japonica seems to have been known in China and Japan in the past. However, Chinese taxonomists do not recognize R. japonica at present. In Japan, only two populations are known, and although these populations flower every year, seed reproduction has not been observed. In this study, we aimed to reveal the phylogenetic relationships and levels of genetic diversity of R. japonica. A haplotype network based on two chloroplast DNA regions (trnL-trnF and rps16) showed that the sequences of R. japonica were distinguishable by three or four sites of indels from the most closely related species, R. chingii, consistent with the separate species status of R. japonica. An analysis of genetic diversity using twelve microsatellite loci showed that all of the ramets of R. japonica collected from two geographically isolated populations had an identical multilocus genotype, including identical heterozygous genotypes at six loci. This result indicated asexual origin of all sampled ramets. This study also suggests that the absence of sexual reproduction of R. japonica is explained by self-incompatibility combined with only a single genet remaining in the R. japonica populations.

Genetic diversity and population structure of Nuphar submersa (Nymphaeaceae), a critically endangered aquatic plant endemic to Japan, and implications for its conservation

Nuphar submersa (Nymphaeaceae) is a critically endangered freshwater macrophyte indigenous to central Japan, with only four small extant populations represented across its entire range. We investigated the genotypic and genetic diversity as well as the genetic structure of all extant individuals of N. submersa based on analysis of 15 microsatellite loci. Among 278 individual ramets, 52 multilocus genotypes were detected: 30 genotypes in Nikko City (NIK), 18 in Nasukarasuyama City (NAS), 3 in Mooka City (MOK), and 1 in Sakura City (SAK). The average number of alleles per locus ranged from 1.20 to 1.93, whereas the observed and expected heterozygosities ranged from 0.11 to 0.33 and from 0.10 to 0.24, respectively. With the exception of SAK, all populations contained multiple clones, but our results indicated low levels of within-population genetic diversity. The populations NIK and NAS comprised few large or middle-sized genets and many small genets. The populations NIK and NAS were suggested to comprise large old, old fragmented, and/or young small genets resulting from seedling establishment. All four populations were differentiated, and gene flow between the populations was restricted (average level of gene flow (Nm) = 0.122, G' _ST  = 0.639). Of the total genetic diversity, 67.20 and 9.13% were attributable to inter- and intra-population diversity, respectively. STRUCTURE analysis revealed two or three well-differentiated groups of populations. Cluster I comprised one population (NIK) and cluster II comprised the remaining populations at K = 2. The populations NIK, NAS, and the remaining populations were assigned to clusters I, II, and III, respectively, at K = 3. For conservation practices, we recommend that each cluster be regarded as a different management unit. We further suggest that artificial gene flow among MOK and SAK populations is an appropriate option, whereas NIK should not be reinforced with genotypes from the remaining populations.

Microsatellite markers for Nuphar japonica (Nymphaeaceae), an aquatic plant in the agricultural ecosystem of Japan

PREMISE OF THE STUDY

Nuphar species (Nymphaeaceae) are representative aquatic plants in irrigation ponds in Japanese agricultural ecosystems. We developed 15 polymorphic microsatellite markers for N. japonica and confirmed their utility for its close relatives N. oguraensis var. akiensis and N. ×saijoensis, which originated from natural hybridization between N. japonica and N. oguraensis.

METHODS AND RESULTS

Genetic variation was characterized in 15 polymorphic loci in three populations of N. japonica. The average number of alleles per locus was 3.47 (range = 2-9; n = 32), and the average expected heterozygosity per locus was 0.84 (range = 0.5-1.0); 11 loci were amplified in N. oguraensis var. akiensis and 15 in N. ×saijoensis.

CONCLUSIONS

The polymorphic microsatellite markers developed in this study will be useful for investigating the levels of genetic diversity within remnant populations of Nuphar taxa and could provide a valuable tool for conservation genetics of these taxa.

Vertical Structure of Phyllosphere Fungal Communities in a Tropical Forest in Thailand Uncovered by High-Throughput Sequencing

Phyllosphere fungi harbor a tremendous species diversity and play important ecological roles. However, little is known about their distribution patterns within forest ecosystems. We examined how species diversity and community composition of phyllosphere fungi change along a vertical structure in a tropical forest in Thailand. Fungal communities in 144 leaf samples from 19 vertical layers (1.28-34.4 m above ground) of 73 plant individuals (27 species) were investigated by metabarcoding analysis using Ion Torrent sequencing. In total, 1,524 fungal operational taxonomic units (OTUs) were detected among 890,710 reads obtained from the 144 leaf samples. Taxonomically diverse fungi belonging to as many as 24 orders of Ascomycota and 21 orders of Basidiomycota were detected, most of which inhabited limited parts of the lowest layers closest to the forest floor. Species diversity of phyllosphere fungi was the highest in the lowest layers closest to the forest floor, decreased with increasing height, and lowest in the canopy; 742 and 55 fungal OTUs were detected at the lowest and highest layer, respectively. On the layers close to the forest floor, phyllosphere fungal communities were mainly composed of low frequency OTUs and largely differentiated among plant individuals. Conversely, in the canopy, fungal communities consisted of similar OTUs across plant individuals, and as many as 86.1%-92.7% of the OTUs found in the canopy (≥22 m above ground) were also distributed in the lower layers. Overall, our study showed the variability of phyllosphere fungal communities along the vertical gradient of plant vegetation and environmental conditions, suggesting the significance of biotic and abiotic variation for the species diversity of phyllosphere fungi.

Complex pollination of a tropical Asian rainforest canopy tree by flower-feeding thrips and thrips-feeding predators

PREMISE OF THE STUDY

In tropical rainforests of Southeast Asia, a highly fecund thrips (Thrips spp.) responds rapidly to the mass flowering at multiple-year intervals characteristic of certain species such as the canopy tree studied here, Shorea acuminata, by feeding on flower resources. However, past DNA analyses of pollen adherent to thrips bodies revealed that the thrips promoted a very high level of self-pollination. Here, we identified the pollinator that contributes to cross-pollination and discuss ways that the pollination system has adapted to mass flowering.

METHODS

By comparing the patterns of floral visitation and levels of genetic diversity in adherent pollen loads among floral visitors, we evaluated the contribution of each flower visitor to pollination.

KEY RESULTS

The big-eyed bug, Geocoris sp., a major thrips predator, was an inadvertent pollinator, and importantly contributed to cross-pollination. The total outcross pollen adhering to thrips was approximately 30% that on the big-eyed bugs. Similarly, 63% of alleles examined in S. acuminata seeds and seedlings occurred in pollen adhering to big-eyed bugs; about 30% was shared with pollen from thrips.

CONCLUSIONS

During mass flowering, big-eyed bugs likely travel among flowering S. acuminata trees, attracted by the abundant thrips. Floral visitation patterns of big-eyed bugs vs. other insects suggest that these bugs can maintain their population size between flowering by preying upon another thrips (Haplothrips sp.) that inhabits stipules of S. acuminata throughout the year and quickly respond to mass flowering. Thus, thrips and big-eyed bugs are essential components in the pollination of S. acuminata.

Isolation and characterization of microsatellite loci for the endangered wetland plant Adenophora palustris (Campanulaceae)

PREMISE OF THE STUDY

Adenophora palustris (Campanulaceae) is an endangered wetland plant species in Japan. Although it is widely distributed in East Asia, only six extant populations are known in Japan, with fewer than 1000 individuals in total. We developed 15 microsatellite markers for this species and confirmed their utility for the closely related species A. triphylla var. japonica.

METHODS AND RESULTS

Ten polymorphic loci were characterized for genetic variation within three populations of A. palustris. The number of alleles per locus ranged from four to 15, with an average of 9.3; the expected heterozygosity ranged from 0.48 to 0.89, with an average of 0.74. Nine loci were successfully amplified in A. triphylla var. japonica, and three of these loci showed polymorphism.

CONCLUSIONS

These markers are useful for investigating genetic diversity and gene flow within and among remnant populations of A. palustris in Japan, and the results will provide crucial information for conservation.