Identification of an EMS-induced causal mutation in a gene required for boron-mediated root development by low-coverage genome re-sequencing in Arabidopsis

Next-generation sequencing (NGS) technologies enable the rapid production of an enormous quantity of sequence data. These powerful new technologies allow the identification of mutations by whole-genome sequencing. However, most reported NGS-based mapping methods, which are based on bulked segregant analysis, are costly and laborious. To address these limitations, we designed a versatile NGS-based mapping method that consists of a combination of low- to medium-coverage multiplex SOLiD (Sequencing by Oligonucleotide Ligation and Detection) and classical genetic rough mapping. Using only low to medium coverage reduces the SOLiD sequencing costs and, since just 10 to 20 mutant F₂ plants are required for rough mapping, the operation is simple enough to handle in a laboratory with limited space and funding. As a proof of principle, we successfully applied this method to identify the CTR1, which is involved in boron-mediated root development, from among a population of high boron requiring Arabidopsis thaliana mutants. Our work demonstrates that this NGS-based mapping method is a moderately priced and versatile method that can readily be applied to other model organisms.     

Exploiting the Extraordinary Genetic Polymorphism of Ciona for Developmental Genetics with Whole Genome Sequencing

Studies in tunicates such as Ciona have revealed new insights into the evolutionary origins of chordate development. Ciona populations are characterized by high levels of natural genetic variation, between 1 and 5%. This variation has provided abundant material for forward genetic studies. In the current study, we make use of deep sequencing and homozygosity mapping to map spontaneous mutations in outbred populations. With this method we have mapped two spontaneous developmental mutants. In Ciona intestinalis we mapped a short-tail mutation with strong phenotypic similarity to a previously identified mutant in the related species Ciona savignyi. Our bioinformatic approach mapped the mutation to a narrow interval containing a single mutated gene, α-laminin3,4,5, which is the gene previously implicated in C. savignyi. In addition, we mapped a novel genetic mutation disrupting neural tube closure in C. savignyi to a T-type Ca²⁺ channel gene. The high efficiency and unprecedented mapping resolution of our study is a powerful advantage for developmental genetics in Ciona, and may find application in other outbred species.     

Independent FLC Mutations as Causes of Flowering-Time Variation in Arabidopsis thaliana and Capsella rubella

Capsella rubella is an inbreeding annual forb closely related to Arabidopsis thaliana, a model species widely used for studying natural variation in adaptive traits such as flowering time. Although mutations in dozens of genes can affect flowering of A. thaliana in the laboratory, only a handful of such genes vary in natural populations. Chief among these are FRIGIDA (FRI) and FLOWERING LOCUS C (FLC). Common and rare FRI mutations along with rare FLC mutations explain a large fraction of flowering-time variation in A. thaliana. Here we document flowering time under different conditions in 20 C. rubella accessions from across the species’ range. Similar to A. thaliana, vernalization, long photoperiods and elevated ambient temperature generally promote flowering. In this collection of C. rubella accessions, we did not find any obvious loss-of-function FRI alleles. Using mapping-by-sequencing with two strains that have contrasting flowering behaviors, we identified a splice-site mutation in FLC as the likely cause of early flowering in accession 1408. However, other similarly early C. rubella accessions did not share this mutation. We conclude that the genetic basis of flowering-time variation in C. rubella is complex, despite this very young species having undergone an extreme genetic bottleneck when it split from C. grandiflora a few tens of thousands of years ago.     

Single nucleotide polymorphisms generated by genotyping by sequencing to characterize genome-wide diversity,linkage disequilibrium,and selective sweeps in cultivated watermelon

Background

A large single nucleotide polymorphism (SNP) dataset was used to analyze genome-wide diversity in a diverse collection of watermelon cultivars representing globally cultivated, watermelon genetic diversity. The marker density required for conducting successful association mapping depends on the extent of linkage disequilibrium (LD) within a population. Use of genotyping by sequencing reveals large numbers of SNPs that in turn generate opportunities in genome-wide association mapping and marker-assisted selection, even in crops such as watermelon for which few genomic resources are available. In this paper, we used genome-wide genetic diversity to study LD, selective sweeps, and pairwise F_ST distributions among worldwide cultivated watermelons to track signals of domestication.

Results

We examined 183 Citrullus lanatus var. lanatus accessions representing domesticated watermelon and generated a set of 11,485 SNP markers using genotyping by sequencing. With a diverse panel of worldwide cultivated watermelons, we identified a set of 5,254 SNPs with a minor allele frequency of ≥ 0.05, distributed across the genome. All ancestries were traced to Africa and an admixture of various ancestries constituted secondary gene pools across various continents. A sliding window analysis using pairwise F_ST values was used to resolve selective sweeps. We identified strong selection on chromosomes 3 and 9 that might have contributed to the domestication process. Pairwise analysis of adjacent SNPs within a chromosome as well as within a haplotype allowed us to estimate genome-wide LD decay. LD was also detected within individual genes on various chromosomes. Principal component and ancestry analyses were used to account for population structure in a genome-wide association study. We further mapped important genes for soluble solid content using a mixed linear model.

Conclusions

Information concerning the SNP resources, population structure, and LD developed in this study will help in identifying agronomically important candidate genes from the genomic regions underlying selection and for mapping quantitative trait loci using a genome-wide association study in sweet watermelon.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-767) contains supplementary material, which is available to authorized users.     

An ultra-dense library resource for rapid deconvolution of mutations that cause phenotypes in Escherichia coli

With the wide availability of whole-genome sequencing (WGS), genetic mapping has become the rate-limiting step, inhibiting unbiased forward genetics in even the most tractable model organisms. We introduce a rapid deconvolution resource and method for untagged causative mutations after mutagenesis, screens, and WGS in Escherichia coli. We created Deconvoluter—ordered libraries with selectable insertions every 50 kb in the E. coli genome. The Deconvoluter method uses these for replacement of untagged mutations in the genome using a phage-P1-based gene-replacement strategy. We validate the Deconvoluter resource by deconvolution of 17 of 17 phenotype-altering mutations from a screen of N-ethyl-N-nitrosourea-induced mutants. The Deconvoluter resource permits rapid unbiased screens and gene/function identification and will enable exploration of functions of essential genes and undiscovered genes/sites/alleles not represented in existing deletion collections. This resource for unbiased forward-genetic screens with mapping-by-sequencing (‘forward genomics’) demonstrates a strategy that could similarly enable rapid screens in many other microbes.     

Caenorhabditis elegans: A Genetic Guide to Parasitic Nematode Biology

The advent of parasite genome sequencing projects, as well as an increase in biology-directed gene discovery, promises to reveal genes encoding many of the key molecules required for nematode-host interactions. However, distinguishing parasitism genes from those merely required for nematode viability remains a substantial challenge. Although this will ultimately require a functional test in the host or parasite, the free-living nematode Caenorhabditis elegans can be exploited as a heterologous system to determine function of candidate parasitism genes. Studies of C. elegans also have revealed genetic networks, such as the dauer pathway, that may also be important adaptations for parasitism. As a more directed means of identifying parasitism traits, we developed classical genetics for Heterodera glycines and have used this approach to map genes conferring host resistance-breaking phenotypes. It is likely that the C. elegans and H. glycines genomes will be at least partially syntenic, thus permitting predictive physical mapping of H. glycines genes of interest.     

Fast forward genetics to identify mutations causing a high light tolerant phenotype in Chlamydomonas reinhardtii by whole-genome-sequencing

Background

High light tolerance of microalgae is a desired phenotype for efficient cultivation in large scale production systems under fluctuating outdoor conditions. Outdoor cultivation requires the use of either wild-type or non-GMO derived mutant strains due to safety concerns. The identification and molecular characterization of such mutants derived from untagged forward genetics approaches was limited previously by the tedious and time-consuming methods involving techniques such as classical meiotic mapping. The combination of mapping with next generation sequencing technologies offers alternative strategies to identify genes involved in high light adaptation in untagged mutants.

Results

We used the model alga Chlamydomonas reinhardtii in a non-GMO mutation strategy without any preceding crossing step or pooled progeny to identify genes involved in the regulatory processes of high light adaptation. To generate high light tolerant mutants, wildtype cells were mutagenized only to a low extent, followed by a stringent selection. We performed whole-genome sequencing of two independent mutants hit1 and hit2 and the parental wildtype. The availability of a reference genome sequence and the removal of shared bakground variants between the wildtype strain and each mutant, enabled us to identify two single nucleotide polymorphisms within the same gene Cre02.g085050, hereafter called LRS1 (putative Light Response Signaling protein 1). These two independent single amino acid exchanges are both located in the putative WD40 propeller domain of the corresponding protein LRS1. Both mutants exhibited an increased rate of non-photochemical-quenching (NPQ) and an improved resistance against chemically induced reactive oxygen species. In silico analyses revealed homology of LRS1 to the photoregulatory protein COP1 in plants.

Conclusions

In this work we identified the nuclear encoded gene LRS1 as an essential factor for high light adaptation in C. reinhardtii. The causative random mutation within this gene was identified by a rapid and efficient method, avoiding any preceding crossing step, meiotic mapping, or pooled progeny. Our results open up new insights into mechanisms of high light adaptation in microalgae and at the same time provide a simplified strategy for non-GMO forward genetics, a crucial precondition that could result in the identification of key factors for economically relevant biological processes within algae.     

PHYTOCHROME C Is an Essential Light Receptor for Photoperiodic Flowering in the Temperate Grass,Brachypodium distachyon

We show that in the temperate grass, Brachypodium distachyon, PHYTOCHROME C (PHYC), is necessary for photoperiodic flowering. In loss-of-function phyC mutants, flowering is extremely delayed in inductive photoperiods. PHYC was identified as the causative locus by utilizing a mapping by sequencing pipeline (Cloudmap) optimized for identification of induced mutations in Brachypodium. In phyC mutants the expression of Brachypodium homologs of key flowering time genes in the photoperiod pathway such as GIGANTEA (GI), PHOTOPERIOD 1 (PPD1/PRR37), CONSTANS (CO), and florigen/FT are greatly attenuated. PHYC also controls the day-length dependence of leaf size as the effect of day length on leaf size is abolished in phyC mutants. The control of genes upstream of florigen production by PHYC was likely to have been a key feature of the evolution of a long-day flowering response in temperate pooid grasses.     

A novel termini analysis theory using HTS data alone for the identification of Enterococcus phage EF4-like genome termini

Background

Enterococcus faecalis and Enterococcus faecium are typical enterococcal bacterial pathogens. Antibiotic resistance means that the identification of novel E. faecalis and E. faecium phages against antibiotic-resistant Enterococcus have an important impact on public health. In this study, the E. faecalis phage IME-EF4, E. faecium phage IME-EFm1, and both their hosts were antibiotic resistant. To characterize the genome termini of these two phages, a termini analysis theory was developed to provide a wealth of terminal sequence information directly, using only high-throughput sequencing (HTS) read frequency statistics.

Results

The complete genome sequences of phages IME-EF4 and IME-EFm1 were determined, and our termini analysis theory was used to determine the genome termini of these two phages. Results showed 9 bp 3′ protruding cohesive ends in both IME-EF4 and IME-EFm1 genomes by analyzing frequencies of HTS reads. For the positive strands of their genomes, the 9 nt 3′ protruding cohesive ends are 5′-TCATCACCG-3′ (IME-EF4) and 5′-GGGTCAGCG-3′ (IME-EFm1). Further experiments confirmed these results. These experiments included mega-primer polymerase chain reaction sequencing, terminal run-off sequencing, and adaptor ligation followed by run-off sequencing.

Conclusion

Using this termini analysis theory, the termini of two newly isolated antibiotic-resistant Enterococcus phages, IME-EF4 and IME-EFm1, were identified as the byproduct of HTS. Molecular biology experiments confirmed the identification. Because it does not require time-consuming wet lab termini analysis experiments, the termini analysis theory is a fast and easy means of identifying phage DNA genome termini using HTS read frequency statistics alone. It may aid understanding of phage DNA packaging.     

Simple Method for Fluorescence DNA In Situ Hybridization to Squashed Chromosomes

DNA in situ hybridization (DNA ISH) is a commonly used method for mapping sequences to specific chromosome regions. This approach is particularly effective at mapping highly repetitive sequences to heterochromatic regions, where computational approaches face prohibitive challenges. Here we describe a streamlined protocol for DNA ISH that circumvents formamide washes that are standard steps in other DNA ISH protocols. Our protocol is optimized for hybridization with short single strand DNA probes that carry fluorescent dyes, which effectively mark repetitive DNA sequences within heterochromatic chromosomal regions across a number of different insect tissue types. However, applications may be extended to use with larger probes and visualization of single copy (non-repetitive) DNA sequences. We demonstrate this method by mapping several different repetitive sequences to squashed chromosomes from Drosophila melanogaster neural cells and Nasonia vitripennis spermatocytes. We show hybridization patterns for both small, commercially synthesized probes and for a larger probe for comparison. This procedure uses simple laboratory supplies and reagents, and is ideal for investigators who have little experience with performing DNA ISH.     

Genetic mapping of two components of reproductive isolation between two sibling species of moths,Ostrinia nubilalis and O. scapulalis

We report the quantitative trait loci (QTL) mapping of reproductive isolation traits between Ostrinia nubilalis (the European corn borer) and its sibling species O. scapulalis (the Adzuki bean borer), focusing on two traits: mating isolation (mi) and pheromone production (Pher). Four genetic maps were generated from two backcross families, with two maps (one chromosomal map and one linkage map) per backcross. We located 165–323 AFLP markers on these four maps, resulting in the identification of 27–31 linkage groups, depending on the map considered. No-choice mating experiments with the offspring of each backcross led to the detection of at least two QTLs for mi in different linkage groups. QTLs underlying Pher were located in a third linkage group. The Z heterochromosome was identified by a specific marker (Tpi) and did not carry any of these QTLs. Finally, we considered the global divergence between the two sibling species, distortions of segregation throughout the genome, and the location and effect of mi and Pher QTLs in light of the known candidate genes for reproductive isolation within the genus Ostrinia and, more broadly, in phytophagous insects.     

Polyphasic characterisation of Burkholderia cepacia complex species isolated from children with cystic fibrosis

Cystic fibrosis (CF) patients with Burkholderia cepacia complex (Bcc) pulmonary infections have high morbidity and mortality. The aim of this study was to compare different methods for identification of Bcc species isolated from paediatric CF patients. Oropharyngeal swabs from children with CF were used to obtain isolates of Bcc samples to evaluate six different tests for strain identification. Conventional (CPT) and automatised (APT) phenotypic tests, polymerase chain reaction (PCR)-recA, restriction fragment length polymorphism-recA, recAsequencing, and matrix-assisted laser desorption ionization-time of flight (MALDI-TOF) were applied. Bacterial isolates were also tested for antimicrobial susceptibility. PCR-recA analysis showed that 36 out of the 54 isolates were Bcc. Kappa index data indicated almost perfect agreement between CPT and APT, CPT and PCR-recA, and APT and PCR-recA to identify Bcc, and MALDI-TOF and recAsequencing to identify Bcc species. The recAsequencing data and the MALDI-TOF data agreed in 97.2% of the isolates. Based on recA sequencing, the most common species identified were Burkholderia cenocepacia IIIA (33.4%),Burkholderia vietnamiensis (30.6%), B. cenocepaciaIIIB (27.8%), Burkholderia multivorans (5.5%), and B. cepacia (2.7%). MALDI-TOF proved to be a useful tool for identification of Bcc species obtained from CF patients, although it was not able to identify B. cenocepacia subtypes.     

Editor's Choice: Polygenic Molecular Architecture Underlying Non-Sexual Cell Aggregation in Budding Yeast

Cell aggregation in unicellular organisms, induced by either cell non-sexual adhesion to yield flocs and biofilm, or pheromone-driving sexual conjugation is of great significance in cellular stress response, medicine, and brewing industries. Most current literatures have focused on one form of cell aggregation termed flocculation and its major molecular determinants, the flocculation (FLO) family genes. Here, we implemented a map-based approach for dissecting the molecular basis of non-sexual cell aggregation in Saccharomyces cerevisiae. Genome-wide mapping has identified four major quantitative trait loci (QTL) underlying nature variation in the cell aggregation phenotype. High-resolution mapping following up with knockout and allele replacement experiments resolved the QTL into the underlying genes (AMN1, RGA1, FLO1, and FLO8) or even into the causative nucleotide. Genetic variation in the QTL genes can explain up to 46% of phenotypic variation of this trait. Of these genes, AMN1 plays the leading role, differing from the FLO family members, in regulating expression of cell clumping phenotype through inducing cell segregation defect. These findings provide novel insights into the molecular mechanism of how cell aggregation is regulated in budding yeast, and the data will be directly implicated to understand the molecular basis and evolutionary implications of cell aggregation in other fungus species.     

Fluorescent in situ Hybridization on Mitotic Chromosomes of Mosquitoes

Fluorescent in situ hybridization (FISH) is a technique routinely used by many laboratories to determine the chromosomal position of DNA and RNA probes. One important application of this method is the development of high-quality physical maps useful for improving the genome assemblies for various organisms. The natural banding pattern of polytene and mitotic chromosomes provides guidance for the precise ordering and orientation of the genomic supercontigs. Among the three mosquito genera, namely Anopheles, Aedes, and Culex, a well-established chromosome-based mapping technique has been developed only for Anopheles, whose members possess readable polytene chromosomes ¹. As a result of genome mapping efforts, 88% of the An. gambiae genome has been placed to precise chromosome positions ^2,3 . Two other mosquito genera, Aedes and Culex, have poorly polytenized chromosomes because of significant overrepresentation of transposable elements in their genomes ^{4, 5, 6}. Only 31 and 9% of the genomic supercontings have been assigned without order or orientation to chromosomes of Ae. aegypti ⁷ and Cx. quinquefasciatus ⁸, respectively. Mitotic chromosome preparation for these two species had previously been limited to brain ganglia and cell lines. However, chromosome slides prepared from the brain ganglia of mosquitoes usually contain low numbers of metaphase plates ⁹. Also, although a FISH technique has been developed for mitotic chromosomes from a cell line of Ae. aegypti ¹⁰, the accumulation of multiple chromosomal rearrangements in cell line chromosomes ¹¹ makes them useless for genome mapping. Here we describe a simple, robust technique for obtaining high-quality mitotic chromosome preparations from imaginal discs (IDs) of 4^th instar larvae which can be used for all three genera of mosquitoes. A standard FISH protocol ¹² is optimized for using BAC clones of genomic DNA as a probe on mitotic chromosomes of Ae. aegypti and Cx. quinquefasciatus, and for utilizing an intergenic spacer (IGS) region of ribosomal DNA (rDNA) as a probe on An. gambiae chromosomes. In addition to physical mapping, the developed technique can be applied to population cytogenetics and chromosome taxonomy/systematics of mosquitoes and other insect groups.