Chikako Matsuba

The utility of QTL-linked markers to detect selective sweeps in natural populations — a case study of the EDA gene and a linked marker in threespine stickleback

Sequence polymorphisms in coding genes and variability in quantitative trait loci (QTL)‐linked markers can be used to uncover the evolutionary mechanisms of traits involved in adaptive processes. We studied sequence variation in the EDA gene and allelic variation in 18 microsatellites — one of which (Gac4174) is linked with the EDA QTL — in low, partially and completely plated morphs from eight threespine stickleback European populations. The results agree with previous studies in that EDA polymorphism is closely related to plate number variation: EDA sequences grouped populations into low and completely plated morphs, whereas microsatellites failed to do so. Furthermore, partially plated fish were heterozygous with respect to the distinctive EDA alleles for completely and low plated morphs, indicating that completely plated morph alleles are not entirely dominant in controlling the expression of lateral plate number. An examination of population differentiation in plate number with quantitative genetic methods revealed that the degree of differentiation exceeded that expected from genetic drift alone (QST > FST). Our results support the adaptive genetic differentiation of plate morphs and the view that distinctive EDA gene polymorphism occurs in similar sites across the distribution range of this species. Yet, allele frequency differentiation in the Gac4174 microsatellite locus, informative in experimental crosses for plate number variation, did not differ from that of neutral markers and, was therefore unable to detect the signature of natural selection responsible for population divergence.

Sex reversal and primary sex ratios in the common frog (Rana temporaria)

Sex reversal has been suggested to have profound implications for the evolution of sex chromosomes and population dynamics in ectotherms. Occasional sex reversal of genetic males has been hypothesized to prevent the evolutionary decay of nonrecombining Y chromosomes caused by the accumulation of deleterious mutations. At the same time, sex reversals can have a negative effect on population growth rate. Here, we studied phenotypic and genotypic sex in the common frog (Rana temporaria) in a subarctic environment, where strongly female‐biased sex ratios have raised the possibility of frequent sex reversals. We developed two novel sex‐linked microsatellite markers for the species and used them with a third, existing marker and a Bayesian modelling approach to study the occurrence of sex reversal and to determine primary sex ratios in egg clutches. Our results show that a significant proportion (0.09, 95% credible interval: 0.04–0.18) of adults that were genetically female expressed the male phenotype, but there was no evidence of sex reversal of genetic males that is required for counteracting the degeneration of Y chromosome. The primary sex ratios were mostly equal, but three clutches consisted only of genetic females and three others had a significant female bias. Reproduction of the sex‐reversed genetic females appears to create all‐female clutches potentially skewing the population level adult sex‐ratio consistent with field observations. However, based on a simulation model, such a bias is expected to be small and transient and thus does not fully explain the observed female‐bias in the field.

Temperature, stress and spontaneous mutation in Caenorhabditis briggsae and Caenorhabditis elegans.

Mutation rate often increases with environmental temperature, but establishing causality is complicated. Asymmetry between physiological stress and deviation from the optimal temperature means that temperature and stress are often confounded. We allowed mutations to accumulate in two species of Caenorhabditis for approximately 100 generations at 18°C and for approximately 165 generations at 26°C; 26°C is stressful for Caenorhabditis elegans but not for Caenorhabditis briggsae. We report mutation rates at a set of microsatellite loci and estimates of the per-generation decay of fitness (ΔMw), the genomic mutation rate for fitness (U) and the average effect of a new mutation (E[a]), assayed at both temperatures. In C. elegans, the microsatellite mutation rate is significantly greater at 26°C than at 18°C whereas in C. briggsae there is only a slight, non-significant increase in mutation rate at 26°C, consistent with stress-dependent mutation in C. elegans. The fitness data from both species qualitatively reinforce the microsatellite results. The fitness results of C. elegans are potentially complicated by selection but also suggest temperature-dependent mutation; the difference between the two species suggests that physiological stress plays a significant role in the mutational process.

RECOMBINATION RATE BETWEEN SEX CHROMOSOMES DEPENDS ON PHENOTYPIC SEX IN THE COMMON FROG

We show that the recombination rate between the sex chromosomes is controlled by phenotypic, rather than genotypic, sex in sex‐reversed common frogs. This supports the recent hypothesis that in ectothermic vertebrates sex reversal can prevent the progressive accumulation of mutations to nonrecombining Y chromosomes and their subsequent evolutionary decay.

Isolation and characterization of 145 polymorphic microsatellite loci for the common frog (Rana temporaria).

We describe primers and polymerase chain reaction conditions to amplify 145 di‐, tri‐ and tetranucleotide microsatellite loci from the common frog (Rana temporaria), a species commonly used as a model in ecological and evolutionary research. Primers were tested on 46 individuals from two Fennoscandian populations and yielded an average of six to nine alleles per locus (range = 1–30) depending on the population. Average observed heterozygosities in the two populations were 0.16 (range = 0–0.91) and 0.36 (range = 0–1).

Sequence Variation in the Melanocortin-1 Receptor Gene (Mc1r) Does Not Explain Variation in the Degree of Melanism in a Widespread Amphibian

Variation in nucleotide sequence of the melanocortin-1 receptor gene (Mc1r) is associated with melanism in several mammalian, avian and reptilian species, but no attempts have been made to understand the genetic underpinnings of melanism in amphibians. We isolated the complete coding sequence (945 bps) of Mc1r from the common frog (Rana temporaria) and compared the predicted amino acid sequence with that of fish, reptiles, birds and mammals. We investigated associations between nucleotide substitutions and the level of dorsal melanism among 28 individuals from two populations with pronounced differences in melanism. According to our results, the transmembrane regions of Mc1r are conserved across vertebrates. In the population comparison, we only found five nucleotide sites with synonymous substitutions; none is being associated with the level of melanism. Our results suggest that either other genes or regulatory regions outside the coding sequence of Mc1r are responsible for expression of melanism in R. temporaria.

Invariance (?) of Mutational Parameters for Relative Fitness Over 400 Generations of Mutation Accumulation in Caenorhabditis elegans

Evidence is accumulating that individuals in poor physiologic condition may accumulate mutational damage faster than individuals in good condition. If poor condition results from pre-existing deleterious mutations, the result is “fitness-dependent mutation rate,” which has interesting theoretical implications. Here we report a study in which 10 mutation accumulation (MA) lines of the nematode Caenorhabditis elegans that had previously accumulated mutations for 250 generations under relaxed selection were expanded into sets of “second-order” MA lines and allowed to accumulate mutations for an additional 150 generations. The 10 lines were chosen on the basis of the relative change in fitness over the first 250 generations of MA, five high-fitness lines and five low-fitness lines. On average, the mutational properties (per-generation change in mean relative fitness, mutational variance, and Bateman-Mukai estimates of genomic mutation rate and average mutational effect) of the high-fitness and low-fitness did not differ significantly, and averaged over all lines, the point estimates were extremely close to those of the first-order MA experiment after 200 generations of MA. However, several nonsignificant trends indicate that low-fitness lines may in fact be more likely to suffer mutational damage than high-fitness lines.

The Red Death Meets the Abdominal Bristle: Polygenic Mutation for Susceptibility to a Bacterial Pathogen in Caenorhabditis elegans

Understanding the genetic basis of susceptibility to pathogens is an important goal of medicine and of evolutionary biology. A key first step toward understanding the genetics and evolution of any phenotypic trait is characterizing the role of mutation. However, the rate at which mutation introduces genetic variance for pathogen susceptibility in any organism is essentially unknown. Here, we quantify the per‐generation input of genetic variance by mutation (VM) for susceptibility of Caenorhabditis elegans to the pathogenic bacterium Pseudomonas aeruginosa (defined as the median time of death, LT50). VM for LT50 is slightly less than VM for a variety of life‐history and morphological traits in this strain of C. elegans, but is well within the range of reported values in a variety of organisms. Mean LT50 did not change significantly over 250 generations of mutation accumulation. Comparison of VM to the standing genetic variance (VG) implies a strength of selection against new mutations of a few tenths of a percent. These results suggest that the substantial standing genetic variation for susceptibility of C. elegans to P. aeruginosa can be explained by polygenic mutation coupled with purifying selection.

Tempo, mode, and fitness effects of mutation in Caenorhabditis elegans over 400 generations of minimal selection

The mutational process varies at many levels, from within genomes to among taxa. Many mechanisms have been linked to variation in mutation, but understanding of the evolution of the mutational process is rudimentary. Physiological condition is often implicated as a source of variation in microbial mutation rate and may contribute to mutation rate variation in multicellular organisms. Deleterious mutations are a ubiquitous source of variation in condition. We test the hypothesis that the mutational process depends on the underlying mutation load in two groups of Caenorhabditis elegans mutation accumulation (MA) lines that differ in their starting mutation loads. “First-Order MA” (O1MA) lines maintained under minimal selection for ~250 generations were divided into high-fitness and low-fitness groups and sets of “second-order MA” (O2MA) lines derived from each O1MA line were maintained for ~150 additional generations. Genomes of 48 O2MA lines and their O1MA progenitors were sequenced. There is significant variation among O2MA lines in base-substitution rate (μbs), but no effect of initial fitness, whereas the indel rate is greater in high-fitness O2MA lines. Overall, μbs is positively correlated with recombination and proximity to short tandem repeats and negatively correlated with 1 Kb GC content. However, multiple logistic regression shows mutability is sufficiently predicted by the three-nucleotide motif. ~90% of the variance in standing nucleotide variation is explained by mutability. Total mutation rate increased in the O2MA lines, as predicted by the “drift barrier” model of mutation rate evolution. These data, combined with experimental estimates of fitness, suggest that epistasis is synergistic.Abstract The mutational process varies at many levels, from within genomes to among taxa. Many mechanisms have been linked to variation in mutation, but understanding of the evolution of the mutational process is rudimentary. Physiological condition is often implicated as a source of variation in microbial mutation rate and may contribute to mutation rate variation in multicellular organisms. Deleterious mutations are a ubiquitous source of variation in condition. We test the hypothesis that the mutational process depends on the underlying mutation load in two groups of Caenorhabditis elegans mutation accumulation (MA) lines that differ in their starting mutation loads. “First-Order MA” (O1MA) lines maintained under minimal selection for ∼250 generations were divided into high-fitness and low-fitness groups and sets of “second-order MA” (O2MA) lines derived from each O1MA line were maintained for ∼150 additional generations. Genomes of 48 O2MA lines and their progenitors were sequenced. There is significant variation among O2MA lines in base-substitution rate (µbs), but no effect of initial fitness, whereas the indel rate is greater in high-fitness O2MA lines. Overall, µbs is positively correlated with recombination and proximity to short tandem repeats and negatively correlated with 10 bp and 1 Kb GC content. However, probability of mutation is well-predicted by the three-nucleotide motif. ∼90% of the variance in standing nucleotide variation is explained by mutability. Total mutation rate increased in the O2MA lines, as predicted by the “drift barrier” model of mutation rate evolution. These data, combined with experimental estimates of fitness, suggest that epistasis is synergistic.

The Epigenomic Landscape of Pituitary Adenomas Reveals Specific Alterations and Differentiates Among Acromegaly, Cushing's Disease and Endocrine-Inactive Subtypes

Purpose: Pituitary adenomas are one of the most common benign neoplasms of the central nervous system. Although emerging evidence suggests roles for both genetic and epigenetic factors in tumorigenesis, the degree to which these factors contribute to disease remains poorly understood. Experimental Design: A multiplatform analysis was performed to identify the genomic and epigenomic underpinnings of disease among the three major subtypes of surgically resected pituitary adenomas in 48 patients: growth hormone (GH)–secreting (n = 17), adrenocorticotropic hormone (ACTH)–secreting (n = 13, including 3 silent-ACTH adenomas), and endocrine-inactive (n = 18). Whole-exome sequencing was used to profile the somatic mutational landscape, whole-transcriptome sequencing was used to identify disease-specific patterns of gene expression, and array-based DNA methylation profiling was used to examine genome-wide patterns of DNA methylation. Results: Recurrent single-nucleotide and small indel somatic mutations were infrequent among the three adenoma subtypes. However, somatic copy-number alterations (SCNA) were identified in all three pituitary adenoma subtypes. Methylation analysis revealed adenoma subtype-specific DNA methylation profiles, with GH-secreting adenomas being dominated by hypomethylated sites. Likewise, gene-expression patterns revealed adenoma subtype-specific profiles. Integrating DNA methylation and gene-expression data revealed that hypomethylation of promoter regions are related with increased expression of GH1 and SSTR5 genes in GH-secreting adenomas and POMC gene in ACTH-secreting adenomas. Finally, multispectral IHC staining of immune-related proteins showed abundant expression of PD-L1 among all three adenoma subtypes. Conclusions: Taken together, these data stress the contribution of epigenomic alterations to disease-specific etiology among adenoma subtypes and highlight potential targets for future immunotherapy-based treatments. This article reveals novel insights into the epigenomics underlying pituitary adenomas and highlights how differences in epigenomic states are related to important transcriptome alterations that define adenoma subtypes. Clin Cancer Res; 24(17); 4126–36. ©2018 AACR.