Targeted exome sequencing of unselected heavy‐ion beam‐irradiated populations reveals less‐biased mutation characteristics in the rice genome

Heavy‐ion beams have been widely utilized as a novel and effective mutagen for mutation breeding in diverse plant species, but the induced mutation spectrum is not fully understood at the genome scale. We describe the development of a multiplexed and cost‐efficient whole‐exome sequencing procedure in rice, and its application to characterize an unselected population of heavy‐ion beam‐induced mutations. The bioinformatics pipeline identified single‐nucleotide mutations as well as small and large (>63 kb) insertions and deletions, and showed good agreement with the results obtained with conventional polymerase chain reaction (PCR) and sequencing analyses. We applied the procedure to analyze the mutation spectrum induced by heavy‐ion beams at the population level. In total, 165 individual M₂ lines derived from six irradiation conditions as well as eight pools from non‐irradiated ‘Nipponbare’ controls were sequenced using the newly established target exome sequencing procedure. The characteristics and distribution of carbon‐ion beam‐induced mutations were analyzed in the absence of bias introduced by visual mutant selections. The average (±SE) number of mutations within the target exon regions was 9.06 ± 0.37 induced by 150 Gy irradiation of dry seeds. The mutation frequency changed in parallel to the irradiation dose when dry seeds were irradiated. The total number of mutations detected by sequencing unselected M₂ lines was correlated with the conventional mutation frequency determined by the occurrence of morphological mutants. Therefore, mutation frequency may be a good indicator for sequencing‐based determination of the optimal irradiation condition for induction of mutations.     

Genome‐wide genetic variation discovery in Chinese Taihu pig breeds using next generation sequencing

The Chinese Taihu pig breeds are an invaluable component of the world's pig genetic resources, and they are the most prolific breeds of swine in the world. In this study, the genomes of 252 pigs of the six indigenous breeds in the Taihu Lake region were sequenced using the genotyping by genome reducing and sequencing approach. A total of 950 million good reads were obtained using an Illumina Hiseq2000 at an average depth of 13× (for SNP calling) and an average coverage of 2.3%. In total, 122 632 indels, 31 444 insertions, 44 056 deletions and 455 CNVs (copy number variants) were identified in the genomes of the pigs. Approximately 2.3% of these genetic markers were mapped to gene exon regions, and 25% were in QTL regions related to economically important traits. The KEGG pathway or GO enrichment analyses revealed that genetic variants assumed to be large‐effect mutations were significantly overrepresented in 22 SNP, 56 indel, 26 insertion, 28 deletion and three CNV gene sets. A total of 343 breed‐specific SNPs were also identified in the six Chinese indigenous pigs. The findings from this study can contribute to future investigations of the genetic diversity, population structure, positive selection signals and molecular evolutionary history of these pigs at the genome level and can serve as a valuable reference for improving the breeding and cultivation of these pigs.     

A macaque's-eye view of human insertions and deletions: differences in mechanisms

Insertions and deletions (indels) cause numerous genetic diseases and lead to pronounced evolutionary differences among genomes. The macaque sequences provide an opportunity to gain insights into the mechanisms generating these mutations on a genome-wide scale by establishing the polarity of indels occurring in the human lineage since its divergence from the chimpanzee. Here we apply novel regression techniques and multiscale analyses to demonstrate an extensive regional indel rate variation stemming from local fluctuations in divergence, GC content, male and female recombination rates, proximity to telomeres, and other genomic factors. We find that both replication and, surprisingly, recombination are significantly associated with the occurrence of small indels. Intriguingly, the relative inputs of replication versus recombination differ between insertions and deletions, thus the two types of mutations are likely guided in part by distinct mechanisms. Namely, insertions are more strongly associated with factors linked to recombination, while deletions are mostly associated with replication-related features. Indel as a term misleadingly groups the two types of mutations together by their effect on a sequence alignment. However, here we establish that the correct identification of a small gap as an insertion or a deletion (by use of an outgroup) is crucial to determining its mechanism of origin. In addition to providing novel insights into insertion and deletion mutagenesis, these results will assist in gap penalty modeling and eventually lead to more reliable genomic alignments.     

Different mutational function of low‐ and high‐linear energy transfer heavy‐ion irradiation demonstrated by whole‐genome resequencing of Arabidopsis mutants

Heavy‐ion irradiation is a powerful mutagen that possesses high linear energy transfer (LET). Several studies have indicated that the value of LET affects DNA lesion formation in several ways, including the efficiency and the density of double‐stranded break induction along the particle path. We assumed that the mutation type can be altered by selecting an appropriate LET value. Here, we quantitatively demonstrate differences in the mutation type induced by irradiation with two representative ions, Ar ions (LET: 290 keV μm^?1) and C ions (LET: 30.0 keV μm^?1), by whole‐genome resequencing of the Arabidopsis mutants produced by these irradiations. Ar ions caused chromosomal rearrangements or large deletions (≥100 bp) more frequently than C ions, with 10.2 and 2.3 per mutant genome under Ar‐ and C‐ion irradiation, respectively. Conversely, C ions induced more single‐base substitutions and small indels (<100 bp) than Ar ions, with 28.1 and 56.9 per mutant genome under Ar‐ and C‐ion irradiation, respectively. Moreover, the rearrangements induced by Ar‐ion irradiation were more complex than those induced by C‐ion irradiation, and tended to accompany single base substitutions or small indels located close by. In conjunction with the detection of causative genes through high‐throughput sequencing, selective irradiation by beams with different effects will be a powerful tool for forward genetics as well as studies on chromosomal rearrangements.     

Comprehensive identification of mutations induced by heavy‐ion beam irradiation in Arabidopsis thaliana

Heavy‐ion beams are widely used for mutation breeding and molecular biology. Although the mutagenic effects of heavy‐ion beam irradiation have been characterized by sequence analysis of some restricted chromosomal regions or loci, there have been no evaluations at the whole‐genome level or of the detailed genomic rearrangements in the mutant genomes. In this study, using array comparative genomic hybridization (array‐CGH) and resequencing, we comprehensively characterized the mutations in Arabidopsis thaliana genomes irradiated with Ar or Fe ions. We subsequently used this information to investigate the mutagenic effects of the heavy‐ion beams. Array‐CGH demonstrated that the average number of deleted areas per genome were 1.9 and 3.7 following Ar‐ion and Fe‐ion irradiation, respectively, with deletion sizes ranging from 149 to 602 180 bp; 81% of the deletions were accompanied by genomic rearrangements. To provide a further detailed analysis, the genomes of the mutants induced by Ar‐ion beam irradiation were resequenced, and total mutations, including base substitutions, duplications, in/dels, inversions, and translocations, were detected using three algorithms. All three resequenced mutants had genomic rearrangements. Of the 22 DNA fragments that contributed to the rearrangements, 19 fragments were responsible for the intrachromosomal rearrangements, and multiple rearrangements were formed in the localized regions of the chromosomes. The interchromosomal rearrangements were detected in the multiply rearranged regions. These results indicate that the heavy‐ion beams led to clustered DNA damage in the chromosome, and that they have great potential to induce complicated intrachromosomal rearrangements. Heavy‐ion beams will prove useful as unique mutagens for plant breeding and the establishment of mutant lines.     

Genome‐wide survey of artificial mutations induced by ethyl methanesulfonate and gamma rays in tomato

Genome‐wide mutations induced by ethyl methanesulfonate (EMS) and gamma irradiation in the tomato Micro‐Tom genome were identified by a whole‐genome shotgun sequencing analysis to estimate the spectrum and distribution of whole‐genome DNA mutations and the frequency of deleterious mutations. A total of ~370 Gb of paired‐end reads for four EMS‐induced mutants and three gamma‐ray‐irradiated lines as well as a wild‐type line were obtained by next‐generation sequencing technology. Using bioinformatics analyses, we identified 5920 induced single nucleotide variations and insertion/deletion (indel) mutations. The predominant mutations in the EMS mutants were C/G to T/A transitions, while in the gamma‐ray mutants, C/G to T/A transitions, A/T to T/A transversions, A/T to G/C transitions and deletion mutations were equally common. Biases in the base composition flanking mutations differed between the mutagenesis types. Regarding the effects of the mutations on gene function, >90% of the mutations were located in intergenic regions, and only 0.2% were deleterious. In addition, we detected 1 140 687 spontaneous single nucleotide polymorphisms and indel polymorphisms in wild‐type Micro‐Tom lines. We also found copy number variation, deletions and insertions of chromosomal segments in both the mutant and wild‐type lines. The results provide helpful information not only for mutation research, but also for mutant screening methodology with reverse‐genetic approaches.     

Zinc finger nuclease‐mediated targeting of multiple transgenes to an endogenous soybean genomic locus via non‐homologous end joining

Emerging genome editing technologies hold great promise for the improvement of agricultural crops. Several related genome editing methods currently in development utilize engineered, sequence‐specific endonucleases to generate DNA double strand breaks (DSBs) at user‐specified genomic loci. These DSBs subsequently result in small insertions/deletions (indels), base substitutions or incorporation of exogenous donor sequences at the target site, depending on the application. Targeted mutagenesis in soybean (Glycine max) via non‐homologous end joining (NHEJ)‐mediated repair of such DSBs has been previously demonstrated with multiple nucleases, as has homology‐directed repair (HDR)‐mediated integration of a single transgene into target endogenous soybean loci using CRISPR/Cas9. Here we report targeted integration of multiple transgenes into a single soybean locus using a zinc finger nuclease (ZFN). First, we demonstrate targeted integration of biolistically delivered DNA via either HDR or NHEJ to the FATTY ACID DESATURASE 2‐1a (FAD2‐1a) locus of embryogenic cells in tissue culture. We then describe ZFN‐ and NHEJ‐mediated, targeted integration of two different multigene donors to the FAD2‐1a locus of immature embryos. The largest donor delivered was 16.2 kb, carried four transgenes, and was successfully transmitted to T₁ progeny of mature targeted plants obtained via somatic embryogenesis. The insertions in most plants with a targeted, 7.1 kb, NHEJ‐integrated donor were perfect or near‐perfect, demonstrating that NHEJ is a viable alternative to HDR for gene targeting in soybean. Taken together, these results show that ZFNs can be used to generate fertile transgenic soybean plants with NHEJ‐mediated targeted insertions of multigene donors at an endogenous genomic locus.     

Genetic variation and genetic structure of five Chinese indigenous pig populations in Jiangsu Province revealed by sequencing data

In this study, we investigated the genetic variants, including SNPs and indels (short insertions or deletions, less than 50 bp in length), in the genomes and genetic structures of five pig populations (in the northern Taihu Lake region, Jiangsu Province) using the genotyping by genome reducing and sequencing (GGRS) approach. A total of 581 million good reads with an average depth of 11× and an average coverage of 2.16% were used to call variants. In general, 202 106 SNPs and 34 415 indels were obtained, of which 2690 SNPs and 224 indels were capable of inducing protein‐coding changes. The genes containing these variants were extracted for functional annotation. The results of gene enrichment analysis revealed that the SNPs under investigation may be associated with reproduction, disease resistance, meat quality and adipose tissue traits, whereas the indels were associated mainly with adipose tissue and disease. Analysis of the genetic structure showed that each population displayed comparable, large differentiations from the others, indicating their uniqueness. In conclusion, the results of our study provide the first genomic overview of the genetic variants and population structures of five Chinese indigenous pig populations.     

SHORT INDELS ARE SUBJECT TO INSERTION‐BIASED GENE CONVERSION

Recombination between homologous loci is accompanied by formation of heteroduplexes. Repairing mismatches in heteroduplexes often leads to single nucleotide substitutions in a process known as gene conversion. Gene conversion was shown to be GC‐biased in different organisms; that is, a W(A or T)→S(G or C) substitution is more likely in this process than a S→W substitution. Here, we show that the insertion/deletion ratio for short noncoding indels that reach fixation between species is positively correlated with the recombination rate in Drosophila melanogaster, Homo sapiens, and Saccharomyces cerevisiae. This correlation is both due to an increase of the fixation rate of insertions and decrease of the fixation rate of deletions in regions of high recombination. Whole‐genome data on indel polymorphism and divergence in D. melanogaster rule out mutation biases and selection as the cause of this trend, pointing to insertion‐biased gene conversion as the most likely explanation. The bias toward insertions is the strongest for single‐nucleotide indels, and decreases with indel length. In regions of high recombination rate this bias leads to an up to ～5‐fold excess of fixed short insertions over deletions, and substantially affects the evolution of DNA segments.     

A chromosome-scale genome assembly of Malus domestica,a multi-stress resistant apple variety

Hanfu apple is the main cultivar grown in the cool areas of Northeast, Northwest, and North China. Here, we proposed a chromosome-level Hanfu genome assembly using PacBio, Illumina and Hi-C sequencing data. The total contig length was 628.99 Mb, with scaffold and contig N50 sizes of 36.18 Mb and 1.25 Mb, respectively. The Hanfu genome had a total of 39,617 genes, of which we predicted the function for 38,816. Evolutionary analysis showed that Hanfu may have undergone a γ-event, a recent whole-genome duplication. A comparative analysis was conducted on the genomes of Hanfu and homozygous triploid HFTH1, which were cultured using the anthers of diploid Hanfu apples. Three variants were identified, including 2,155,184 single nucleotide polymorphisms (SNPs), 413,108 insertions/deletions (indels), and 7,587 structural variants (SVs).This high-quality genome will provide a reference for the genetic improvement of apples and the breeding of more varieties with high resistance and high quality.     

The majority of recent short DNA insertions in the human genome are tandem duplications

Nucleotide substitutions, insertions, and deletions constitute the principal molecular mechanisms generating genetic variation on small length scales. In contrast to substitutions, the nature of short DNA insertions and deletions (indels) is far less understood. With the recent availability of whole-genome multiple alignments between human and other primates, detailed investigations on indel characteristics and origin have come within reach. Here, we show that the majority of short (1-100 bp) DNA insertions in the human lineage are tandem duplications of directly adjacent sequence segments with conserved polarity. Indels in microsatellites comprise only a small fraction. The underlying molecular processes generating indels do not necessarily rely on the presence of preexisting duplicates, as would be expected for unequal crossing over, as well as replication slippage. Instead, our findings point toward a mechanism that preferentially occurs in the male germline and is not recombination-mediated. Surprisingly, nonframeshifting tandem duplications and deletions in coding regions still occur at approximately 50% of their genomic background rates. As is already well established in the context of gene and segmental duplications, our results demonstrate that duplications are also likely to constitute the predominant process for rapid generation of new genetic material and function on smaller scales.     

Robust estimation of survival and contribution of captive‐bred Mallards Anas platyrhynchos to a wild population in a large‐scale release programme

The survival of captive‐bred individuals from release into the wild to their first breeding season is crucial to assess the success of reintroduction or translocation programmes, and to assess their potential impact of wild populations. However, assessing the survival of captive‐bred individuals following their release is often complicated by immediate dispersal once in the wild. Here, we apply Lindberg's robust design model, a method that incorporates emigration from the study site, to obtain true estimates of survival of captive‐bred Mallards Anas platyrhynchos, a common duck species released on a large scale in Europe since the 1970s. Overall survival rate from release in July until the onset of the next breeding season in April was low (0.18 ± 0.07 se) and equivalent to half the first‐year survival of local wild Mallards. Higher overall detectability and temporary emigration during the hunting period revealed movements in response to hunting pressure. Such low survival of released Mallards during their first year may help prevent large‐scale genetic mixing with the wild population. Nevertheless, by combining our results with regional waterfowl counts, we estimated that a minimum of 34% of the Mallards in the region were of captive origin at the onset of the breeding season. Although most released birds quickly die, restocking for hunting may be of sufficient magnitude to affect the wild population through genetic homogenization or loss of local adaptation. Robust design protocols allow for the estimation of true survival estimates by controlling for permanent and temporary emigration and may require only a moderate increase in fieldwork effort.     

Genetic variation in Miscanthus × giganteus and the importance of estimating genetic distance thresholds for differentiating clones

Miscanthus × giganteus (Mxg) is an important bioenergy feedstock crop, however, genetic diversity among legacy cultivars may be severely constrained. Only one introduction from Japan to Denmark of this sterile, triploid, vegetatively propagated crop was recorded in the 1930s. We sought to determine if the Mxg cultivars in North America were all synonyms, and if they were derived from the European introduction. We used 64 nuclear and five chloroplast simple sequence repeat (SSR) markers to estimate genetic similarity for 27 Mxg accessions from North America, and compared them with six accessions from Europe, including the species’ type‐specimen. A subset of accessions was also evaluated by restriction‐site associated DNA sequencing (RAD‐seq). In addition, we assessed the potential of new crosses to increase Mxg genetic diversity by comparing eight new triploid Mxg progeny grown from seed, along with samples of the parental species M. sacchariflorus and M. sinensis. Estimates of genotyping error rates were essential for distinguishing between experimental error and true genotypic differences among accessions. Given differences in estimated error rates and costs per marker for SSRs and RAD‐seq, the former is currently more cost‐effective for determining if two accessions are genetically identical. We concluded that all of the Mxg legacy cultivars were derived via vegetative propagation from a single genet. In contrast with the Mxg legacy cultivars, genetic similarity to the type‐specimen of eight new triploid Mxg progeny ranged from 0.46 to 0.56. Though genetic diversity among the Mxg legacy cultivars is critically low, new crosses can provide much‐needed variation to growers.     

Genetic changes in a novel breeding population of Brassica napus synthesized from hundreds of crosses between B. rapa and B. carinata

Introgression of genomic variation between and within related crop species is a significant evolutionary approach for population differentiation, genome reorganization and trait improvement. Using the Illumina Infinium Brassica 60K SNP array, we investigated genomic changes in a panel of advanced generation new‐type Brassica napus breeding lines developed from hundreds of interspecific crosses between 122 Brassica rapa and 74 Brassica carinata accessions, and compared them with representative accessions of their three parental species. The new‐type B. napus population presented rich genetic diversity and abundant novel genomic alterations, consisting of introgressions from B. rapa and B. carinata, novel allelic combinations, reconstructed linkage disequilibrium patterns and haplotype blocks, and frequent deletions and duplications (nonrandomly distributed), particularly in the C subgenome. After a much shorter, but very intensive, selection history compared to traditional B. napus, a total of 15 genomic regions with strong selective sweeps and 112 genomic regions with putative signals of selective sweeps were identified. Some of these regions were associated with important agronomic traits that were selected for during the breeding process, while others were potentially associated with restoration of genome stability and fertility after interspecific hybridization. Our results demonstrate how a novel method for population‐based crop genetic improvement can lead to rapid adaptation, restoration of genome stability and positive responses to artificial selection.     

A strong deletion bias in nonallelic gene conversion

Gene conversion is the unidirectional transfer of genetic information between orthologous (allelic) or paralogous (nonallelic) genomic segments. Though a number of studies have examined nucleotide replacements, little is known about length difference mutations produced by gene conversion. Here, we investigate insertions and deletions produced by nonallelic gene conversion in 338 Drosophila and 10,149 primate paralogs. Using a direct phylogenetic approach, we identify 179 insertions and 614 deletions in Drosophila paralogs, and 132 insertions and 455 deletions in primate paralogs. Thus, nonallelic gene conversion is strongly deletion-biased in both lineages, with almost 3.5 times as many conversion-induced deletions as insertions. In primates, the deletion bias is considerably stronger for long indels and, in both lineages, the per-site rate of gene conversion is orders of magnitudes higher than that of ordinary mutation. Due to this high rate, deletion-biased nonallelic gene conversion plays a key role in genome size evolution, leading to the cooperative shrinkage and eventual disappearance of selectively neutral paralogs.     

DNA indels in coding regions reveal selective constraints on protein evolution in the human lineage

Background  

Insertions and deletions of DNA segments (indels) are together with substitutions the major mutational processes that generate genetic variation. Here we focus on recent DNA insertions and deletions in protein coding regions of the human genome to investigate selective constraints on indels in protein evolution.     

Next‐generation sequencing (NGS)‐based identification of induced mutations in a doubly mutagenized tomato (Solanum lycopersicum) population

The identification of mutations in targeted genes has been significantly simplified by the advent of TILLING (Targeting Induced Local Lesions In Genomes), speeding up the functional genomic analysis of animals and plants. Next‐generation sequencing (NGS) is gradually replacing classical TILLING for mutation detection, as it allows the analysis of a large number of amplicons in short durations. The NGS approach was used to identify mutations in a population of Solanum lycopersicum (tomato) that was doubly mutagenized by ethylmethane sulphonate (EMS). Twenty‐five genes belonging to carotenoids and folate metabolism were PCR‐amplified and screened to identify potentially beneficial alleles. To augment efficiency, the 600‐bp amplicons were directly sequenced in a non‐overlapping manner in Illumina MiSeq, obviating the need for a fragmentation step before library preparation. A comparison of the different pooling depths revealed that heterozygous mutations could be identified up to 128‐fold pooling. An evaluation of six different software programs (camba , crisp , gatk unified genotyper , lofreq , snver and vipr ) revealed that no software program was robust enough to predict mutations with high fidelity. Among these, crisp and camba predicted mutations with lower false discovery rates. The false positives were largely eliminated by considering only mutations commonly predicted by two different software programs. The screening of 23.47 Mb of tomato genome yielded 75 predicted mutations, 64 of which were confirmed by Sanger sequencing with an average mutation density of 1/367 Kb. Our results indicate that NGS combined with multiple variant detection tools can reduce false positives and significantly speed up the mutation discovery rate.     

Gamma irradiation‐induced disease resistance of pear (Pyrus pyrifolia “Niitaka”) against Penicillium expansum

In this study, the effects of gamma irradiation on the resistance of pear fruit against Penicillium expansum, the causal agent of blue mould disease, were investigated. A low dose of gamma irradiation for 14 days increased the disease resistance and firmness of pear fruits. Remarkably, exposure to 200 Gy of gamma irradiation significantly maintained fruit firmness, markedly reduced disease incidence and enhanced the activity of defence‐related enzymes (e.g., β‐1,3‐glucanase, phenylalanine ammonia lyase, peroxidase and polyphenol oxidase) and expression of pathogenesis‐related (PR) genes (e.g., PR‐1, PR‐3 and PR‐4). Therefore, the gamma irradiation‐induced resistance against P. expansum involves both metabolic changes and the induction of expression of defence‐related genes. In addition, scanning electron microscopic analysis revealed that gamma irradiation significantly inhibits the growth of P. expansum. These results suggest that exposure of mature harvested pear fruits to artificial gamma irradiation confers fungal disease resistance; therefore, gamma irradiation represents an important strategy for controlling postharvest diseases in pear fruit.