Protein intake from plant sources and from all the animal food items considered were included in the multivariable model that was also adjusted for total caloric intake and percentage of energy from saturated fat, polyunsaturated fat, monounsaturated fat, and trans-fat (all continuous), multivitamin use (yes or no), smoking status (never, past, or current [1-14, and ≥15 cigarettes/d]), pack-years of smoking (in women, ≤15, 16-25, 26-45, and ≥46; in men, <10, 11-24, 25-44, and ≥45), body mass index (calculated as weight in kilograms divided by height in meters squared; <23.0, 23.0-24.9, 25.0-26.9, 27.0-29.9, 30.0-34.9, and ≥35), physical activity (quintiles), alcohol consumption (in women, 0, 0.1-5.0, 5.1-15.0, and >15.0 g/d; in men, 0, 0.1-10.0, 10.1-20.0, and >20.0 g/d), history of hypertension diagnosis (yes or no), glycemic index (in quintiles), and intake of whole grains, total fiber, fruits, and vegetables (all in quintiles). CVD indicates cardiovascular disease; HR, hazard ratio. Error bars indicate 95% CIs.
eMethods. Assessment and Statistical Analysis
eTable 1. Risk for All-Cause Mortality According to Percentage of Energy From Animal and Plant Protein Intake in the Nurses’ Health Study and Health Professionals Follow-up Study
eTable 2. Hazard Ratio (95% CI) for All-Cause Mortality According to Percentage of Energy From Animal and Plant Protein Intake by Duration of Follow-up
eTable 3. Risk for All-Cause Mortality According to Percentage of Energy From Animal and Plant Protein Intake With Adjustment for Propensity Score
eTable 4. Association of Animal and Plant Protein Intake With All-Cause Mortality According to History of Diabetes
eTable 5. Median Protein Intake (% of Energy) of Various Food Origins and Their Correlations
eTable 6. Hazard Ratio (95% CI) of All-Cause and Cause-Specific Mortality According to Percentage of Energy Intake From Protein of Various Animal Foods With Mutual Adjustment for Each Other
eFigure 1. Median Intake and Interquartile Range of Animal and Plant Protein in Women (A) and Men (B) During Follow-up
eFigure 2. Association of Animal and Plant Protein Intake With All-Cause Mortality According to Age and Lifestyle Factors
Song M, Fung TT, Hu FB, Willett WC, Longo VD, Chan AT, Giovannucci EL. Association of Animal and Plant Protein Intake With All-Cause and Cause-Specific Mortality. JAMA Intern Med. 2016;176(10):1453-1463. doi:10.1001/jamainternmed.2016.4182
Defining what represents a macronutritionally balanced diet remains an open question and a high priority in nutrition research. Although the amount of protein may have specific effects, from a broader dietary perspective, the choice of protein sources will inevitably influence other components of diet and may be a critical determinant for the health outcome.
To examine the associations of animal and plant protein intake with the risk for mortality.
Design, Setting, and Participants
This prospective cohort study of US health care professionals included 131 342 participants from the Nurses’ Health Study (1980 to end of follow-up on June 1, 2012) and Health Professionals Follow-up Study (1986 to end of follow-up on January 31, 2012). Animal and plant protein intake was assessed by regularly updated validated food frequency questionnaires. Data were analyzed from June 20, 2014, to January 18, 2016.
Main Outcomes and Measures
Hazard ratios (HRs) for all-cause and cause-specific mortality.
Of the 131 342 participants, 85 013 were women (64.7%) and 46 329 were men (35.3%) (mean [SD] age, 49  years). The median protein intake, as assessed by percentage of energy, was 14% for animal protein (5th-95th percentile, 9%-22%) and 4% for plant protein (5th-95th percentile, 2%-6%). After adjusting for major lifestyle and dietary risk factors, animal protein intake was not associated with all-cause mortality (HR, 1.02 per 10% energy increment; 95% CI, 0.98-1.05; P for trend = .33) but was associated with higher cardiovascular mortality (HR, 1.08 per 10% energy increment; 95% CI, 1.01-1.16; P for trend = .04). Plant protein was associated with lower all-cause mortality (HR, 0.90 per 3% energy increment; 95% CI, 0.86-0.95; P for trend < .001) and cardiovascular mortality (HR, 0.88 per 3% energy increment; 95% CI, 0.80-0.97; P for trend = .007). These associations were confined to participants with at least 1 unhealthy lifestyle factor based on smoking, heavy alcohol intake, overweight or obesity, and physical inactivity, but not evident among those without any of these risk factors. Replacing animal protein of various origins with plant protein was associated with lower mortality. In particular, the HRs for all-cause mortality were 0.66 (95% CI, 0.59-0.75) when 3% of energy from plant protein was substituted for an equivalent amount of protein from processed red meat, 0.88 (95% CI, 0.84-0.92) from unprocessed red meat, and 0.81 (95% CI, 0.75-0.88) from egg.
Conclusions and Relevance
High animal protein intake was positively associated with cardiovascular mortality and high plant protein intake was inversely associated with all-cause and cardiovascular mortality, especially among individuals with at least 1 lifestyle risk factor. Substitution of plant protein for animal protein, especially that from processed red meat, was associated with lower mortality, suggesting the importance of protein source.
Defining what represents a macronutritionally balanced diet remains an open question and a high priority in nutrition research.1,2Quiz Ref ID In short-term randomized clinical trials, substitution of protein for carbohydrate has been shown to favor weight management, decrease blood pressure, and improve cardiometabolic biomarkers, including blood lipid and lipoprotein profiles and glycemic regulation.3- 5 These beneficial effects are partly dependent on weight loss and possibly owing to the enhanced postprandial satiety and energy expenditure when exchanging protein for carbohydrate.6 Therefore, high-protein and low-carbohydrate diets have been promoted for weight loss and health improvement. Although the amount and type of protein may have specific effects,7 such as insulinlike growth factor 1 levels,8 from a broader dietary perspective, the choice of protein sources will inevitably influence other components of diet, including macronutrients, micronutrients, and phytochemicals, that can in turn influence health outcomes. Therefore, taking into account food sources is critical to better understand the health effect of protein intake and fine-tune dietary recommendations.
To date, data examining protein sources in relation to mortality are sparse. Although no association was found between animal or plant protein and all-cause mortality in a cohort of postmenopausal women, substitution of plant protein for animal protein was associated with lower mortality due to cardiovascular disease (CVD).9 A positive association between animal protein and mortality was also found in the other study using the National Health and Nutrition Examination Survey.8 Nevertheless, these data are far from conclusive owing to several limitations, including the relatively small sample size, single assessment of diet at baseline, and lack of data on detailed food sources of animal and plant protein. Therefore, we used data from 2 large US cohort studies with repeated measures of diet and up to 32 years of follow-up to prospectively examine animal protein vs plant protein in relation to the risk for all-cause and cause-specific mortality and to perform an isocaloric substitution analysis for a variety of food sources of protein.
Question What is the association of the source of protein intake with mortality in US adults?
Findings In this cohort study, high intake of animal protein was positively associated with mortality, with the inverse true for high intake of plant protein, especially among individuals with at least 1 lifestyle risk factor. Replacement of animal protein with plant protein was associated with lower mortality, suggesting the importance of protein source.
Meaning Public health recommendations should focus on improvement of protein sources.
The Nurses’ Health Study (NHS)10 included 121 700 US registered female nurses who were aged 30 to 55 years in 1976; for this study, data were collected from 1980 to June 1, 2012. The Health Professionals Follow-up Study (HPFS)11 included 51 529 US male health care professionals who were aged 40 to 75 years in 1986; data were collected from 1986 to January 31, 2012. Details of the 2 cohorts have been described elsewhere.10,11 Briefly, follow-up questionnaires were administered at baseline enrollment and every 2 years thereafter to collect lifestyle and medical information. Dietary intake was assessed by the food frequency questionnaires (FFQs) every 4 years. The follow-up rates were 95.4% in the NHS and 95.9% in the HPFS. The study protocol was approved by the institutional review board at the Brigham and Women’s Hospital and the Harvard T. H. Chan School of Public Health. All participants provided written informed consent.
Among participants who returned baseline questionnaires, we excluded those who had a history of cancer (except nonmelanoma skin cancer), CVD, or diabetes at baseline, left more than 10 items blank on the baseline FFQ in the NHS and more than 70 items blank in the HPFS, or reported implausible energy intake levels (<500 or >3500 kcal/d for women, or <800 or >4200 kcal/d for men). After exclusions, 85 013 women and 46 329 men were available for the analysis.
In each FFQ, participants were asked how often, on average, they consumed a standardized portion size of each food during the previous year. The mean daily nutrient intake was calculated by multiplying the consumption frequency of each food item by its nutrient content and then summing across all foods. Animal and plant protein intake was expressed as a percentage of total energy consumption. Major sources of animal protein included processed and unprocessed red meat, poultry, dairy products, fish, and egg. Major food contributors to plant protein included bread, cereals, pasta, nuts, beans, and legumes. We derived protein intake from processed red meat by summing the products between intake frequency (servings per day) and the protein content (grams per serving) for various processed red meats (ie, bacon, beef or pork hot dogs, salami, bologna or other processed meat sandwiches, other processed meats [eg, sausage, kielbasa]). Similar calculations were performed for protein intake from unprocessed red meat, poultry, fish, egg, and dairy. Food frequency questionnaires have demonstrated good validity in assessing protein intake. The Spearman correlation coefficient of intake assessed by the FFQs and 7-day dietary record was 0.56 for animal protein and 0.66 for plant protein,12 as detailed in eMethods in the Supplement.
We identified deaths from state statistics records, the National Death Index, next of kin, and the postal system. Using these methods, we were able to ascertain more than 96% of the deaths in each cohort.13 Cause of death was identified from death certificates or review of medical records by physicians. For this analysis, we assessed all-cause mortality and deaths due to CVD (International Classification of Diseases, Eighth Revision, codes 390-458), cancer (International Classification of Diseases, Eighth Revision, codes 140 to 207), and other causes.
Data were analyzed from June 20, 2014, to January 18, 2016. We calculated person-time of follow-up for each participant from the age in months at the return date of the baseline FFQ (1980 for the NHS and 1986 for the HPFS) until the age in months at the date of death, loss to follow-up, or end of follow-up (June 1, 2012, for the NHS and January 31, 2012, for the HPFS), whichever came first. We used time-varying Cox proportional hazards regression models with age as the time scale to estimate the hazard ratio (HR) and 95% CI for mortality associated with animal and plant protein intake.
To reduce random within-person variation and to best represent long-term dietary intake, we calculated the cumulative mean protein intake from our repeated FFQs.14 We stopped updating dietary information when a participant reported a diagnosis of cancer (except nonmelanoma skin cancer), diabetes, stroke, coronary heart disease, or angina, because these conditions may lead to dietary change.15
We used a nutrient density model with adjustment for total energy intake and the percentage of energy from various fats (saturated, polyunsaturated, monounsaturated, and trans-fat).16 Thus, the coefficient for animal and plant protein reflects the substitution effect of an equal amount of energy from protein for carbohydrate. In the multivariable analysis, we further adjusted for several potential dietary and lifestyle confounding factors, including multivitamin use, smoking status, pack-years of smoking, body mass index, physical activity, alcohol consumption, history of hypertension diagnosis, glycemic index, and intake of whole grains, total fiber, fruits, and vegetables. To address the possibility of residual confounding, we further adjusted for a propensity score that reflected associations of protein consumption with potential confounding covariates.17 Details about covariate assessment and propensity score analysis are provided in the eMethods in the Supplement.
We performed stratified analyses by age and lifestyle factors and evaluated the interaction via a likelihood ratio test. To minimize the confounding effect and test for potential modification by an overall lifestyle pattern, we further performed a stratified analysis according to a priori–defined healthy lifestyle pattern, as characterized by never smoking or ever smoking for fewer than 5 pack-years, never or moderate alcohol intake (<14 g/d in women and <28 g/d in men), body mass index (calculated as weight in kilograms divided by height in meters squared) of at least 18.5 and less than 25.0, and physical activity of at least 150 min/wk at a moderate level or at least 75 min/wk at a vigorous level (equivalent to ≥7.5 metabolic equivalent h/wk) as recommended.18 Likewise, given the previous report that protein intake was associated with a higher risk for diabetes-related mortality,8 we examined the protein-mortality association according to the history of diabetes.
Finally, we estimated the effect of substituting 3% of energy from plant protein for an equivalent amount of animal protein from various sources, including processed and unprocessed red meat, poultry, fish, egg, and dairy, by simultaneously including these protein items as continuous variables in the multivariable model. The HRs and 95% CIs for the isoprotein substitution effect were derived from the difference between the regression coefficients, variance, and covariance.19
The analyses were first conducted in each cohort separately, and because no appreciable difference was detected by cohort (eTable 1 in the Supplement), we then conducted the pooled analysis using the sex-stratified Cox proportional hazards regression model in the combined data set. More details about statistical analysis are provided in the eMethods in the Supplement.
Quiz Ref IDIn the 2 cohorts with 3 540 791 person-years of follow-up, we documented 36 115 deaths, of which 8851 were due to CVD, 13 159 were due to cancer, and 14 105 were due to other causes.Quiz Ref ID Participants’ median intake, as assessed by percentage of energy, was 14% (5th-95th percentile, 9%-22%) for animal protein and 4% (5th-95th percentile, 2%-6%) for plant protein. Animal protein intake decreased, whereas plant protein intake increased over time throughout follow-up (eFigure 1 in the Supplement). Table 1 shows the basic characteristics of participants according to protein intake. Compared with participants who consumed no more than 10% of energy from animal protein, those consuming more than 18% were slightly heavier and less physically active and consumed more fats (especially saturated fat) and less fiber and plant foods. In contrast, those with higher plant protein intake demonstrated a clustering of positive health behaviors and had a substantially healthier diet than those with lower plant protein consumption.
As shown in Table 2, higher intake of animal protein was associated with higher CVD mortality. After adjusting for major lifestyle and dietary risk factors, the HR per 10% increment of animal protein intake from total energy intake was 1.02 (95% CI, 0.98-1.05; P for trend = .33) for all-cause mortality and 1.08 (95% CI, 1.01-1.16; P for trend = .04) for CVD mortality. Quiz Ref IDIn contrast, a higher level of plant protein intake was associated with lower mortality, with the multivariable HR per 3% increment of total energy intake of 0.90 (95% CI, 0.86-0.95; P for trend < .001) for all-cause mortality and 0.88 (95% CI, 0.80-0.97; P for trend = .007) for CVD mortality. The associations did not differ by duration of follow-up (eTable 2 in the Supplement). We did not detect any statistically significant nonlinear relationship between protein intake and mortality by spline analysis (data not shown). The results remained largely unchanged when we adjusted for a propensity score that predicted protein intake levels (eTable 3 in the Supplement).
The increased mortality associated with higher animal protein intake was more pronounced among obese participants (P for interaction = .008) and those with heavy alcohol intake (P for interaction = .06) (eFigure 2 in the Supplement). The association between higher plant protein intake and lower mortality was stronger among participants who were 65 years or younger or older than 80 years, currently smoked, consumed at least 14 g/d of alcohol, were overweight or obese, and were physically inactive (P for interaction ≤.02 for all).
Because most of the statistically significant associations were seen among participants with an unhealthy lifestyle, we further divided participants into healthy- and unhealthy-lifestyle groups according to a priori–defined criteria. Table 3 shows the basic characteristics of the 2 groups. Participants in the healthy-lifestyle group demonstrated slightly more homogeneous distributions in health behaviors than those in the unhealthy-lifestyle group. At similar amounts of protein intake, protein sources differed between the 2 groups. Compared with the healthy-lifestyle group, the unhealthy-lifestyle group with similar animal protein intake consumed more unprocessed and processed red meat, eggs, and high-fat dairy products, but less chicken, fish, and low-fat dairy products. At similar amounts of plant protein intake, the unhealthy-lifestyle group consumed less fiber, fruit, vegetables, and whole grains than the healthy-lifestyle group.
Table 4 shows the associations of protein intake and mortality in the 2 groups. The positive association with all-cause mortality for animal protein intake and the inverse association for plant protein intake were restricted to the unhealthy-lifestyle group (P for interaction <.001), although the association with animal protein intake did not reach statistical significance. In the unhealthy-lifestyle group, the multivariable HR per 10% increment of animal protein was 1.03 (95% CI, 0.99-1.07; P for trend = .16) and the HR per 3% increment of plant protein was 0.90 (95% CI, 0.85-0.95; P for trend < .001). Similar results were observed for CVD mortality. When stratified by history of diabetes, the positive association with all-cause mortality for animal protein intake and the inverse association for plant protein intake appeared to be stronger among participants with diabetes than those without diabetes (P for interaction = .06 and P for interaction .02, respectively; eTable 4 in the Supplement).
Finally, we examined the substitution association of different protein sources with mortality. The mean protein intake from various foods and their correlations are shown in eTable 5 in the Supplement, and their individual associations with mortality are summarized in eTable 6 in the Supplement. Protein intake from processed red meat was strongly associated with mortality, whereas no association was found for protein from fish or poultry. The Figure presents the HRs for mortality with substitution of 3% energy from plant protein for the same amount of animal protein from different food sources. The HRs for all-cause mortality were 0.66 (95% CI, 0.59-0.75) when 3% of energy from plant protein was substituted for an equivalent amount of protein from processed red meat; 0.88 (95% CI, 0.84-0.92), from unprocessed red meat; 0.94 (95% CI, 0.90-0.99), from poultry; 0.94 (95% CI, 0.89-0.99), from fish; 0.81 (95% CI, 0.75-0.88), from egg; and 0.92 (95% CI, 0.87-0.96), from dairy. The substitution associations were generally stronger for death due to CVD and other causes than those due to cancer, except substitution of egg, for which substitution of 3% energy plant protein was associated with 17% lower mortality due to cancer (95% CI, 7%-27%).
After adjusting for other dietary and lifestyle factors, animal protein intake was associated with a higher risk for CVD mortality, whereas higher plant protein intake was associated with lower all-cause and cardiovascular mortality. However, in the stratified analysis, these associations were confined to participants with at least 1 lifestyle risk factor. Moreover, we observed that substitution of plant protein for animal protein from a variety of food sources, particularly processed red meat, was associated with a lower risk for mortality, suggesting that the protein source is important for long-term health.
Although short-term randomized clinical trials have shown a beneficial effect of high protein intake,3,4,20,21 the long-term health consequences of protein intake remain controversial.8,9,22- 25 In a randomized clinical trial with a 2-year intervention, 4 calorie-restricted diets with different macronutrient compositions did not show a difference in the effects on weight loss or on improvement of lipid profiles and insulin levels.26 When protein is substituted for other macronutrients, the dietary source of protein appears to be a critical determinant of the outcome.
To our knowledge, only 2 cohort studies8,9 have examined animal and plant protein intake in relation to mortality. In the Iowa Women’s Health Study,9 although neither animal nor plant protein was associated with all-cause mortality, an inverse association was found between plant protein and CVD mortality, and substituting plant protein for animal protein was associated with a substantially lower CVD mortality. In a recent report from the National Health and Nutrition Examination Survey III,8 higher protein intake was related to an increased risk for all-cause mortality among participants younger than 65 years. However, when animal protein intake was controlled for, this association was eliminated, suggesting that animal protein was responsible for the effect of higher protein intake, if any, on increased mortality. Although a direct comparison of these studies is difficult, given the variation in the study methods,27 these data together with our current findings support the importance of protein sources for the long-term health outcome and suggest that plants constitute a preferred protein source compared with animal foods.
Indeed, unlike animal protein, plant protein has not been associated with increased insulinlike growth factor 1 levels28,29 and has been linked to lower blood pressure,30- 32 reduced low-density lipoprotein levels,32- 34 and improved insulin sensitivity.35 Substitution of plant protein for animal protein has been related to a lower incidence of CVD36- 39 and type 2 diabetes.40- 42 Moreover, although a high intake of red meat, particularly processed red meat, has been associated with increased mortality in a recent meta-analysis of 13 cohort studies,43 high consumption of nuts, a major contributor to plant protein, has been associated lower CVD and all-cause mortality.44 These results underscore the importance of protein sources for risk assessment and suggest that other components in protein-rich foods (eg, sodium,45 nitrates, and nitrites46 in processed red meat), in addition to protein per se, may have a critical health effect.
Interestingly, in this study, we found that the association of animal and plant protein with mortality varied by lifestyle factors, and any statistically significant protein-mortality associations were restricted to participants with at least 1 of the unhealthy behaviors, including smoking, heavy alcohol intake, overweight or obesity, and physical inactivity. Several reasons may explain these findings. First, given the remaining variation of health behaviors across protein intake categories in the unhealthy-lifestyle group, residual confounding from lifestyle factors may contribute to the observed protein-mortality associations. However, our results are robust to adjustment for a wide spectrum of potential confounders and the propensity score. Second, our results suggest that the adverse effects of high animal protein intake and beneficial effects of plant protein may be enhanced by other unhealthy lifestyle choices and become evident among the subgroup of individuals with these behaviors who may already have had some underlying inflammatory or metabolic disorders. Finally, as shown in Table 3, participants with a similar intake and with and without a healthy lifestyle demonstrated distinct profiles of protein sources. Those with unhealthy lifestyles consumed more processed and unprocessed red meat, whereas the healthy-lifestyle group consumed more fish and chicken as animal protein sources, suggesting that different protein sources, at least in part, contributed to the observed variation in the protein-mortality associations according to lifestyle factors. This hypothesis is supported by our substitution analysis results. Although substituting plant foods for various animal foods was associated with a lower mortality, red meat, especially processed red meat, showed a much stronger association than fish and poultry, which themselves were not associated with mortality (eTable 6 in the Supplement). In fact, protein from certain fish, such as cod, has been suggested to improve the lipid profile, glycemic control, and insulin sensitivity.35,47,48
The strengths of the present study included the large sample size, repeated dietary assessments, and high follow-up rate of the 2 well-established cohorts for up to 32 years. Moreover, we collected detailed data on a wide spectrum of lifestyle factors that allowed for rigorous confounding adjustment and subgroup analysis. In addition, to facilitate public health recommendations, we calculated protein intake according to food sources and assessed the substitution effect for protein of various origins.
Quiz Ref IDA limitation of the study is the moderately higher protein consumption (median, 19% of calories) in our study population compared with the general US population (15%-16%),49,50 thus limiting our ability to assess the effect of the very low end of intake. Furthermore, as an observational study, residual confounding could not be excluded. However, our results are robust to the multivariable adjustment and propensity score analysis, and any confounding effect may have been minimized in our stratified analysis according to lifestyle profile.
Although higher intake of animal protein was associated with higher cardiovascular mortality and higher intake of plant protein was associated with lower mortality, these associations were confined to participants with at least 1 lifestyle risk factor. Substitution of plant protein for animal protein, especially from processed red meat, may confer a substantial health benefit. Therefore, public health recommendations should focus on improvement of protein sources.
Correction: This article was corrected on October 3, 2016, to fix the Abstract and text.
Corresponding Author: Mingyang Song, MD, ScD, Clinical and Translational Epidemiology Unit, Division of Gastroenterology, Massachusetts General Hospital and Harvard Medical School Bartlett Hall Extension, Room 906, 55 Fruit St, Boston, MA 02114 (email@example.com).
Accepted for Publication: June 7, 2016.
Published Online: August 1, 2016. doi:10.1001/jamainternmed.2016.4182
Author Contributions: Drs Song and Giovannucci had full access to all the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis.
Study concept and design: Song, Hu, Willett, Longo, Giovannucci.
Acquisition, analysis, or interpretation of data: Song, Fung, Hu, Chan, Giovannucci.
Drafting of the manuscript: Song, Chan.
Critical revision of the manuscript for important intellectual content: All authors.
Statistical analysis: Song, Fung, Giovannucci.
Obtained funding: Chan, Giovannucci.
Administrative, technical, or material support: Hu, Willett, Chan.
Study supervision: Chan, Giovannucci.
Conflict of Interest Disclosures: None reported.
Funding/Support: This study was supported by the grants UM1 CA186107, P01 CA87969, and UM1 CA167552 from the National Institutes of Health.
Role of the Funder/Sponsor: The funding sources had no role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication.
Additional Information: Dr Song is a trainee of the Harvard Transdisciplinary Research Center on Energetics and Cancer.
Additional Contributions: We thank the participants and staff of the Nurses’ Health Study and the Health Professionals Follow-up Study for their valuable contributions and the following state cancer registries for their help: Alabama, Arizona, Arkansas, California, Colorado, Connecticut, Delaware, Florida, Georgia, Idaho, Illinois, Indiana, Iowa, Kentucky, Louisiana, Maine, Maryland, Massachusetts, Mississippi, Nebraska, New Hampshire, New Jersey, New York, North Carolina, North Dakota, Ohio, Oklahoma, Oregon, Pennsylvania, Rhode Island, South Carolina, Tennessee, Texas, Virgina, Washington, and Wyoming.