Incorporating Spinal Muscular Atrophy Analysis by Next-Generation Sequencing into a Comprehensive Multigene Panel for Neuromuscular Disorders.

Background: Spinal muscular atrophy (SMA) is traditionally molecularly diagnosed by multiplex ligation-dependent probe amplification or quantitative polymerase chain reaction (qPCR). SMA analyses are not routinely incorporated into gene panel analyses for individuals with suspected SMA or broader neuromuscular indications. Aim: We sought to determine whether a next-generation sequencing (NGS) approach that integrates SMA analyses into a multigene neuromuscular disorders panel could detect undiagnosed SMA. Materials and Methods: Sequence and copy number variants of the SMN1/SMN2 genes were simultaneously analyzed in samples from 5304 unselected individuals referred for testing using an NGS-based 122-gene neuromuscular panel. This diagnostic approach was validated using DNA from 68 individuals who had been previously diagnosed with SMA via quantitative PCR for SMN1/SMN2. Results: Homozygous loss of SMN1 was detected in 47 unselected individuals. Heterozygous loss of SMN1 was detected in 118 individuals; 8 had an indeterminate variant in "SMN1 or SMN2" that supported an SMA diagnosis but required additional disambiguation. Of the remaining SMA carriers, 44 had pathogenic variants in other genes. Concordance rates between NGS and qPCR were 100% and 93% for SMN1 and SMN2 copy numbers, respectively. Where there was disagreement, phenotypes were more consistent with the SMN2 results from NGS. Conclusion: Integrating NGS-based SMA testing into a multigene neuromuscular panel allows a single assay to diagnose SMA while comprehensively assessing the spectrum of variants that can occur in individuals with broad differential diagnoses or nonspecific/overlapping neuromuscular features.

Acrophiarin (antibiotic S31794/F-1) from Penicillium arenicola shares biosynthetic features with both Aspergillus- and Leotiomycete-type echinocandins.

The antifungal echinocandin lipopeptide, acrophiarin, was circumscribed in a patent in 1979. We confirmed that the producing strain NRRL 8095 is Penicillium arenicola and other strains of P. arenicola produced acrophiarin and acrophiarin analogues. Genome sequencing of NRRL 8095 identified the acrophiarin gene cluster. Penicillium arenicola and echinocandin-producing Aspergillus species belong to the family Aspergillaceae of the Eurotiomycetes, but several features of acrophiarin and its gene cluster suggest a closer relationship with echinocandins from Leotiomycete fungi. These features include hydroxy-glutamine in the peptide core instead of a serine or threonine residue, the inclusion of a non-heme iron, α-ketoglutarate-dependent oxygenase for hydroxylation of the C3 of the glutamine, and a thioesterase. In addition, P. arenicola bears similarity to Leotiomycete echinocandin-producing species because it exhibits self-resistance to exogenous echinocandins. Phylogenetic analysis of the genes of the echinocandin biosynthetic family indicated that most of the predicted proteins of acrophiarin gene cluster exhibited higher similarity to the predicted proteins of the pneumocandin gene cluster of the Leotiomycete Glarea lozoyensis than to those of the echinocandin B gene cluster from A. pachycristatus. The fellutamide gene cluster and related gene clusters are recognized as relatives of the echinocandins. Inclusion of the acrophiarin gene cluster into a comprehensive phylogenetic analysis of echinocandin gene clusters indicated the divergent evolutionary lineages of echinocandin gene clusters are descendants from a common ancestral progenitor. The minimal 10-gene cluster may have undergone multiple gene acquisitions or losses and possibly horizontal gene transfer after the ancestral separation of the two lineages.

Heterozygote Advantage Is a Common Outcome of Adaptation in Saccharomyces cerevisiae.

Adaptation in diploids is predicted to proceed via mutations that are at least partially dominant in fitness. Recently, we argued that many adaptive mutations might also be commonly overdominant in fitness. Natural (directional) selection acting on overdominant mutations should drive them into the population but then, instead of bringing them to fixation, should maintain them as balanced polymorphisms via heterozygote advantage. If true, this would make adaptive evolution in sexual diploids differ drastically from that of haploids. The validity of this prediction has not yet been tested experimentally. Here, we performed four replicate evolutionary experiments with diploid yeast populations (Saccharomyces cerevisiae) growing in glucose-limited continuous cultures. We sequenced 24 evolved clones and identified initial adaptive mutations in all four chemostats. The first adaptive mutations in all four chemostats were three copy number variations, all of which proved to be overdominant in fitness. The fact that fitness overdominant mutations were always the first step in independent adaptive walks supports the prediction that heterozygote advantage can arise as a common outcome of directional selection in diploids and demonstrates that overdominance of de novo adaptive mutations in diploids is not rare.

The Valley-of-Death: reciprocal sign epistasis constrains adaptive trajectories in a constant, nutrient limiting environment.

The fitness landscape is a powerful metaphor for describing the relationship between genotype and phenotype for a population under selection. However, empirical data as to the topography of fitness landscapes are limited, owing to difficulties in measuring fitness for large numbers of genotypes under any condition. We previously reported a case of reciprocal sign epistasis (RSE), where two mutations individually increased yeast fitness in a glucose-limited environment, but reduced fitness when combined, suggesting the existence of two peaks on the fitness landscape. We sought to determine whether a ridge connected these peaks so that populations founded by one mutant could reach the peak created by the other, avoiding the low-fitness "Valley-of-Death" between them. Sequencing clones after 250 generations of further evolution provided no evidence for such a ridge, but did reveal many presumptive beneficial mutations, adding to a growing body of evidence that clonal interference pervades evolving microbial populations.

Phenotypic and genotypic convergences are influenced by historical contingency and environment in yeast.

Different organisms have independently and recurrently evolved similar phenotypic traits at different points throughout history. This phenotypic convergence may be caused by genotypic convergence and in addition, constrained by historical contingency. To investigate how convergence may be driven by selection in a particular environment and constrained by history, we analyzed nine life-history traits and four metabolic traits during an experimental evolution of six yeast strains in four different environments. In each of the environments, the population converged toward a different multivariate phenotype. However, the evolution of most traits, including fitness components, was constrained by history. Phenotypic convergence was partly associated with the selection of mutations in genes involved in the same pathway. By further investigating the convergence in one gene, BMH1, mutated in 20% of the evolved populations, we show that both the history and the environment influenced the types of mutations (missense/nonsense), their location within the gene itself, as well as their effects on multiple traits. However, these effects could not be easily predicted from ancestors' phylogeny or past selection. Combined, our data highlight the role of pleiotropy and epistasis in shaping a rugged fitness landscape.

Whole genome, whole population sequencing reveals that loss of signaling networks is the major adaptive strategy in a constant environment.

Molecular signaling networks are ubiquitous across life and likely evolved to allow organisms to sense and respond to environmental change in dynamic environments. Few examples exist regarding the dispensability of signaling networks, and it remains unclear whether they are an essential feature of a highly adapted biological system. Here, we show that signaling network function carries a fitness cost in yeast evolving in a constant environment. We performed whole-genome, whole-population Illumina sequencing on replicate evolution experiments and find the major theme of adaptive evolution in a constant environment is the disruption of signaling networks responsible for regulating the response to environmental perturbations. Over half of all identified mutations occurred in three major signaling networks that regulate growth control: glucose signaling, Ras/cAMP/PKA and HOG. This results in a loss of environmental sensitivity that is reproducible across experiments. However, adaptive clones show reduced viability under starvation conditions, demonstrating an evolutionary tradeoff. These mutations are beneficial in an environment with a constant and predictable nutrient supply, likely because they result in constitutive growth, but reduce fitness in an environment where nutrient supply is not constant. Our results are a clear example of the myopic nature of evolution: a loss of environmental sensitivity in a constant environment is adaptive in the short term, but maladaptive should the environment change.

Analysis of the Saccharomyces cerevisiae pan-genome reveals a pool of copy number variants distributed in diverse yeast strains from differing industrial environments.

Although the budding yeast Saccharomyces cerevisiae is arguably one of the most well-studied organisms on earth, the genome-wide variation within this species--i.e., its "pan-genome"--has been less explored. We created a multispecies microarray platform containing probes covering the genomes of several Saccharomyces species: S. cerevisiae, including regions not found in the standard laboratory S288c strain, as well as the mitochondrial and 2-µm circle genomes-plus S. paradoxus, S. mikatae, S. kudriavzevii, S. uvarum, S. kluyveri, and S. castellii. We performed array-Comparative Genomic Hybridization (aCGH) on 83 different S. cerevisiae strains collected across a wide range of habitats; of these, 69 were commercial wine strains, while the remaining 14 were from a diverse set of other industrial and natural environments. We observed interspecific hybridization events, introgression events, and pervasive copy number variation (CNV) in all but a few of the strains. These CNVs were distributed throughout the strains such that they did not produce any clear phylogeny, suggesting extensive mating in both industrial and wild strains. To validate our results and to determine whether apparently similar introgressions and CNVs were identical by descent or recurrent, we also performed whole-genome sequencing on nine of these strains. These data may help pinpoint genomic regions involved in adaptation to different industrial milieus, as well as shed light on the course of domestication of S. cerevisiae.

Reciprocal sign epistasis between frequently experimentally evolved adaptive mutations causes a rugged fitness landscape.

The fitness landscape captures the relationship between genotype and evolutionary fitness and is a pervasive metaphor used to describe the possible evolutionary trajectories of adaptation. However, little is known about the actual shape of fitness landscapes, including whether valleys of low fitness create local fitness optima, acting as barriers to adaptive change. Here we provide evidence of a rugged molecular fitness landscape arising during an evolution experiment in an asexual population of Saccharomyces cerevisiae. We identify the mutations that arose during the evolution using whole-genome sequencing and use competitive fitness assays to describe the mutations individually responsible for adaptation. In addition, we find that a fitness valley between two adaptive mutations in the genes MTH1 and HXT6/HXT7 is caused by reciprocal sign epistasis, where the fitness cost of the double mutant prohibits the two mutations from being selected in the same genetic background. The constraint enforced by reciprocal sign epistasis causes the mutations to remain mutually exclusive during the experiment, even though adaptive mutations in these two genes occur several times in independent lineages during the experiment. Our results show that epistasis plays a key role during adaptation and that inter-genic interactions can act as barriers between adaptive solutions. These results also provide a new interpretation on the classic Dobzhansky-Muller model of reproductive isolation and display some surprising parallels with mutations in genes often associated with tumors.

Transient genotype-by-environment interactions following environmental shock provide a source of expression variation for essential genes.

Understanding complex genotype-by-environment interactions (GEIs) is crucial for understanding phenotypic variation. An important factor often overlooked in GEI studies is time. We measured the contribution of GEIs to expression variation in four nonlaboratory Saccharomyces cerevisiae strains responding dynamically to a 25 degrees -37 degrees heat shock. GEI was a major force explaining expression variation, affecting 55% of the genes analyzed. Importantly, almost half of these expression patterns showed GEI influence only during the transition between environments, but not in acclimated cells. This class reveals a genotype-by-environment-by-time interaction that affected expression of a large fraction of yeast genes. Strikingly, although transcripts subject to persistent GEI effects were enriched for nonessential genes with upstream TATA elements, those displaying transient GEIs were enriched for essential genes regardless of TATA regulation. Genes subject to persistent GEI influences showed relaxed constraint on acclimated gene expression compared to the average yeast gene, whereas genes restricted to transient GEIs did not. We propose that transient GEI during the transition between environments provides a previously unappreciated source of expression variation, particularly for essential genes.

Variations in stress sensitivity and genomic expression in diverse S. cerevisiae isolates.

Interactions between an organism and its environment can significantly influence phenotypic evolution. A first step toward understanding this process is to characterize phenotypic diversity within and between populations. We explored the phenotypic variation in stress sensitivity and genomic expression in a large panel of Saccharomyces strains collected from diverse environments. We measured the sensitivity of 52 strains to 14 environmental conditions, compared genomic expression in 18 strains, and identified gene copy-number variations in six of these isolates. Our results demonstrate a large degree of phenotypic variation in stress sensitivity and gene expression. Analysis of these datasets reveals relationships between strains from similar niches, suggests common and unique features of yeast habitats, and implicates genes whose variable expression is linked to stress resistance. Using a simple metric to suggest cases of selection, we found that strains collected from oak exudates are phenotypically more similar than expected based on their genetic diversity, while sake and vineyard isolates display more diverse phenotypes than expected under a neutral model. We also show that the laboratory strain S288c is phenotypically distinct from all of the other strains studied here, in terms of stress sensitivity, gene expression, Ty copy number, mitochondrial content, and gene-dosage control. These results highlight the value of understanding the genetic basis of phenotypic variation and raise caution about using laboratory strains for comparative genomics.