PubMed Health⌕ Search

Biomedical subjects

Nicholas G Martin

Publications and source records attributed to Nicholas G Martin.

At least 19 recordsLinked to original sources

Robust inference and correlates from genetic associations with personality.

Personality traits describe stable differences in how people think, feel and behave, and how they interact with and experience their social and physical environments1,2. Many questions remain unanswered about associations between DNA and personality traits, such as their robustness, their generalizability and the biological and social pathways through which they act. Here we meta-analyse data across 46 cohorts comprising 611,037 to 1.14 million participants with European-like and African-like genomes for genome-wide association studies (GWAS) of the Big Five personality traits (extraversion, agreeableness, conscientiousness, neuroticism and openness to experience), and data from up to 50,725 participants for within-family GWAS. We identify 1,260 lead genetic variants associated with personality, including 824 novel variants3. Common genetic variants explain a moderate 4.8-9.3% of the variance in measures of each trait, and 9.3-13.3% among instruments with typical measurement reliability. Genetic associations with personality are highly consistent but not identical across geography, reporter (self versus close other), age group and measurement instrument, and we find minimal spousal assortment for personality in recent history. In contrast to many other social and behavioural traits4,5, within-family GWAS and polygenic index analyses indicate that genetic associations with personality are minimally confounded by the shared family environment. Polygenic prediction, genetic correlation and Mendelian randomization analyses indicate that personality traits have widespread, potentially causal associations with consequential behaviours and life outcomes. Overall, we find that the genetic architecture of personality is robustly generalizable, minimally confounded and widely relevant to human experience.

Journal Article↗

Unravelling sex differences in the genetic architecture of anxiety.

BACKGROUND: Anxiety disorders show striking sex differences in prevalence, symptoms, and clinical characteristics, shaping how they manifest and are experienced. METHODS: Here, we report the first sex-specific meta-analysis of genome-wide association studies (GWAS) of anxiety, leveraging two of the largest biobank datasets, UK Biobank and All of Us, comprising 85,042 female cases with 196,789 controls and 36,732 male cases with 136,924 controls. Functional annotation, sex-specific polygenic scores (PGS), and genetic correlations were performed to assess genetic differences and functional implications. RESULTS: In females, 21 lead SNPs were significantly associated with anxiety, compared to five in males. Although the genetic correlation between sexes was high, it was significantly different from one, indicating partially distinct genetic architectures. In addition, both the SNP-based observed and liability-scale heritabilities (assuming a 2:1 female-to-male prevalence ratio) were significantly higher in females. Gene-based tests and functional prioritization identified different genes associated with anxiety in females and males. Moreover, genetic correlation analyses revealed stronger associations of female anxiety with attention-deficit/hyperactivity disorder (ADHD) and body mass index (BMI), whereas male anxiety showed stronger correlations with waist-hip-ratio-adjusted BMI. CONCLUSIONS: While the overall genetic architecture of anxiety is largely shared, our findings reveal distinct sex-specific genetic associations and correlations, highlighting the value of analyzing the sexes separately to uncover genetic signals that may be masked in sex-combined samples.

Female↗

Genetics of major depressive disorder in a homogeneous population with uniform phenotyping.

Harmonized phenotyping and diverse population-specific studies are crucial for advancing gene discovery in psychiatric genetics. We conducted a genome-wide association (GWAS) mega-analysis of DSM-defined lifetime major depressive disorder (MDD) in 64 941 participants (25.7% cases) from the Dutch BIObanks Netherlands Internet Collaboration (BIONIC) consortium. Liability-scale SNP-based heritability was 12.0% (SE = 1.4%) as estimated by LDSC (assuming a lifetime prevalence of 15%) and 26.6% (SE = 1.1%) when estimated by LDAK-REML on individual-level genotype data, indicating substantial common-variant signal in this clinically harmonized sample. The genetic correlation with the latest major depression GWAS from the Psychiatric Genomics Consortium (PGC-MD) was high (rG = 0.89, SE = 0.048). Polygenic scores (PGSs) based on BIONIC predicted depression in UK Biobank, and PGSs derived from PGC-MD predicted MDD in BIONIC, supporting transferability of depression polygenic signal across cohorts and phenotype definitions. Within-family PGS analyses in twins suggested that the observed prediction was not primarily driven by detectable family-level confounding, and twin concordance for MDD increased with polygenic burden. We identified one genome-wide significant locus, indexed by rs3818852 in PALMD, but this finding currently lacks independent replication and should be interpreted cautiously. Finally, genetic correlation and latent causal variable analyses identified multiple traits showing shared or directionally consistent genetic associations with MDD. Together, these findings underscore the value of clinically harmonized phenotyping in regional biobank collaborations for studying the genetic architecture of MDD.

Humans↗

Genome-wide association meta-regression identifies stem cell lineage orchestration as a key driver of acne risk.

Over 85% of the population experience acne at some point in their lives, with its severity spanning a quantitative spectrum, from mild, transient outbreaks to more persistent, severe forms of the condition. Moderate to severe disease poses a substantial global burden arising from both the physical and psychological impacts of this highly visible condition. The analytical approach taken in this study aimed to address the impact of variation in the dichotomisation of acne case control status, driven by ascertainment and study design, on effect size estimates across independent genetic association studies of acne. Through a fixed intercept meta-regression framework, we combined evidence genome-wide for association with acne across studies in which case-control status had been ascertained in different settings, allowing for different severity threshold definitions. Across a combined sample of 73,997 cases and 1,103,940 controls of European, South Asian and African American ancestry we identify genetic variation at 165 genomic loci that influence acne risk. There is evidence for both shared and ancestry specific components to the genetic susceptibility to acne and for sex differences in the magnitude of effect of risk alleles at three loci. We observe that common genetic variation explains 13.4% of acne heritability on the liability scale. Consistent with the hypothesis that genetic risk primarily operates at the level of individual pilosebaceous units, a polygenic score derived from this case-control study of acne susceptibility is associated with both self-reported and clinically assessed acne severity in adolescence, further strengthening the link between genetic risk and disease severity. Prioritisation of causal genes at the identified acne risk loci, provides genetic validation of the targets of established and emerging acne therapies, including retinoid treatments. The identified acne risk loci are enriched for genes encoding downstream effectors of RXRA signalling, including SOX9 and components of the WNT and p53 pathways. Illustrating that the control of stem cell lineage plasticity and cellular fate are important mechanisms through which genetic variation influences acne susceptibility within the pilosebaceous unit.

Journal Article↗

Genome-wide association study of adolescent-onset depression.

Adolescent depression is a heritable psychiatric condition with rising global prevalence and severe long-term outcomes, yet its biological underpinnings remain poorly understood. We conducted the first genome-wide association study of adolescent-onset depression, comprising 102,428 cases (diagnosis or clinical symptom thresholds) and 286,911 controls, including diverse ancestries. Cross-ancestry meta-analysis identified 52 independent variants across 17 loci; European-only analysis found 61 variants at 29 loci, with a SNP-based heritability of 9.8%. Comparative analyses revealed two genes unique to adolescent-onset versus lifetime depression, enriched in neuronal subtypes, and two genes as potential drug repurposing targets. Polygenic scores were associated with adolescent-onset depression across ancestries, persistent depression trajectories, more severe outcomes, as well as reduced cortical volume, surface area and white matter integrity. Genetic correlation and Mendelian randomisation analyses support shared genetic liability and causal links with early puberty and modifiable health and behavioural risk factors. These findings uncover novel genetic loci and refine biological pathways underlying adolescent-onset depression, revealing age-specific mechanisms and early intervention opportunities.

Journal Article↗

Genome-wide copy number variation association study in anorexia nervosa.

This study represents the first large-scale investigation of rare (<1% population frequency) copy number variants (CNVs) in anorexia nervosa (AN). Large, rare CNVs are reported to be causally associated with anthropometric traits, neurodevelopmental disorders, and schizophrenia, yet their role in the genetic basis of AN is unclear. Using genome-wide association study (GWAS) array data from the Anorexia Nervosa Genetics Initiative (ANGI), which included 7414 AN case and 5044 controls, we investigated the association of 67 well-established syndromic CNVs and 178 pleiotropic disease-risk dosage-sensitive CNVs with AN. To identify novel CNV regions (CNVRs) that increase the risk of AN, we conducted genome-wide association studies with a focus on rare CNV-breakpoints (CNV-GWAS). We found no net enrichment of rare CNVs, either deletions or duplications, in AN, and none of the well-established syndromic or pleiotropic CNVs had a significant association with AN status. However, the CNV-GWAS found 21 nominally associated CNVRs that contribute to AN risk, covering protein-coding genes implicated in synaptic function, metabolic/mitochondrial factors, and lipid characteristics, like the CD36 (7q21.11) gene, which transports long-chain fatty acids into cells. CNVRs intersecting genes previously related to neurodevelopmental traits include deletions of NRXN1 intron 5 (2p16.3), IMMP2L (7q31.1), and PTPRD (9p23). Overall, given that our study is well powered to detect the CNV burden level reported for schizophrenia, we can conclude that rare CNVs have a limited role in the etiology of AN, as reported for bipolar disorder. Our nominal associations for the 21 discovered CNVRs are consistent with AN being a metabo-psychiatric trait, as demonstrated by the common genetic architecture of AN, and we provide association results to allow for replication in future research.

Humans↗

1q21.1 distal copy number variants are associated with cerebral and cognitive alterations in humans.

Low-frequency 1q21.1 distal deletion and duplication copy number variant (CNV) carriers are predisposed to multiple neurodevelopmental disorders, including schizophrenia, autism and intellectual disability. Human carriers display a high prevalence of micro- and macrocephaly in deletion and duplication carriers, respectively. The underlying brain structural diversity remains largely unknown. We systematically called CNVs in 38 cohorts from the large-scale ENIGMA-CNV collaboration and the UK Biobank and identified 28 1q21.1 distal deletion and 22 duplication carriers and 37,088 non-carriers (48% male) derived from 15 distinct magnetic resonance imaging scanner sites. With standardized methods, we compared subcortical and cortical brain measures (all) and cognitive performance (UK Biobank only) between carrier groups also testing for mediation of brain structure on cognition. We identified positive dosage effects of copy number on intracranial volume (ICV) and total cortical surface area, with the largest effects in frontal and cingulate cortices, and negative dosage effects on caudate and hippocampal volumes. The carriers displayed distinct cognitive deficit profiles in cognitive tasks from the UK Biobank with intermediate decreases in duplication carriers and somewhat larger in deletion carriers-the latter potentially mediated by ICV or cortical surface area. These results shed light on pathobiological mechanisms of neurodevelopmental disorders, by demonstrating gene dose effect on specific brain structures and effect on cognitive function.

Brain↗

The genetics of voting: an Australian twin study.

Previously we and others have shown evidence for genetic influences on political attitudes and sociodemographic indicators (Martin 1987; Posner et al. 1996; Truett et al. 1992; Eaves et al. 1999). However, the nature of the relationship between political attitudes, social indictors and voting behavior has not been investigated. While heritability estimates for social and political attitudes have been reported in previous research, the heritability for vote choice has not. Furthermore, if vote choice is heritable, it is unclear whether the heritable component can be accounted for through the genetic influence on related social and political traits, or if there exists a unique genetic component specific to voting behavior. In mailed surveys of adult Australian twins, we asked respondents to indicate their usual voting preference as well as attitudes on contemporary individual political items. When vote choice was dichotomized as Labor versus Conservative, twin correlations were r (mz) = 0.81 (1,661 pairs), and r (dz) = 0.69 (1,727 pairs) consistent with modest genetic influence (a (2) = 0.24). However, multivariate genetic analysis showed no unique genetic contribution to voting preference; rather, the genetic influence in vote choice could be explained by shared genetic influences in perceived social class, church attendance and certain key political attitude items.

Adult↗

The ongoing adaptive evolution of ASPM and Microcephalin is not explained by increased intelligence.

Recent studies have made great strides towards identifying putative genetic events underlying the evolution of the human brain and its emergent cognitive capacities. One of the most intriguing findings is the recurrent identification of adaptive evolution in genes associated with primary microcephaly, a developmental disorder characterized by severe reduction in brain size and intelligence, reminiscent of the early hominid condition. This has led to the hypothesis that the adaptive evolution of these genes has contributed to the emergence of modern human cognition. As with other candidate loci, however, this hypothesis remains speculative due to the current lack of methodologies for characterizing the evolutionary function of these genes in humans. Two primary microcephaly genes, ASPM and Microcephalin, have been implicated not only in the adaptive evolution of the lineage leading to humans, but in ongoing selective sweeps in modern humans as well. The presence of both the putatively adaptive and neutral alleles at these loci provides a unique opportunity for using normal trait variation within humans to test the hypothesis that the recent selective sweeps are driven by an advantage in cognitive abilities. Here, we report a large-scale association study between the adaptive alleles of these genes and normal variation in several measures of IQ. Five independent samples were used, totaling 2393 subjects, including both family-based and population-based datasets. Our overall findings do not support a detectable association between the recent adaptive evolution of either ASPM or Microcephalin and changes in IQ. As we enter the post-genomic era, with the number of candidate loci underlying human evolution growing rapidly, our findings highlight the importance of direct experimental validation in elucidating their evolutionary role in shaping the human phenotype.

Adolescent↗

Genetic influences on exercise participation in 37,051 twin pairs from seven countries.

BACKGROUND: A sedentary lifestyle remains a major threat to health in contemporary societies. To get more insight in the relative contribution of genetic and environmental influences on individual differences in exercise participation, twin samples from seven countries participating in the GenomEUtwin project were used. METHODOLOGY: Self-reported data on leisure time exercise behavior from Australia, Denmark, Finland, Norway, The Netherlands, Sweden and United Kingdom were used to create a comparable index of exercise participation in each country (60 minutes weekly at a minimum intensity of four metabolic equivalents). PRINCIPAL FINDINGS: Modest geographical variation in exercise participation was revealed in 85,198 subjects, aged 19-40 years. Modeling of monozygotic and dizygotic twin resemblance showed that genetic effects play an important role in explaining individual differences in exercise participation in each country. Shared environmental effects played no role except for Norwegian males. Heritability of exercise participation in males and females was similar and ranged from 48% to 71% (excluding Norwegian males). CONCLUSIONS: Genetic variation is important in individual exercise behavior and may involve genes influencing the acute mood effects of exercise, high exercise ability, high weight loss ability, and personality. This collaborative study suggests that attempts to find genes influencing exercise participation can pool exercise data across multiple countries and different instruments.

Adult↗

Replicated effects of sex and genotype on gene expression in human lymphoblastoid cell lines.

The expression level for 15,887 transcripts in lymphoblastoid cell lines from 19 monozygotic twin pairs (10 male, 9 female) were analysed for the effects of genotype and sex. On an average, the effect of twin pairs explained 31% of the variance in normalized gene expression levels, consistent with previous broad sense heritability estimates. The effect of sex on gene expression levels was most noticeable on the X chromosome, which contained 15 of the 20 significantly differentially expressed genes. A high concordance was observed between the sex difference test statistics and surveys of genes escaping X chromosome inactivation. Notably, several autosomal genes showed significant differences in gene expression between the sexes despite much of the cellular environment differences being effectively removed in the cell lines. A publicly available gene expression data set from the CEPH families was used to validate the results. The heritability of gene expression levels as estimated from the two data sets showed a highly significant positive correlation, particularly when both estimates were close to one and thus had the smallest standard error. There was a large concordance between the genes significantly differentially expressed between the sexes in the two data sets. Analysis of the variability of probe binding intensities within a probe set indicated that results are robust to the possible presence of polymorphisms in the target sequences.

Adult↗

Significant evidence of one or more susceptibility loci for endometriosis with near-Mendelian inheritance on chromosome 7p13-15.

BACKGROUND: Endometriosis is a common disease with a heritable component. The collaborative International Endogene Study consists of two data sets (Oxford and Australia) comprising 1176 families with multiple affected. The aim was to investigate whether the apparent concentration of cases in a proportion of families could be explained by one or more rare variants with (near-)Mendelian autosomal inheritance. METHODS AND RESULTS: Linkage analyses (aimed at finding chromosomal regions harbouring disease-predisposing genes) were conducted in families with three or more affected (Oxford: n = 52; Australia: n = 196). In the Oxford data set, a non-parametric linkage score (Kong & Cox (K&C) Log of ODds (LOD)) of 3.52 was observed on chromosome 7p (genome-wide significance P = 0.011). A parametric MOD score (equal to maximum LOD maximized over 357 possible inheritance models) of 3.89 was found at 65.72 cM (D7S510) for a dominant model with reduced penetrance. After including the Australian data set, the non-parametric K&C LOD of the combined data set was 1.46 at 57.3 cM; the parametric analysis found an MOD score of 3.30 at D7S484 (empirical significance: P = 0.035) for a recessive model with high penetrance. Critical recombinant analysis narrowed the probable region of linkage down to overlapping 6.4 Mb and 11 Mb intervals containing 48 and 96 genes, respectively. CONCLUSIONS: This is the first report to suggest that there may be one or more high-penetrance susceptibility loci for endometriosis with (near-)Mendelian inheritance.

Australia↗

Novel genes identified in a high-density genome wide association study for nicotine dependence.

Tobacco use is a leading contributor to disability and death worldwide, and genetic factors contribute in part to the development of nicotine dependence. To identify novel genes for which natural variation contributes to the development of nicotine dependence, we performed a comprehensive genome wide association study using nicotine dependent smokers as cases and non-dependent smokers as controls. To allow the efficient, rapid, and cost effective screen of the genome, the study was carried out using a two-stage design. In the first stage, genotyping of over 2.4 million single nucleotide polymorphisms (SNPs) was completed in case and control pools. In the second stage, we selected SNPs for individual genotyping based on the most significant allele frequency differences between cases and controls from the pooled results. Individual genotyping was performed in 1050 cases and 879 controls using 31 960 selected SNPs. The primary analysis, a logistic regression model with covariates of age, gender, genotype and gender by genotype interaction, identified 35 SNPs with P-values less than 10(-4) (minimum P-value 1.53 x 10(-6)). Although none of the individual findings is statistically significant after correcting for multiple tests, additional statistical analyses support the existence of true findings in this group. Our study nominates several novel genes, such as Neurexin 1 (NRXN1), in the development of nicotine dependence while also identifying a known candidate gene, the beta3 nicotinic cholinergic receptor. This work anticipates the future directions of large-scale genome wide association studies with state-of-the-art methodological approaches and sharing of data with the scientific community.

Case-Control Studies↗

Cholinergic nicotinic receptor genes implicated in a nicotine dependence association study targeting 348 candidate genes with 3713 SNPs.

Nicotine dependence is one of the world's leading causes of preventable death. To discover genetic variants that influence risk for nicotine dependence, we targeted over 300 candidate genes and analyzed 3713 single nucleotide polymorphisms (SNPs) in 1050 cases and 879 controls. The Fagerström test for nicotine dependence (FTND) was used to assess dependence, in which cases were required to have an FTND of 4 or more. The control criterion was strict: control subjects must have smoked at least 100 cigarettes in their lifetimes and had an FTND of 0 during the heaviest period of smoking. After correcting for multiple testing by controlling the false discovery rate, several cholinergic nicotinic receptor genes dominated the top signals. The strongest association was from an SNP representing CHRNB3, the beta3 nicotinic receptor subunit gene (P = 9.4 x 10(-5)). Biologically, the most compelling evidence for a risk variant came from a non-synonymous SNP in the alpha5 nicotinic receptor subunit gene CHRNA5 (P = 6.4 x 10(-4)). This SNP exhibited evidence of a recessive mode of inheritance, resulting in individuals having a 2-fold increase in risk of developing nicotine dependence once exposed to cigarette smoking. Other genes among the top signals were KCNJ6 and GABRA4. This study represents one of the most powerful and extensive studies of nicotine dependence to date and has found novel risk loci that require confirmation by replication studies.

Adult↗

Replication of reported linkages for dyslexia and spelling and suggestive evidence for novel regions on chromosomes 4 and 17.

We report the first genome-wide linkage analysis for reading and spelling in a sample of 403 families of twins, aged between 12 and 25 years taken from the normal population and unselected for reading ability. These traits showed heritabilities of 0.52-0.73, and support for linkage exceeded replication levels (lod > 1.44) of seven of the 11 linkages reported in dyslexic samples, namely: 2q22.3, 3p12-q13, 6q11.2, 7q32, 15q21.1, 18p21, and Xq27.3. For five of these (2q22.3, 6q11.2, 7q32, 18p21, and Xq27), this study provides the first independent replication. 1p34-36 and 2p15-16 received some support, with lods of 1.2 and 0.83, respectively, whereas two regions received little support (6p23-21.3 and 11p15.5). This study also identified two novel linkages at 4p15.33-16.1 and 17p13.3, which received suggestive support (max. lod 2.08 and 1.99, respectively).

Chromosome Mapping↗

A genome-wide scan for naevus count: linkage to CDKN2A and to other chromosome regions.

High numbers of melanocytic naevi (moles), and mutations in the p16 gene (CDKN2A), are two strong risk factors for cutaneous malignant melanoma. We have previously reported linkage of mole count to the CDKN2A locus. Here, we report genome-wide scans for mole counts (differentiated into flat, raised and atypical subtypes) with a total of 796 microsatellite markers for 424 families with 1024 twins and siblings, plus genotypes for 690 parents. Inclusion of 221 pairs of MZ twins enabled separation of shared environmental and polygenic influences, so placing an upper limit to estimates of QTL variance. Maximum likelihood multipoint variance component methods were used to assess linkage of naevus count. Sex, age, body surface area, skin colour, hair colour, sunburn and facial freckles were included as covariates. Peak linkage of flat mole count was to regions on chromosomes 2, 9, 8 and 17 with lod scores 2.95, 2.95, 2.50 and 2.15, respectively. The support for linkage to the CDKN2A gene region (9p21) increased to 3.42 when additional fine mapping markers were added. For raised mole count, there was suggestive evidence of linkage in our sample to chromosome 16 (lod=1.87), and for atypical mole count on chromosomes 1, 6 and X with lod scores of 2.20, 2.00 and 2.00, respectively. The multivariate linkage peaks generally match those from individual trait analyses, with the exception of a new peak on chromosome 4 (point-wise empirical P-value=0.001). We replicate our earlier finding of linkage to CDKN2A and discovering linkage to several novel regions that may also influence risk of the development of malignant melanoma.

Adolescent↗

HLA and genomewide allele sharing in dizygotic twins.

Gametic selection during fertilization or the effects of specific genotypes on the viability of embryos may cause a skewed transmission of chromosomes to surviving offspring. A recent analysis of transmission distortion in humans reported significant excess sharing among full siblings. Dizygotic (DZ) twin pairs are a special case of the simultaneous survival of two genotypes, and there have been reports of DZ pairs with excess allele sharing around the HLA locus, a candidate locus for embryo survival. We performed an allele-sharing study of 1,592 DZ twin pairs from two independent Australian cohorts, of which 1,561 pairs were informative for linkage on chromosome 6. We also analyzed allele sharing in 336 DZ twin pairs from The Netherlands. We found no evidence of excess allele sharing, either at the HLA locus or in the rest of the genome. In contrast, we found evidence of a small but significant (P=.003 for the Australian sample) genomewide deficit in the proportion of two alleles shared identical by descent among DZ twin pairs. We reconciled conflicting evidence in the literature for excess genomewide allele sharing by performing a simulation study that shows how undetected genotyping errors can lead to an apparent deficit or excess of allele sharing among sibling pairs, dependent on whether parental genotypes are known. Our results imply that gene-mapping studies based on affected sibling pairs that include DZ pairs will not suffer from false-positive results due to loci involved in embryo survival.

Adolescent↗

Cohort trends in prevalence and spousal concordance for smoking.

Evaluation of the pattern of smoking prevalence over time could provide insight for disease prevention and tobacco control policy. Examination of assortative mating is important in assessing the relative contribution of genetic and environmental influences in behavior genetics. A significant spousal concordance for smoking could imply assortative mating behavior and have potential influence on the risk for offspring smoking. Using intergenerational data from large scale adult twin panels in Australia, the present study examined the patterns of prevalence and spousal concordance of lifetime smoking across a wide range of birth cohorts. Data were collected for more than 16,000 twins, their spouses, parents, and parents-reported grandparents. The measurement models were used to calculate prevalence and spousal concordance for smoking. Prevalence of lifetime smoking for females increased significantly over successive cohorts, and was higher for twins (0.28-0.49) than for their parents or grandparents (0.06-0.42). Smoking prevalence was still higher in males than in females for the twins (0.49-0.69), and parental or grandparental (0.62-0.77) generations, but the gender differences largely decreased in younger cohorts. Moderate but significant spousal correlations were found across cohorts and generations, which were higher for twins (0.39-0.55) than for their parents and grandparents (0.19-0.49). Despite using a simplified smoking assessment, findings in this study demonstrated changes of lifetime smoking patterns across birth cohorts and the presence of assortative mating for smoking behavior.

Australia↗