||Natsuki Moriguchi, Kentaro Uchiyama, Ryutaro Miyagi, Etsuko Moritsuka, Aya Takahashi, Koichiro Tamura, Yoshihiko Tsumura, Kosuke M. Teshima, Hidenori Tachida, Junko Kusumi, Inferring the demographic history of Japanese cedar, Cryptomeria japonica, using amplicon sequencing, Heredity, 10.1038/s41437-019-0198-y, 2019.01, The evolution of a species depends on multiple forces, such as demography and natural selection. To understand the trajectory and driving forces of evolution of a target species, it is first necessary to uncover that species’ population history, such as past and present population sizes, subdivision and gene flow, by using appropriate genetic markers. Cryptomeria japonica is a long-lived monoecious conifer species that is distributed in Japan. There are two main lines (omote-sugi and ura-sugi), which are distinguished by apparent differences in morphological traits that may have contributed to their local adaptation. The evolution of these morphological traits seems to be related to past climatic changes in East Asia, but no precise estimate is available for the divergence time of these two lines and the subsequent population dynamics in this species. Here, we analyzed the nucleotide variations at 120 nuclear genes in 94 individuals by using amplicon sequencing in combination with high-throughput sequencing technologies. Our analysis indicated that the population on Yakushima Island, the southern distribution limit of C. japonica in Japan, diverged from the other populations 0.85 million years ago (MYA). The divergence time of the other populations on mainland Japan was estimated to be 0.32 MYA suggesting that the divergence of omote-sugi and ura-sugi might have occurred before the last glacial maximum. Although we found modest levels of gene flow between the present populations, the long-term isolation and environmental heterogeneity caused by climatic changes might have contributed to the differentiation of the lines and their local adaptation..
||Junko Kusumi, Motoshi Ichinose, Masaru Iizuka, Effects of gene duplication, epistasis, recombination and gene conversion on the fixation time of compensatory mutations, Journal of Theoretical Biology, 10.1016/j.jtbi.2019.02.001, 467, 134-141, 2019.04, Gene duplication is one of the major mechanisms of molecular evolution. Gene duplication enables copies of a gene to accumulate mutations through functional redundancy. If a gene encodes a specific protein that interacts with other proteins, RNA, or DNA, the relaxation of selective constraints caused by gene duplication might contribute to the fixation of compensatory mutations that occur at the interacting sites. In this study, we investigate the effect of gene duplication, epistasis among the duplicated copies and gene conversion on the fixation time of compensatory mutations by extending the original model of compensatory evolution proposed by Kimura in 1985. Our simulation results reveal that the time to fixation of compensatory mutations can be decreased remarkably by gene duplication if one of the duplicated loci can completely mask the deleterious effects of a mutation that occurs at the other locus. Conversely, the fixation time can be increased by gene duplication if such functional compensation is weak. We also show that the combination of the degree of functional compensation and the rate of gene conversion between duplicate loci have contrasting effects on the time to fixation of compensatory mutations..
||Yuka Ikezaki, Yoshihisa Suyama, Beth A. Middleton, Yoshihiko Tsumura, Kosuke Teshima, Hidenori Tachida, Junko Kusumi, Inferences of population structure and demographic history for Taxodium distichum, a coniferous tree in North America, based on amplicon sequencing analysis, American journal of Botany, 10.3732/ajb.1600046, 103, 11, 2016.12, PREMISE OF THE STUDY: Studies of natural genetic variation can elucidate the genetic basis of phenotypic variation and the past population structure of species. Our study species, Taxodium distichum, is a unique conifer that inhabits the ood plains and swamps of North America. Morphological and eco- logical di erences in two varieties, T. distichum var. distichum (bald cypress) and T. distichum var. imbricarium (pond cypress), are well known, but little is known about the level of genetic di erentiation between the varieties and the demographic history of local populations.
METHODS: We analyzed nucleotide polymorphisms at 47 nuclear loci from 96 individuals collected from the Mississippi River Alluvial Valley (MRAV), and Gulf Coastal populations in Texas, Louisiana, and Florida using high-throughput DNA sequencing. Standard population genetic statistics were calculated, and demographic parameters were estimated using a composite-likelihood approach.
KEY RESULTS: Taxodium distichum in North America can be divided into at least three genetic groups, bald cypress in the MRAV and Texas, bald cypress in Florida, and pond cypress in Florida. The levels of genetic di erentiation among the groups were low but signi cant. Several loci showed the signatures of positive selection, which might be responsible for local adaptation or varietal di erentiation.
CONCLUSIONS: Bald cypress was genetically di erentiated into two geographical groups, and the boundary was located between the MRAV and Florida. This di erentiation could be explained by population expansion from east to west. Despite the overlap of the two varieties’ ranges, they were genetically di erentiated in Florida. The estimated demographic parameters suggested that pond cypress split from bald cypress during the late Miocene..
||Junko Kusumi, Motoshi Ichinose, Masasuke Takefu, Robert Piskol, Wolfgang Stephan, Masaru Iizuka, A model of compensatory molecular evolution involving multiple sites in RNA molecules, Journal of Theoretical Biology, doi:10.1016/j.jtbi.2015.10.008, 388, 96-107, 2016.01, Consider two sites under compensatory fitness interaction, such as a Watson–Crick base pair in an RNA helix or two interacting residues in a protein. A mutation at any one of these two sites may reduce the fitness of an individual. However, fitness may be restored by the occurrence of a second mutation at the other site. Kimura modeled this process using a two-locus haploid fitness scheme with two alleles at each locus. He predicted that compensatory evolution following this model is very rare unless selection against the deleterious single mutations is weak and linkage between the interacting sites is tight. Here we investigate the question whether the rate of compensatory evolution increases if we take the context of the two directly interacting sites into account. By “context”, we mean the effect of neighboring sites in an RNA helix. Interaction between the focal pair of sites under consideration and the context may lead to so-called indirect compensation. Thus, extending Kimura׳s classical model of compensatory evolution, we study the effects of both direct and indirect compensation on the rate of compensatory evolution. It is shown that the effects of indirect compensation are very strong. We find that recombination does not slow down the rate of compensatory evolution as predicted by the classical model. Instead, compensatory substitutions may be relatively frequent, even if linkage between the focal interacting sites is loose, selection against deleterious mutations is strong, and mutation rate is low. We compare our theoretical results with data on RNA secondary structures from vertebrate introns..
||Junko Kusumi, Yoshihiko Tsumura, Hidenori Tachida, Evolutionary rate variation in two conifer species, Taxodium distichum (L.) Rich. var. distichum (baldcypress) and Cryptomeria japonica (Thunb. ex L.f.) D. Don (Sugi, Japanese cedar), Genes & Genetic Systems, http://doi.org/10.1266/ggs.14-00079, 90, 305-315, 2015.12, With the advance of sequencing technologies, large-scale data of expressed sequence tags and full-length cDNA sequences have been reported for several conifer species. Comparative analyses of evolutionary rates among diverse taxa provide insights into taxon-specific molecular evolutionary features and into the origin of variation in evolutionary rates within genomes and between species. Here, we estimated evolutionary rates in two conifer species, Taxodium distichum and Cryptomeria japonica, to illuminate the molecular evolutionary features of these species, using hundreds of genes and employing Chamaecyparis obtusa as an outgroup. Our results show that the mutation rates based on synonymous substitution rates (dS) of T. distichum and C. japonica are approximately 0.67×10-9 and 0.59×10-9/site/year, respectively, which are 15-25 times lower than those of annual angiosperms. We found a significant positive correlation between dS and GC3. This implies that a local mutation bias, such as context dependency of the mutation bias, exists within the genomes of T. distichum and C. japonica, and/or that selection acts on synonymous sites in these species. In addition, the means of the ratios of synonymous to nonsynonymous substitution rate in the two species are almost the same, suggesting that the average intensity of functional constraint is constant between the lineages. Finally, we tested the possibility of positive selection based on the site model, and detected one candidate gene for positive selection..
||Junko Kusumi, Zhi-Hui Su, Isolation and characterization of 15 polymorphic microsatellite markers for the fig-pollinating wasp, Blastophaga nipponica (Hymenoptera: Agaonidae), Applied Entomology and Zoology, 10.1007/s13355-014-0267-x, 2014.04, We developed microsatellite markers for the fig-pollinating wasp Blastophaga nipponica Grandi using a dual-suppression-PCR technique. Twenty-one candidates of microsatellite loci were obtained, of which 15 yielded scorable patterns. The degree of polymorphism for the 15 loci was further characterized using summary statistics describing the genetic variation in 60 individuals from three natural populations in Japan. All 15 loci were polymorphic and yielded 2–27 alleles per locus. Overall observed heterozygosity (HO) and expected heterozygosity(HE) were 0.465 and 0.631, respectively. As expected,based on the inbreeding tendency of this species, the mean inbreeding coefficient (FIS) was high (= 0.255). These markers will contribute to studies on the population structure of this species..
||Motoshi Ichinose, Masaru Iizuka, Junko Kusumi, Masasuke Takefu, Models of compensatory molecular evolution: Effects of back mutation, JOURNAL OF THEORETICAL BIOLOGY, 10.1016/j.jtbi.2013.01.011, 323, 1-10, 2013.01, Compensatory mutations are individually deleterious but appropriate combinations of mutants are harmless. For several models of compensatory molecular evolution, we consider the effects of back mutation. It is shown that the effects of back mutation on the rate of compensatory molecular evolution are weak. Further we estimate the values of selection parameter of deleterious single mutants for the models of compensatory molecular evolution both with and without back mutation using sequence data of folded RNA molecules and compare them with previous results..
||Ayako Tanaka & Masato Ohtani & Yoshihisa Suyama &
Nobuyuki Inomata & Yoshihiko Tsumura &
Beth A. Middleton & Hidenori Tachida & Junko Kusumi, Population genetic structure of a widespread coniferous tree, Taxodium distichum [L.] Rich. (Cupressaceae), in the Mississippi River Alluvial Valley and Florida, Tree Genetics & Genomes, 2012.03, Studies of genetic variation can elucidate the structure
of present and past populations aswell as the genetic basis
of the phenotypic variability of species. Taxodium distichum is
a coniferous tree dominant in lowland river flood plains and
swamps of the southeastern USA which exhibits morphological
variability and adaption to stressful habitats. This study
provides a survey of the Mississippi River Alluvial Valley
(MAV) and Florida to elucidate their population structure and
the extent of genetic differentiation between the two regions
and sympatric varieties, including bald cypress (var. distichum)
and pond cypress (var. imbricatum). We determined
the genotypes of 12 simple sequence repeat loci totaling 444
adult individuals from 18 natural populations. Bayesian clustering
analysis revealed high levels of differentiation between
the MAV and the Florida regions. Within the MAV region,
there was a significant correlation between genetic and geographical
distances. In addition, we found that there was
almost no genetic differentiation between the varieties. Most
genetic variation was found within individuals (76.73 %),
1.67 % among individuals within population, 15.36 % among
populations within the regions, and 9.23 % between regions
within the variety. Our results suggest that (1) the populations
of theMAVand the Florida regions are divided into two major
genetic groups, which might originate from different glacial
refugia, and (2) the patterns of genetic differentiation and
phenotypic differentiation were not parallel in this species..
||Junko Kusumi; Hiroshi Azuma; Hsy-Yu Tzeng; Lien-Siang Chou; Yan-Qiong Peng;
Keiko Nakamura; Zhi-Hui Su, Phylogenetic analyses suggest a hybrid origin of the figs (Moraceae: Ficus) that are endemic to
the Ogasawara (Bonin) Islands, Japan, Molecular Phylogenetics and Evolution, http://dx.doi.org/10.1016/j.bbr.2011.03.031, 63, １, 168-179, 2012.01, The Ogasawara Islands are oceanic islands and harbor a unique endemic flora. There are
three fig species (Ficus boninsimae, F. nishimurae and F. iidaiana) endemic to the Ogasawara Islands,
and these species have been considered to be closely related to Ficus erecta, and to have diverged
within the islands. However, this hypothesis remains uncertain. To investigate this issue, we assessed
the phylogenetic relationships of the Ogasawara figs and their close relatives occurring in Japan,
Taiwan and South China based on six plastid genome regions, nuclear ITS region and two nuclear
genes. The plastid genome-based tree indicated a close relationship between the Ogasawara figs and F.
erecta, whereas some of the nuclear gene-based trees suggested this relationship was not so close. In
addition, the phylogenetic analyses of the pollinating wasps associated with these fig species based on
the nuclear 28S rRNA and mitochondrial cytB genes suggested that the fig-pollinating wasps of F.
erecta are not sister to those of the Ogasawara figs. These results suggest the occurrence of an early
hybridization event(s) in the lineage leading to the Ogasawara figs..
||Junko Kusumi, Zidong, Li, Tomoyuki Kado, Yoshihiko Tsumura, Beth A. Middleton and Hidenori Tachida, Multilocus patterns of nucleotide polymorphism and demographic change in Taxodium distichum (Cupressaceae) in the lower Mississippi River alluvial valley, American journal of botany, 10.3732/ajb.1000082, 97, 11, 1848-1857, 2010.11.
||Junko Kusumi、Aya Sato、Hidenori Tachida, Relaxation of Functional Constraint on Light-Independent Protochlorophyllide Oxidoreductase in Thuja, Molecular Biology and Evolution, 10.1093/molbev/msj097, 23, 5, 941-948, Relaxation of Functional Constraint on Light-Independent Protochlorophyllide Oxidoreductase in Thuja, 2006.05.
||Junko Kusumi and Hidenori Tachida, Compositional properties of green-plant plastid genomes., Journal of Molecular evolution, 10.1007/s00239-004-0086-8, 60, 4, 417-425, 2005.04.
||Junko Kusumi, Yoshihiko Tsumura, Hiroshi Yoshimaru and Hidenori Tachida, Molecular evolution of nuclear genes in Cupressaceae, a group of conifer trees., Molecular biology and evolution, 19, 5, 736-747, 2002.10.