CooVar: co-occurring variant analyzer

The α-arrestin SUP-13/ARRD-15 promotes isoform turnover of actin-interacting protein 1 in Caenorhabditis elegans striated muscle

Precise arrangement of actin, myosin, and other regulatory components in a sarcomeric pattern is critical for producing contractile forces in striated muscles. Actin-interacting protein 1 (AIP1), also known as WD-repeat protein 1 (WDR1), is one of essential factors that regulate sarcomeric assembly of actin filaments. In the nematode Caenorhabditis elegans, mutation in unc-78, encoding one of the two AIP1 isoforms, causes severe disorganization of sarcomeric actin filaments and near paralysis, but mutation in sup-13 suppresses the unc-78-mutant phenotypes to restore nearly normal sarcomeric actin organization and worm motility. Here, we identified that sup-13 is a nonsense allele of arrd-15 encoding an α-arrestin. The sup-13/arrd-15 mutation suppressed the phenotypes of unc-78 null mutant but required aipl-1 that encodes a second AIP1 isoform. aipl-1 was normally expressed highly in embryos and downregulated in mature muscle. However, in the sup-13/arrd-15 mutant, the AIPL-1 protein was maintained at high levels in adult muscle to compensate for the absence of the UNC-78 protein. The sup-13/arrd-15 mutation caused accumulation of ubiquitinated AIPL-1 protein, suggesting that a normal function of sup-13/arrd-15 is to enhance degradation of ubiquitinated AIPL-1, thereby promoting transition of AIP1 isoforms from AIPL-1 to UNC-78 in developing muscle. These results suggest that α-arrestin is a novel factor to promote isoform turnover by enhancing protein degradation.

Genetic factors inherited from both diploid parents interact to affect genome stability and fertility in resynthesized allotetraploid Brassica napus

Established allopolyploids are known to be genomically stable and fertile. However, in contrast, most newly resynthesized allopolyploids are infertile and meiotically unstable. Identifying the genetic factors responsible for genome stability in newly formed allopolyploid is key to understanding how 2 genomes come together to form a species. One hypothesis is that established allopolyploids may have inherited specific alleles from their diploid progenitors which conferred meiotic stability. Resynthesized Brassica napus lines are often unstable and infertile, unlike B. napus cultivars. We tested this hypothesis by characterizing 41 resynthesized B. napus lines produced by crosses between 8 Brassica rapa and 8 Brassica oleracea lines for copy number variation resulting from nonhomologous recombination events and fertility. We resequenced 8 B. rapa and 5 B. oleracea parent accessions and analyzed 19 resynthesized lines for allelic variation in a list of meiosis gene homologs. SNP genotyping was performed using the Illumina Infinium Brassica 60K array for 3 individuals per line. Self-pollinated seed set and genome stability (number of copy number variants) were significantly affected by the interaction between both B. rapa and B. oleracea parental genotypes. We identified 13 putative meiosis gene candidates which were significantly associated with frequency of copy number variants and which contained putatively harmful mutations in meiosis gene haplotypes for further investigation. Our results support the hypothesis that allelic variants inherited from parental genotypes affect genome stability and fertility in resynthesized rapeseed.

Botrytis cinerea Loss and Restoration of Virulence during In Vitro Culture Follows Flux in Global DNA Methylation

Pathogenic fungi can lose virulence after protracted periods of culture, but little is known of the underlying mechanisms. Here, we present the first analysis of DNA methylation flux at a single-base resolution for the plant pathogen B. cinerea and identify differentially methylated genes/genomic regions associated with virulence erosion during in vitro culture. Cultures were maintained for eight months, with subcultures and virulence testing every month. Methylation-sensitive amplified polymorphisms were performed at monthly intervals to characterise global changes to the pathogen’s genome during culture and also on DNA from mycelium inoculated onto Arabidopsis thaliana after eight months in culture. Characterisation of culture-induced epialleles was assessed by whole-genome re-sequencing and whole-genome bisulfite sequencing. Virulence declined with time in culture and recovered after inoculation on A. thaliana. Variation detected by methylation-sensitive amplified polymorphisms followed virulence changes during culture. Whole-genome (bisulfite) sequencing showed marked changes in global and local methylation during culture but no significant genetic changes. We imply that virulence is a non-essential plastic character that is at least partly modified by the changing levels of DNA methylation during culture. We hypothesise that changing DNA methylation during culture may be responsible for the high virulence/low virulence transition in B. cinerea and speculate that this may offer fresh opportunities to control pathogen virulence.

Whole genome sequencing facilitates intragenic variant interpretation following modifier screening in C. elegans

Background

Intragenic modifiers (in-phase, second-site variants) are known to have dramatic effects on clinical outcomes, affecting disease attributes such as severity or age of onset. However, despite their clinical importance, the focus of many genetic screens in model systems is on the discovery of extragenic variants, with many labs still relying upon more traditional methods to identify modifiers. However, traditional methods such as PCR and Sanger sequencing can be time-intensive and do not permit a thorough understanding of the intragenic modifier effects in the context of non-isogenic genomic backgrounds.

Results

Here, we apply high throughput approaches to identify and understand intragenic modifiers using Caenorhabditis elegans. Specifically, we applied whole genome sequencing (WGS) to a mutagen-induced forward genetic screen to identify intragenic suppressors of a temperature-sensitive zyg-1(it25) allele in C. elegans. ZYG-1 is a polo kinase that is important for centriole function and cell divisions, and mutations that truncate its human orthologue, PLK4, have been associated with microcephaly. Combining WGS and CRISPR/Cas9, we rapidly identify intragenic modifiers, show that these variants are distributed non-randomly throughout zyg-1 and that genomic context plays an important role on phenotypic outcomes.

Conclusions

Ultimately, our work shows that WGS facilitates high-throughput identification of intragenic modifiers in clinically relevant genes by reducing hands-on research time and overall costs and by allowing thorough understanding of the intragenic phenotypic effects in the context of different genetic backgrounds.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08142-8.

An inferred functional impact map of genetic variants in rice

Interpreting the functional impacts of genetic variants (GVs) is an important challenge for functional genomic studies in crops and next-generation breeding. Previous studies in rice (Oryza sativa) have focused mainly on the identification of GVs, whereas systematic functional annotation of GVs has not yet been performed. Here, we present a functional impact map of GVs in rice. We curated haplotype information for 17 397 026 GVs from sequencing data of 4726 rice accessions. We quantitatively evaluated the effects of missense mutations in coding regions in each haplotype based on the conservation of amino acid residues and obtained the effects of 918 848 non-redundant missense GVs. Furthermore, we generated high-quality chromatin accessibility (CA) data from six representative rice tissues and used these data to train deep convolutional neural network models to predict the impacts of 5 067 405 GVs for CA in regulatory regions. We characterized the functional properties and tissue specificity of the GV effects and found that large-effect GVs in coding and regulatory regions may be subject to selection in different directions. Finally, we demonstrated how the functional impact map could be used to prioritize causal variants in mapping populations. This impact map will be a useful resource for accelerating gene cloning and functional studies in rice, and can be freely queried in RiceVarMap V2.0 (http://ricevarmap.ncpgr.cn).

Functional effects of protein variants

Genetic and other variations frequently affect protein functions. Scientific articles can contain confusing descriptions about which function or property is affected, and in many cases the statements are pure speculation without any experimental evidence. To clarify functional effects of protein variations of genetic or non-genetic origin, a systematic conceptualisation and framework are introduced. This framework describes protein functional effects on abundance, activity, specificity and affinity, along with countermeasures, which allow cells, tissues and organisms to tolerate, avoid, repair, attenuate or resist (TARAR) the effects. Effects on abundance discussed include gene dosage, restricted expression, mis-localisation and degradation. Enzymopathies, effects on kinetics, allostery and regulation of protein activity are subtopics for the effects of variants on activity. Variation outcomes on specificity and affinity comprise promiscuity, specificity, affinity and moonlighting. TARAR mechanisms redress variations with active and passive processes including chaperones, redundancy, robustness, canalisation and metabolic and signalling rewiring. A framework for pragmatic protein function analysis and presentation is introduced. All of the mechanisms and effects are described along with representative examples, most often in relation to diseases. In addition, protein function is discussed from evolutionary point of view. Application of the presented framework facilitates unambiguous, detailed and specific description of functional effects and their systematic study.

Adaptive Evolution of Geobacter sulfurreducens in Coculture with Pseudomonas aeruginosa

Geobacter and Pseudomonas spp. cohabit many of the same environments, where Geobacter spp. often dominate. Both bacteria are capable of extracellular electron transfer (EET) and play important roles in biogeochemical cycling. Although they recently in 2017 were demonstrated to undergo direct interspecies electron transfer (DIET) with one another, the genetic evolution of this syntrophic interaction has not been examined. Here, we use whole-genome sequencing of the cocultures before and after adaptive evolution to determine whether genetic selection is occurring. We also probe their interaction on a temporal level and determine whether their interaction dynamics change over the course of adaptive evolution. This study brings to light the multifaceted nature of interactions between just two microorganisms within a controlled environment and will aid in improving metabolic models of microbial communities comprising these two bacteria.

Interactions between microorganisms in mixed communities are highly complex, being either syntrophic, neutral, predatory, or competitive. Evolutionary changes can occur in the interaction dynamics between community members as they adapt to coexistence. Here, we report that the syntrophic interaction between Geobacter sulfurreducens and Pseudomonas aeruginosa coculture change in their dynamics over evolutionary time. Specifically, Geobacter sp. dominance increases with adaptation within the cocultures, as determined through quantitative PCR and fluorescence in situ hybridization. This suggests a transition from syntrophy to competition and demonstrates the rapid adaptive capacity of Geobacter spp. to dominate in cocultures with P. aeruginosa. Early in coculture establishment, two single-nucleotide variants in the G. sulfurreducensfabI and tetR genes emerged that were strongly selected for throughout coculture evolution with P. aeruginosa phenazine wild-type and phenazine-deficient mutants. Sequential window acquisition of all theoretical spectra-mass spectrometry (SWATH-MS) proteomics revealed that the tetR variant cooccurred with the upregulation of an adenylate cyclase transporter, CyaE, and a resistance-nodulation-division (RND) efflux pump notably known for antibiotic efflux. To determine whether antibiotic production was driving the increased expression of the multidrug efflux pump, we tested Pseudomonas-derived phenazine-1-carboxylic acid (PHZ-1-CA) for its potential to inhibit Geobacter growth and drive selection of the tetR and fabI genetic variants. Despite its inhibitory properties, PHZ-1-CA did not drive variant selection, indicating that other antibiotics may drive overexpression of the efflux pump and CyaE or that a novel role exists for these proteins in the context of this interaction.

A-to-I RNA Editing Uncovers Hidden Signals of Adaptive Genome Evolution in Animals

In animals, the most common type of RNA editing is the deamination of adenosines (A) into inosines (I). Because inosines basepair with cytosines (C), they are interpreted as guanosines (G) by the cellular machinery and genomically encoded G alleles at edited sites mimic the function of edited RNAs. The contribution of this hardwiring effect on genome evolution remains obscure. We looked for population genomics signatures of adaptive evolution associated with A-to-I RNA edited sites in humans and Drosophila melanogaster. We found that single nucleotide polymorphisms at edited sites occur 3 (humans) to 15 times (Drosophila) more often than at unedited sites, the nucleotide G is virtually the unique alternative allele at edited sites and G alleles segregate at higher frequency at edited sites than at unedited sites. Our study reveals that a significant fraction of coding synonymous and nonsynonymous as well as silent and intergenic A-to-I RNA editing sites are likely adaptive in the distantly related human and Drosophila lineages.

Parallel and nonparallel genomic responses contribute to herbicide resistance in Ipomoea purpurea, a common agricultural weed

The repeated evolution of herbicide resistance has been cited as an example of genetic parallelism, wherein separate species or genetic lineages utilize the same genetic solution in response to selection. However, most studies that investigate the genetic basis of herbicide resistance examine the potential for changes in the protein targeted by the herbicide rather than considering genome-wide changes. We used a population genomics screen and targeted exome re-sequencing to uncover the potential genetic basis of glyphosate resistance in the common morning glory, Ipomoea purpurea, and to determine if genetic parallelism underlies the repeated evolution of resistance across replicate resistant populations. We found no evidence for changes in 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS), glyphosate's target protein, that were associated with resistance, and instead identified five genomic regions that showed evidence of selection. Within these regions, genes involved in herbicide detoxification-cytochrome P450s, ABC transporters, and glycosyltransferases-are enriched and exhibit signs of selective sweeps. One region under selection shows parallel changes across all assayed resistant populations whereas other regions exhibit signs of divergence. Thus, while it appears that the physiological mechanism of resistance in this species is likely the same among resistant populations, we find patterns of both similar and divergent selection across separate resistant populations at particular loci.

Enhanced Growth of Pilin-Deficient Geobacter sulfurreducens Mutants in Carbon Poor and Electron Donor Limiting Conditions

Geobacter sulfurreducens pili enable extracellular electron transfer and play a role in secretion of c-type cytochromes such as OmcZ. PilA-deficient mutants of G. sulfurreducens have previously been shown to accumulate cytochromes within their membranes. This cytochrome retaining phenotype allowed for enhanced growth of PilA-deficient mutants in electron donor and carbon-limited conditions where formate and fumarate are provided as the sole electron donor and acceptor with no supplementary carbon source. Conversely, wild-type G. sulfurreducens, which has normal secretion of cytochromes, has comparative limited growth in these conditions. This growth is further impeded for OmcZ-deficient and OmcS-deficient mutants. A PilB-deficient mutant which prevents pilin production but allows for secretion of OmcZ had moderate growth in these conditions, indicating a role for cytochrome localization to enabling survival in the electron donor and carbon-limited conditions. To determine which pathways enhanced growth using formate, Sequential Window Acquisition of all Theoretical Mass Spectra mass spectrometry (SWATH-MS) proteomics of formate adapted PilA-deficient mutants and acetate grown wild type was performed. PilA-deficient mutants had an overall decrease in tricarboxylic acid (TCA) cycle enzymes and significant upregulation of electron transport chain associated proteins including many c-type cytochromes and [NiFe]-hydrogenases. Whole genome sequencing of the mutants shows strong convergent evolution and emergence of genetic subpopulations during adaptation to growth on formate. The results described here suggest a role for membrane constrained c-type cytochromes to the enhancement of survival and growth in electron donor and carbon-limited conditions.

Kettin, the large actin-binding protein with multiple immunoglobulin domains, is essential for sarcomeric actin assembly and larval development in Caenorhabditis elegans

Among many essential genes in the nematode Caenorhabditis elegans, let-330 is located on the left arm of chromosome V and was identified as the largest target of a mutagen in this region. However, let-330 gene has not been characterized at the molecular level. Here, we report that two sequenced let-330 alleles are nonsense mutations of ketn-1, a previously characterized gene encoding kettin. Kettin is a large actin-binding protein of 472 kDa with 31 immunoglobulin domains and is expressed in muscle cells in C. elegans. let-330/ketn-1 mutants are homozygous lethal at the first larval stage with mild defects in body elongation. These mutants have severe defects in sarcomeric actin and myosin assembly in striated muscle. However, α-actinin and vinculin, which are components of the dense bodies anchoring actin to the membranes, were not significantly disorganized by let-330/ketn-1 mutation. Kettin localizes to embryonic myofibrils before α-actinin is expressed, and α-actinin deficiency does not affect kettin localization in larval muscle. Depletion of vinculin minimally affects kettin localization but significantly reduces colocalization of actin with kettin in embryonic muscle cells. These results indicate that kettin is an essential protein for sarcomeric assembly of actin filaments in muscle cells.

Inherited allelic variants and novel karyotype changes influence fertility and genome stability in Brassica allohexaploids

Synthetic allohexaploid Brassica hybrids (2n = AABBCC) do not exist naturally, but can be synthesized by crosses between diploid and/or allotetraploid Brassica species. Using these hybrids, we aimed to identify how novel allohexaploids restore fertility and normal meiosis after formation. Chromosome inheritance, genome structure, fertility and meiotic behaviour were assessed in three segregating allohexaploid populations derived from the cross (B. napus × B. carinata) × B. juncea using a combination of molecular marker genotyping, phenotyping and cytogenetics. Plants with unbalanced A-C translocations in one direction (where a C-genome chromosome fragment replaces an A-genome fragment) but not the other (where an A-genome fragment replaces a C-genome fragment) showed significantly reduced fertility across all populations. Genomic regions associated with fertility contained several meiosis genes with putatively causal mutations inherited from the parents (copies of SCC2 in the A genome, PAIR1/PRD3, PRD1 and ATK1/KATA/KIN14a in the B genome, and MSH2 and SMC1/TITAN8 in the C genome). Reduced seed fertility associated with the loss of chromosome fragments from only one subgenome following homoeologous exchanges could comprise a mechanism for biased genome fractionation in allopolyploids. Pre-existing meiosis gene variants present in allotetraploid parents may help to stabilize meiosis in novel allohexaploids.

Systematic identification and annotation of multiple-variant compound effects at transcription factor binding sites in human genome

Understanding the functional effects of genetic variants is crucial in modern genomics and genetics. Transcription factor binding sites (TFBSs) are one of the most important cis-regulatory elements. While multiple tools have been developed to assess functional effects of genetic variants at TFBSs, they usually assume that each variant works in isolation and neglect the potential "interference" among multiple variants within the same TFBS. In this study, we presented COPE-TFBS (Context-Oriented Predictor for variant Effect on Transcription Factor Binding Site), a novel method that considers sequence context to accurately predict variant effects on TFBSs. We systematically re-analyzed the sequencing data from both the 1000 Genomes Project and the Genotype-Tissue Expression (GTEx) Project via COPE-TFBS, and identified numbers of novel TFBSs, transformed TFBSs and discordantly annotated TFBSs resulting from multiple variants, further highlighting the necessity of sequence context in accurately annotating genetic variants. COPE-TFBS is freely available for academic use at http://cope.cbi.pku.edu.cn/.

Genomic Identification and Functional Characterization of Essential Genes in Caenorhabditis elegans

Using combined genetic mapping, Illumina sequencing, bioinformatics analyses, and experimental validation, we identified 60 essential genes from 104 lethal mutations in two genomic regions of Caenorhabditis elegans totaling ∼14 Mb on chromosome III(mid) and chromosome V(left). Five of the 60 genes had not previously been shown to have lethal phenotypes by RNA interference depletion. By analyzing the regions around the lethal missense mutations, we identified four putative new protein functional domains. Furthermore, functional characterization of the identified essential genes shows that most are enzymes, including helicases, tRNA synthetases, and kinases in addition to ribosomal proteins. Gene Ontology analysis indicated that essential genes often encode for enzymes that conduct nucleic acid binding activities during fundamental processes, such as intracellular DNA replication, transcription, and translation. Analysis of essential gene shows that they have fewer paralogs, encode proteins that are in protein interaction hubs, and are highly expressed relative to nonessential genes. All these essential gene traits in C. elegans are consistent with those of human disease genes. Most human orthologs (90%) of the essential genes in this study are related to human diseases. Therefore, functional characterization of essential genes underlines their importance as proxies for understanding the biological functions of human disease genes.

Accurately annotate compound effects of genetic variants using a context-sensitive framework

In genomics, effectively identifying the biological effects of genetic variants is crucial. Current methods handle each variant independently, assuming that each variant acts in a context-free manner. However, variants within the same gene may interfere with each other, producing combinational (compound) rather than individual effects. In this work, we introduce COPE, a gene-centric variant annotation tool that integrates the entire sequential context in evaluating the functional effects of intra-genic variants. Applying COPE to the 1000 Genomes dataset, we identified numerous cases of multiple-variant compound effects that frequently led to false-positive and false-negative loss-of-function calls by conventional variant-centric tools. Specifically, 64 disease-causing mutations were identified to be rescued in a specific genomic context, thus potentially contributing to the buffering effects for highly penetrant deleterious mutations. COPE is freely available for academic use at http://cope.cbi.pku.edu.cn.

Whole-Genome Sequence and Variant Analysis of W303, a Widely-Used Strain of Saccharomyces cerevisiae

The yeast Saccharomyces cerevisiae has emerged as a superior model organism. Selection of distinct laboratory strains of S. cerevisiae with unique phenotypic properties, such as superior mating or sporulation efficiencies, has facilitated advancements in research. W303 is one such laboratory strain that is closely related to the first completely sequenced yeast strain, S288C. In this work, we provide a high-quality, annotated genome sequence for W303 for utilization in comparative analyses and genome-wide studies. Approximately 9500 variations exist between S288C and W303, affecting the protein sequences of ∼700 genes. A listing of the polymorphisms and divergent genes is provided for researchers interested in identifying the genetic basis for phenotypic differences between W303 and S288C. Several divergent functional gene families were identified, including flocculation and sporulation genes, likely representing selection for desirable laboratory phenotypes. Interestingly, remnants of ancestor wine strains were found on several chromosomes. Finally, as a test of the utility of the high-quality reference genome, variant mapping revealed more accurate identification of accumulated mutations in passaged mismatch repair-defective strains.

Targeted deep sequencing of flowering regulators in Brassica napus reveals extensive copy number variation

Gene copy number variation (CNV) is increasingly implicated in control of complex trait networks, particularly in polyploid plants like rapeseed (Brassica napus L.) with an evolutionary history of genome restructuring. Here we performed sequence capture to assay nucleotide variation and CNV in a panel of central flowering time regulatory genes across a species-wide diversity set of 280 B. napus accessions. The genes were chosen based on prior knowledge from Arabidopsis thaliana and related Brassica species. Target enrichment was performed using the Agilent SureSelect technology, followed by Illumina sequencing. A bait (probe) pool was developed based on results of a preliminary experiment with representatives from different B. napus morphotypes. A very high mean target coverage of ~670x allowed reliable calling of CNV, single nucleotide polymorphisms (SNPs) and insertion-deletion (InDel) polymorphisms. No accession exhibited no CNV, and at least one homolog of every gene we investigated showed CNV in some accessions. Some CNV appear more often in specific morphotypes, indicating a role in diversification.

Towards a whole-genome sequence for rye (Secale cereale L.)

We report on a whole-genome draft sequence of rye (Secale cereale L.). Rye is a diploid Triticeae species closely related to wheat and barley, and an important crop for food and feed in Central and Eastern Europe. Through whole-genome shotgun sequencing of the 7.9-Gbp genome of the winter rye inbred line Lo7 we obtained a de novo assembly represented by 1.29 million scaffolds covering a total length of 2.8 Gbp. Our reference sequence represents nearly the entire low-copy portion of the rye genome. This genome assembly was used to predict 27 784 rye gene models based on homology to sequenced grass genomes. Through resequencing of 10 rye inbred lines and one accession of the wild relative S. vavilovii, we discovered more than 90 million single nucleotide variants and short insertions/deletions in the rye genome. From these variants, we developed the high-density Rye600k genotyping array with 600 843 markers, which enabled anchoring the sequence contigs along a high-density genetic map and establishing a synteny-based virtual gene order. Genotyping data were used to characterize the diversity of rye breeding pools and genetic resources, and to obtain a genome-wide map of selection signals differentiating the divergent gene pools. This rye whole-genome sequence closes a gap in Triticeae genome research, and will be highly valuable for comparative genomics, functional studies and genome-based breeding in rye.

Post-polyploidisation morphotype diversification associates with gene copy number variation

Genetic models for polyploid crop adaptation provide important information relevant for future breeding prospects. A well-suited model is Brassica napus, a recent allopolyploid closely related to Arabidopsis thaliana. Flowering time is a major adaptation trait determining life cycle synchronization with the environment. Here we unravel natural genetic variation in B. napus flowering time regulators and investigate associations with evolutionary diversification into different life cycle morphotypes. Deep sequencing of 35 flowering regulators was performed in 280 diverse B. napus genotypes. High sequencing depth enabled high-quality calling of single-nucleotide polymorphisms (SNPs), insertion-deletions (InDels) and copy number variants (CNVs). By combining these data with genotyping data from the Brassica 60 K Illumina® Infinium SNP array, we performed a genome-wide marker distribution analysis across the 4 ecogeographical morphotypes. Twelve haplotypes, including Bna.FLC.A10, Bna.VIN3.A02 and the Bna.FT promoter on C02_random, were diagnostic for the diversification of winter and spring types. The subspecies split between oilseed/kale (B. napus ssp. napus) and swedes/rutabagas (B. napus ssp. napobrassica) was defined by 13 haplotypes, including genomic rearrangements encompassing copies of Bna.FLC, Bna.PHYA and Bna.GA3ox1. De novo variation in copies of important flowering-time genes in B. napus arose during allopolyploidisation, enabling sub-functionalisation that allowed different morphotypes to appropriately fine-tune their lifecycle.

Flowering Time Gene Variation in Brassica Species Shows Evolutionary Principles

Flowering time genes have a strong influence on successful reproduction and life cycle adaptation. However, their regulation is highly complex and only well understood in diploid model systems. For crops with a polyploid background from the genus Brassica, data on flowering time gene variation are scarce, although indispensable for modern breeding techniques like marker-assisted breeding. We have deep-sequenced all paralogs of 35 Arabidopsis thaliana flowering regulators using Sequence Capture followed by Illumina sequencing in two selected accessions of the vegetable species Brassica rapa and Brassica oleracea, respectively. Using these data, we were able to call SNPs, InDels and copy number variations (CNVs) for genes from the total flowering time network including central flowering regulators, but also genes from the vernalisation pathway, the photoperiod pathway, temperature regulation, the circadian clock and the downstream effectors. Comparing the results to a complementary data set from the allotetraploid species Brassica napus, we detected rearrangements in B. napus which probably occurred early after the allopolyploidisation event. Those data are both a valuable resource for flowering time research in those vegetable species, as well as a contribution to speciation genetics.

Degenerin channel activation causes caspase-mediated protein degradation and mitochondrial dysfunction in adult C. elegans muscle

BACKGROUND

Declines in skeletal muscle structure and function are found in various clinical populations, but the intramuscular proteolytic pathways that govern declines in these individuals remain relatively poorly understood. The nematode Caenorhabditis elegans has been developed into a model for identifying and understanding these pathways. Recently, it was reported that UNC-105/degenerin channel activation produced muscle protein degradation via an unknown mechanism.

METHODS

Generation of transgenic and double mutant C. elegans, RNAi, and drug treatments were utilized to assess molecular events governing protein degradation. Western blots were used to measure protein content. Cationic dyes and adenosine triphosphate (ATP) production assays were utilized to measure mitochondrial function.

RESULTS

unc-105 gain-of-function mutants display aberrant muscle protein degradation and a movement defect; both are reduced in intragenic revertants and in let-2 mutants that gate the hyperactive UNC-105 channel. Degradation is not suppressed by interventions suppressing proteasome-mediated, autophagy-mediated, or calpain-mediated degradation nor by suppressors of degenerin-induced neurodegeneration. Protein degradation, but not the movement defect, is decreased by treatment with caspase inhibitors or RNAi against ced-3 or ced-4. Adult unc-105 muscles display a time-dependent fragmentation of the mitochondrial reticulum that is associated with impaired mitochondrial membrane potential and that correlates with decreased rates of maximal ATP production. Reduced levels of CED-4, which is sufficient to activate CED-3 in vitro, are observed in unc-105 mitochondrial isolations.

CONCLUSIONS

Constitutive cationic influx into muscle appears to cause caspase degradation of cytosolic proteins as the result of mitochondrial dysfunction, which may be relevant to ageing and sarcopenia.

Comprehensive genome and epigenome characterization of CHO cells in response to evolutionary pressures and over time

The most striking characteristic of CHO cells is their adaptability, which enables efficient production of proteins as well as growth under a variety of culture conditions, but also results in genomic and phenotypic instability. To investigate the relative contribution of genomic and epigenetic modifications towards phenotype evolution, comprehensive genome and epigenome data are presented for six related CHO cell lines, both in response to perturbations (different culture conditions and media as well as selection of a specific phenotype with increased transient productivity) and in steady state (prolonged time in culture under constant conditions). Clear transitions were observed in DNA‐methylation patterns upon each perturbation, while few changes occurred over time under constant conditions. Only minor DNA‐methylation changes were observed between exponential and stationary growth phase; however, throughout a batch culture the histone modification pattern underwent continuous adaptation. Variation in genome sequence between the six cell lines on the level of SNPs, InDels, and structural variants is high, both upon perturbation and under constant conditions over time. The here presented comprehensive resource may open the door to improved control and manipulation of gene expression during industrial bioprocesses based on epigenetic mechanisms. Biotechnol. Bioeng. 2016;113: 2241–2253. © 2016 The Authors. Biotechnology and Bioengineering Published by Wiley Periodicals, Inc.

Species-wide genome sequence and nucleotide polymorphisms from the model allopolyploid plant Brassica napus

Brassica napus (oilseed rape, canola) is one of the world's most important sources of vegetable oil for human nutrition and biofuel, and also a model species for studies investigating the evolutionary consequences of polyploidisation. Strong bottlenecks during its recent origin from interspecific hybridisation, and subsequently through intensive artificial selection, have severely depleted the genetic diversity available for breeding. On the other hand, high-throughput genome profiling technologies today provide unprecedented scope to identify, characterise and utilise genetic diversity in primary and secondary crop gene pools. Such methods also enable implementation of genomic selection strategies to accelerate breeding progress. The key prerequisite is availability of high-quality sequence data and identification of high-quality, genome-wide sequence polymorphisms representing relevant gene pools. We present comprehensive genome resequencing data from a panel of 52 highly diverse natural and synthetic B. napus accessions, along with a stringently selected panel of 4.3 million high-confidence, genome-wide SNPs. The data is of great interest for genomics-assisted breeding and for evolutionary studies on the origins and consequences in allopolyploidisation in plants.

Bulk Segregant Analysis Reveals the Genetic Basis of a Natural Trait Variation in Fission Yeast

Although the fission yeast Schizosaccharomyces pombe is a well-established model organism, studies of natural trait variations in this species remain limited. To assess the feasibility of segregant-pool-based mapping of phenotype-causing genes in natural strains of fission yeast, we investigated the cause of a maltose utilization defect (Mal(-)) of the S. pombe strain CBS5557 (originally known as Schizosaccharomyces malidevorans). Analyzing the genome sequence of CBS5557 revealed 955 nonconservative missense substitutions, and 61 potential loss-of-function variants including 47 frameshift indels, 13 early stop codons, and 1 splice site mutation. As a side benefit, our analysis confirmed 146 sequence errors in the reference genome and improved annotations of 27 genes. We applied bulk segregant analysis to map the causal locus of the Mal(-) phenotype. Through sequencing the segregant pools derived from a cross between CBS5557 and the laboratory strain, we located the locus to within a 2.23-Mb chromosome I inversion found in most S. pombe isolates including CBS5557. To map genes within the inversion region that occupies 18% of the genome, we created a laboratory strain containing the same inversion. Analyzing segregants from a cross between CBS5557 and the inversion-containing laboratory strain narrowed down the locus to a 200-kb interval and led us to identify agl1, which suffers a 5-bp deletion in CBS5557, as the causal gene. Interestingly, loss of agl1 through a 34-kb deletion underlies the Mal(-) phenotype of another S. pombe strain CGMCC2.1628. This work adapts and validates the bulk segregant analysis method for uncovering trait-gene relationship in natural fission yeast strains.

Targeted Sequencing Reveals Large-Scale Sequence Polymorphism in Maize Candidate Genes for Biomass Production and Composition

A major goal of maize genomic research is to identify sequence polymorphisms responsible for phenotypic variation in traits of economic importance. Large-scale detection of sequence variation is critical for linking genes, or genomic regions, to phenotypes. However, due to its size and complexity, it remains expensive to generate whole genome sequences of sufficient coverage for divergent maize lines, even with access to next generation sequencing (NGS) technology. Because methods involving reduction of genome complexity, such as genotyping-by-sequencing (GBS), assess only a limited fraction of sequence variation, targeted sequencing of selected genomic loci offers an attractive alternative. We therefore designed a sequence capture assay to target 29 Mb genomic regions and surveyed a total of 4,648 genes possibly affecting biomass production in 21 diverse inbred maize lines (7 flints, 14 dents). Captured and enriched genomic DNA was sequenced using the 454 NGS platform to 19.6-fold average depth coverage, and a broad evaluation of read alignment and variant calling methods was performed to select optimal procedures for variant discovery. Sequence alignment with the B73 reference and de novo assembly identified 383,145 putative single nucleotide polymorphisms (SNPs), of which 42,685 were non-synonymous alterations and 7,139 caused frameshifts. Presence/absence variation (PAV) of genes was also detected. We found that substantial sequence variation exists among genomic regions targeted in this study, which was particularly evident within coding regions. This diversification has the potential to broaden functional diversity and generate phenotypic variation that may lead to new adaptations and the modification of important agronomic traits. Further, annotated SNPs identified here will serve as useful genetic tools and as candidates in searches for phenotype-altering DNA variation. In summary, we demonstrated that sequencing of captured DNA is a powerful approach for variant discovery in maize genes.

Spectrum of variations in dog-1/FANCJ and mdf-1/MAD1 defective Caenorhabditis elegans strains after long-term propagation

Background

Whole and partial chromosome losses or gains and structural chromosome changes are hallmarks of human tumors. Guanine-rich DNA, which has a potential to form a G-quadruplex (G₄) structure, is particularly vulnerable to changes. In Caenorhabditis elegans, faithful transmission of G-rich DNA is ensured by the DOG-1/FANCJ deadbox helicase.

Results

To identify a spectrum of mutations, after long-term propagation, we combined whole genome sequencing (WGS) and oligonucleotide array Comparative Genomic Hybridization (oaCGH) analysis of a C. elegans strain that was propagated, in the absence of DOG-1 and MDF-1/MAD1, for a total of 470 generations, with samples taken for long term storage (by freezing) in generations 170 and 270. We compared the genomes of F₁₇₀ and F₄₇₀ strains and identified 94 substitutions, 17 InDels, 3 duplications, and 139 deletions larger than 20 bp. These homozygous variants were predicted to impact 101 protein-coding genes. Phenotypic analysis of this strain revealed remarkable fitness recovery indicating that mutations, which have accumulated in the strain, are not only tolerated but also cooperate to achieve long-term population survival in the absence of DOG-1 and MDF-1. Furthermore, deletions larger than 20 bp were the only variants that frequently occurred in G-rich DNA. We showed that 126 of the possible 954 predicted monoG/C tracts, larger than 14 bp, were deleted in unc-46 mdf-1 such-4; dog-1 F₄₇₀ (JNC170).

Conclusions

Here, we identified variants that accumulated in C. elegans’ genome after long-term propagation in the absence of DOG-1 and MDF-1. We showed that DNA sequences, with G₄-forming potential, are vulnerable to deletion-formation in this genetic background.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1402-y) contains supplementary material, which is available to authorized users.

Mutational landscape of intrahepatic cholangiocarcinoma

Intrahepatic cholangiocarcinoma (ICC) is a fatal primary liver cancer (PLC) that affects 5-10% of all PLCs. Here we sequence tumour and matching control sample pairs of a large cohort of 103 ICC patients in China, resulting in the identification of an ICC-specific somatic mutational signature that is associated with liver inflammation, fibrosis and cirrhosis. We further uncover 25 significantly mutated genes including eight potential driver genes (TP53, KRAS, IDH1, PTEN, ARID1A, EPPK1, ECE2 and FYN). We find that TP53-defective ICC patients are more likely to be HBsAg-seropositive, whereas mutations in the oncogene KRAS are nearly exclusively found in HBsAg-seronegative ICC patients. Three pathways (Ras/phosphatidylinositol-4,5-bisphosphate 3-kinase signalling, p53/cell cycle signalling and transforming growth factor-β/Smad signalling), genes important for epigenetic regulation and oxidative phosphorylation are substantially affected in ICC. We reveal mutations in this study that may be valuable for designing further studies, better diagnosis and effective therapies.

The genome of the chicken DT40 bursal lymphoma cell line

The chicken DT40 cell line is a widely used model system in the study of multiple cellular processes due to the efficiency of homologous gene targeting. The cell line was derived from a bursal lymphoma induced by avian leukosis virus infection. In this study we characterized the genome of the cell line using whole genome shotgun sequencing and single nucleotide polymorphism array hybridization. The results indicate that wild-type DT40 has a relatively normal karyotype, except for whole chromosome copy number gains, and no karyotype variability within stocks. In a comparison to two domestic chicken genomes and the Gallus gallus reference genome, we found no unique mutational processes shaping the DT40 genome except for a mild increase in insertion and deletion events, particularly deletions at tandem repeats. We mapped coding sequence mutations that are unique to the DT40 genome; mutations inactivating the PIK3R1 and ATRX genes likely contributed to the oncogenic transformation. In addition to a known avian leukosis virus integration in the MYC gene, we detected further integration sites that are likely to de-regulate gene expression. The new findings support the hypothesis that DT40 is a typical transformed cell line with a relatively intact genome; therefore, it is well-suited to the role of a model system for DNA repair and related processes. The sequence data generated by this study, including a searchable de novo genome assembly and annotated lists of mutated genes, will support future research using this cell line.

High-throughput capturing and characterization of mutations in essential genes of Caenorhabditis elegans

BACKGROUND

Essential genes are critical for the development of all organisms and are associated with many human diseases. These genes have been a difficult category to study prior to the availability of balanced lethal strains. Despite the power of targeted mutagenesis, there are limitations in identifying mutations in essential genes. In this paper, we describe the identification of coding regions for essential genes mutated using forward genetic screens in Caenorhabditis elegans. The lethal mutations described here were isolated and maintained by a wild-type allele on a rescuing duplication.

RESULTS

We applied whole genome sequencing to identify the causative molecular lesion resulting in lethality in existing C. elegans mutant strains. These strains are balanced and can be easily maintained for subsequent characterization. Our method can be effectively used to analyze mutations in a large number of essential genes. We describe here the identification of 64 essential genes in a region of chromosome I covered by the duplication sDp2. Of these, 42 are nonsense mutations, six are splice signal mutations, one deletion, and 15 are non-synonymous mutations. Many of the essential genes in this region function in cell cycle, transcriptional regulation, and RNA processing.

CONCLUSIONS

The essential genes identified here are represented by mutant strains, many of which have more than one mutant allele. The genetic resource can be utilized to further our understanding of essential gene function and will be applicable to the study of C. elegans development, conserved cellular function, and ultimately lead to improved human health.

Genome-wide variations in a natural isolate of the nematode Caenorhabditis elegans

Background

Increasing genetic and phenotypic differences found among natural isolates of C. elegans have encouraged researchers to explore the natural variation of this nematode species.

Results

Here we report on the identification of genomic differences between the reference strain N2 and the Hawaiian strain CB4856, one of the most genetically distant strains from N2. To identify both small- and large-scale genomic variations (GVs), we have sequenced the CB4856 genome using both Roche 454 (~400 bps single reads) and Illumina GA DNA sequencing methods (101 bps paired-end reads). Compared to previously described variants (available in WormBase), our effort uncovered twice as many single nucleotide variants (SNVs) and increased the number of small InDels almost 20-fold. Moreover, we identified and validated large insertions, most of which range from 150 bps to 1.2 kb in length in the CB4856 strain. Identified GVs had a widespread impact on protein-coding sequences, including 585 single-copy genes that have associated severe phenotypes of reduced viability in RNAi and genetics studies. Sixty of these genes are homologs of human genes associated with diseases. Furthermore, our work confirms previously identified GVs associated with differences in behavioural and biological traits between the N2 and CB4856 strains.

Conclusions

The identified GVs provide a rich resource for future studies that aim to explain the genetic basis for other trait differences between the N2 and CB4856 strains.