Research Institutue of Evolutionary Biology

Butterflies are a diverse and charismatic insect group that are thought to have evolved with plants and dispersed throughout the world in response to key geological events. However, these hypotheses have not been extensively tested because a comprehensive phylogenetic framework and datasets for butterfly larval hosts and global distributions are lacking. We sequenced 391 genes from nearly 2,300 butterfly species, sampled from 90 countries and 28 specimen collections, to reconstruct a new phylogenomic tree of butterflies representing 92% of all genera. Our phylogeny has strong support for nearly all nodes and demonstrates that at least 36 butterfly tribes require reclassification. Divergence time analyses imply an origin ~100 million years ago for butterflies and indicate that all but one family were present before the K/Pg extinction event. We aggregated larval host datasets and global distribution records and found that butterflies are likely to have first fed on Fabaceae and originated in what is now the Americas. Soon after the Cretaceous Thermal Maximum, butterflies crossed Beringia and diversified in the Palaeotropics. Our results also reveal that most butterfly species are specialists that feed on only one larval host plant family. However, generalist butterflies that consume two or more plant families usually feed on closely related plants.

In this study, we analysed the mitochondrial DNA D-loop region of Japanese native chickens to clarify their phylogenetic relationships, possible maternal origin and routes of introduction into Japan. Seven haplogroups (Types A-G) were identified. Types A-C were observed in Jidori, Shokoku and related breeds. However, Type C was absent in Shokoku, which was introduced from China, while most Indonesian native chickens were included in the Type C haplogroup. Types D-G were observed in Shamo and related breeds. Type E had a close genetic relationship with Chinese native chickens. Our results indicate that some breeds were not introduced into Japan as suggested in conventional literature, based on low nucleotide diversity of certain chicken breeds. Sequences originating from China and Korea could be clearly distinguished from those originating from Southeast Asia. In each group, domestic chickens were divided into the Jidori-Shokoku and Shamo groups. These results indicate that Chinese and Korean chickens were derived from Southeast Asia. Following the domestication of red junglefowl, a non-game type chicken was developed, and it spread to China. A game type chicken was developed in each area. Both non-game and game chickens formed the foundation of Japanese native chickens.

Blocking characteristics of diamond junction field-effect transistors were evaluated at room temperature (RT) and 200 <sup xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">°</sup> C. A high source-drain bias (breakdown voltage) of 566 V was recorded at RT, whereas it increased to 608 V at 200 <sup xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">°</sup> C. The positive temperature coefficient of the breakdown voltage indicates the avalanche breakdown of the device. We found that the breakdown occurred at the drain edge of the p-n junction between p-channel and the n <sup xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">+</sup> -gates. All four devices measured in this letter showed a maximum gate-drain bias over 500 V at RT and 600 V at 200 <sup xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">°</sup> C.

Populations are shaped by their history. It is crucial to interpret population structure in an evolutionary context. Pairwise FST measures population structure, whereas population-specific FST measures deviation from the ancestral population. To understand the current population structure and a population's history of range expansion, we propose a representation method that overlays population-specific FST estimates on a sampling location map, and on an unrooted neighbor-joining tree and a multi-dimensional scaling plot inferred from a pairwise FST distance matrix. We examined the usefulness of our procedure using simulations that mimicked population colonization from an ancestral population and by analyzing published human, Atlantic cod, and wild poplar data. Our results demonstrated that population-specific FST values identify the source population and trace the evolutionary history of its derived populations. Conversely, pairwise FST values represent the current population structure. By integrating the results of both estimators, we obtained a new picture of the population structure that incorporates evolutionary history. The generalized least squares estimate of genome-wide population-specific FST indicated that the wild poplar population expanded its distribution to the north, where daylight hours are long in summer, to coastal areas with abundant rainfall, and to the south where summers are dry. Genomic data highlight the power of the bias-corrected moment estimators of FST, whether global, pairwise, or population-specific, that provide unbiased estimates of FST. All FST moment estimators described in this paper have reasonable processing times and are useful in population genomics studies.

Renaming schizophrenia is a potential strategy to reduce the stigma attached to people with schizophrenia. However, the overall associations between renaming schizophrenia and stigma-related outcomes have not been fully elucidated. We conducted a systematic review of studies that empirically examined the outcomes between new or alternative terms and old or existing terms for schizophrenia. We searched for relevant articles in eight bibliographic databases, conducted a Google search, examined reference lists, and contacted relevant experts. We found a total of 2601 reference records, and 23 articles were included in this review. Overall, in countries where schizophrenia has been renamed, the name changes may be associated with improvements in adults' attitudes toward people with schizophrenia, and with increased diagnosis announcement. However, studies conducted in countries where schizophrenia has not been renamed report inconsistent findings. In addition, renaming may not influence portrayals of schizophrenia in the media. Most studies included in our review had a risk of bias in their methodology, and we employed a vote-counting method to synthesize study results; therefore, the impacts of renaming are still inconclusive. Future studies are needed to address the following issues: use of univariate descriptive statistics, adjustment for confounding variables, use of reliable measures, and employing a question that addresses the image of split or multiple personalities. Evidence is limited regarding the associations between renaming and stigma experienced by both people with schizophrenia and their families (e.g., perceived stigma, self-stigma, discrimination experience, and burden). Further research in these populations is needed to confirm the effects of renaming schizophrenia.

Humans divide themselves up into groups based on a shared cultural identity and common descent. Culturally inherited differences in dress, language, and institutions are often used as symbolic markers of the boundaries of these ethnic groups. Relatively little is known about the function of such ethnic groups, and why ethnic diversity is high in some regions yet lower in others. In this paper, we demonstrate how investigating the spatial distribution of ethnolinguistic groups can reveal the factors that affect the origin and maintenance of human ethnic group diversity. Here we describe the use of a Geographic Information System to construct a large database that integrates information about languages with a number of environmental, ecological, and ethnographic variables. Using these data on the spatial distribution of ethnolinguistic groups, we employ a hierarchical linear modeling approach to test a variety of hypotheses concerning the function of such groups. Despite revealing intriguing spatial patterns such as the latitudinal gradient in ethnolinguistic diversity, previous analyses suggested that the direct effects of environmental variables on the distribution of ethnolinguistic groups were in fact quite small. Here we show that the strength of the relationship between ethnolinguistic area and environmental variables is stronger in societies whose primary mode of subsistence is foraging. We then go on to demonstrate this same finding using the estimated native distributions of ethnolinguistic group in the Americas and Australia. In particular, Net Primary Productivity is shown to be a good predictor of the area covered by ethnolinguistic groups in foragers but not in agriculturalists. This provides support for the idea that the factors affecting ethnic diversity have changed in a systematic way with changes in subsistence strategies and social organization. We highlight future avenues for spatially explicit investigations of the evolution of ethnic diversity, and suggest that the evolutionary ecological approach adopted here may provide important insights into processes affecting ethnic diversity in the modern world.

OBJECTIVE: Disruptions in the dopamine system have been observed in psychiatric disorders. Since dopamine is mainly produced in the ventral tegmental area (VTA), elucidating the differences in the VTA neural network across psychiatric disorders would facilitate a greater understanding of the pathophysiological mechanisms underlying these disorders. However, no study has compared VTA-seed-based functional connectivity across psychiatric disorders. Therefore, we conducted a resting-state functional magnetic resonance imaging (rs-fMRI) study to perform a seed-based fMRI analysis, using the VTA as a seed. METHODS: We included participants with major depressive disorder (MDD; n = 45), schizophrenia (n = 32), and bipolar disorder (BPD; n = 30), along with healthy control participants (n = 46) who were matched for age, gender, and handedness. RESULTS: The results showed that patients with MDD and BPD had altered VTA-related connectivity in the superior frontal gyrus, frontal pole regions, hippocampus, cerebellum, and posterior cingulate cortex. Some of these differences in connectivity were also found between affective disorders and schizophrenia; however, there were no differences between the schizophrenia and control groups. Connectivity between the VTA and the hippocampus was correlated with positive symptoms in the schizophrenia group. The connectivity was not associated with medication dose, and the results remained significant after controlling for dose. CONCLUSIONS: The results suggest that altered brain functional connectivity related to VTA networks could be associated with the distinctive pathophysiologies of psychiatric disorders, especially affective disorders.

There exists diversity among individuals in difficulty controlling body weight. Body weight control, or obesity prevention, requires cognitive control over ingestive behavior, which may account for the diverse ability of body weight control. The caudate nuclei, especially the dorsal area, have been shown to play critical roles in ingestive behaviors, which significantly influences body weight control. However, the practice of body weight control is dependent on the body weight status, because the current obesity status determines the need for body weight control. To elucidate the underlying neural mechanism that accounts for individual differences in obesity prevention, we aimed to isolate functional caudate connectivity responsible for the underlying tendency of obesity prevention, independent of the current obesity status, using resting state fMRI data, body mass index (BMI), and assessment of ingestive behavior from 185 individuals from the NKI-Rockland sample. The underlying tendency of obesity prevention was estimated from BMI and behavioral and cognitive components of food intake. Functional connectivities between the caudate head and the whole brain were tested as a function of the estimated tendency in a voxel-wise manner. The bilateral precuneus showed inverse association between its connectivity to the caudate and the estimated tendency. Caudate-precuneus connectivity may have significant implications to understanding personal differences that accounts for the success in body weight control.

A series of emotional events successively occur in temporal context. The present study investigated how physiological and psychological responses are modulated by emotional context. Skin conductance response (SCR), heart rate, corrugator activity, zygomatic activity, and subjective feelings during emotional picture viewing were measured. To create an emotional context, a unpleasant or pleasant picture was preceded by three types of pictures, i.e., unpleasant, pleasant, and neutral pictures, resulting in six pairings. The results showed that viewing an unpleasant picture attenuated pleasant feelings induced by the following pleasant picture. On the other hand, preceding pleasant pictures decreased SCR to the following pictures. The effects of contextual modulation on emotional responses might be due to the informative function of pre-existing feelings; unpleasant feelings signal a threatening environment, whereas pleasant feelings signal a benign environment. With respect to facial muscle activities, viewing a pleasant picture decreased corrugator activity in response to the preceding picture. These findings suggest several types of contextual modulation effects on psychological, autonomic, and somatic responses to emotional stimuli.

It has been hypothesized that a higher genetic risk of bipolar disorder (BD) is associated with greater creativity. Given the clinical importance of bipolar disorder and the importance of creativity to human society and cultural development, it is essential to reveal their associations and the neural basis of the genetic risk of bipolar disorder to gain insight into its etiology. However, despite the previous demonstration of the associations of polygenic risk score (PRS) of BD and creative jobs, the associations of BD-PRS and creativity measured by the divergent thinking (CMDT) and regional gray matter volume (rGMV) as well as regional white matter volume (rWMV) have not been investigated. Using psychological analyses and whole-brain voxel-by-voxel analyses, we examined these potential associations in 1558 young, typically developing adult students. After adjusting for confounding variables and multiple comparisons, a greater BD-PRS was associated with a greater total CMDT fluency score, and a significant relationship was found in fluency subscores. A greater BD-PRS was also associated with lower total mood disturbance. Neuroimaging analyses revealed that the BD-PRS was associated with greater rGMV in the right inferior frontal gyrus, which is a consistently affected area in BD, as well as a greater rWMV in the left middle frontal gyrus, which has been suggested to play a central role in the increased creativity associated with the risk of BD with creativity. These findings suggest a relationship between the genetic risk of BD and CMDT and prefrontal cortical structures among young educated individuals.

OBJECTIVES: Hybridization between rhesus (Macaca mulatta) and long-tailed (M. fascicularis) macaques has become a focal point of interest. The majority of such studies have evaluated their genetics, but not their morphological characters. MATERIALS AND METHODS: We analyzed morphological characters of eight free-ranging populations of Indochinese rhesus and long-tailed macaques distributed at the proposed hybrid zone (15.75-21.58° N) in comparison with one population each of Chinese and Indian-derived rhesus macaques and three populations of Sundaic long-tailed macaques. RESULTS: Chinese and Indian rhesus macaques had a heavier body mass, longer crown-rump length, shorter relative facial length and relative tail length, and a greater contrast of reddish and yellowish dorsal pelage color than the Sundaic long-tailed macaques for which the latter three parameters could be used to visually discriminate between the two species. Although the morphological characters of Indochinese rhesus and long-tailed macaques were intermediate between the Chinese/Indian rhesus and Sundaic long-tailed macaques, they were more similar to their respective conspecifics. The species-specific characters of a shorter tail (<70%) and a bipartite pelage color pattern were retained in the Indochinese rhesus macaques while the longer tail (>90%) and no bipartite pattern was found in the Indochinese long-tailed macaques. No morphological cline was observed across the species and the variations were abrupt to some extent. DISCUSSION: The hybridization between rhesus and long-tailed macaques may be results of multiple contacts and isolations over a long period of time, thus their evolutionary history should not be drawn solely by genetic or morphological analysis.

Increasing evidence shows that the midbrain dopaminergic system is involved in various functions. However, details of the role of the midbrain dopaminergic system in these functions are still to be determined in humans. Considering that the ventral tegmental area (VTA) and substantia nigra (SN) in the midbrain are the primary dopamine producers, creating reliable anatomical templates of the VTA and SN through neuroimaging studies would be useful for achieving a detailed understanding of this dopaminergic system. Although VTA and SN anatomical templates have been created, no specific templates exist for the Asian population. Thus, we conducted anatomical and resting-state functional magnetic resonance imaging (rs-fMRI) studies to create VTA and SN templates for the Asian population. First, a neuromelanin-sensitive MRI technique was used to visualize the VTA and SN, and then individual hand-drawn VTA and SN regions of interests (ROIs) were traced on a small sample of neuromelanin-sensitive MRIs (dataset 1). Second, individual hand-drawn VTA and SN ROIs were normalized to create normalized VTA and SN templates for the Asian population. Third, a seed-based functional connectivity analysis was performed on rs-fMRI data using hand-drawn ROIs to calculate neural networks of VTA and SN in dataset 1. Fourth, a seed-based functional connectivity analysis was performed using VTA and SN seeds that were created based on normalized templates from dataset 1. Subsequently, a seed-based functional connectivity analysis was performed using VTA and SN seeds in another, larger sample (dataset 2) to assess whether neural networks of VTA or SN seeds from dataset 1 would be replicated in dataset 2. The Asian VTA template was smaller and located in a more posterior and inferior part of the midbrain compared to the published VTA template, while the Asian SN template, relative to the published SN template, did not differ in size but was located in the more inferior part of the midbrain. The neural networks of the VTA and SN seeds in dataset 1 were replicated in dataset 2. Altogether, our normalized template of the VTA and SN could be used for measuring fMRI activities related to the VTA and SN in the Asian population.

Obesity causes a wide range of systemic diseases and is associated with mood and anxiety disorders. It is also associated with dopaminergic reward system function. However, the relationships between microstructural properties of the dopaminergic system and body mass index (BMI) have not been investigated. In this study, we investigated the associations of BMI with mean diffusivity (MD), diffusion tensor imaging measure in areas of the dopaminergic system (MDDS) in 435 healthy young adults with mild obesity and without obesity (BMI < 40). We detected the association between greater BMI and lower MD of the right globus pallidus and the right putamen. These results suggest that the property of the dopaminergic system is associated with BMI among young adults with mild obesity and without obesity.

Problem-solving can be facilitated with instructions or hints, which provide information about given problems. The proper amount of instruction that should be provided for learners is controversial. Research shows that tasks with intermediate difficulty induce the largest sense of accomplishment (SA), leading to an intrinsic motivation for learning. To investigate the effect of instructions, we prepared three instruction levels (No hint, Indirect hint, and Direct hint) for the same insight-problem types. We hypothesized that indirect instructions impose intermediate difficulty for each individual, thereby inducing the greatest SA per person. Based on previous neuroimaging studies that showed involvement of the bilateral caudate in learning and motivation, we expected SA to be processed in this reward system. We recruited twenty-one participants, and investigated neural activations during problem solving by functional magnetic resonance imaging (fMRI). We confirmed that the Indirect hint, which imposed intermediate difficulty, induced the largest SA among the three instruction types. Using fMRI, we showed that activations in the bilateral caudate and anterior cingulate cortex (ACC) were significantly modulated by SA. In the bilateral caudate, the indirect hint induced the largest activation, while the ACC seemed to reflect the difference between correct and incorrect trials. Importantly, such activation pattern was independent of notations (number or letter). Our results indicate that SA is represented in the reward system, and that the Indirect instruction effectively induces such sensation.

Abstract. Pupil diameter change is indicative of emotional processing. Most previous findings regarding pupillary response and emotion have reported that the pupil enlarges in response to the presentation of emotional perceptual stimuli (e.g., visual images) within several seconds. It is considered that such stimuli activate the sympathetic nervous system, leading to pupil dilation. In order to examine the effects of emotions similar to daily emotional experiences of mood, the present study examined pupil diameter changes and their relationships with subjective emotional changes while recalling a topic of stressful interpersonal events in daily life. The data of 20 university students (11 males, M age = 20.36 ± 2.38 years; 9 females, M age = 22.33 ± 3.57) were analyzed. In the experimental task, participants were instructed to recall their memories concerning the topic through instructions and questions presented on a monitor, which proceeded at their own pace, through a key press. Subsequently, after baseline and instruction periods, participants were instructed to freely recall their memories. They were then asked to respond silently to a series of questions concerning the freely recalled memories. In the analysis, we compared the pupil diameters between these different periods and observed that pupil diameters significantly decreased during the response period relative to the free recall or baseline periods. Furthermore, pupil constrictions during the response period were negatively correlated with increases in negative affect scale scores. Pupil constriction, which is indicative of decreased arousal level and parasympathetic activation, was presumably caused by multiple factors including less cognitive difficulty and a relatively long experimental task period. As the result of a less tonic mode in the response period, the attention of participants might be more successfully focused on ongoing tasks, which might lead to optimal performance in recalling memories, possibly leading to correlations between pupil diameter and negative emotional changes.

Repeated gene targeting by recombinase-mediated cassette exchange (RMCE) is an efficient tool for the study of gene function and regulation because of the high predictability and repeatability of gene expression. We have developed the site-directed integration (SDI) vector system for Agrobacterium-mediated transformation to precisely integrate a single copy of a desired gene into a predefined chromosomal locus in the absence of any coexpressed selection marker gene (Nanto et al. Plant Biotechnol J 3:203-214, 2005; Nanto and Ebinuma Transgenic Res 17:337-344, 2008; Nanto et al. Plant Cell Rep 28:777-785, 2009; Ebinuma and Nanto (2009) Marker-free targeted transformation, in Molecular techniques in crop improvement (2nd Edition). (Jain, S. M. and Brar, D. S. eds.), Springer Netherlands, pp. 527-543; Ebinuma and Nanto in preparation). The SDI vector system consists of a target vector to introduce the target cassette and an exchange vector to reintroduce the exchange cassette for gene replacement. We describe the molecular design and experimental protocol that can efficiently enrich RMCE events through the removal of randomly integrated copies and select clean marker-free targeted transgenic plants by using a negative marker.

Abstract During the process of range expansion, populations encounter a variety of environments. They respond to the local environments by modifying their mutually interacting traits. Common approaches of landscape analysis include first focusing on the genes that undergo diversifying selection or directional selection in response to environmental variation. To understand the whole history of populations, it is ideal to capture the history of their range expansion with reference to the series of surrounding environments and to infer the multitrait coadaptation. To this end, we propose a complementary approach; it is an exploratory analysis using up‐to‐date methods that integrate population genetic features and features of selection on multiple traits. First, we conduct correspondence analysis of site frequency spectra, traits, and environments with auxiliary information of population‐specific fixation index ( F ST ). This visualizes the structure and the ages of populations and helps infer the history of range expansion, encountered environmental changes, and selection on multiple traits. Next, we further investigate the inferred history using an admixture graph that describes the population split and admixture. Finally, principal component analysis of the selection on edge‐by‐trait (SET) matrix identifies multitrait coadaptation and the associated edges of the admixture graph. We introduce a newly defined factor loadings of environmental variables in order to identify the environmental factors that caused the coadaptation. A numerical simulation of one‐dimensional stepping‐stone population expansion showed that the exploratory analysis reconstructed the pattern of the environmental selection that was missed by analysis of individual traits. Analysis of a public dataset of natural populations of black cottonwood in northwestern America identified the first principal component (PC) coadaptation of photosynthesis‐ vs growth‐related traits responding to the geographical clines of temperature and day length. The second PC coadaptation of volume‐related traits suggested that soil condition was a limiting factor for aboveground environmental selection.

Based on emergence records from trap nests and naturally occurring nests, we determined the brood parasitism behaviour of Irenangelus punctipleuris (Pompilidae: Ceropalinae) for the first time and that of Dipogon (Nipponodipogon) nagasei (Pompilidae: Pepsinae) as a complement to a previous study. The former species was found to parasitise Auplopus carbonarius and A. pygialis (Pompilidae: Pepsinae), while the latter species was found to parasitise Di. (Deuteragenia) sperconsus (Pepsinae). When using the well-established system of categories of brood parasitism in Pompilidae, Irenangelus corresponds to type A brood parasitism ('highwayman') and Nipponodipogon corresponds to type B brood parasitism ('sneak thief'). In comparison with brood parasitism in birds, we propose new concepts, namely 'generalists' and 'specialists' of brood parasitism, which correspond to type A and type B brood parasites, respectively. We discuss the differences in behavioural patterns and morphological features between the two types. We also confirm the proposed concepts by their molecular phylogenetic positions.

Abstract Our motor system uses sensory feedback to keep desired performance. From this view, motor fluctuation is not simply ‘noise’ inevitably caused in the nervous system but would play a role in generating variations to explore better outcomes via sensory feedback. Vocalization system offers a good model for studying such sensory-motor interactions since we regulate vocalization by hearing our own voice. This behavior is typically observed as compensatory responses in vocalized pitch, or fundamental frequency ( f o ), when artificial f o shifts were induced in the auditory feedback. However, the relationship between adaptive regulation and motor exploration in vocalization has remained unclear. Here we investigated behavioral variability in spontaneous vocal f o and compensatory responses against f o shifts in the feedback, and demonstrated that larger spontaneous fluctuation correlates with greater compensation in vocal f o . This correlation was found in slow components (≤ 5 Hz) of the spontaneous fluctuation but not in fast components (between 6 and 30 Hz), and the slow one was amplified during the compensatory responses. Furthermore, the compensatory ratio was reduced when large f o shifts were applied to the auditory feedback, as if reflecting the range of motor exploration. All these findings consistently suggest the functional role of motor variability in the exploration of better vocal outcomes.

Abstract Genetic stock identification (GSI) is a major management tool of Pacific salmon ( Oncorhynchus Spp.) that has provided rich genetic baseline data of allozymes, microsatellites, and single nucleotide polymorphisms (SNPs) across the Pacific Rim. Here, we analyzed published data sets for adult chum salmon ( Oncorhynchus keta ), namely 10 microsatellites, 53 SNPs, and a mitochondrial DNA locus (mtDNA3, control region and NADH-3 combined) from 495 locations in the same distribution range ( n = 61,813). TreeMix analysis of the microsatellite loci identified the highest level of genetic drift towards Japanese/Korean populations and suggested two admixture events from Japan/Korea to Russia and the Alaskan Peninsula. The SNPs had been purposively collected from rapidly evolving genes to increase the power of GSI. The highest expected heterozygosity was observed in Japanese/Korean populations for microsatellites, whereas it was highest in Western Alaskan populations for SNPs, reflecting the SNP discovery process. By regressing the SNP population structures on those of the microsatellites, we estimated the selection on the SNP loci according to deviations from the predicted structures. Specifically, we matched the sampling locations of the SNPs with those of the microsatellites according to geographical information and performed regression analyses of SNP allele frequencies on the two coordinates of multi-dimensional scaling (MDS) of matched locations obtained from microsatellite pairwise F ST values. The MDS first axis indicated a latitudinal cline in American and Russian populations, whereas the second axis found a differentiation of Japanese/Korean populations. The top five outlier SNPs were mtDNA3 (combined locus of the control region and NADH-3), U502241 (unknown), GnRH373, ras1362, and TCP178, which were consistently identified by principal component analysis. We summarized the functions of the 53 nuclear SNPs and mtDNA3 locus by referring to a gene database system and discussed the functions of the outlier SNPs and fitness of chum salmon.