Complete sequence of heterogenous-composition mitochondrial genome (Brassica napus) and its exogenous source
- Juan Wang†1,
- Jinjin Jiang†1,
- Xiaoming Li1,
- Aimin Li2,
- Yongtai Zhang2,
- Rongzhan Guan3 and
- Youping Wang1Email author
© Wang et al.; licensee BioMed Central Ltd. 2012
Received: 24 July 2012
Accepted: 27 November 2012
Published: 28 November 2012
Unlike maternal inheritance of mitochondria in sexual reproduction, somatic hybrids follow no obvious pattern. The introgressed segment orf138 from the mitochondrial genome of radish (Raphanus sativus) to its counterpart in rapeseed (Brassica napus) demonstrates that this inheritance mode derives from the cytoplasm of both parents. Sequencing of the complete mitochondrial genome of five species from Brassica family allowed the prediction of other extraneous sources of the cybrids from the radish parent, and the determination of their mitochondrial rearrangement.
We obtained the complete mitochondrial genome of Ogura-cms-cybrid (oguC) rapeseed. To date, this is the first time that a heterogeneously composed mitochondrial genome was sequenced. The 258,473 bp master circle constituted of 33 protein-coding genes, 3 rRNA sequences, and 23 tRNA sequences. This mitotype noticeably holds two copies of atp9 and is devoid of cox 2-2. Relative to nap mitochondrial genome, 40 point mutations were scattered in the 23 protein-coding genes. atp6 even has an abnormal start locus whereas tatC has an abnormal end locus. The rearrangement of the 22 syntenic regions that comprised 80.11% of the genome was influenced by short repeats. A pair of large repeats (9731 bp) was responsible for the multipartite structure. Nine unique regions were detected when compared with other published Brassica mitochondrial genome sequences. We also found six homologous chloroplast segments (Brassica napus).
The mitochondrial genome of oguC is quite divergent from nap and pol, which are more similar with each other. We analyzed the unique regions of every genome of the Brassica family, and found that very few segments were specific for these six mitotypes, especially cam, jun, and ole, which have no specific segments at all. Therefore, we conclude that the most specific regions of oguC possibly came from radish. Compared with the chloroplast genome, six identical regions were found in the seven mitochondrial genomes, which show that the Brassica family has a stable chloroplast-derived source.
The major function of the mitochondria (mt), as a semiautonomous organelle, in plant growth and development is to provide energy through oxidative phosphorylation . In different to the small mt genome of animals (~16 kb), plants have longer mtDNA ranging from 200kb to 2000kb [2, 3]. To date, several mt genomes from fertile and sterile plant species have been sequenced, including Arabidopsis thaliana, Oryza sativa[5–7], Beta vulgaris[8, 9], Zea mays[10, 11], Nicotiana tabacum, Triticum aestivum[13, 14], and five species from the Brassica genus, i.e., B. napus (pol, nap), B. rapa (cam), B. oleracea (ole), B. juncea (jun), and B. carinata (car) [15–17]. The sequencing results indicated that apart from ribosomal protein genes, protein-coding genes are also relatively conserved both in nucleotide sequence and in number. However, the non-coding sequences are quite inconsistent among species, and even within the same species. The presence of large and short repeats is responsible for the dynamic multipartite structures, reorganization, and recombination .
In higher plants, mitochondrial inheritance usually follows the maternal origin during sexual hybridization. However, much more complicated modes are detected in somatic hybridization, wherein mt genome inheritance is derived from either or both biparents . In the latter pattern, part of the mt genome, including cytoplasmic male sterility (CMS) genes, can be transferred from the donor parent to the receptor parent, and the introgressed segment experiences extensive rearrangement and recombination with the mtDNA of the receptor one. Orf138, originally identified in radish, was transferred successfully to various species, including Arabidopsis, B. napus, and B. oleracea by somatic hybridization [19–25].
CMS genes have a defect in the production of functional pollen. Generally, genes associated with CMS genes are located in the periphery of certain known mitochondrial genes and are cotranscribed with them . T-urf13 (orf115) was the first identified aberrant gene in the Texas (T)-cytoplasm of maize, which encodes a 13 kDa membrane-spanning polypeptide that depolarizes the mitochondria and leads to cell death [27–29]. In the BT (Boro II)-type CMS line of rice, orf79 was cotranscribed with the atp6 gene forming a 2.0 kb transcript . The expressed protein contains a predicted transmembrane domain .
In the Brassica genus, the complete mt genomes of five species are sequenced, coupled with the basic feature of published CMS genes, which allows the detection of the extraneous source from donor parent (radish) of somatic hybrids.
Genome size and nucleotide sequence in the genic region
SNP in protein-coding genes of mtDNA between oguC and nap
Reconstruction of the nap-CMS cybrid mitochondrial genome
Reorganization of the mitochondrial genome
Unique region of the genome
Unique region found in oguC
included in the large repeats
part of orf101-4
Homology with rapeseed chloroplast genome
Homologous segments to chloroplast of rapeseed found in oguC mitotype
ORFs and predicted CMS-related chimeric ORF in this genome
We detected 39 ORFs in this genome, with the shortest size equal to 303 bp, which summed to 7.41% of the mitotype. Of the 39 ORFs, 23 (similarity ≥ 99%) were shared in one or more Brassica genomes, which were remotely related to CMS. However, they are likely functional genes as these later-discovered genes ccm, orf25 (atp4), and orfB (atp8) [37, 38]. Of the remaining 16 unique ORFs in the oguC mitochondrial genome, 8 ORFs were totally not matched and 8 were partly identical to those present in the other six mitotypes. Five, including orf138, which is the oguC-related CMS gene, out of eight non-matched ORFs were completely or partially situated in the unique regions (Table 2). Three common ORFs and three unique ORFs were located in the cp-derived domains (Table 3). Among them, orf210 and orf344-1 were highly similar to 2 segments of the beta subunit of RNA polymerase, which were wholly situated in the chloroplast genome with a length of 1072 amino acids. When the intact nucleotide sequence of beta subunit of RNA polymerase was aligned with oguC mtDNA, H1 was found to be a truncated RNA polymerase beta subunit with 97% similarity. Some point mutations and indels resulted in the production of these two ORFs. Similar to orf344-1, a truncated ribulose-1, 5-bisphosphate carboxylase/oxygenase large subunit (rbcL) gene from cp genome evolved into orf313. Based on the three features of CMS related genes, namely, unique to the given mitotype, membrane-spanning domains and near to the functional genes [26, 39], orf138, which encodes a 19 kDa transmembrane protein that showed toxicity to bacterial growth, can be suppressed by the nuclear Rfo locus [40–43].
When nap mtDNA was taken as the control, 22 syntenic regions were detected in total. Estimating the number of recombination events was difficult because of the many syntenic regions. However, relative to two more similar mitotypes (pol and nap) that have 13 syntenic region with the same analysis criterion (length ≥ 1000 bp, similarity ≥ 95%) , it showed complex reconstruction. oguC mitotype must have undergone complicated changes and evolutionary events when the cytoplasm of two cells contacted each other.
The composition of mitochondrial genomes of 7 mitotypes from Brassica family*
Genome size (bp)
Unique region (%)
cp-derived sequences (%)
Large repeat (bp)
Short repeat (%)
Using one genome of the seven mitotypes as the control to find the unique regions for every mitotype, the percentage ranged from 0% in three mitotypes (cam, jun, and ole) to 8.60% in oguC. Both nap and car contained three shorter specific segments constituting 0.74% and 1.57% of these two genomes, respectively. A 620 bp unique segment located in pol contributed 0.28% to its total genome. In terms of the percentage, at least 7% of the exogenous sequences from radish mtDNA coexist with the oguC mtDNA (Table 4). When searched against the NCBI databases using those specific segments, similar alignments to that of oguC were obtained, some of which resembled those in Arabidopsis thaliana.
We also predicted the cp-derived sequence, which was intriguing because the seven mitotypes were blasted for the identical six segments with identities more than 95% (Table 4). However, because of the large copy of R1, five cp-derived segments had two copies in ole. From the cp-derived data, we found that Brassica species have stable sources of chloroplast sequence.
This study finished mtDNA sequencing of a Ogura-cms-cybrid (oguC), which derived from somatic fusion between Brassica napus and sterile radish. By contrast to one or more of six other Brassica lines, we reasonably speculated that tatC gene and 2 unique regions, U3 and U7, must be introgressed from radish. In addition, the rearrangement mediated by large and short repeats between these two parental mtDNAs extensively existed. With regard to the evolution of this integrated CMS mtDNA, more data need to be known.
Plant material and mitochondrial genome extraction
Seed of Brassica napus (ogu-CMS cybrid, oguC) was kindly provided by Norddeutsche Pflanzenzucht, Hans-Georg Lembke KG (Germany). The etiolated one-week-old Brassica napus seedlings were prepared; the mitochondria and mtDNA extraction were performed following previously published methods . To satisfy the requirements for 454 sequences, the minimum criterion for sample concentration was 50 ng/μl and total amount was equal to at least 20 μg.
A shotgun library that includes short and long paired end libraries were constructed simultaneously, which was followed by emulsion-based clonal amplification (emPCR) for DNA library bead enrichment. Finally, a genome sequencer FLX operation was conducted and the system output was derived. The contigs were joined by PCR sequencing. For the oguC genome, high quality read number, high quality bases, average read length, and sequencing depth were 8387, 3,913,351 bp, 469.3 bp, 15.2X, respectively. For SNP analysis and unique regions in oguC resequence was done.
The genes scattered in this genome were annotated using the Blast service of NCBI. tRNAscan  and ORF finder (http://www.ncbi.nlm.nih.gov/gorf/gorf.html) were used to identify the tRNA sequences and potential ORFs, respectively. The unique regions of seven genomes were dug out with MultiPipMaker . BlastN was used to discover large repeats. Short repeats were detected using commercial software developed by Shanghai Majorbio Bio-pharm Biotechnology Company (China). The accession numbers of the mtDNA are listed in Table 4.
This work was supported by the NSFC project (30971812, 31171581), a Project Funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions, Program of International S &T Cooperation of China (1021) and the project of Jiangsu Province (BRA2010141, CXLX11_0998).
- Nathalie P, Michael H, Luigi P, Ferdinando P: The growing family of mitochondrial carriers in Arabidopsis. Trends Plant Sci. 2004, 9: 138-146. 10.1016/j.tplants.2004.01.007.View ArticleGoogle Scholar
- Oda K, Yamato K, Ohta E, Nakamura Y, Takemura M, Nozato N, Akashi K, Kanegae T, Ogura Y, Kohchi T: Gene organization deduced from the complete sequence of livewort Marchantia polymorpha mitochondrial DNA primitive form of plant mitochondrial genome. Mol Biol. 1992, 223: 1-7. 10.1016/0022-2836(92)90708-R.View ArticleGoogle Scholar
- Gray MW, Lang BL, Cedergren R, Golding GB, Lemieux C, Sankoff D, Turmel M, Brossard N, Delage E, Littlejohn TG, et al: Genome structure and gene content in protist mitochondrial DNAs. Nucleic Acids Res. 1998, 26: 865-878. 10.1093/nar/26.4.865.PubMed CentralView ArticlePubMedGoogle Scholar
- Unseld M, Marienfeld JR, Brandt P, Brandt P, Brennicke A: The complete genome of Arabidopsis thaliana contains 57 genes in 366,924 nucleotides. Nat Genet. 1997, 15: 57-61.View ArticlePubMedGoogle Scholar
- Notsu Y, Masood S, Nishikawa T, Kubo N, Akiduki G, Nakazono M, Hirai A, Kadowaki K: The complete sequence of the rice (Oryza sativa L.) mitochondrial genome: frequent DNA sequence acquisition and loss during the evolution of flowering plants. Mol Genet Genomics. 2002, 268: 434-445. 10.1007/s00438-002-0767-1.View ArticlePubMedGoogle Scholar
- Tian X, Zheng J, Hu S, Yu J: The rice mitochondrial genomes and their variations. Plant Physiol. 2006, 140: 401-410. 10.1104/pp.105.070060.PubMed CentralView ArticlePubMedGoogle Scholar
- Fujii S, Kazama T, Yamada M, Toriyama K: Discovery of global genomic re-organization based on comparison of two newly sequenced rice mitochondrial genomes with cytoplasmic male sterility-related genes. BMC Genomics. 2010, 11: 209-10.1186/1471-2164-11-209.PubMed CentralView ArticlePubMedGoogle Scholar
- Satoh M, Kubo T, Nishizaw S, Estiati A, Itchoda N, Mikami T: The cytoplasmic male-sterile type and normal type mitochondrial genomes of sugar beet share the same complement of genes of known function but differ in the content of expressed ORFs. Mol Genet Genomics. 2004, 272: 247-256. 10.1007/s00438-004-1058-9.View ArticlePubMedGoogle Scholar
- Kubo T, Nishizawa S, Sugawara A, Itchoda N, Estiati A, Mikami T: The complete nucleotide sequence of the mitochondrial genome of sugar beet (Beta vulgaris L.) reveals a novel gene for tRNA(Cys) (GCA). Nucleic Acids Res. 2000, 28: 2571-2576. 10.1093/nar/28.13.2571.PubMed CentralView ArticlePubMedGoogle Scholar
- Clifton SW, Minx P, Fauron CM, Gibson M, Allen JO, Sun H, Thompson M, Barbazuk WB, Kanuganti S, Tayloe C, et al: Sequence and comparative analysis of the maize NB mitochondrial genome. Plant Physiol. 2004, 136: 486-3503.View ArticleGoogle Scholar
- Allen JO, Fauron CM, Minx P, Roark L, Oddiraju S, Lin GN, Meyer L, Sun H, Kim K, Wang C, et al: Comparison among two fertile and three male-sterile mitochondrial genomes of maize. Genetics. 2007, 177: 1173-1192. 10.1534/genetics.107.073312.PubMed CentralView ArticlePubMedGoogle Scholar
- Sugiyama Y, Watase Y, Nagase M, Makita N, Yagura S, Hirai A, Sugiura M: The complete nucleotide sequence and multipartite organization of the tobacco mitochondrial genome: comparative analysis of mitochondrial genomes in higher plants. Mol Genet Genomics. 2005, 272: 603-615. 10.1007/s00438-004-1075-8.View ArticlePubMedGoogle Scholar
- Ogihara Y, Yamazaki Y, Murai K, Kanno A, Terachi T, Shiina T, Miyashita N, Nasuda S, Nakamura C, Mori N, et al: Structural dynamics of cereal mitochondrial genomes as revealed by complete nucleotide sequencing of the wheat mitochondrial genome. Nucleic Acids Res. 2005, 33: 6235-6250. 10.1093/nar/gki925.PubMed CentralView ArticlePubMedGoogle Scholar
- Liu HT, Cui P, Zhan KH, Qiang L, Zhuo GY, Guo XL, Ding F, Yang WL, Liu DC, et al: Comparative analysis of mitochondrial genomes between a wheat K-type cytoplasmic male sterility (CMS) line and its maintainer line. BMC Genomics. 2011, 12: 163-10.1186/1471-2164-12-163.PubMed CentralView ArticlePubMedGoogle Scholar
- Chen JM, Guan RZ, Chang SX, Du TQ, Zhang HS, Xing H: Substoichiometrically different mitotypes coexist in mitochondrial genomes of Brassica napus L. PLoS One. 2011, 6: 1-8.Google Scholar
- Handa H: The complete nucleotide sequence and RNA editing content of the mitochondrial genome of rapeseed (Brassica napus L.): comparative analysis of the mitochondrial genomes of rapeseed and Arabidopsis thaliana. Nucleic Acids Res. 2003, 31: 5907-5916. 10.1093/nar/gkg795.PubMed CentralView ArticlePubMedGoogle Scholar
- Chang S, Yang TT, Du TQ, Huang YQ, Chen JM, Yan LY, He LB, Guan RZ: Mitochondrial genome sequencing helps show the evolutionary mechanism of mitochondrial genome formation in Brassica. BMC Genomics. 2011, 12: 479-10.1186/1471-2164-12-479.View ArticleGoogle Scholar
- Liu JH, Xu XY, Deng XX: Intergeneric somatic hybridization and its application to crop genetic improvement. Plant Cell, Tiss Org Cult. 2005, 82: 19-44. 10.1007/s11240-004-6015-0.View ArticleGoogle Scholar
- Yamagishi H, Glimelius K: Somatic hybrids between Arabidopsis thaliana and cytoplasmic male-sterile radish (Raphanus sativus). Plant Cell Rep. 2003, 22: 52-58. 10.1007/s00299-003-0655-0.View ArticlePubMedGoogle Scholar
- Sakai T, Imamura J: Intergeneric transfer of cytoplasmic male sterility between Raphanus sativus (CMS line) and Brassica napus through cytoplast-protoplast fusion. Theor Appl Genet. 1990, 80: 421-427.View ArticlePubMedGoogle Scholar
- Jourdan PS, Earle ED, Mutschler MA: Synthesis of male sterile, triazine-resistant Brassica napus by somatic hybridization between cytoplasmic male sterile B. oleracea and atrazine-resistant B. campestris. Theor Appl Genet. 1989, 78: 445-455.PubMedGoogle Scholar
- Pelletier G, Primard C, Vedel F, Chetrit P, Remy R, Rousselle , Renard M: Intergeneric cytoplasmic hybridization in cruciferae by protoplast fusion. Mol Gen Genet. 1983, 191: 244-250. 10.1007/BF00334821.View ArticleGoogle Scholar
- Bannerot H, Boulidard L, Cauderon Y, Tempe J: Transfer of cytoplasmic male sterility from Raphanus sativus to Brassica oleracea. Proc. Eucarpia Meet. Cruciferae. 1974, 25: 52-54.Google Scholar
- Krishnasamy S, Makaroff CA: Characterization of the radish mitochondrial orfB locus: Possible relationship with male sterility in Ogura radish. Curr Genet. 1993, 24: 156-163. 10.1007/BF00324680.View ArticlePubMedGoogle Scholar
- Bonhomme S, Budar F, Lancelin D, Small I, Defrance MC, Pelletier G: Sequence and transcript analysis of the Nco2.5 Ogura-specific fragment correlated with cytoplasmic male sterility in Brassica cybrids. Mol Gen Genet. 1992, 235: 340-348. 10.1007/BF00279379.View ArticlePubMedGoogle Scholar
- Jing B, Heng S, Tong D, Wan Z, Fu T, Tu J, Ma C, Yi B, Wen J, Shen J: A male sterility-associated cytotoxic protein ORF288 in Brassica juncea caused aborted pollen development. J Exp Bot. 2012, 63: 1285-1295. 10.1093/jxb/err355.PubMed CentralView ArticlePubMedGoogle Scholar
- Dewey RE, Timothy DH, Levings CS: A mitochondrial protein associated with cytoplasmic male sterility in the T cytoplasm of maize. Proc Natl Acad Sci, USA. 1987, 84: 5374-5378. 10.1073/pnas.84.15.5374.PubMed CentralView ArticlePubMedGoogle Scholar
- Wise RP, Pring DR, Gengenbach BG: Mutation to male fertility and toxin insensitivity in Texas (T)-cytoplasm maize is associated with a frameshift in a mitochondrial open reading frame. Proc Natl Acad Sci, USA. 1987, 84: 2858-2862. 10.1073/pnas.84.9.2858.PubMed CentralView ArticlePubMedGoogle Scholar
- Wise RP, Broson CR, Schnable PS, Homer HT: The genetics, pathology, and molecular biology of T-cytoplasm male sterility in maize. Adv Agro. 1999, 65: 79-130.View ArticleGoogle Scholar
- Akagi H, Nakamura A, Yokozeki-Misono Y, Inagaki A, Takahashi H, Mori K, Fujimura T: Positional cloning of the rice Rf-1 gene, a restorer of BT-type cytoplasmic male sterility that encodes a mitochondria-targeting PPR protein. Theor Appl Genet. 2004, 108: 1449-1457. 10.1007/s00122-004-1591-2.View ArticlePubMedGoogle Scholar
- Wang ZH, Zou Y, Li X, et al: Cytoplasmic male sterility of rice with boro II cytoplasm is caused by a cytotoxic peptide and is restored by two related PPR motif genes via distinct modes of mRNA silencing. Plant Cell. 2006, 18: 676-687. 10.1105/tpc.105.038240.PubMed CentralView ArticlePubMedGoogle Scholar
- Palmer JD, Shields CR: Tripartite structure of the Brassica campestris mitochondrial genome. Nature. 1984, 307: 437-440. 10.1038/307437a0.View ArticleGoogle Scholar
- Satoh M, Nemoto Y, Kawano S, Nagata T, Hirokawa H, Kuroiwa T: Organization of heterogeneous mitochondrial DNA molecules in mitochondrial nuclei of cultured tobacco cells. Protoplasma. 1993, 175: 112-120. 10.1007/BF01385008.View ArticleGoogle Scholar
- Andre C, Levy A, Walbot V: Small repeated sequences and the structure of plant mitochondrial genomes. Trends Genet. 1992, 8: 128-132.PubMedGoogle Scholar
- Oshima M, Kikuchi R, Imamura J, Handa H: Origin of the CMS gene locus in rapeseed cybrid mitochondria: active and inactive recombination produces the complex CMS gene region in the mitochondrial genomes of Brassicaceae. Genes Genet Syst. 2010, 85: 311-318. 10.1266/ggs.85.311.View ArticlePubMedGoogle Scholar
- Nakazono M, Nishiwaki S, Tsutsumi N: A chloroplast-derived sequence is utilized as a source of promoter sequences for the gene for subunit 9 of NADH dehydrogenase (nad9) in rice mitochondria. Mol Gen Genet. 1996, 252: 371-378.PubMedGoogle Scholar
- Handa H, Bonnard HG, Grienenberger JM: The rapeseed mitochondrial gene encoding a homologue of the bacterial protein Cc11 is divided into independently transcribed reading frames. Mol Gen Genet. 1996, 252: 292-302. 10.1007/BF02173775.View ArticlePubMedGoogle Scholar
- Heazlewood JL, Whelan J, Millar AH: The products of the mitochondrial orf25 and orfB genes are FO components in the plant F1FO ATP synthase. FEBS Lett. 2003, 540: 201-205. 10.1016/S0014-5793(03)00264-3.View ArticlePubMedGoogle Scholar
- Jing B, Heng S, Tong D, Wan Z, Fu T, Tu J, Ma C, Yi B, Wen J, Shen J: A male sterility-associated cytotoxic protein ORF288 in Brassica juncea causes aborted pollen development. J Exp Bot. 63: 1285-1295.
- Duroc Y, Gaillard C, Hiard S, Defrance MC, Pelletier G, Budar F: Biochemical and functional characterization of ORF138, a mitochondrial protein responsible for Ogura cytoplasmic male sterility in Brassicaceae. Biochimie. 2005, 87: 1089-1100. 10.1016/j.biochi.2005.05.009.View ArticlePubMedGoogle Scholar
- Desloire S, Gherbi H, Laloui W, Marhadour S, Clouet V, Cattolico L, Falentin C, Giancola S, Renard M, Budar F, et al: Identification of the fertility restoration locus, Rfo, in radish, as a member of the pentatricopeptide-repeat protein family. EMBO Rep. 2003, 4: 588-594. 10.1038/sj.embor.embor848.PubMed CentralView ArticlePubMedGoogle Scholar
- Giancola S, Marhadour S, Desloire S, Clouet V, Falentin-Guyomarc’h H, Laloui W, Falentin C, Pelletier G, Renard M, Bendahmane A, et al: Characterization of a radish introgression carrying the Ogura fertility restorer gene Rfo in rapeseed, using the Arabidopsis genome sequence and radish genetic mapping. Theor Appl Genet. 2003, 107: 1442-1451. 10.1007/s00122-003-1381-2.View ArticlePubMedGoogle Scholar
- Uyttewaal M, Arnal N, Quadrado M, Martin-Canadell A, Vrielynck N, Hiard S, Gherbi H, Bendahmane A, Budar F, Mireau H: Characterization of Raphanus sativus pentatricopeptide repeat proteins encoded by the fertility restorer locus for ogura cytoplasmic male sterility. Plant Cell. 2008, 20: 3331-3345. 10.1105/tpc.107.057208.PubMed CentralView ArticlePubMedGoogle Scholar
- Schattner P, Brooks AN, Lowe TM: The tRNAscan-SE, snoscan and snoGPS wed servers for the detection of tRNAs and snoRNAs. Nucl Acids Res. 2005, 33: W686-W689. 10.1093/nar/gki366.PubMed CentralView ArticlePubMedGoogle Scholar
- Schwartz S, Zhang Z, Frazer KA, Smit A, Riemer C, Bouck J, Gibbs R, Hardison R, Miller W: PipMaker-a web server for aligning two genomic DNA sequences. Genome Res. 2000, 10: 577-586. 10.1101/gr.10.4.577.PubMed CentralView ArticlePubMedGoogle Scholar