1
|
Liu H, Wang L, Guo Z, Xu Q, Fan W, Xu Y, Hu J, Zhang Y, Tang J, Xie M, Zhou Z, Hou S. Genome-wide association and selective sweep analyses reveal genetic loci for FCR of egg production traits in ducks. Genet Sel Evol 2021; 53:98. [PMID: 34930109 PMCID: PMC8690979 DOI: 10.1186/s12711-021-00684-5] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2020] [Accepted: 11/10/2021] [Indexed: 12/30/2022] Open
Abstract
Background As a major economic trait in poultry, egg production efficiency attracts widespread interest in breeding and production. However, limited information is available about the underlying genetic architecture of egg production traits in ducks. In this paper, we analyzed six egg production-related traits in 352 F2 ducks derived from reciprocal crosses between mallard and Pekin ducks. Results Feed conversation ratio (FCR) was positively correlated with feed intake but negatively correlated with egg-related traits, including egg weight and egg production, both phenotypically and genetically. Estimates of pedigree-based heritability were higher than 0.2 for all traits investigated, except hip-width. Based on whole-genome sequencing data, we conducted genome-wide association studies to identify genomic regions associated with these traits. In total, 11 genomic regions were associated with FCR. No genomic regions were identified as significantly associated with hip-width, total feed intake, average daily feed intake, and total egg production. Analysis of selective sweeps between mallard and Pekin ducks confirmed three of these genomic regions on chromosomes 13, 3 and 6. Within these three regions, variants in candidate genes that were in linkage disequilibrium with the GWAS leader single nucleotide polymorphisms (SNPs) (Chr13:2,196,728, P = 7.05 × 10–14; Chr3:76,991,524, P = 1.06 × 10–12; Chr6:20,356,803, P = 1.14 × 10–10) were detected. Thus, we identified 31 potential candidate genes associated with FCR, among which the strongest candidates are those that are highly expressed in tissues involved in reproduction and nervous system functions of ducks: CNTN4, CRBR, GPR63, KLHL32, FHL5, TRNT1, MANEA, NDUFAF4, and SCD. Conclusions For the first time, we report the identification of genomic regions that are associated with FCR in ducks and our results illustrate the genomic changes that occurred during their domestication and are involved in egg production efficiency. Supplementary Information The online version contains supplementary material available at 10.1186/s12711-021-00684-5.
Collapse
Affiliation(s)
- Hehe Liu
- State Key Laboratory of Animal Nutrition; Key Laboratory of Animal (Poultry) Genetics Breeding and Reproduction, Ministry of Agriculture and Rural Affairs, Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, 100193, China.,Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, College of Animal Science and Technology, Sichuan Agricultural University, Chengdu, 613000, China
| | - Lei Wang
- Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, College of Animal Science and Technology, Sichuan Agricultural University, Chengdu, 613000, China
| | - Zhanbao Guo
- State Key Laboratory of Animal Nutrition; Key Laboratory of Animal (Poultry) Genetics Breeding and Reproduction, Ministry of Agriculture and Rural Affairs, Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, 100193, China
| | - Qian Xu
- Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, College of Animal Science and Technology, Sichuan Agricultural University, Chengdu, 613000, China
| | - Wenlei Fan
- State Key Laboratory of Animal Nutrition; Key Laboratory of Animal (Poultry) Genetics Breeding and Reproduction, Ministry of Agriculture and Rural Affairs, Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, 100193, China
| | - Yaxi Xu
- State Key Laboratory of Animal Nutrition; Key Laboratory of Animal (Poultry) Genetics Breeding and Reproduction, Ministry of Agriculture and Rural Affairs, Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, 100193, China
| | - Jian Hu
- State Key Laboratory of Animal Nutrition; Key Laboratory of Animal (Poultry) Genetics Breeding and Reproduction, Ministry of Agriculture and Rural Affairs, Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, 100193, China
| | - Yunsheng Zhang
- State Key Laboratory of Animal Nutrition; Key Laboratory of Animal (Poultry) Genetics Breeding and Reproduction, Ministry of Agriculture and Rural Affairs, Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, 100193, China
| | - Jing Tang
- State Key Laboratory of Animal Nutrition; Key Laboratory of Animal (Poultry) Genetics Breeding and Reproduction, Ministry of Agriculture and Rural Affairs, Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, 100193, China
| | - Ming Xie
- State Key Laboratory of Animal Nutrition; Key Laboratory of Animal (Poultry) Genetics Breeding and Reproduction, Ministry of Agriculture and Rural Affairs, Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, 100193, China
| | - Zhengkui Zhou
- State Key Laboratory of Animal Nutrition; Key Laboratory of Animal (Poultry) Genetics Breeding and Reproduction, Ministry of Agriculture and Rural Affairs, Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, 100193, China
| | - Shuisheng Hou
- State Key Laboratory of Animal Nutrition; Key Laboratory of Animal (Poultry) Genetics Breeding and Reproduction, Ministry of Agriculture and Rural Affairs, Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, 100193, China.
| |
Collapse
|
2
|
Johnsson M, Henriksen R, Wright D. The neural crest cell hypothesis: no unified explanation for domestication. Genetics 2021; 219:iyab097. [PMID: 34849908 PMCID: PMC8633120 DOI: 10.1093/genetics/iyab097] [Citation(s) in RCA: 18] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [MESH Headings] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2021] [Accepted: 02/11/2021] [Indexed: 12/03/2022] Open
Affiliation(s)
- Martin Johnsson
- Department of Animal Breeding and Genetics, Swedish University of Agricultural Sciences, Uppsala 750 07, Sweden
| | - Rie Henriksen
- IFM Biology, Linköping University, Linköping 58183, Sweden
| | - Dominic Wright
- IFM Biology, Linköping University, Linköping 58183, Sweden
| |
Collapse
|
3
|
Wang L, Guo J, Xi Y, Ma S, Li Y, He H, Wang J, Han C, Bai L, Mustafa A, Liu H, Li L. Understanding the Genetic Domestication History of the Jianchang Duck by Genotyping and Sequencing of Genomic Genes Under Selection. G3 (BETHESDA, MD.) 2020; 10:1469-1476. [PMID: 32165372 PMCID: PMC7202016 DOI: 10.1534/g3.119.400893] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/22/2019] [Accepted: 03/01/2020] [Indexed: 12/11/2022]
Abstract
The Jianchang duck is mainly distributed in Southwest China, and has the characteristics of fast growth rate and strong abilities in lipid deposition in the liver. In order to investigate the effects of domestication process on formation of the unique characteristics of Jianchang duck, the whole genome of sixteen individuals and three pooling of Jianchang duck were re-sequenced, and genome data of 70 mallards and 83 domestic ducks from thirteen different places in China were obtained from NCBI. The population stratification and evolution analysis showed gene exchanges existed between the Jianchang and other domestic duck populations, as well as Jianchang ducks and mallards. Genomic comparison between mallards and Jianchang ducks showed genes, including CNTN1, CHRNA9, and SHANK2, which is involved in brain and nerve development, experienced strong positive selection in the process of Jianchang duck domestication. The genomic comparison between Jianchang and domestic duck populations showed that HSD17B12 and ESM1, which affect lipid metabolism, experienced strong positive selection during the domestication process. FST analysis among populations of Jianchang duck with different plumage colors indicated that MITF was related to the phenotype of a white feather, while MC1R was related to the phenotype of hemp feather. Our results provided a base for the domestication process of Jianchang duck and the genomic genes for unique traits.
Collapse
Affiliation(s)
- Lei Wang
- Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, College of Animal Science and Technology, Sichuan Agricultural University, Chengdu, P.R. China
| | - Jiazhong Guo
- Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, College of Animal Science and Technology, Sichuan Agricultural University, Chengdu, P.R. China
| | - Yang Xi
- Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, College of Animal Science and Technology, Sichuan Agricultural University, Chengdu, P.R. China
| | - Shengchao Ma
- Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, College of Animal Science and Technology, Sichuan Agricultural University, Chengdu, P.R. China
| | - Yanying Li
- Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, College of Animal Science and Technology, Sichuan Agricultural University, Chengdu, P.R. China
| | - Hua He
- Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, College of Animal Science and Technology, Sichuan Agricultural University, Chengdu, P.R. China
| | - Jiwen Wang
- Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, College of Animal Science and Technology, Sichuan Agricultural University, Chengdu, P.R. China
| | - Chunchun Han
- Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, College of Animal Science and Technology, Sichuan Agricultural University, Chengdu, P.R. China
| | - Lili Bai
- Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, College of Animal Science and Technology, Sichuan Agricultural University, Chengdu, P.R. China
| | - Ahsan Mustafa
- Institute of Animal Nutrition, Key Laboratory for Animal Disease-Resistance Nutrition of China, Ministry of Education, Sichuan Agricultural University, Chengdu, P.R. China
| | - Hehe Liu
- Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, College of Animal Science and Technology, Sichuan Agricultural University, Chengdu, P.R. China
| | - Liang Li
- Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, College of Animal Science and Technology, Sichuan Agricultural University, Chengdu, P.R. China
| |
Collapse
|
4
|
Salas CA, Yopak KE, Warrington RE, Hart NS, Potter IC, Collin SP. Ontogenetic shifts in brain scaling reflect behavioral changes in the life cycle of the pouched lamprey Geotria australis. Front Neurosci 2015; 9:251. [PMID: 26283894 PMCID: PMC4517384 DOI: 10.3389/fnins.2015.00251] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/17/2015] [Accepted: 07/03/2015] [Indexed: 12/11/2022] Open
Abstract
Very few studies have described brain scaling in vertebrates throughout ontogeny and none in lampreys, one of the two surviving groups of the early agnathan (jawless) stage in vertebrate evolution. The life cycle of anadromous parasitic lampreys comprises two divergent trophic phases, firstly filter-feeding as larvae in freshwater and secondly parasitism as adults in the sea, with the transition marked by a radical metamorphosis. We characterized the growth of the brain during the life cycle of the pouched lamprey Geotria australis, an anadromous parasitic lamprey, focusing on the scaling between brain and body during ontogeny and testing the hypothesis that the vast transitions in behavior and environment are reflected in differences in the scaling and relative size of the major brain subdivisions throughout life. The body and brain mass and the volume of six brain structures of G. australis, representing six points of the life cycle, were recorded, ranging from the early larval stage to the final stage of spawning and death. Brain mass does not increase linearly with body mass during the ontogeny of G. australis. During metamorphosis, brain mass increases markedly, even though the body mass does not increase, reflecting an overall growth of the brain, with particularly large increases in the volume of the optic tectum and other visual areas of the brain and, to a lesser extent, the olfactory bulbs. These results are consistent with the conclusions that ammocoetes rely predominantly on non-visual and chemosensory signals, while adults rely on both visual and olfactory cues.
Collapse
Affiliation(s)
- Carlos A Salas
- Neuroecology Group, School of Animal Biology and UWA Oceans Institute, The University of Western Australia Crawley, WA, Australia
| | - Kara E Yopak
- Neuroecology Group, School of Animal Biology and UWA Oceans Institute, The University of Western Australia Crawley, WA, Australia
| | - Rachael E Warrington
- Neuroecology Group, School of Animal Biology and UWA Oceans Institute, The University of Western Australia Crawley, WA, Australia
| | - Nathan S Hart
- Neuroecology Group, School of Animal Biology and UWA Oceans Institute, The University of Western Australia Crawley, WA, Australia
| | - Ian C Potter
- Centre for Fish and Fisheries Research, School of Veterinary and Life Sciences, Murdoch University Murdoch, WA, Australia
| | - Shaun P Collin
- Neuroecology Group, School of Animal Biology and UWA Oceans Institute, The University of Western Australia Crawley, WA, Australia
| |
Collapse
|
5
|
Kawabe S, Matsuda S, Tsunekawa N, Endo H. Ontogenetic Shape Change in the Chicken Brain: Implications for Paleontology. PLoS One 2015; 10:e0129939. [PMID: 26053849 PMCID: PMC4460028 DOI: 10.1371/journal.pone.0129939] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/13/2015] [Accepted: 05/14/2015] [Indexed: 11/19/2022] Open
Abstract
Paleontologists have investigated brain morphology of extinct birds with little information on post-hatching changes in avian brain morphology. Without the knowledge of ontogenesis, assessing brain morphology in fossil taxa could lead to misinterpretation of the phylogeny or neurosensory development of extinct species. Hence, it is imperative to determine how avian brain morphology changes during post-hatching growth. In this study, chicken brain shape was compared at various developmental stages using three-dimensional (3D) geometric morphometric analysis and the growth rate of brain regions was evaluated to explore post-hatching morphological changes. Microscopic MRI (μMRI) was used to acquire in vivo data from living and post-mortem chicken brains. The telencephalon rotates caudoventrally during growth. This change in shape leads to a relative caudodorsal rotation of the cerebellum and myelencephalon. In addition, all brain regions elongate rostrocaudally and this leads to a more slender brain shape. The growth rates of each brain region were constant and the slopes from the growth formula were parallel. The dominant pattern of ontogenetic shape change corresponded with interspecific shape changes due to increasing brain size. That is, the interspecific and ontogenetic changes in brain shape due to increased size have similar patterns. Although the shape of the brain and each brain region changed considerably, the volume ratio of each brain region did not change. This suggests that the brain can change its shape after completing functional differentiation of the brain regions. Moreover, these results show that consideration of ontogenetic changes in brain shape is necessary for an accurate assessment of brain morphology in paleontological studies.
Collapse
Affiliation(s)
- Soichiro Kawabe
- Gifu Prefectural Museum, Gifu, Japan
- The University Museum, The University of Tokyo, Tokyo, Japan
| | - Seiji Matsuda
- Department of Anatomy and Embryology, School of Medicine, Ehime University, Ehime, Japan
| | - Naoki Tsunekawa
- Department of Veterinary Anatomy, The University of Tokyo, Tokyo, Japan
| | - Hideki Endo
- The University Museum, The University of Tokyo, Tokyo, Japan
| |
Collapse
|
6
|
van Ballegooijen M, Rutter CM, Knudsen AB, Zauber AG, Savarino JE, Lansdorp-Vogelaar I, Boer R, Feuer EJ, Habbema JDF, Kuntz KM. Clarifying differences in natural history between models of screening: the case of colorectal cancer. Med Decis Making 2011; 31:540-9. [PMID: 21673187 DOI: 10.1177/0272989x11408915] [Citation(s) in RCA: 37] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
BACKGROUND Microsimulation models are important decision support tools for screening. However, their complexity makes them difficult to understand and limits realization of their full potential. Therefore, it is important to develop documentation that clarifies their structure and assumptions. The authors demonstrate this problem and explore a solution for natural history using 3 independently developed colorectal cancer screening models. METHODS The authors first project effectiveness and cost-effectiveness of colonoscopy screening for the 3 models (CRC-SPIN, SimCRC, and MISCAN). Next, they provide a conventional presentation of each model, including information on structure and parameter values. Finally, they report the simulated reduction in clinical cancer incidence following a one-time complete removal of adenomas and preclinical cancers for each model. They call this new measure the maximum clinical incidence reduction (MCLIR). RESULTS Projected effectiveness varies widely across models. For example, estimated mortality reduction for colonoscopy screening every 10 years from age 50 to 80 years, with surveillance in adenoma patients, ranges from 65% to 92%. Given only conventional information, it is difficult to explain these differences, largely because differences in structure make parameter values incomparable. In contrast, the MCLIR clearly shows the impact of model differences on the key feature of natural history, the dwell time of preclinical disease. Dwell times vary from 8 to 25 years across models and help explain differences in projected screening effectiveness. CONCLUSIONS The authors propose a new measure, the MCLIR, which summarizes the implications for predicted screening effectiveness of differences in natural history assumptions. Including the MCLIR in the standard description of a screening model would improve the transparency of these models.
Collapse
Affiliation(s)
| | | | - Amy B Knudsen
- Institute for Technology Assessment, Massachusetts General Hospital, Boston (ABK)
| | - Ann G Zauber
- Department of Epidemiology and Biostatistics, Memorial SloanKettering Cancer Center, New York, New York (AGZ)
| | | | - Iris Lansdorp-Vogelaar
- Department of Public Health, Erasmus MC, Rotterdam, the Netherlands (MvB, IL-V, RB, JDFH)
| | - Rob Boer
- Department of Public Health, Erasmus MC, Rotterdam, the Netherlands (MvB, IL-V, RB, JDFH)
| | - Eric J Feuer
- Division of Cancer Control and Population Sciences, National Cancer Institute, Washington, DC (EJF)
| | - J Dik F Habbema
- Department of Public Health, Erasmus MC, Rotterdam, the Netherlands (MvB, IL-V, RB, JDFH)
| | - Karen M Kuntz
- Division of Health Policy and Management, School of Public Health, University of Minnesota, Minneapolis (KMK)
| |
Collapse
|
7
|
Rutter CM, Savarino JE. An evidence-based microsimulation model for colorectal cancer: validation and application. Cancer Epidemiol Biomarkers Prev 2010; 19:1992-2002. [PMID: 20647403 DOI: 10.1158/1055-9965.epi-09-0954] [Citation(s) in RCA: 39] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023] Open
Abstract
BACKGROUND The Colorectal Cancer Simulated Population model for Incidence and Natural history (CRC-SPIN) is a new microsimulation model for the natural history of colorectal cancer that can be used for comparative effectiveness studies of colorectal cancer screening modalities. METHODS CRC-SPIN simulates individual event histories associated with colorectal cancer, based on the adenoma-carcinoma sequence: adenoma initiation and growth, development of preclinical invasive colorectal cancer, development of clinically detectable colorectal cancer, death from colorectal cancer, and death from other causes. We present the CRC-SPIN structure and parameters, data used for model calibration, and model validation. We also provide basic model outputs to further describe CRC-SPIN, including annual transition probabilities between various disease states and dwell times. We conclude with a simple application that predicts the impact of a one-time colonoscopy at age 50 on the incidence of colorectal cancer assuming three different operating characteristics for colonoscopy. RESULTS CRC-SPIN provides good prediction of both the calibration and the validation data. Using CRC-SPIN, we predict that a one-time colonoscopy greatly reduces colorectal cancer incidence over the subsequent 35 years. CONCLUSIONS CRC-SPIN is a valuable new tool for combining expert opinion with observational and experimental results to predict the comparative effectiveness of alternative colorectal cancer screening modalities. IMPACT Microsimulation models such as CRC-SPIN can serve as a bridge between screening and treatment studies and health policy decisions by predicting the comparative effectiveness of different interventions. As such, it is critical to publish model descriptions that provide insight into underlying assumptions along with validation studies showing model performance.
Collapse
Affiliation(s)
- Carolyn M Rutter
- Group Health Research Institute, 1630 Minor Avenue, Seattle, WA 98101, USA.
| | | |
Collapse
|
8
|
Abstract
Microsimulation models that describe disease processes synthesize information from multiple sources and can be used to estimate the effects of screening and treatment on cancer incidence and mortality at a population level. These models are characterized by simulation of individual event histories for an idealized population of interest. Microsimulation models are complex and invariably include parameters that are not well informed by existing data. Therefore, a key component of model development is the choice of parameter values. Microsimulation model parameter values are selected to reproduce expected or known results though the process of model calibration. Calibration may be done by perturbing model parameters one at a time or by using a search algorithm. As an alternative, we propose a Bayesian method to calibrate microsimulation models that uses Markov chain Monte Carlo. We show that this approach converges to the target distribution and use a simulation study to demonstrate its finite-sample performance. Although computationally intensive, this approach has several advantages over previously proposed methods, including the use of statistical criteria to select parameter values, simultaneous calibration of multiple parameters to multiple data sources, incorporation of information via prior distributions, description of parameter identifiability, and the ability to obtain interval estimates of model parameters. We develop a microsimulation model for colorectal cancer and use our proposed method to calibrate model parameters. The microsimulation model provides a good fit to the calibration data. We find evidence that some parameters are identified primarily through prior distributions. Our results underscore the need to incorporate multiple sources of variability (i.e., due to calibration data, unknown parameters, and estimated parameters and predicted values) when calibrating and applying microsimulation models.
Collapse
Affiliation(s)
- Carolyn M. Rutter
- Carolyn M. Rutter is Senior Investigator, Group Health Center for Health Studies, Seattle, WA 98101 () and Affiliate Professor, Departments of Biostatistics and Health Services, University of Washington, WA 98195. Diana L. Miglioretti is Senior Investigator, Group Health Center for Health Studies, Seattle, WA 98101 and Affiliate Associate Professor, Department of Biostatistics, University of Washington, WA 98195. James E. Savarino is Programmer, Group Health Center for Health Studies, Seattle, WA 98101
| | - Diana L. Miglioretti
- Carolyn M. Rutter is Senior Investigator, Group Health Center for Health Studies, Seattle, WA 98101 () and Affiliate Professor, Departments of Biostatistics and Health Services, University of Washington, WA 98195. Diana L. Miglioretti is Senior Investigator, Group Health Center for Health Studies, Seattle, WA 98101 and Affiliate Associate Professor, Department of Biostatistics, University of Washington, WA 98195. James E. Savarino is Programmer, Group Health Center for Health Studies, Seattle, WA 98101
| | - James E. Savarino
- Carolyn M. Rutter is Senior Investigator, Group Health Center for Health Studies, Seattle, WA 98101 () and Affiliate Professor, Departments of Biostatistics and Health Services, University of Washington, WA 98195. Diana L. Miglioretti is Senior Investigator, Group Health Center for Health Studies, Seattle, WA 98101 and Affiliate Associate Professor, Department of Biostatistics, University of Washington, WA 98195. James E. Savarino is Programmer, Group Health Center for Health Studies, Seattle, WA 98101
| |
Collapse
|