PMERGE: Computational filtering of paralogous sequences from RAD‐seq data

Restriction‐site associated DNA sequencing (RAD‐seq) can identify and score thousands of genetic markers from a group of samples for population‐genetics studies. One challenge of de novo RAD‐seq analysis is to distinguish paralogous sequence variants (PSVs) from true single‐nucleotide polymorphisms (SNPs) associated with orthologous loci. In the absence of a reference genome, it is difficult to differentiate true SNPs from PSVs, and their impact on downstream analysis remains unclear. Here, we introduce a network‐based approach, PMERGE that connects fragments based on their DNA sequence similarity to identify probable PSVs. Applying our method to de novo RAD‐seq data from 150 Atlantic salmon (Salmo salar) samples collected from 15 locations across the Southern Newfoundland coast allowed the identification of 87% of total PSVs identified through alignment to the Atlantic salmon genome. Removal of these paralogs altered the inferred population structure, highlighting the potential impact of filtering in RAD‐seq analysis. PMERGE is also applied to a green crab (Carcinus maenas) data set consisting of 242 samples from 11 different locations and was successfully able to identify and remove the majority of paralogous loci (62%). The PMERGE software can be run as part of the widely used Stacks analysis package.     

New view of population genetics of zooplankton: RAD‐seq analysis reveals population structure of the North Atlantic planktonic copepod Centropages typicus

Detection of population genetic structure of zooplankton at medium‐to‐small spatial scales in the absence of physical barriers has remained challenging and controversial. The large population sizes and high rates of gene flow characteristic of zooplankton have made resolution of geographical differentiation very difficult, especially when using few genetic markers and assuming equilibrium conditions. Next‐generation sequencing now allows simultaneous sampling of hundreds to thousands of genetic markers; new analytical approaches allow studies under nonequilibrium conditions and directional migration. Samples of the North Atlantic Ocean planktonic copepod, Centropages typicus, were analysed using restriction site‐associated DNA (RAD) sequencing on a PROTON platform. Although prior studies revealed no genetic differentiation of populations across the geographical range of the species, analysis of RAD tags showed significant structure across the North Atlantic Ocean. We also compared the likelihood for models of connectivity among NW Atlantic populations under various directional flow scenarios that replicate oceanographic conditions of the sampled domain. High‐density marker sampling with RAD sequencing markedly outperformed other technical and analytical approaches in detection of population genetic structure and characterization of connectivity of this high geneflow zooplankton species.     

Population genomic analyses from low‐coverage RAD‐Seq data: a case study on the non‐model cucurbit bottle gourd

Restriction site‐associated DNA sequencing (RAD‐Seq), a next‐generation sequencing‐based genome ‘complexity reduction’ protocol, has been useful in population genomics in species with a reference genome. However, the application of this protocol to natural populations of genomically underinvestigated species, particularly under low‐to‐medium sequencing depth, has not been well justified. In this study, a Bayesian method was developed for calling genotypes from an F₂ population of bottle gourd [Lagenaria siceraria (Mol.) Standl.] to construct a high‐density genetic map. Low‐depth genome shotgun sequencing allowed the assembly of scaffolds/contigs comprising approximately 50% of the estimated genome, of which 922 were anchored for identifying syntenic regions between species. RAD‐Seq genotyping of a natural population comprising 80 accessions identified 3226 single nuclear polymorphisms (SNPs), based on which two sub‐gene pools were suggested for association with fruit shape. The two sub‐gene pools were moderately differentiated, as reflected by the Hudson's F_ST value of 0.14, and they represent regions on LG7 with strikingly elevated F_ST values. Seven‐fold reduction in heterozygosity and two times increase in LD (r²) were observed in the same region for the round‐fruited sub‐gene pool. Outlier test suggested the locus LX3405 on LG7 to be a candidate site under selection. Comparative genomic analysis revealed that the cucumber genome region syntenic to the high F_ST island on LG7 harbors an ortholog of the tomato fruit shape gene OVATE. Our results point to a bright future of applying RAD‐Seq to population genomic studies for non‐model species even under low‐to‐medium sequencing efforts. The genomic resources provide valuable information for cucurbit genome research.     

Fine‐scale population structure of Collichtys lucidus populations inferred from microsatellite markers

The spinyhead croaker Collichthys lucidus (Richardson) is a small sciaenid species distributed along the inshore waters of northwestern Pacific Ocean, and now has been listed as Key Protected Commercial Sources of Aquatic Animals and Plants in China. To delineate stock boundaries and inform conservation policy for its management, samples were collected from eight locations across the Chinese coastal waters and analyzed at nine microsatellite loci. C. lucidus populations showed low genetic diversity (expected heterozygosity = 0.445–0.542; observed heterozygosity = 0.392–0.539; Polymorphism Information Content = 0.268–0.684). Strong genetic fdifferentiation (F_st = 0.065–0.510, all significant after Bonferroni correction) among all populations and high levels of self‐recruitment (89.2%–91.5%) were observed, which suggested limited genetic exchange for this species. Clustering results of discriminant analysis of principal components and STRUCTURE found strong support for obvious genetic clusters (populations FZ, XM and SZ vs. populations SH, YRE, ZS, WZ and ND). The results of the present study not only supported the phylogeographic pattern of north‐south differentiation, but also suggested that C. lucidus populations may be predominantly sustained by self‐replenishment rather than by recruitment from distant populations.     

Development and preliminary evaluation of a genomewide single nucleotide polymorphisms resource generated by RAD‐seq for the small yellow croaker (Larimichthys polyactis)

Recent advances in high‐throughput sequencing technologies have offered the possibility to generate genomewide sequence data to delineate previously unidentified genetic structure, obtain more accurate estimates of demographic parameters and to evaluate potential adaptive divergence. Here, we identified 27 556 single nucleotide polymorphisms for the small yellow croaker (Larimichthys polyactis) using restriction‐site‐associated DNA (RAD) sequencing of 24 individuals from two populations. Significant sources of genetic variation were identified, with an average nucleotide diversity (π) of 0.00105 ± 0.000425 across individuals, and long‐term effective population size was thus estimated to range between 26 172 and 261 716. According to the results, no differentiation between the two populations was detected based on the SNP data set of top quality score per contig or neutral loci. However, the two analysed populations were highly differentiated based on SNP data set of both top F_ST value per contig and the outlier SNPs. Moreover, local adaptation was highlighted by an F_ST‐based outlier tests implemented in LOSITAN and a total of 538 potentially locally selected SNPs were identified. blast2go annotation of contigs containing the outlier SNPs yielded hits for 37 (66%) of 56 significant blastx matches. Candidate genes for local adaptation constituted a wide array of biological functions, including cellular response to oxidative stress, actin filament binding, ion transmembrane transport and synapse assembly. The generated SNP resources in this study provided a valuable tool for future population genetics and genomics studies of L. polyactis.     

Identification of sex‐specific molecular markers using restriction site‐associated DNA sequencing

A major barrier to evolutionary studies of sex determination and sex chromosomes has been a lack of information on the types of sex‐determining mechanisms that occur among different species. This is particularly problematic in groups where most species lack visually heteromorphic sex chromosomes, such as fish, amphibians and reptiles, because cytogenetic analyses will fail to identify the sex chromosomes in these species. We describe the use of restriction site‐associated DNA (RAD) sequencing, or RAD‐seq, to identify sex‐specific molecular markers and subsequently determine whether a species has male or female heterogamety. To test the accuracy of this technique, we examined the lizard Anolis carolinensis. We performed RAD‐seq on seven male and ten female A. carolinensis and found one male‐specific molecular marker. Anolis carolinensis has previously been shown to possess male heterogamety and the recently published A. carolinensis genome facilitated the characterization of the sex‐specific RAD‐seq marker. We validated the male specificity of the new marker using PCR on additional individuals and also found that it is conserved in some other Anolis species. We discuss the utility of using RAD‐seq to identify sex‐determining mechanisms in other species with cryptic or homomorphic sex chromosomes and the implications for the evolution of male heterogamety in Anolis.     

Transatlantic secondary contact in Atlantic Salmon,comparing microsatellites,a single nucleotide polymorphism array and restriction‐site associated DNA sequencing for the resolution of complex spatial structure

Identification of discrete and unique assemblages of individuals or populations is central to the management of exploited species. Advances in population genomics provide new opportunities for re‐evaluating existing conservation units but comparisons among approaches remain rare. We compare the utility of RAD‐seq, a single nucleotide polymorphism (SNP) array and a microsatellite panel to resolve spatial structuring under a scenario of possible trans‐Atlantic secondary contact in a threatened Atlantic Salmon, Salmo salar, population in southern Newfoundland. Bayesian clustering indentified two large groups subdividing the existing conservation unit and multivariate analyses indicated significant similarity in spatial structuring among the three data sets. mtDNA alleles diagnostic for European ancestry displayed increased frequency in southeastern Newfoundland and were correlated with spatial structure in all marker types. Evidence consistent with introgression among these two groups was present in both SNP data sets but not the microsatellite data. Asymmetry in the degree of introgression was also apparent in SNP data sets with evidence of gene flow towards the east or European type. This work highlights the utility of RAD‐seq based approaches for the resolution of complex spatial patterns, resolves a region of trans‐Atlantic secondary contact in Atlantic Salmon in Newfoundland and demonstrates the utility of multiple marker comparisons in identifying dynamics of introgression.     

Strengths and potential pitfalls of hay transfer for ecological restoration revealed by RAD‐seq analysis in floodplain Arabis species

Achieving high intraspecific genetic diversity is a critical goal in ecological restoration as it increases the adaptive potential and long‐term resilience of populations. Thus, we investigated genetic diversity within and between pristine sites in a fossil floodplain and compared it to sites restored by hay transfer between 1997 and 2014. RAD‐seq genotyping revealed that the stenoecious floodplain species Arabis nemorensis is co‐occurring with individuals that, based on ploidy, ITS‐sequencing and morphology, probably belong to the close relative Arabis sagittata, which has a documented preference for dry calcareous grasslands but has not been reported in floodplain meadows. We show that hay transfer maintains genetic diversity for both species. Additionally, in A. sagittata, transfer from multiple genetically isolated pristine sites resulted in restored sites with increased diversity and admixed local genotypes. In A. nemorensis, transfer did not create novel admixture dynamics because genetic diversity between pristine sites was less differentiated. Thus, the effects of hay transfer on genetic diversity also depend on the genetic make‐up of the donor communities of each species, especially when local material is mixed. Our results demonstrate the efficiency of hay transfer for habitat restoration and emphasize the importance of prerestoration characterization of microgeographic patterns of intraspecific diversity of the community to guarantee that restoration practices reach their goal, that is maximize the adaptive potential of the entire restored plant community. Overlooking these patterns may alter the balance between species in the community. Additionally, our comparison of summary statistics obtained from de novo‐ and reference‐based RAD‐seq pipelines shows that the genomic impact of restoration can be reliably monitored in species lacking prior genomic knowledge.     

Population genomics reveals multiple drivers of population differentiation in a sex‐role‐reversed pipefish

A major goal of molecular ecology is to identify the causes of genetic and phenotypic differentiation among populations. Population genomics is suitably poised to tackle these key questions by diagnosing the evolutionary mechanisms driving divergence in nature. Here, we set out to investigate the evolutionary processes underlying population differentiation in the Gulf pipefish, Syngnathus scovelli. We sampled approximately 50 fish from each of 12 populations distributed from the Gulf coast of Texas to the Atlantic coast of Florida and performed restriction‐site‐associated DNA sequencing to identify SNPs throughout the genome. After imposing quality and stringency filters, we selected a panel of 6348 SNPs present in all 12 populations, 1753 of which were not physically linked. We identified a genome‐wide pattern of isolation by distance, in addition to a more substantial genetic break separating populations in the Gulf of Mexico from those in the Atlantic. We also used several divergence outlier approaches and tests for genotype–environment correlations to identify 400 SNPs putatively involved in local adaptation. Patterns of phenotypic differentiation and variation diverged from the overall genomic pattern, suggesting that selection, phenotypic plasticity or demographic factors may be shaping phenotypes in distinct populations. Overall, our results suggest that population divergence is driven by a variety of factors in S. scovelli, including neutral processes and selection on multiple traits.     

Genetic diversity of Chinese cattle revealed by Y‐SNP and Y‐STR markers

With its vast territory and complex natural environment, China boasts rich cattle genetic resources. To gain the further insight into the genetic diversity and paternal origins of Chinese cattle, we analyzed the polymorphism of Y‐SNPs (UTY19 and ZFY10) and Y‐STRs (INRA189 and BM861) in 34 Chinese cattle breeds/populations, including 606 males representative of 24 cattle breeds/populations collected in this study as well as previously published data for 302 bulls. Combined genotypic data identified 14 Y‐chromosome haplotypes that represented three haplogroups. Y2‐104‐158 and Y2‐102‐158 were the most common taurine haplotypes detected mainly in northern and central China, whereas the indicine haplotype Y3‐88‐156 predominates in southern China. Haplotypes Y2‐108‐158, Y2‐110‐158, Y2‐112‐158 and Y3‐92‐156 were private to Chinese cattle. The population structure revealed by multidimensional scaling analysis differentiated Tibetan cattle from the other three groups of cattle. Analysis of molecular variance showed that the majority of the genetic variation was explained by the genetic differences among groups. Overall, our study indicates that Chinese cattle retain high paternal diversity (H = 0.607 ± 0.016) and probably much of the original lineages that derived from the domestication center in the Near East without strong admixture from commercial cattle carrying Y1 haplotypes.     

RAD‐sequencing for estimating genomic relatedness matrix‐based heritability in the wild: A case study in roe deer

Estimating the evolutionary potential of quantitative traits and reliably predicting responses to selection in wild populations are important challenges in evolutionary biology. The genomic revolution has opened up opportunities for measuring relatedness among individuals with precision, enabling pedigree‐free estimation of trait heritabilities in wild populations. However, until now, most quantitative genetic studies based on a genomic relatedness matrix (GRM) have focused on long‐term monitored populations for which traditional pedigrees were also available, and have often had access to knowledge of genome sequence and variability. Here, we investigated the potential of RAD‐sequencing for estimating heritability in a free‐ranging roe deer (Capreolous capreolus) population for which no prior genomic resources were available. We propose a step‐by‐step analytical framework to optimize the quality and quantity of the genomic data and explore the impact of the single nucleotide polymorphism (SNP) calling and filtering processes on the GRM structure and GRM‐based heritability estimates. As expected, our results show that sequence coverage strongly affects the number of recovered loci, the genotyping error rate and the amount of missing data. Ultimately, this had little effect on heritability estimates and their standard errors, provided that the GRM was built from a minimum number of loci (above 7,000). Genomic relatedness matrix‐based heritability estimates thus appear robust to a moderate level of genotyping errors in the SNP data set. We also showed that quality filters, such as the removal of low‐frequency variants, affect the relatedness structure of the GRM, generating lower h² estimates. Our work illustrates the huge potential of RAD‐sequencing for estimating GRM‐based heritability in virtually any natural population.     

Population genetic structure of a common host predicts the spread of white‐nose syndrome,an emerging infectious disease in bats

Landscape complexity influences patterns of animal dispersal, which in turn may affect both gene flow and the spread of pathogens. White‐nose syndrome (WNS) is an introduced fungal disease that has spread rapidly throughout eastern North America, causing massive mortality in bat populations. We tested for a relationship between the population genetic structure of the most common host, the little brown myotis (Myotis lucifugus), and the geographic spread of WNS to date by evaluating logistic regression models of WNS risk among hibernating colonies in eastern North America. We hypothesized that risk of WNS to susceptible host colonies should increase with both geographic proximity and genetic similarity, reflecting historical connectivity, to infected colonies. Consistent with this hypothesis, inclusion of genetic distance between infected and susceptible colonies significantly improved models of disease spread, capturing heterogeneity in the spatial expansion of WNS despite low levels of genetic differentiation among eastern populations. Expanding our genetic analysis to the continental range of little brown myotis reveals strongly contrasting patterns of population structure between eastern and western North America. Genetic structure increases markedly moving westward into the northern Great Plains, beyond the current distribution of WNS. In western North America, genetic differentiation of geographically proximate populations often exceeds levels observed across the entire eastern region, suggesting infrequent and/or locally restricted dispersal, and thus relatively limited opportunities for pathogen introduction in western North America. Taken together, our analyses suggest a possibly slower future rate of spread of the WNS pathogen, at least as mediated by little brown myotis.     

The genetic component of human longevity: New insights from the analysis of pathway‐based SNP‐SNP interactions

In human longevity studies, single nucleotide polymorphism (SNP) analysis identified a large number of genetic variants with small effects, yet not easily replicable in different populations. New insights may come from the combined analysis of different SNPs, especially when grouped by metabolic pathway. We applied this approach to study the joint effect on longevity of SNPs belonging to three candidate pathways, the insulin/insulin‐like growth factor signalling (IIS), DNA repair and pro/antioxidant. We analysed data from 1,058 tagging SNPs in 140 genes, collected in 1825 subjects (1,089 unrelated nonagenarians from the Danish 1905 Birth Cohort Study and 736 Danish controls aged 46–55 years) for evaluating synergic interactions by SNPsyn. Synergies were further tested by the multidimensional reduction (MDR) approach, both intra‐ and interpathways. The best combinations (FDR<0.0001) resulted those encompassing IGF1R‐rs12437963 and PTPN1‐rs6067484, TP53‐rs2078486 and ERCC2‐rs50871, TXNRD1‐rs17202060 and TP53‐rs2078486, the latter two supporting a central role of TP53 in mediating the concerted activation of the DNA repair and pro‐antioxidant pathways in human longevity. Results were consistently replicated with both approaches, as well as a significant effect on longevity was found for the GHSR gene, which also interacts with partners belonging to both IIS and DNA repair pathways (PAPPA, PTPN1, PARK7, MRE11A). The combination GHSR‐MREA11, positively associated with longevity by MDR, was further found influencing longitudinal survival in nonagenarian females (p = .026). Results here presented highlight the validity of SNP‐SNP interactions analyses for investigating the genetics of human longevity, confirming previously identified markers but also pointing to novel genes as central nodes of additional networks involved in human longevity.     

Extending RAD tag analysis to microbial ecology: a comparison between MultiLocus Sequence Typing and 2b‐RAD to investigate Listeria monocytogenes genetic structure

The advent of next‐generation sequencing (NGS) has dramatically changed bacterial typing technologies, increasing our ability to differentiate bacterial isolates. Despite it is now possible to sequence a bacterial genome in a few days and at reasonable costs, most genetic analyses do not require whole‐genome sequencing, which also remains impractical for large population samples due to the cost of individual library preparation and bioinformatics. More traditional sequencing approaches, however, such as MultiLocus Sequence Typing (mlst ) are quite laborious and time‐consuming, especially for large‐scale analyses. In this study, a genotyping approach based on restriction site‐associated (RAD) tag sequencing, 2b‐RAD, was applied to characterize Listeria monocytogenes strains. To verify the feasibility of the method, an in silico analysis was performed on 30 available complete genomes. For the same set of strains, in silico mlst analysis was conducted as well. Subsequently, 2b‐RAD and mlst analyses were experimentally carried out on 58 isolates collected from food samples or food‐processing sites. The obtained results demonstrate that 2b‐RAD predicts mlst types and often provides more detailed information on population structure than mlst . Moreover, the majority of variants differentiating identical sequence type isolates mapped against accessory fragments, thus providing additional information to characterize strains. Although mlst still represents a reliable typing method, large‐scale studies on molecular epidemiology and public health, as well as bacterial phylogenetics, population genetics and biosafety could benefit of a low cost and fast turnaround time approach such as the 2b‐RAD analysis proposed here.     

Population structure of a vector‐borne plant parasite

Parasites are among the most diverse groups of life on Earth, yet complex natural histories often preclude studies of their speciation processes. The biology of parasitic plants facilitates in situ collection of data on both genetic structure and the mechanisms responsible for that structure. Here, we studied the role of mating, dispersal and establishment in host race formation of a parasitic plant. We investigated the population genetics of a vector‐borne desert mistletoe (Phoradendron californicum) across two legume host tree species (Senegalia greggii and Prosopis velutina) in the Sonoran desert using microsatellites. Consistent with host race formation, we found strong host‐associated genetic structure in sympatry, little genetic variation due to geographic site and weak isolation by distance. We hypothesize that genetic differentiation results from differences in the timing of mistletoe flowering by host species, as we found initial flowering date of individual mistletoes correlated with genetic ancestry. Hybrids with intermediate ancestry were detected genetically. Individuals likely resulting from recent, successful establishment events following dispersal between the host species were detected at frequencies similar to hybrids between host races. Therefore, barriers to gene flow between the host races may have been stronger at mating than at dispersal. We also found higher inbreeding and within‐host individual relatedness values for mistletoes on the more rare and isolated host species (S. greggii). Our study spanned spatial scales to address how interactions with both vectors and hosts influence parasitic plant structure with implications for parasite virulence evolution and speciation.     

Population genetics and sympatric divergence of the freshwater gudgeon,Gobiobotia filifer,in the Yangtze River inferred from mitochondrial DNA

The ecosystem and Pleistocene glaciations play important roles in population demography. The freshwater gudgeon, Gobiobotia filifer, is an endemic benthic fish in the Yangtze River and is a good model for ecological and evolutionary studies. This study aimed to decode the population structure of G. filifer in the Yangtze River and reveal whether divergence occurred before or after population radiation. A total of 292 specimens from eight locations in the upper and middle reaches of the Yangtze River were collected from 2014 to 2016 and analyzed via mitochondrial DNA Cyt b gene sequencing. A moderately high level of genetic diversity was found without structures among the population. However, phylogenetic and network topology showed two distinct haplotype groups, and each group contained a similar proportion of individuals from all sampled sites. This suggested the existence of two genetically divergent source populations in G. filifer. We deduced that a secondary contact of distinct glacial refugia was the main factor creating sympatric populations of G. filifer, and climate improvement promoted population expansion and colonization.     

Population structure and linkage disequilibrium in a large collection of Fusarium oxysporum strains analysed through iPBS markers

In the current study, 160 pathogenic strains of Fusarium oxysporum collected from tomato, eggplant and pepper were studied. Eighteen inter‐primer binding site (iPBS)‐retrotransposon primers were used, and these primers generated 205 scorable polymorphic bands. The number of polymorphic bands per primer varied between 9 and 19, with a mean of 11 bands per primer. The highest polymorphism information content (PIC) value was determined as 0.27, and the lowest was 0.05. The unweighted pair‐group method with arithmetic averages (UPGMA) dendrogram including a heat map revealed that the 160 pathogenic strains of F. oxysporum were divided into two main clusters. The first cluster mainly included F. oxysporum f. sp. capsici (FOC) and F. oxysporum f. sp. melongenae (FOMG) isolates. The second cluster mainly comprised F. oxysporum f. sp. lycopersici (FOL) and F. oxysporum f. sp. radicis lycopersici (FORL) isolates. The highest percentage of loci in significant linkage disequilibrium (LD) was detected for FOL, whereas the lowest level of LD was found for FOC, and 95.2%, 99.4%, 99.1% and 97.4% of the relative kinship estimates were less than 0.4 for FOL, FOMG, FORL and FOC, respectively. LD differences were detected among formae speciales, and LD was higher in FOL as compare to FOC species. The findings of this study confirm that iPBS‐retrotransposon markers are highly polymorphic at the intraspecific level in Fusarium spp.     

Survival and abundance of short‐finned pilot whales in the archipelago of Madeira,NE Atlantic

Estimates of population parameters for the short‐finned pilot whale, Globicephala macrorhynchus, are scarce in literature, contributing to an International Union for Conservation of Nature (IUCN) status of Data Deficient. In this study, photo‐identification data collected over 7 yr from Madeira were used to estimate for the first time survivorship, capture probability, and abundance in this species using mark‐recapture methodology. The Cormack‐Jolly‐Seber model estimated that the adult island‐associated (i.e., resident and regular visitor) whales had a constant survival rate of 0.960 (95% CI: 0.853–0.990) and an annual capture probability varying between 0.372 (CI: 0.178–0.619) and 0.843 (CI: 0.619–0.947). A parameterization of the Jolly‐Seber model estimated that 140 island‐associated whales (CI: 131–151) used the area throughout the course of the study. Based on a closed population model, the most precise (lower CV) annual estimate of the total number of pilot whales using the southern and eastern waters of Madeira (~900 km²) in a 3 mo period covering summer/autumn was 334 animals (CI: 260–437). No trend was observed. Despite including biases, the approach used in this study provided plausible estimates of population parameters, which can contribute to the regional conservation strategies.     

A single‐nucleotide polymorphism of miR‐146a and psoriasis: an association and functional study

Epidermal growth factor receptor (EGFR), which is overexpressed in psoriatic lesions, has been proven to contribute to the hyperproliferation of keratinocytes in psoriasis. Single nucleotide polymorphisms (SNPs) involved in miRNAs that can regulate the expression of EGFR could potentially influence the development of psoriasis. The present study investigated the association between a functional SNP of rs2910164 in miR‐146a and the risk of psoriasis in the Chinese Han population. A total of 521 Han Chinese patients with psoriasis and 582 healthy controls were recruited in this study. The miR‐146a rs2910164 SNP was genotyped by polymerase chain reaction‐restriction fragment length polymorphism. Overall, a significantly increased risk of psoriasis was associated with the rs2910164 miR‐146a CG and GG genotypes (adjusted OR, 1.38; 95% CI, 1.06–1.80). Furthermore, the rs2910164G allele in miR‐146a attenuated its inhibitory regulation on the expression of EGFR as well as the proliferation of human keratinocytes, and lowered the level of miR‐146a in the psoriatic lesions. These findings indicate that the rs2910164G allele in miR‐146a weakens its suppression on the proliferation of keratinocytes probably through the decreased inhibition of the target gene, EGFR, which may account for the increased risk of psoriasis in this study population.