Cytogenetic analysis of quinoa chromosomes using nanoscale imaging and spectroscopy techniques

Chromosome-level genome assemblies for two quinoa inbred lines from northern and southern highlands of Altiplano where quinoa originated

Quinoa is emerging as a key seed crop for global food security due to its ability to grow in marginal environments and its excellent nutritional properties. Because quinoa is partially allogamous, we have developed quinoa inbred lines necessary for molecular genetic analysis. Our comprehensive genomic analysis showed that the quinoa inbred lines fall into three genetic subpopulations: northern highland, southern highland, and lowland. Lowland and highland quinoa are the same species, but have very different genotypes and phenotypes. Lowland quinoa has relatively small grains and a darker grain color, and is widely tested and grown around the world. In contrast, the white, large-grained highland quinoa is grown in the Andean highlands, including the region where quinoa originated, and is exported worldwide as high-quality quinoa. Recently, we have shown that viral vectors can be used to regulate endogenous genes in quinoa, paving the way for functional genomics to reveal the diversity of quinoa. However, although a high-quality assembly has recently been reported for a lowland quinoa line, genomic resources of the quality required for functional genomics are not available for highland quinoa lines. Here we present high-quality chromosome-level genome assemblies for two highland inbred quinoa lines, J075 representing the northern highland line and J100 representing the southern highland line, using PacBio HiFi sequencing and dpMIG-seq. In addition, we demonstrate the importance of verifying and correcting reference-based scaffold assembly with other approaches such as linkage maps. The assembled genome sizes of J075 and J100 are 1.29 and 1.32 Gb, with contigs N50 of 66.3 and 12.6 Mb, and scaffold N50 of 71.2 and 70.6 Mb, respectively, comprising 18 pseudochromosomes. The repetitive sequences of J075 and J100 represent 72.6% and 71.5% of the genome, the majority of which are long terminal repeats, representing 44.0% and 42.7% of the genome, respectively. The de novo assembled genomes of J075 and J100 were predicted to contain 65,303 and 64,945 protein-coding genes, respectively. The high quality genomes of these highland quinoa lines will facilitate quinoa functional genomics research on quinoa and contribute to the identification of key genes involved in environmental adaptation and quinoa domestication.

Genomic basis of seed colour in quinoa inferred from variant patterns using extreme gradient boosting

Quinoa is an agriculturally important crop species originally domesticated in the Andes of central South America. One of its most important phenotypic traits is seed colour. Seed colour variation is determined by contrasting abundance of betalains, a class of strong antioxidant and free radicals scavenging colour pigments only found in plants of the order Caryophyllales. However, the genetic basis for these pigments in seeds remains to be identified. Here we demonstrate the application of machine learning (extreme gradient boosting) to identify genetic variants predictive of seed colour. We show that extreme gradient boosting outperforms the classical genome-wide association approach. We provide re-sequencing and phenotypic data for 156 South American quinoa accessions and identify candidate genes potentially controlling betalain content in quinoa seeds. Genes identified include novel cytochrome P450 genes and known members of the betalain synthesis pathway, as well as genes annotated as being involved in seed development. Our work showcases the power of modern machine learning methods to extract biologically meaningful information from large sequencing data sets.

The genotype-dependent phenotypic landscape of quinoa in salt tolerance and key growth traits

Cultivation of quinoa (Chenopodium quinoa), an annual pseudocereal crop that originated in the Andes, is spreading globally. Because quinoa is highly nutritious and resistant to multiple abiotic stresses, it is emerging as a valuable crop to provide food and nutrition security worldwide. However, molecular analyses have been hindered by the genetic heterogeneity resulting from partial outcrossing. In this study, we generated 136 inbred quinoa lines as a basis for the molecular identification and characterization of gene functions in quinoa through genotyping and phenotyping. Following genotyping-by-sequencing analysis of the inbred lines, we selected 5,753 single-nucleotide polymorphisms (SNPs) in the quinoa genome. Based on these SNPs, we show that our quinoa inbred lines fall into three genetic sub-populations. Moreover, we measured phenotypes, such as salt tolerance and key growth traits in the inbred quinoa lines and generated a heatmap that provides a succinct overview of the genotype–phenotype relationship between inbred quinoa lines. We also demonstrate that, in contrast to northern highland lines, most lowland and southern highland lines can germinate even under high salinity conditions. These findings provide a basis for the molecular elucidation and genetic improvement of quinoa and improve our understanding of the evolutionary process underlying quinoa domestication.

Virus-Mediated Transient Expression Techniques Enable Functional Genomics Studies and Modulations of Betalain Biosynthesis and Plant Height in Quinoa

Quinoa (Chenopodium quinoa), native to the Andean region of South America, has been recognized as a potentially important crop in terms of global food and nutrition security since it can thrive in harsh environments and has an excellent nutritional profile. Even though challenges of analyzing the complex and heterogeneous allotetraploid genome of quinoa have recently been overcome, with the whole genome-sequencing of quinoa and the creation of genotyped inbred lines, the lack of technology to analyze gene function in planta is a major limiting factor in quinoa research. Here, we demonstrate that two virus-mediated transient expression techniques, virus-induced gene silencing (VIGS) and virus-mediated overexpression (VOX), can be used in quinoa. We show that apple latent spherical virus (ALSV) can induce gene silencing of quinoa phytoene desaturase (CqPDS1) in a broad range of quinoa inbred lines derived from the northern and southern highland and lowland sub-populations. In addition, we show that ALSV can be used as a VOX vector in roots. Our data also indicate that silencing a quinoa 3,4-dihydroxyphenylalanine 4,5-dioxygenase gene (CqDODA1) or a cytochrome P450 enzyme gene (CqCYP76AD1) inhibits betalain production and that knockdown of a reduced-height gene homolog (CqRHT1) causes an overgrowth phenotype in quinoa. Moreover, we show that ALSV can be transmitted to the progeny of quinoa plants. Thus, our findings enable functional genomics in quinoa, ushering in a new era of quinoa research.

Quantitative Distribution of DNA, RNA, Histone and Proteins Other than Histone in Mammalian Cells, Nuclei and a Chromosome at High Resolution Observed by Scanning Transmission Soft X-Ray Microscopy (STXM)

Soft X-ray microscopy was applied to study the quantitative distribution of DNA, RNA, histone, and proteins other than histone (represented by BSA) in mammalian cells, apoptotic nuclei, and a chromosome at spatial resolutions of 100 to 400 nm. The relative distribution of closely related molecules, such as DNA and RNA, was discriminated by the singular value decomposition (SVD) method using aXis2000 software. Quantities of nucleic acids and proteins were evaluated using characteristic absorption properties due to the 1s–π * transition of N=C in nucleic acids and amide in proteins, respectively, in the absorption spectra at the nitrogen K absorption edge. The results showed that DNA and histone were located in the nucleus. By contrast, RNA was clearly discriminated and found mainly in the cytoplasm. Interestingly, in a chromosome image, DNA and histone were found in the center, surrounded by RNA and proteins other than histone. The amount of DNA in the chromosome was estimated to be 0.73 pg, and the content of RNA, histone, and proteins other than histone, relative to DNA, was 0.48, 0.28, and 4.04, respectively. The method we present in this study could be a powerful approach for the quantitative molecular mapping of biological samples at high resolution.

Discrimination of DNA and RNA distribution in a mammalian cell by scanning transmission soft X-ray microscopy

Soft X-ray spectromicroscopy was applied to study the distribution of DNA and RNA in a mammalian cell at the spatial resolution of 400 nm. The relative distribution of DNA and RNA was examined by the SVD (singular value decomposition) method in aXis2000 program using combined full spectra of DNA and RNA at the absorption edge regions of carbon, nitrogen and oxygen. The absorption of nucleic acid was evaluated using 1s-π* transitions in the NEXAFS spectra at the nitrogen K absorption edge and distributed to DNA and RNA according to the relative level obtained above. The present results revealed the usefulness of the SVD method to discriminate closely related molecules such as DNA and RNA.

Development of novel InDel markers and genetic diversity in Chenopodium quinoa through whole-genome re-sequencing

Background

Quinoa (Chenopodium quinoa Willd.) is a balanced nutritional crop, but its breeding improvement has been limited by the lack of information on its genetics and genomics. Therefore, it is necessary to obtain knowledge on genomic variation, population structure, and genetic diversity and to develop novel Insertion/Deletion (InDel) markers for quinoa by whole-genome re-sequencing.

Results

We re-sequenced 11 quinoa accessions and obtained a coverage depth between approximately 7× to 23× the quinoa genome. Based on the 1453-megabase (Mb) assembly from the reference accession Riobamba, 8,441,022 filtered bi-allelic single nucleotide polymorphisms (SNPs) and 842,783 filtered InDels were identified, with an estimated SNP and InDel density of 5.81 and 0.58 per kilobase (kb). From the genomic InDel variations, 85 dimorphic InDel markers were newly developed and validated. Together with the 62 simple sequence repeat (SSR) markers reported, a total of 147 markers were used for genotyping the 129 quinoa accessions. Molecular grouping analysis showed classification into two major groups, the Andean highland (composed of the northern and southern highland subgroups) and Chilean coastal, based on combined STRUCTURE, phylogenetic tree and PCA (Principle Component Analysis) analyses. Further analysis of the genetic diversity exhibited a decreasing tendency from the Chilean coast group to the Andean highland group, and the gene flow between subgroups was more frequent than that between the two subgroups and the Chilean coastal group. The majority of the variations (approximately 70%) were found through an analysis of molecular variation (AMOVA) due to the diversity between the groups. This was congruent with the observation of a highly significant F_ST value (0.705) between the groups, demonstrating significant genetic differentiation between the Andean highland type of quinoa and the Chilean coastal type. Moreover, a core set of 16 quinoa germplasms that capture all 362 alleles was selected using a simulated annealing method.

Conclusions

The large number of SNPs and InDels identified in this study demonstrated that the quinoa genome is enriched with genomic variations. Genetic population structure, genetic core germplasms and dimorphic InDel markers are useful resources for genetic analysis and quinoa breeding.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4093-8) contains supplementary material, which is available to authorized users.

Draft genome sequence of an inbred line of Chenopodium quinoa, an allotetraploid crop with great environmental adaptability and outstanding nutritional properties

Chenopodium quinoa Willd. (quinoa) originated from the Andean region of South America, and is a pseudocereal crop of the Amaranthaceae family. Quinoa is emerging as an important crop with the potential to contribute to food security worldwide and is considered to be an optimal food source for astronauts, due to its outstanding nutritional profile and ability to tolerate stressful environments. Furthermore, plant pathologists use quinoa as a representative diagnostic host to identify virus species. However, molecular analysis of quinoa is limited by its genetic heterogeneity due to outcrossing and its genome complexity derived from allotetraploidy. To overcome these obstacles, we established the inbred and standard quinoa accession Kd that enables rigorous molecular analysis, and presented the draft genome sequence of Kd, using an optimized combination of high-throughput next generation sequencing on the Illumina Hiseq 2500 and PacBio RS II sequencers. The de novo genome assembly contained 25 k scaffolds consisting of 1 Gbp with N50 length of 86 kbp. Based on these data, we constructed the free-access Quinoa Genome DataBase (QGDB). Thus, these findings provide insights into the mechanisms underlying agronomically important traits of quinoa and the effect of allotetraploidy on genome evolution.

Discrimination Between Cervical Cancer Cells and Normal Cervical Cells Based on Longitudinal Elasticity Using Atomic Force Microscopy

The mechanical properties of cells are considered promising biomarkers for the early diagnosis of cancer. Recently, atomic force microscopy (AFM)-based nanoindentation technology has been utilized for the examination of cell cortex mechanics in order to distinguish malignant cells from normal cells. However, few attempts to evaluate the biomechanical properties of cells have focused on the quantification of the non-homogeneous longitudinal elasticity of cellular structures. In the present study, we applied a variation of the method of Carl and Schillers to investigate the differences between longitudinal elasticity of human cervical squamous carcinoma cells (CaSki) and normal cervical epithelial cells (CRL2614) using AFM. The results reveal a three-layer heterogeneous structure in the probing volume of both cell types studied. CaSki cells exhibited a lower whole-cell stiffness and a softer nuclei zone compared to the normal counterpart cells. Moreover, a better differentiated cytoskeleton was found in the inner cytoplasm/nuclei zone of the normal CRL2614 cells, whereas a deeper cytoskeletal distribution was observed in the probing volume of the cancerous counterparts. The sensitive cortical panel of CaSki cells, with a modulus of 0.35~0.47 kPa, was located at 237~225 nm; in normal cells, the elasticity was 1.20~1.32 kPa at 113~128 nm. The present improved method may be validated using the conventional Hertz-Sneddon method, which is widely reported in the literature. In conclusion, our results enable the quantification of the heterogeneous longitudinal elasticity of cancer cells, in particular the correlation with the corresponding depth. Preliminary results indicate that our method may potentially be applied to improve the detection of cancerous cells and provide insights into the pathophysiology of the disease.

Introduction of soft X-ray spectromicroscopy as an advanced technique for plant biopolymers research

Soft X-ray absorption spectroscopy coupled with nano-scale microscopy has been widely used in material science, environmental science, and physical sciences. In this work, the advantages of soft X-ray absorption spectromicroscopy for plant biopolymer research were demonstrated by determining the chemical sensitivity of the technique to identify common plant biopolymers and to map the distributions of biopolymers in plant samples. The chemical sensitivity of soft X-ray spectroscopy to study biopolymers was determined by recording the spectra of common plant biopolymers using soft X-ray and Fourier Transform mid Infrared (FT-IR) spectroscopy techniques. The soft X-ray spectra of lignin, cellulose, and polygalacturonic acid have distinct spectral features. However, there were no distinct differences between cellulose and hemicellulose spectra. Mid infrared spectra of all biopolymers were unique and there were differences between the spectra of water soluble and insoluble xylans. The advantage of nano-scale spatial resolution exploited using soft X-ray spectromicroscopy for plant biopolymer research was demonstrated by mapping plant cell wall biopolymers in a lentil stem section and compared with the FT-IR spectromicroscopy data from the same sample. The soft X-ray spectromicroscopy enables mapping of biopolymers at the sub-cellular (~30 nm) resolution whereas, the limited spatial resolution in the micron scale range in the FT-IR spectromicroscopy made it difficult to identify the localized distribution of biopolymers. The advantages and limitations of soft X-ray and FT-IR spectromicroscopy techniques for biopolymer research are also discussed.

Overview of nanoscale NEXAFS performed with soft X-ray microscopes

Today, in material science nanoscale structures are becoming more and more important. Not only for the further miniaturization of semiconductor devices like carbon nanotube based transistors, but also for newly developed efficient energy storage devices, gas sensors or catalytic systems nanoscale and functionalized materials have to be analysed. Therefore, analytical tools like near-edge X-ray absorption fine structure (NEXAFS) spectroscopy has to be applied on single nanostructures. Scanning transmission X-ray microscopes (STXM) as well as full-field transmission X-ray microscopes (TXM) allow the required spatial resolution to study individual nanostructures. In the soft X-ray energy range only STXM was used so far for NEXAFS studies. Due to its unique setup, the TXM operated by the Helmholtz-Zentrum Berlin (HZB) at the electron storage ring BESSY II is the first one in the soft X-ray range which can be used for NEXAFS spectroscopy studies which will be shown in this review. Here we will give an overview of the different microscopes used for NEXAFS studies and describe their advantages and disadvantages for different samples.