Genome-wide association mapping and candidate gene analysis for seed shape in soybean (Glycine max)
Xue Zhao A , Wenjing Li A , Xiaoyue Zhao A , Jinyang Wang A , Zhiyang Liu B , Yingpeng Han A C and Wenbin Li A C
+ Author Affiliations
- Author Affiliations
A Key Laboratory of Soybean Biology in Chinese Ministry of Education (Northeastern Key Laboratory of Soybean Biology and Genetics & Breeding in Chinese Ministry of Agriculture), Northeast Agricultural University, Harbin, China, 150030.
B School of the Grain and Economic Crops, Harbin Academy of Agricultural Science, Harbin, China, 150028.
C Corresponding authors. Email: hyp234286@aliyun.com; wenbinli@neau.edu.cn
Crop and Pasture Science 70(8) 684-693 https://doi.org/10.1071/CP19028
Submitted: 19 January 2019 Accepted: 2 July 2019 Published: 27 August 2019
Abstract
Seed shape (SS) of soybean (Glycine max (L.) Merr.) is an important morphological trait that significantly affects the quality of marketable seed. Study of the genetic architecture of SS is important and basic to soybean molecular breeding. In the present study, a natural soybean population of 202 diverse accessions mainly from China was used to analyse the genetic basis of SS via genome-wide association analysis (GWAS), which was based on single-nucleotide polymorphisms (SNP) generated by specific-locus amplified fragment sequencing method. In total, 27 335 SNPs were finally identified with minor allele frequencies >5%. By using GWAS, 14 quantitative trait nucleotides (QTNs) were identified to be associated with seed length, 13 with seed width and 21 with seed thickness in four tested environments. Among these QTNs, 21 QTNs overlapped or were located in the linked genomic regions of the reported quantitative trait loci related to SS or seed weight; and the other 27 QTNs were novel loci for SS. Ten QTNs showed environmental stability and were detected under at least two environments. In total, 83 genes were predicted in the 200-kbp flanking region of six stable QTNs that could be detected under >three environments. Gene-based association analysis was performed by using 38 accessions of diverse SS; 778 SNPs were found in the 83 genes based on 38 accessions, and 270 SNPs from 41 genes were found significantly associated with SS. Twenty-eight genes were environmentally stable and/or pleiotropic in controlling two or more SS-related traits at the same time. The identified loci along with the candidate genes could be of great value for studying the molecular mechanisms underlying SS and improving the potential seed yield of soybean.
Additional keywords: candidate genes, molecular assisted breeding, SLAF-seq.
References
Bradbury PJ, Zhang Z, Kroon DE, Casstevens TM, Ramdoss Y, Buckler ES (2007) TASSEL: software for association mapping of complex traits in diverse samples. Bioinformatics 23, 2633–2635.| TASSEL: software for association mapping of complex traits in diverse samples.Crossref | GoogleScholarGoogle Scholar | 17586829PubMed |
Cai H, Morishima H (2002) QTL clusters reflect character associations in wild and cultivated rice. Theoretical and Applied Genetics 104, 1217–1228.
| QTL clusters reflect character associations in wild and cultivated rice.Crossref | GoogleScholarGoogle Scholar |
Cober ER, Voldeng HD, Fregeau-Reid JA (1997) Heritability of seed shape and seed size in soybean. Crop Science 37, 1767–1769.
| Heritability of seed shape and seed size in soybean.Crossref | GoogleScholarGoogle Scholar |
Du W, Wang M, Fu S (2009) Mapping QTLs for seed yield and drought susceptibility index in soybean (Glycine max L.) across different environments. Journal of Genetics and Genomics 36, 721–731.
| Mapping QTLs for seed yield and drought susceptibility index in soybean (Glycine max L.) across different environments.Crossref | GoogleScholarGoogle Scholar | 20129399PubMed |
Fan C, Xiong Y, Mao H, Lu T, Han B, Xu C, Li X, Zhang Q (2006) GS3, a major QTL for grain length and weight and minor QTL for grain width and thickness in rice, encodes a putative transmembrane protein. Theoretical and Applied Genetics 112, 1164–1171.
| GS3, a major QTL for grain length and weight and minor QTL for grain width and thickness in rice, encodes a putative transmembrane protein.Crossref | GoogleScholarGoogle Scholar | 16453132PubMed |
Fatokun CA, Menancio-Hautea DI, Danesh D, Young ND (1992) Evidence for orthologous seed weight genes in cowpea and mung bean based on RFLP mapping. Genetics 132, 841–846.
Han Y, Li D, Zhu D, Li H, Li X, Teng W, Li W (2012) QTL analysis of soybean seed weight across multi-genetic backgrounds and environments. Theoretical and Applied Genetics 125, 671–683.
| QTL analysis of soybean seed weight across multi-genetic backgrounds and environments.Crossref | GoogleScholarGoogle Scholar | 22481120PubMed |
Hayes AJ, Ma GR, Buss GR, Saghai Maroof MA (2000) Molecular marker mapping of Rsv4, a gene conferring resistance to all known strains of soybean mosaic virus. Crop Science 40, 1434–1437.
| Molecular marker mapping of Rsv4, a gene conferring resistance to all known strains of soybean mosaic virus.Crossref | GoogleScholarGoogle Scholar |
Hoeck JA, Fehr WR, Shoemaker RC, Welke GA, Johnson SL, Cianzio SR (2003) Molecular marker analysis of soybean seed size. Crop Science 43, 68–74.
| Molecular marker analysis of soybean seed size.Crossref | GoogleScholarGoogle Scholar |
Hu ZB, Zhang HR, Kan GZ, Ma DY, Zhang D, Shi GX, Hong DL, Zhang GZ, Yu DY (2013) Determination of the genetic architecture of seed size and shape via linkage and association analysis in soybean (Glycine max L. Merr.). Genetica 141, 247–254.
| Determination of the genetic architecture of seed size and shape via linkage and association analysis in soybean (Glycine max L. Merr.).Crossref | GoogleScholarGoogle Scholar |
Jofuku KD, Omidyar PK, Gee Z, Okamuro JK (2005) Control of seed mass and seed yield by the floral homeotic gene APETALA2. Proceedings of the National Academy of Sciences of the United States of America 102, 3117–3122.
| Control of seed mass and seed yield by the floral homeotic gene APETALA2.Crossref | GoogleScholarGoogle Scholar | 15708974PubMed |
Jun TH, Freewalt K, Michel AP, Mian R (2014) Identification of novel QTL for leaf traits in soybean. Plant Breeding 133, 61–66.
| Identification of novel QTL for leaf traits in soybean.Crossref | GoogleScholarGoogle Scholar |
Kato S, Sayama T, Fujii K, Yumoto S, Kono Y, Hwang TY, Kikuchi A, Takada Y, Yu T, Tatsuhiko S, Masao I (2014) A major and stable QTL associated with seed weight in soybean across multiple environments and genetic backgrounds. Theoretical & Applied Genetics 127, 1365–1374.
| A major and stable QTL associated with seed weight in soybean across multiple environments and genetic backgrounds.Crossref | GoogleScholarGoogle Scholar |
Kuroda Y, Kaga A, Tomooka N, Yano H, Takada Y, Kato S, Vaughan D (2013) QTL affecting fitness of hybrids between wild and cultivated soybeans in experimental fields. Ecology and Evolution 3, 2150–2168.
| QTL affecting fitness of hybrids between wild and cultivated soybeans in experimental fields.Crossref | GoogleScholarGoogle Scholar | 23919159PubMed |
Li H, Durbin R (2009) Fast and accurate short read alignment with Burrows–Wheeler transform. Bioinformatics 25, 1754–1760.
Li RQ, Yu C, Li YR, Lam TW, Yiu SM, Kristiansen K, Wang J (2009a) SOAP2: an improved ultrafast tool for short read alignment. Bioinformatics 25, 1966–1967.
| SOAP2: an improved ultrafast tool for short read alignment.Crossref | GoogleScholarGoogle Scholar |
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R (2009b) The sequence alignment/map format and SAMtools. Bioinformatics 25, 2078–2079.
| The sequence alignment/map format and SAMtools.Crossref | GoogleScholarGoogle Scholar | 19505943PubMed |
Li X, Yan W, Agrama H, Jia L, Shen X, Jackson A, Moldenhauer K, Yeater K, McClung A, Wu D (2011a) Mapping QTLs for improving grain yield using the USDA rice mini-core collection. Planta 234, 347–361.
| Mapping QTLs for improving grain yield using the USDA rice mini-core collection.Crossref | GoogleScholarGoogle Scholar | 21479810PubMed |
Li Y, Fan C, Xiong Y, Jiang Y, Luo L, Sun L, Shao D, Xu C, Li X, Xiao J, He Y, Zhang Q (2011b) Natural variation in GS5 plays an important role in regulating grain size and yield in rice. Nature Genetics 43, 1266–1269.
| Natural variation in GS5 plays an important role in regulating grain size and yield in rice.Crossref | GoogleScholarGoogle Scholar | 22019783PubMed |
Liang Q, Cheng X, Mei M, Yan X, Liao H (2010) QTL analysis of root traits as related to phosphorus efficiency in soybean. Annals of Botany 106, 223
| QTL analysis of root traits as related to phosphorus efficiency in soybean.Crossref | GoogleScholarGoogle Scholar | 20472699PubMed |
Lipka AE, Tian F, Wang Q, Peiffer J, Li M, Bradbury P, Gore M, Buckler E, Zhang Z (2012) GAPIT: genome association and prediction integrated tool. Bioinformatics 28, 2397–2399.
| GAPIT: genome association and prediction integrated tool.Crossref | GoogleScholarGoogle Scholar | 22796960PubMed |
Liu J, Van Eck J, Cong B, Tanksley SD (2002) A new class of regulatory genes underlying the cause of pear-shaped tomato fruit. Proceedings of the National Academy of Sciences of the United States of America 99, 13302–13306.
| A new class of regulatory genes underlying the cause of pear-shaped tomato fruit.Crossref | GoogleScholarGoogle Scholar | 12242331PubMed |
Lu X, Xiong Q, Cheng T, Li QT, Liu XL, Bi YD, Li W, Zhang WK, Ma B, Lai YC (2017) A PP2C-1 allele underlying a quantitative trait locus enhances soybean 100-seed weight. Molecular Plant 10, 670–684.
| A PP2C-1 allele underlying a quantitative trait locus enhances soybean 100-seed weight.Crossref | GoogleScholarGoogle Scholar | 28363587PubMed |
Mansur LM, Orf JH, Chase K, Jarvik T, Cregan PB, Lark KG (1996) Genetic mapping of agronomic traits using recombinant inbred lines of soybean. Crop Science 36, 1327–1336.
| Genetic mapping of agronomic traits using recombinant inbred lines of soybean.Crossref | GoogleScholarGoogle Scholar |
Mao H, Sun S, Yao J, Wang C, Yu S, Xu C, Li X, Zhang Q (2010) Linking differential domain functions of the GS3 protein to natural variation of grain size in rice. Proceedings of the National Academy of Sciences of the United States of America 107, 19579–19584.
| Linking differential domain functions of the GS3 protein to natural variation of grain size in rice.Crossref | GoogleScholarGoogle Scholar | 20974950PubMed |
Mian MAR, Bailey MA, Tamulonis JP, Shipe ER, Carter TE, Parrott WA, Ashley DA, Hussey RS, Boerma HR (1996) Molecular markers associated with seed weight in two soybean populations. Theoretical and Applied Genetics 93, 1011–1016.
| Molecular markers associated with seed weight in two soybean populations.Crossref | GoogleScholarGoogle Scholar |
Nelson RL, Wang P (1989) Variation and evaluation of seed shape in soybean. Crop Science 29, 147–150.
| Variation and evaluation of seed shape in soybean.Crossref | GoogleScholarGoogle Scholar |
Niu Y, Xu Y, Liu XF, Yang SX, Wei SP, Xie FT, Zhang YM (2013) Association mapping for seed size and shape traits in soybean cultivars. Molecular Breeding 31, 785–794.
| Association mapping for seed size and shape traits in soybean cultivars.Crossref | GoogleScholarGoogle Scholar |
Ohto MA, Fischer RL, Goldberg RB, Nakamura K, Harada JJ (2005) Control of seed mass by APETALA2. Proceedings of the National Academy of Sciences of the United States of America 102, 3123–3128.
| Control of seed mass by APETALA2.Crossref | GoogleScholarGoogle Scholar | 15708976PubMed |
Orf JH, Chase K, Jarvik T, Mansur LM, Cregan PB, Adler FR, Lark KG (1999) Genetics of soybean agronomic traits: I. Comparison of three related recombinant inbred populations. Crop Science 39, 1642–1651.
| Genetics of soybean agronomic traits: I. Comparison of three related recombinant inbred populations.Crossref | GoogleScholarGoogle Scholar |
Palomeque L, Li-Jun L, Li W, Hedges B, Cober E, Rajcan I (2009) QTL in mega-environments: II. Agronomic trait QTL co-localized with seed yield QTL detected in a population derived from a cross of high-yielding adapted × high-yielding exotic soybean lines Theoretical and Applied Genetics 119, 429–436.
| QTL in mega-environments: II. Agronomic trait QTL co-localized with seed yield QTL detected in a population derived from a cross of high-yielding adapted × high-yielding exotic soybean linesCrossref | GoogleScholarGoogle Scholar | 19462149PubMed |
Pathan SM, Vuong T, Clark K, Lee J, Shannon J, Roberts C, Ellersieck M, Burton J, Burton J, Cregan P, Hyten D, Nguyen H, Sleper D (2013) Genetic mapping and confirmation of quantitative trait loci for seed protein and oil contents and seed weight in soybean. Crop Science 53, 765–774.
| Genetic mapping and confirmation of quantitative trait loci for seed protein and oil contents and seed weight in soybean.Crossref | GoogleScholarGoogle Scholar |
Qin H, Guo W, Zhang Y, Zhang T (2008) QTL mapping of yield and fiber traits based on a four-way cross population in Gossypium hirsutum L. Theoretical and Applied Genetics 117, 883–894.
| QTL mapping of yield and fiber traits based on a four-way cross population in Gossypium hirsutum L.Crossref | GoogleScholarGoogle Scholar | 18604518PubMed |
Quarrie SA, Quarrie SP, Radosevic R, Rancic D, Kaminska A, Barnes JD, Leverington M, Ceoloni C, Dodig D (2006) Dissecting a wheat QTL for yield present in a range of environments: from the QTL to candidate genes. Journal of Experimental Botany 57, 2627–2637.
| Dissecting a wheat QTL for yield present in a range of environments: from the QTL to candidate genes.Crossref | GoogleScholarGoogle Scholar | 16831847PubMed |
Rector BG, All JN, Parrott WA, Boerma HR (1999) Quantitative trait loci for antixenosis resistance to corn earworm in soybean. Crop Science 39, 531–538.
| Quantitative trait loci for antixenosis resistance to corn earworm in soybean.Crossref | GoogleScholarGoogle Scholar |
Salas P, Oyarzo-Llaipen JC, Wang D, Chase K, Mansur L (2006) Genetic mapping of seed shape in three populations of recombinant inbred lines of soybean (Glycine max L. Merr.). Theoretical and Applied Genetics 113, 1459–1466.
| Genetic mapping of seed shape in three populations of recombinant inbred lines of soybean (Glycine max L. Merr.).Crossref | GoogleScholarGoogle Scholar | 17036219PubMed |
Schmutz J, Cannon SB, Schlueter J, Ma J, Mitros T, Nelson W, Hyten DL, Song Q, Thelen JJ, Cheng J, Xu D, Helllsten U, May GD, Yu Y, Sakurai T, Umezawa T, Bhattacharyya M, Sandhu D, Valliyodan B, Lindquist E, Peto M, Grant D, Shu S, Goodstein D, Barry K, Futrell-Griggs M, Abernathy B, Du J, Tian Z, Zhu L, Gill N, Joshi T, Libault M, Sethuraman A, Zhang XC, Shinozaki K, Nguyen HT, Wing RA, Cregan P, Specht J, Grimwood J, Rokhsar D, Stacey G, Shoemaker RC, Jackson SA (2010) Genome sequence of the palaeopolyploid soybean. Nature 463, 178–183.
| Genome sequence of the palaeopolyploid soybean.Crossref | GoogleScholarGoogle Scholar | 20075913PubMed |
Schruff MC, Spielman M, Tiwari S, Adams S, Fenby N, Scott RJ (2006) The AUXIN RESPONSE FACTOR 2 gene of Arabidopsis links auxin signalling, cell division, and the size of seeds and other organs. Development 133, 251–261.
| The AUXIN RESPONSE FACTOR 2 gene of Arabidopsis links auxin signalling, cell division, and the size of seeds and other organs.Crossref | GoogleScholarGoogle Scholar | 16339187PubMed |
Sebolt AM, Shoemaker RC, Diers BW (2000) Analysis of a quantitative trait locus allele from wild soybean that increases seed protein concentration in soybean. Crop Science 40, 1438–1444.
| Analysis of a quantitative trait locus allele from wild soybean that increases seed protein concentration in soybean.Crossref | GoogleScholarGoogle Scholar |
Shappley ZW, Jenkins JN, Zhu J, McCarty JC (1998) Quantitative trait loci associated with yield and fiber traits of upland cotton. Journal of Cotton Science 4, 153–163.
Song X, Huang W, Shi M, Zhu M, Lin H (2007) A QTL for rice grain width and weight encodes a previously unknown RING-type E3 ubiquitin ligase. Nature Genetics 39, 623–630.
| A QTL for rice grain width and weight encodes a previously unknown RING-type E3 ubiquitin ligase.Crossref | GoogleScholarGoogle Scholar | 17417637PubMed |
Stoll M, Kwitek-Black AE, Cowley AW, Harris EL, Harrap SB, Krieger JE, Printz MP, Provoost AP, Sassard J, Jacob HJ (2000) New target regions for human hypertension via comparative genomics. Genome Research 10, 473–482.
| New target regions for human hypertension via comparative genomics.Crossref | GoogleScholarGoogle Scholar | 10779487PubMed |
Sun X, Shantharaj D, Kang X, Ni M (2010) Transcriptional and hormonal signaling control of Arabidopsis seed development. Current Opinion in Plant Biology 13, 611–620.
| Transcriptional and hormonal signaling control of Arabidopsis seed development.Crossref | GoogleScholarGoogle Scholar | 20875768PubMed |
Sun X, Liu D, Zhang X, Li W, Liu H, Hong W, Jiang C, Guan N, Ma C, Zeng H, Xu C, Song J, Huang L, Zheng H (2013) SLAF-seq: an efficient method of large-scale de novo SNP discovery and genotyping using high-throughput sequencing. PLoS One 8, e58700
| SLAF-seq: an efficient method of large-scale de novo SNP discovery and genotyping using high-throughput sequencing.Crossref | GoogleScholarGoogle Scholar | 24391853PubMed |
Teng W, Han Y, Du Y, Sun D, Qiu L, Sun G, Li W (2009) QTL analyses of seed weight during the development of soybean (Glycine max L. Merr.). Heredity 102, 372
| QTL analyses of seed weight during the development of soybean (Glycine max L. Merr.).Crossref | GoogleScholarGoogle Scholar | 18971958PubMed |
Teng WL, Sui MN, Li W, Wu DP, Zhao X, Li HY, Han YP, Li WB (2018) Identification of quantitative trait loci underlying seed shape in soybean across multiple environments. The Journal of Agricultural Science 156, 3–12.
| Identification of quantitative trait loci underlying seed shape in soybean across multiple environments.Crossref | GoogleScholarGoogle Scholar |
van der Knaap E, Tanksley SD (2001) Identification and character1ization of a novel locus controlling early fruit development in tomato. Theoretical and Applied Genetics 103, 353–358.
| Identification and character1ization of a novel locus controlling early fruit development in tomato.Crossref | GoogleScholarGoogle Scholar |
Wan JM (2006) Perspectives of molecular design breeding in crops. Acta Agronomica Sinica 32, 455–462.
Wang D, Graef GL, Procopiuk AM, Diers BW (2004) Identification of putative QTL that underlie yield in interspecific soybean backcross populations. Theoretical and Applied Genetics 108, 458–467.
| Identification of putative QTL that underlie yield in interspecific soybean backcross populations.Crossref | GoogleScholarGoogle Scholar | 14504749PubMed |
Wang J, Wan X, Crossa J, Crouch J, Weng J, Zhai H, Wan J (2006) QTL mapping of grain length in rice (Oryza sativa L.) using chromosome segment substitution lines. Genetical Research 88, 93–104.
| QTL mapping of grain length in rice (Oryza sativa L.) using chromosome segment substitution lines.Crossref | GoogleScholarGoogle Scholar | 17125584PubMed |
Wang A, Garcia D, Zhang H, Feng K, Chaudhury A, Berger F, Peacock WJ, Dennis ES, Luo M (2010) The VQ motif protein IKU1 regulates endosperm growth and seed size in Arabidopsis. The Plant Journal 63, 670–679.
| The VQ motif protein IKU1 regulates endosperm growth and seed size in Arabidopsis.Crossref | GoogleScholarGoogle Scholar | 20545893PubMed |
Wang S, Wu K, Yuan Q, Liu X, Liu Z, Lin X, Zeng R, Zhu H, Dong G, Qian Q, Zhang G, Fu X (2012) Control of grain size, shape and quality by OsSPL16 in rice. Nature Genetics 44, 950–954.
| Control of grain size, shape and quality by OsSPL16 in rice.Crossref | GoogleScholarGoogle Scholar | 22729225PubMed |
Wilson DO, Jr (1995) Storage of orthodox seeds. In ‘Seed quality: basic mechanisms agricultural implications’. (Ed. AS Basra) pp. 173–208. (Food Products Press: New York)
Xie FT, Niu Y, Zhang J, Bu SH, Zhang HZ, Geng QC, Feng JY, Zhang YM (2014) Fine mapping of quantitative trait loci for seed size traits in soybean. Molecular Breeding 34, 2165–2178.
| Fine mapping of quantitative trait loci for seed size traits in soybean.Crossref | GoogleScholarGoogle Scholar |
Yang K, Moon JK, Jeong N, Chun H, Kang S, Back K, Jeong S (2011) Novel major quantitative trait loci regulating the content of isoflavone in soybean seeds. Genes & Genomics 33, 685–692.
| Novel major quantitative trait loci regulating the content of isoflavone in soybean seeds.Crossref | GoogleScholarGoogle Scholar |
Yang HY, Wang WB, He QY, Xiang SH, Tian D, Zhao TJ, Gai JY (2017) Chromosome segment detection for seed size and shape traits using an improved population of wild soybean chromosome segment substitution lines. Physiology and Molecular Biology of Plants 23, 877–889.
| Chromosome segment detection for seed size and shape traits using an improved population of wild soybean chromosome segment substitution lines.Crossref | GoogleScholarGoogle Scholar |
Yao D, Liu ZZ, Zhang J, Liu S, Qu J, Guan S, Pan L, Wang D, Liu J, Wang P (2015) Analysis of quantitative trait loci for main plant traits in soybean. Genetics & Molecular Research 14, 6101–6109.
| Analysis of quantitative trait loci for main plant traits in soybean.Crossref | GoogleScholarGoogle Scholar |
Zhou Y, Ni M (2010) SHORT HYPOCOTYL UNDER BLUE1 truncations and mutations alter its association with a signaling protein complex in Arabidopsis. The Plant Cell 22, 703–715.
| SHORT HYPOCOTYL UNDER BLUE1 truncations and mutations alter its association with a signaling protein complex in Arabidopsis.Crossref | GoogleScholarGoogle Scholar | 20354198PubMed |
Zhou Y, Zhang X, Kang X, Zhao X, Zhang X, Ni M (2009) SHORT HYPOCOTYL UNDER BLUE1 associates with MINISEED3 and HAIKU2 promoters in vivo to regulate Arabidopsis seed development. The Plant Cell 21, 106–117.
| SHORT HYPOCOTYL UNDER BLUE1 associates with MINISEED3 and HAIKU2 promoters in vivo to regulate Arabidopsis seed development.Crossref | GoogleScholarGoogle Scholar | 19141706PubMed |
Zhou Z, Jiang Y, Wang Z, Gou Z, Lyu J, Li W, Yu Y, Shu L, Zhao Y, Ma Y, Fang C, Shen Y, Liu T, Li C, Li Q, Wu M, Wang M, Wu Y, Dong Y, Wan W, Wang X, Ding Z, Gao Y, Xiang H, Zhu B, Lee SH, Wang W, Tian Z (2015) Resequencing 302 wild and cultivated accessions identifies genes related to domestication and improvement in soybean. Nature Biotechnology 33, 408–414.
| Resequencing 302 wild and cultivated accessions identifies genes related to domestication and improvement in soybean.Crossref | GoogleScholarGoogle Scholar | 25643055PubMed |
