Background: Wheat (Triticum ssp.) is an important food source for humans in many regions around the world. However, the ability to understand and modify gene function for crop improvement is hindered by the lack of available genomic resources. TILLING is a powerful reverse genetics approach that combines chemical mutagenesis with a high-throughput screen for mutations. Wheat is specially well-suited for TILLING due to the high mutation densities tolerated by polyploids, which allow for very efficient screens. Despite this, few TILLING populations are currently available. In addition, current TILLING screening protocols require high-throughput genotyping platforms, limiting their use.

Discovery of rare mutations in populations requires methods, such as TILLING (for Targeting Induced Local Lesions in Genomes), for processing and analyzing many individuals in parallel. Previous TILLING protocols employed enzymatic or physical discrimination of heteroduplexed from homoduplexed target DNA. Using mutant populations of rice (Oryza sativa) and wheat (Triticum durum), we developed a method based on Illumina sequencing of target genes amplified from multidimensionally pooled templates representing 768 individuals per experiment. Parallel processing of sequencing libraries was aided by unique tracer sequences and barcodes allowing flexibility in the number and pooling arrangement of targeted genes, species, and pooling scheme. Sequencing reads were processed and aligned to the reference to identify possible single-nucleotide changes, which were then evaluated for frequency, sequencing quality, intersection pattern in pools, and statistical relevance to produce a Bayesian score with an associated confidence threshold. Discovery was robust both in rice and wheat using either bidimensional or tridimensional pooling schemes. The method compared favorably with other molecular and computational approaches, providing high sensitivity and specificity.

BackgroundSalt stress, a major plant environmental stress, is a critical constraint for rice productivity. Dissecting the genetic loci controlling salt tolerance in rice for improving productivity, especially at the flowering stage, remains challenging. Here, we conducted a genome-wide association study (GWAS) of salt tolerance based on exome sequencing of the Thai rice accessions.ResultsPhotosynthetic parameters and cell membrane stability under salt stress at the flowering stage; and yield-related traits of 104 Thai rice (Oryza sativa L.) accessions belonging to the indica subspecies were evaluated. The rice accessions were subjected to exome sequencing, resulting in 112,565 single nucleotide polymorphisms (SNPs) called with a minor allele frequency of at least 5%. LD decay analysis of the panel indicates that the average LD for SNPs at 20 kb distance from each other was 0.34 (r2), which decayed to its half value (~ 0.17) at around 80 kb. By GWAS performed using mixed linear model, two hundred loci containing 448 SNPs on exons were identified based on the salt susceptibility index of the net photosynthetic rate at day 6 after salt stress; and the number of panicles, filled grains and unfilled grains per plant. One hundred and forty six genes, which accounted for 73% of the identified loci, co-localized with the previously reported salt quantitative trait loci (QTLs). The top four regions that contained a high number of significant SNPs were found on chromosome 8, 12, 1 and 2. While many are novel, their annotation is consistent with potential involvement in plant salt tolerance and in related agronomic traits. These significant SNPs greatly help narrow down the region within these QTLs where the likely underlying candidate genes can be identified.ConclusionsInsight into the contribution of potential genes controlling salt tolerance from this GWAS provides further understanding of salt tolerance mechanisms of rice at the flowering stage, which can help improve yield productivity under salinity via gene cloning and genomic selection.Electronic supplementary materialThe online version of this article (10.1186/s12864-018-5317-2) contains supplementary material, which is available to authorized users.

Targeting Induced Local Lesions in Genomes (TILLING) provides a nontransgenic method for reverse genetics that is widely applicable, even in species where other functional resources are missing or expensive to build. The efficiency of TILLING, however, is greatly facilitated by high mutation density. Species vary in the number of mutations induced by comparable mutagenic treatments, suggesting that genetic background may affect the response. Allopolyploid species have often yielded higher mutation density than diploids. To examine the effect of ploidy, we autotetraploidized the Arabidopsis (Arabidopsis thaliana) ecotype Columbia, whose diploid has been used for TILLING extensively, and mutagenized it with 50 mM ethylmethane sulfonate. While the same treatment sterilized diploid Columbia, the tetraploid M1 plants produced good seed. To determine the mutation density, we searched 528 individuals for induced mutations in 15 genes for which few or no knockout alleles were previously available. We constructed tridimensional pools from the genomic DNA of M2 plants, amplified target DNA, and subjected them to Illumina sequencing. The results were analyzed with an improved version of the mutation detection software CAMBa that accepts any pooling scheme. This small population provided a rich resource with approximately 25 mutations per queried 1.5-kb fragment, including on average four severe missense and 1.3 truncation mutations. The overall mutation density of 19.4 mutations Mb -1 is 4 times that achieved in the corresponding diploid accession, indicating that genomic redundancy engenders tolerance to high mutation density. Polyploidization of diploids will allow the production of small populations, such as less than 2,000, that provide allelic series from knockout to mild loss of function for virtually all genes.

Gene dosage variation and the associated changes in gene expression influence a wide variety of traits, ranging from cancer in humans to yield in plants. It is also expected to affect important traits of ecological and agronomic importance in forest trees, but this variation has not been systematically characterized or exploited. Here we performed a comprehensive scan of the Populus genome for dosage-sensitive loci affecting quantitative trait variation for spring and fall phenology and biomass production. The study population was a large collection of clonally propagated F1 hybrid lines of Populus that saturate the genome 10-fold with deletions and insertions (indels) of known sizes and positions. As a group, the phenotypic means of the indel lines consistently differed from control nonindel lines, with an overall negative effect of both insertions and deletions on all biomass-related traits but more diverse effects and an overall wider phenotypic distribution of the indel lines for the phenology-related traits. We also investigated the correlation between gene dosage at specific chromosomal locations and phenotype, to identify dosage quantitative trait loci (dQTL). Such dQTL were detected for most phenotypes examined, but stronger effect dQTL were identified for the phenology-related traits than for the biomass traits. Our genome-wide screen for dosage sensitivity in a higher eukaryote demonstrates the importance of global genomic balance and the impact of dosage on life history traits.

Gene copy number variation is frequent in plant genomes of various species, but the impact of such gene dosage variation on morphological traits is poorly understood. We used a large population of Populus carrying genomically characterized insertions and deletions across the genome to systematically assay the effect of gene dosage variation on a suite of leaf morphology traits. A systems genetics approach was used to integrate insertion and deletion locations, leaf morphology phenotypes, gene expression, and transcriptional network data, to provide an overview of how gene dosage influences morphology. Dosage-sensitive genomic regions were identified that influenced individual or pleiotropic morphological traits. We also identified cis-expression quantitative trait loci (QTL) within these dosage QTL regions, a subset of which modulated trans-expression QTL as well. Integration of data types within a gene co-expression framework identified co-expressed gene modules that are dosage sensitive, enriched for dosage expression QTL, and associated with morphological traits. Functional description of these modules linked dosage-sensitive morphological variation to specific cellular processes, as well as candidate regulatory genes. Together, these results show that gene dosage variation can influence morphological variation through complex changes in gene expression, and suggest that frequently occurring gene dosage variation has the potential to likewise influence quantitative traits in nature.

The operating thermal power transfer mechanisms in a thermal-wave resonant cavity were explored theoretically and experimentally. Both steady-state ac (thermal-wave) and dc temperature rise were considered, and conduction and radiation heat transfer modes were found to co-exist in the cavity. By introducing controlled variable offset dc resistive heating superposed on the fixed-amplitude thermal-wave oscillation, it was also found that the thermal-diffusivity values of the intracavity gas can vary sensitively as a function of the dc temperature rise within a thin boundary layer adjacent to the cavity thermal source (a metallic Cr–Ni alloy strip). This resulted in the measurement of the temperature dependence of the thermal diffusivity of air. Furthermore, the observed dominance of thermal-wave radiation power transfer in the phase channel of the thermal-wave signal at large cavity lengths allowed the measurement of the absolute infrared emissivity of the thin Cr–Ni strip source material: ε=0.091±0.004.

The chromosome axis plays a crucial role in meiotic recombination. Here, we study the function of ASY1, the Arabidopsis homolog of the yeast chromosome axis associated component Hop1. Specifically, we characterized cross-over (CO) distribution in female and male meiosis by deep sequencing of the progeny of an allelic series of asy1 mutants. Combining data from nearly 1000 individual plants, we find that reduced ASY1 activity leads to genomic instability and sometimes drastic genomic rearrangements. We further observed that COs are less frequent and appear in more distal chromosomal regions in plants with no or reduced ASY1 activity, consistent with previous analyses. However, our sequencing approach revealed that the reduction in CO number is not as dramatic as suggested by cytological analyses. Analysis of double mutants of asy1 with mutants with three other CO factors, MUS81, MSH4 and MSH5 as well as the determination of foci number of the CO regulator MLH1 demonstrates that the majority of the COs in asy1, similar to the situation in the wildtype, largely belong to the class I, which are subject to interference. However, these COs are redistributed in asy1 mutants and typically appear much closer than in the wildtype. Hence, ASY1 plays a key role in CO interference that spaces COs along a chromosome. Conversely, since a large proportion of chromosomes do not receive any CO, we conclude that CO assurance, the process that ensures the obligatory assignment of one CO per chromosome, is also affected in asy1 mutants.Significant statementThe regulation of the number and placement of cross-overs (COs) during meiosis is critical to ensure meiotic fidelity and promote new genetic combinations. Here, we investigated the function of one of the proteins of the chromosome axis, which plays a key role in CO formation: ASY1. Our results show that COs in asy1 mutants are positioned closer to each other than in the wildtype and that, despite a roughly similar number of COs, not every chromosome receives a CO. With this, our results shed light on the mechanisms regulating two important but still poorly understood aspects of meiosis: CO assurance, which safeguards at least one CO per chromosome pair, and CO interference, which prevents two COs from occurring close to each other.