Y chromosomal noncoding RNAs regulate autosomal gene expression via piRNAs in mouse testis.

BACKGROUND: Deciphering the functions of Y chromosome in mammals has been slow owing to the presence of repeats. Some of these repeats transcribe coding RNAs, the roles of which have been studied. Functions of the noncoding transcripts from Y chromosomal repeats however, remain unclear. While a majority of the genes expressed during spermatogenesis are autosomal, mice with different deletions of the long arm of the Y chromosome (Yq) were previously also shown to be characterized by subfertility, sterility and sperm abnormalities, suggesting the presence of effectors of spermatogenesis at this location. Here we report a set of novel noncoding RNAs from mouse Yq and explore their connection to some of the autosomal genes expressed in testis. RESULTS: We describe a set of novel mouse male-specific Y long arm (MSYq)-derived long noncoding (lnc) transcripts, named Pirmy and Pirmy-like RNAs. Pirmy shows a large number of splice variants in testis. We also identified Pirmy-like RNAs present in multiple copies at different loci on mouse Y chromosome. Further, we identified eight differentially expressed autosome-encoded sperm proteins in a mutant mouse strain, XYRIIIqdel (2/3 Yq-deleted). Pirmy and Pirmy-like RNAs have homology to 5'/3'UTRs of these deregulated autosomal genes. Several lines of experiments show that these short homologous stretches correspond to piRNAs. Thus, Pirmy and Pirmy-like RNAs act as templates for several piRNAs. In vitro functional assays reveal putative roles for these piRNAs in regulating autosomal genes. CONCLUSIONS: Our study elucidates a set of autosomal genes that are potentially regulated by MSYq-derived piRNAs in mouse testis. Sperm phenotypes from the Yq-deleted mice seem to be similar to that reported in inter-specific male-sterile hybrids. Taken together, this study provides novel insights into possible role of MSYq-derived ncRNAs in male sterility and speciation.

Reconstructing the demographic history of the Himalayan and adjoining populations.

The rugged topography of the Himalayan region has hindered large-scale human migrations, population admixture and assimilation. Such complexity in geographical structure might have facilitated the existence of several small isolated communities in this region. We have genotyped about 850,000 autosomal markers among 35 individuals belonging to the four major populations inhabiting the Himalaya and adjoining regions. In addition, we have genotyped 794 individuals belonging to 16 ethnic groups from the same region, for uniparental (mitochondrial and Y chromosomal DNA) markers. Our results in the light of various statistical analyses suggest a closer link of the Himalayan and adjoining populations to East Asia than their immediate geographical neighbours in South Asia. Allele frequency-based analyses likely support the existence of a specific ancestry component in the Himalayan and adjoining populations. The admixture time estimate suggests a recent westward migration of populations living to the East of the Himalaya. Furthermore, the uniparental marker analysis among the Himalayan and adjoining populations reveal the presence of East, Southeast and South Asian genetic signatures. Interestingly, we observed an antagonistic association of Y chromosomal haplogroups O3 and D clines with the longitudinal distance. Thus, we summarise that studying the Himalayan and adjoining populations is essential for a comprehensive reconstruction of the human evolutionary and ethnolinguistic history of eastern Eurasia.

Uncoupling of X-linked gene silencing from XIST binding by DICER1 and chromatin modulation on human inactive X chromosome.

In mammals, X-inactivation process is achieved by the cis-spreading of long noncoding Xist RNA over one of the female X chromosomes. The Xist binding accumulates histones H3 methylation and H4 hypoacetylation required for X inactivation that leads to proper dosage compensation of the X-linked genes. Co-transcription of Tsix, an antisense copy of Xist, blocks the Xist coating on the Xi. In mice ES cells, an RNase III enzyme Dicer1 disrupts Xist binding and methylated H3K27me3 accumulation on the Xi. Later, multiple reports opposed these findings raising a question regarding the possible role of Dicer1 in murine X silencing. Here, we show that reduction of DICER1 in human female cells increases XIST transcripts without compromising the binding of the XIST and histone tail modifications on the Xi. Moreover, DICER1-depleted cells show differential upregulation of many human X-linked genes by binding different amounts of acetylated histone predominantly on their active promoter sites. Therefore, X-linked gene silencing, which is thought to be coupled with the accumulation of XIST and heterochromatin markers on Xi can be disrupted in DICER1 depleted human cells. These results suggest that DICER1 has no apparent effect on the recruitment of heterochromatic markers on the Xi but is required for inactivation of differentially regulated genes for the maintenance of proper dosage compensation in differentiated cells.

Genetic evidence for recent population mixture in India.

Most Indian groups descend from a mixture of two genetically divergent populations: Ancestral North Indians (ANI) related to Central Asians, Middle Easterners, Caucasians, and Europeans; and Ancestral South Indians (ASI) not closely related to groups outside the subcontinent. The date of mixture is unknown but has implications for understanding Indian history. We report genome-wide data from 73 groups from the Indian subcontinent and analyze linkage disequilibrium to estimate ANI-ASI mixture dates ranging from about 1,900 to 4,200 years ago. In a subset of groups, 100% of the mixture is consistent with having occurred during this period. These results show that India experienced a demographic transformation several thousand years ago, from a region in which major population mixture was common to one in which mixture even between closely related groups became rare because of a shift to endogamy.

The paternal ancestry of Uttarakhand does not imitate the classical caste system of India.

Although, there have been rigorous research on the Indian caste system by several disciplines, it is still one of the most controversial socioscientific topic. Previous genetic studies on the subcontinent have supported a classical hierarchal sharing of genetic component by various castes of India. In the present study, we have used high-resolution mtDNA and Y chromosomal markers to characterize the genetic structuring of the Uttarakhand populations in the context of neighboring regions. Furthermore, we have tested whether the genetic structuring of caste populations at different social levels of this region, follow the classical chaturvarna system. Interestingly, we found that this region showed a high level of variation for East Eurasian ancestry in both maternal and paternal lines of descent. Moreover, the intrapopulation comparison showed a high level of heterogeneity, likely because of different caste hierarchy, interpolated on asymmetric admixture of populations inhabiting on both sides of the Himalayas.

The light skin allele of SLC24A5 in South Asians and Europeans shares identity by descent.

Skin pigmentation is one of the most variable phenotypic traits in humans. A non-synonymous substitution (rs1426654) in the third exon of SLC24A5 accounts for lighter skin in Europeans but not in East Asians. A previous genome-wide association study carried out in a heterogeneous sample of UK immigrants of South Asian descent suggested that this gene also contributes significantly to skin pigmentation variation among South Asians. In the present study, we have quantitatively assessed skin pigmentation for a largely homogeneous cohort of 1228 individuals from the Southern region of the Indian subcontinent. Our data confirm significant association of rs1426654 SNP with skin pigmentation, explaining about 27% of total phenotypic variation in the cohort studied. Our extensive survey of the polymorphism in 1573 individuals from 54 ethnic populations across the Indian subcontinent reveals wide presence of the derived-A allele, although the frequencies vary substantially among populations. We also show that the geospatial pattern of this allele is complex, but most importantly, reflects strong influence of language, geography and demographic history of the populations. Sequencing 11.74 kb of SLC24A5 in 95 individuals worldwide reveals that the rs1426654-A alleles in South Asian and West Eurasian populations are monophyletic and occur on the background of a common haplotype that is characterized by low genetic diversity. We date the coalescence of the light skin associated allele at 22-28 KYA. Both our sequence and genome-wide genotype data confirm that this gene has been a target for positive selection among Europeans. However, the latter also shows additional evidence of selection in populations of the Middle East, Central Asia, Pakistan and North India but not in South India.

Reconstructing Indian population history.

India has been underrepresented in genome-wide surveys of human variation. We analyse 25 diverse groups in India to provide strong evidence for two ancient populations, genetically divergent, that are ancestral to most Indians today. One, the 'Ancestral North Indians' (ANI), is genetically close to Middle Easterners, Central Asians, and Europeans, whereas the other, the 'Ancestral South Indians' (ASI), is as distinct from ANI and East Asians as they are from each other. By introducing methods that can estimate ancestry without accurate ancestral populations, we show that ANI ancestry ranges from 39-71% in most Indian groups, and is higher in traditionally upper caste and Indo-European speakers. Groups with only ASI ancestry may no longer exist in mainland India. However, the indigenous Andaman Islanders are unique in being ASI-related groups without ANI ancestry. Allele frequency differences between groups in India are larger than in Europe, reflecting strong founder effects whose signatures have been maintained for thousands of years owing to endogamy. We therefore predict that there will be an excess of recessive diseases in India, which should be possible to screen and map genetically.

MBL2 variations and malaria susceptibility in Indian populations.

Human mannose-binding lectin (MBL) encoded by the MBL2 gene is a pattern recognition protein and has been associated with many infectious diseases, including malaria. We sought to investigate the contribution of functional MBL2 gene variations to Plasmodium falciparum malaria in well-defined cases and in matched controls. We resequenced the 8.7 kb of the entire MBL2 gene in 434 individuals clinically classified with malaria from regions of India where malaria is endemic. The study cohort included 176 patients with severe malaria, 101 patients with mild malaria, and 157 ethnically matched asymptomatic individuals. In addition, 830 individuals from 32 socially, linguistically, and geographically diverse endogamous populations of India were investigated for the distribution of functional MBL2 variants. The MBL2 -221C (X) allelic variant is associated with increased risk of malaria (mild malaria odds ratio [OR] = 1.9, corrected P value [P(Corr)] = 0.0036; severe malaria OR = 1.6, P(Corr) = 0.02). The exon1 variants MBL2*B (severe malaria OR = 2.1, P(Corr) = 0.036; mild versus severe malaria OR = 2.5, P(Corr) = 0.039) and MBL2*C (mild versus severe malaria OR = 5.4, P(Corr) = 0.045) increased the odds of having malaria. The exon1 MBL2*D/*B/*C variant increased the risk for severe malaria (OR = 3.4, P(Corr) = 0.000045). The frequencies of low MBL haplotypes were significantly higher in severe malaria (14.2%) compared to mild malaria (7.9%) and asymptomatic (3.8%). The MBL2*LYPA haplotypes confer protection, whereas MBL2*LXPA increases the malaria risk. Our findings in Indian populations demonstrate that MBL2 functional variants are strongly associated with malaria and infection severity.

Shared and unique components of human population structure and genome-wide signals of positive selection in South Asia.

South Asia harbors one of the highest levels genetic diversity in Eurasia, which could be interpreted as a result of its long-term large effective population size and of admixture during its complex demographic history. In contrast to Pakistani populations, populations of Indian origin have been underrepresented in previous genomic scans of positive selection and population structure. Here we report data for more than 600,000 SNP markers genotyped in 142 samples from 30 ethnic groups in India. Combining our results with other available genome-wide data, we show that Indian populations are characterized by two major ancestry components, one of which is spread at comparable frequency and haplotype diversity in populations of South and West Asia and the Caucasus. The second component is more restricted to South Asia and accounts for more than 50% of the ancestry in Indian populations. Haplotype diversity associated with these South Asian ancestry components is significantly higher than that of the components dominating the West Eurasian ancestry palette. Modeling of the observed haplotype diversities suggests that both Indian ancestry components are older than the purported Indo-Aryan invasion 3,500 YBP. Consistent with the results of pairwise genetic distances among world regions, Indians share more ancestry signals with West than with East Eurasians. However, compared to Pakistani populations, a higher proportion of their genes show regionally specific signals of high haplotype homozygosity. Among such candidates of positive selection in India are MSTN and DOK5, both of which have potential implications in lipid metabolism and the etiology of type 2 diabetes.

Indian Siddis: African descendants with Indian admixture.

The Siddis (Afro-Indians) are a tribal population whose members live in coastal Karnataka, Gujarat, and in some parts of Andhra Pradesh. Historical records indicate that the Portuguese brought the Siddis to India from Africa about 300-500 years ago; however, there is little information about their more precise ancestral origins. Here, we perform a genome-wide survey to understand the population history of the Siddis. Using hundreds of thousands of autosomal markers, we show that they have inherited ancestry from Africans, Indians, and possibly Europeans (Portuguese). Additionally, analyses of the uniparental (Y-chromosomal and mitochondrial DNA) markers indicate that the Siddis trace their ancestry to Bantu speakers from sub-Saharan Africa. We estimate that the admixture between the African ancestors of the Siddis and neighboring South Asian groups probably occurred in the past eight generations (∼200 years ago), consistent with historical records.

Herders of Indian and European cattle share their predominant allele for lactase persistence.

Milk consumption and lactose digestion after weaning are exclusively human traits made possible by the continued production of the enzyme lactase in adulthood. Multiple independent mutations in a 100-bp region--part of an enhancer--approximately 14-kb upstream of the LCT gene are associated with this trait in Europeans and pastoralists from Saudi Arabia and Africa. However, a single mutation of purported western Eurasian origin accounts for much of observed lactase persistence outside Africa. Given the high levels of present-day milk consumption in India, together with archaeological and genetic evidence for the independent domestication of cattle in the Indus valley roughly 7,000 years ago, we sought to determine whether lactase persistence has evolved independently in the subcontinent. Here, we present the results of the first comprehensive survey of the LCT enhancer region in south Asia. Having genotyped 2,284 DNA samples from across the Indian subcontinent, we find that the previously described west Eurasian -13910 C>T mutation accounts for nearly all the genetic variation we observed in the 400- to 700-bp LCT regulatory region that we sequenced. Geography is a significant predictor of -13910*T allele frequency, and consistent with other genomic loci, its distribution in India follows a general northwest to southeast declining pattern, although frequencies among certain neighboring populations vary substantially. We confirm that the mutation is identical by descent to the European allele and is associated with the same>1 Mb extended haplotype in both populations.

Association between neuropeptide Y gene polymorphisms and alcohol dependence: a case-control study in two independent populations.

BACKGROUND: Alcohol dependence is a chronic, progressive neurobiological brain disorder. Previous research reported an inverse association between ethanol drinking and cerebral neuropeptide Y (NPY) levels. There are conflicting results of studies on NPY gene polymorphisms in association with alcohol dependence in humans. METHODS: To assess the role of the NPY gene in alcohol dependence, we genotyped three polymorphisms--in a sample of 195 subjects from the Kota population (80 alcohol dependence and 115 controls) and 141 subjects from the Badaga population (80 alcohol dependence and 61 controls). Phenotype was defined based on the DSM-IV criteria. Genotyping was performed using sequencing. Association of the NPY gene with alcohol dependence was tested by using logistic regression and haplotype analyses and linkage disequilibrium. RESULTS: All three polymorphisms were found to be in the Hardy-Weinberg equilibrium in both populations. The results of our study reveal a significant association between G1258A and alcohol dependence in both the Kota and Badaga populations. The linkage disequilibrium between the markers is not strong or significant. Haplotype analysis also did not show significant association between the NPY gene and alcohol dependence. CONCLUSION: These data support the hypothesis that alcohol dependence is influenced by the NPY G1258A polymorphism in Indian populations.

EPHX1 gene polymorphisms in alcohol dependence and their distribution among the Indian populations.

BACKGROUND: The microsomal epoxide hydrolase is a phase II enzyme of the biotransformation. The human epoxide hydrolase 1 (EPHX1) gene lies in the chromosomal region 1q42.1 and exhibits polymorphism. Two single nucleotide polymorphisms (SNPs) have been described in the coding region of the EPHX1 gene that produces two protein variants. SUBJECTS AND METHODS: A total of 604 samples belonging to 13 Indian populations were included in this study. Based on the DSM-IV criteria, 184 individuals from Kota population were classified into alcoholism cases (100) and controls (84). Genotypes of Tyr113His and His139Arg polymorphisms in the EPHX1 gene were determined using PCR and sequencing. Associations were tested using Pearson's χ(2) test and haplotype analyses. RESULTS: We found significant association between EPHX1 gene Tyr113His polymorphism and alcoholism in the Kota population (T vs. C: OR = .615, 95% CI = .399-.949, p = .027; TT vs. CC + CT: OR = .536, 95% CI = .297-.969, p = .038). The very slow activity haplotype CA (113His-139His) was also found to be associated with alcohol dependence (p = .048). Analysis of additional populations demonstrated that the Tyr113His polymorphism significantly deviated from Hardy-Weinberg equilibrium in four populations but only one population deviated for the His139Arg locus. All populations shared the four possible two-site haplotypes. Linkage disequilibrium between these two loci was not significant in any of the population studied. CONCLUSION: EPHX1 gene polymorphisms and haplotypes are associated with an increased risk for alcoholism in the Kota population. This is the first report from India that will serve as a template for future investigations of the prevalence of EPHX1 alleles in association with various clinical entities.

Cytochrome b based genetic differentiation of Indian wild pig (Sus scrofa cristatus) and domestic pig (Sus scrofa domestica) and its use in wildlife forensics.

The Indian wild pig (Sus scrofa cristatus) is a protected species and listed in the Indian Wildlife (Protection) Act, 1972. The wild pig is often hunted illegally and sold in market as meat warranting punishment under law. To avoid confusion in identification of these two subspecies during wildlife forensic examinations, we describe genetic differentiation of Indian wild and domestic pigs using a molecular technique. Analysis of sequence generated from the partial fragment (421bp) of mitochondrial DNA (mtDNA) cytochrome b (Cyt b) gene exhibited unambiguous (>3%) genetic variation between Indian wild and domestic pigs. We observed nine forensically informative nucleotide sequence (FINS) variations between Indian wild and domestic pigs. The overall genetic variation described in this study is helpful in forensic identification of the biological samples of wild and domestic pigs. It also helped in differentiating the Indian wild pig from other wild pig races. This study indicates that domestic pigs in India are not descendent of the Indian wild pig, however; they are closer to the other wild pig races found in Asia and Europe.

Population genetic structure in Indian Austroasiatic speakers: the role of landscape barriers and sex-specific admixture.

The geographic origin and time of dispersal of Austroasiatic (AA) speakers, presently settled in south and southeast Asia, remains disputed. Two rival hypotheses, both assuming a demic component to the language dispersal, have been proposed. The first of these places the origin of Austroasiatic speakers in southeast Asia with a later dispersal to south Asia during the Neolithic, whereas the second hypothesis advocates pre-Neolithic origins and dispersal of this language family from south Asia. To test the two alternative models, this study combines the analysis of uniparentally inherited markers with 610,000 common single nucleotide polymorphism loci from the nuclear genome. Indian AA speakers have high frequencies of Y chromosome haplogroup O2a; our results show that this haplogroup has significantly higher diversity and coalescent time (17-28 thousand years ago) in southeast Asia, strongly supporting the first of the two hypotheses. Nevertheless, the results of principal component and "structure-like" analyses on autosomal loci also show that the population history of AA speakers in India is more complex, being characterized by two ancestral components-one represented in the pattern of Y chromosomal and EDAR results and the other by mitochondrial DNA diversity and genomic structure. We propose that AA speakers in India today are derived from dispersal from southeast Asia, followed by extensive sex-specific admixture with local Indian populations.

High prevalence of Arginine to Glutamine substitution at 98, 141 and 162 positions in Troponin I (TNNI3) associated with hypertrophic cardiomyopathy among Indians.

BACKGROUND: Troponin I (TNNI3) is the inhibitory subunit of the thin filament regulatory complex Troponin, which confers calcium-sensitivity to striated muscle actomyosin ATPase activity. Mutations (2-7%) in this gene had been reported in hypertrophic cardiomyopathy patients (HCM). However, the frequencies of mutations and associated clinical presentation have not been established in cardiomyopathy patients of Indian origin, hence we have undertaken this study. METHODS: We have sequenced all the exons, including the exon-intron boundaries of TNNI3 gene in 101 hypertrophic cardiomyopathy patients (HCM), along with 160 healthy controls, inhabited in the same geographical region of southern India. RESULTS: Our study revealed a total of 16 mutations. Interestingly, we have observed Arginine to Glutamine (R to Q) mutation at 3 positions 98, 141 and 162, exclusively in HCM patients with family history of sudden cardiac death. The novel R98Q was observed in a severe hypertrophic obstructive cardiomyopathy patient (HOCM). The R141Q mutation was observed in two familial cases of severe asymmetric septal hypertrophy (ASH++). The R162Q mutation was observed in a ASH++ patient with mean septal thickness of 29 mm, and have also consists of allelic heterogeneity by means of having one more synonymous (E179E) mutation at g.4797: G → A: in the same exon 7, which replaces a very frequent codon (GAG: 85%) with a rare codon (GAA: 14%). Screening for R162Q mutation in all the available family members revealed its presence in 9 individuals, including 7 with allelic heterogeneity (R162Q and E179E) of which 4 were severely affected. We also found 2 novel SNPs, (g.2653; G → A and g.4003 C → T) exclusively in HCM, and in silico analysis of these SNPs have predicted to cause defect in recognition/binding sites for proteins responsible for proper splicing. CONCLUSION: Our study has provided valuable information regarding the prevalence of TNNI3 mutations in Indian HCM patients and its risk assessment, these will help in genetic counseling and to adopt appropriate treatment strategies.

Dopamine transporter (DAT1) VNTR polymorphism and alcoholism in two culturally different populations of south India.

It is well established that the central dopaminergic reward pathway is likely involved in alcohol intake and the progression of alcohol dependence. Dopamine transporter (DAT1) mediates the active re-uptake of DA from the synapse and is a principal regulator of dopaminergic neurotransmission. The gene for the human DAT1 displays several polymorphisms, including a 40-bp variable number of tandem repeats (VNTR) ranging from 3 to 16 copies in the 3'-untranslated region (UTR) of the gene. To assess the role of this gene in alcoholism, we genotyped the VNTR of DAT1 gene in a sample of 206 subjects from the Kota population (111 alcohol dependence cases and 95 controls) and 142 subjects from Badaga population (81 alcohol dependence cases and 61 controls). Both populations inhabit a similar environmental zone, but have different ethnic histories. Phenotype was defined based on the DSM-IV criteria. Genotyping was performed using PCR and electrophoresis. The association of DAT1 with alcoholism was tested by using the Clump v1.9 program which uses the Monte Carlo method. In both Kota and Badaga populations, the allele A10 was the most frequent allele followed by allele A9. The genotypic distribution is in Hardy-Weinberg equilibrium in both cases and control groups of Kota and Badaga populations. The DAT1 VNTR was significantly associated with alcoholism in Badaga population but not in Kota population. Our results suggest that the A9 allele of the DAT gene is involved in vulnerability to alcoholism, but that these associations are population specific.

Novel alleles of HLA-DQ and -DR loci show association with recurrent miscarriages among South Indian women.

BACKGROUND In this study, recurrent miscarriages (RMs) are defined as loss of two or more clinically detectable pregnancies before 20 weeks of gestation. HLA has been thought to play a role in RM. However, the results of earlier studies on the role of different human leucocyte antigen (HLA) genes were conflicting and inconclusive. In the present study, we investigate HLA genes (HLA-DRA, HLA-DRB1, HLA-DQA1 and HLA-DQB1) in RM couples with unknown etiology and normal couples. METHODS Blood samples from 143 RM couples and 150 control couples were analyzed, firstly to validate previously reported association studies and secondly to explore whether any novel alleles or haplotypes specific to Indian populations can be observed to be associated with RM. HLA typing was carried out by DNA sequencing. RESULTS Results suggest an association of the DQB1*03:03:02 allele with RM (odd ratio = 2.66; p(c) = 0.02; confidence interval = 1.47-4.84). Haplotypes of the DQA1 and DQB1 risk alleles also showed a significant association with RM, albeit not after Bonferroni correction for multiple comparisons. CONCLUSIONS HLA-DQB1 appears to have a strong involvement in the manifestation of RM in this population from South India. The current genetic analysis of RM and control couples not only highlights the genes exhibiting a strong etiological role but also reflects the protective nature of some HLA genes against RM. Nevertheless, most of these alleles/haplotypes were not those that are implicated in RM in other ethnic backgrounds, and hence require further validation in other populations of India, from different ethnic and/or geographic backgrounds.

Germline study of AR gene of Indian women with ovarian failure.

OBJECTIVE: Present study was designed for carrying out the mutational analysis of the entire Androgen receptor (AR) gene including two microsatellite (CAG)n, (GGN)n, promoter region in cases of premature ovarian failure (POF) and primary amenorrhea (PA). DESIGN: Previous reports of AR knockout mice model showed POF phenotype, this draws an attention on the role of AR gene in the aetiology of POF for the case-control association studies in POF samples (n = 133), PA samples (n = 63) and control samples (n = 200). RESULTS: We identified six mutations including four novel mutations, i.e. c.636G > A, c.1885 + 9C > A, c.1948A > G, c.1972C > A, and two previously reported mutations, i.e. c.639G > A, c.2319-78T > G. Repeat length variation was noted in the two microsatellite regions CAG and GGN, located in the coding region of exon 1 at the N-terminal region of the AR gene. The CAG repeat length was homogeneously distributed with the same frequency and no association among all cases and controls. The GGN repeat showed a significant association among the SS and SL allele with p = 0.0231 and p = 0.0476, respectively, among the POF/control samples. CONCLUSIONS: Thus, AR gene mutations may play a role in the genetic cause of POF. Identification of the underlying genetic alteration of the AR gene is important for a proper diagnosis of POF subjects.