Genomic Survey Sequencing and Estimation of Genome Size of Ammopiptanthus mongolicus
WANG Xue,ZHOU Jia-yi,SUN Hui-gai,YU Rui-min,GAO Fei,ZHOU Yi-jun
DOI: https://doi.org/10.13430/j.cnki.jpgr.2018.01.016
2018-01-01
Abstract:Ammopiptanthus mongolicus (Maxim.ex Kom.)Cheng f.,one of the tertiary relict plants,is a rare evergreen broad-leaved shrub distributed in desert region of central Asia.A.mongolicus (Maxim.ex Kom.) Cheng f.grows in environment with extremely high level of drought and freezing stresses,and the plant species can be used as an ideal model for stress tolerance research in plant.Previous studies revealed that A.mongolicus (Maxim.ex Kom.) Cheng f.showed high levels of tolerance to various abiotic stresses.Several genes presumed to play important roles in stress tolerance of A.mongolicus (Maxim.ex Kom.)Cheng f.were cloned and characterized,including AnGobsl,AnPLD,AnMCSD,AnBADH,and AnAFP.However,research in A.mongolicus (Maxim.ex Kom.)Cheng f.at the molecular level was hindered by the lack of genome information.There were several genome size prediction methods,including Feulgen densitometry,field gel electrophoresis,quantitative reverse transcription PCR(qRT-PCR),flow cytometry analysis,and K-mer analysis.Flow cytometry was the standard method for estimating genome size.K-mer analysis not only predicted genome size accurately,but also provideed more information about genome,such as heterozygosity and GC percentage.In the present study,flow cytometry analysis was performed to predict the genome size of A.mongolicus (Maxim.ex Kom.) Cheng f..Genomic survey sequencing of A.mongolicus (Maxim.ex Kom.) Cheng f.was then conducted,and a total of 65 Gb sequencing reads was obtained.K-mer analysis using these sequencing reads was conducted to predict the genome size,heterozygosity and GC content of A.mongolicus (Maxim.ex Kom.) Cheng f..Based on the results of flow cytometry and K-mer analyses,the genome size of A.mongolicus (Maxim.ex Kom.)Cheng f.was estimated to be 770-787 Mb.The genome size of A.mongolicus (Maxim.ex Kom.)Cheng f.was close to that of Cicer arietinum L.,and was smaller than that of Glycine max L..The total reads length of contigs was 0.538 Gb,with N50 of 684 bp,and after further assembly,the N50 of scaffolds was 12.09 Kb with a total reads length of 0.602 Gb.The predicted heterozygosity was 0.0844% and the GC content was 36.51%.The low heterozygosity of A.mongolicus (Maxim.ex Kom.) Cheng f.will facilitate the whole genome sequencing of A.mongolicus(Maxim.ex Kom.)Cheng f.,and subsequent gene annotation and comparative genomics study.Simple sequence repeat(SSR) molecular markers were predicted using the assembled genome sequences and 151858 SSRs were obtained.Of all SSR categories,the dinucleotide repeat unit was the largest category,with a percentage of 56.39%,and the AT/TA components was the dominant dinucleotide repeat unit.Our study reported the genome size prediction of A.mongolicus (Maxim.ex Kom.) Cheng f.for the first time and provided a large number of genomic sequences for further research in A.mongolicus (Maxim.ex Kom.) Cheng f..The large number of SSR molecular markers identified in the present study will promote the study of genome mapping,evolutionary biology,and population genetics in A.mongolicus (Maxim.ex Kom.) Cheng f..