Histone H3.3G34-Mutant Interneuron Progenitors Co-opt PDGFRA for Gliomagenesis.

Histone H3.3 glycine 34 to arginine/valine (G34R/V) mutations drive deadly gliomas and show exquisite regional and temporal specificity, suggesting a developmental context permissive to their effects. Here we show that 50% of G34R/V tumors (n = 95) bear activating PDGFRA mutations that display strong selection pressure at recurrence. Although considered gliomas, G34R/V tumors actually arise in GSX2/DLX-expressing interneuron progenitors, where G34R/V mutations impair neuronal differentiation. The lineage of origin may facilitate PDGFRA co-option through a chromatin loop connecting PDGFRA to GSX2 regulatory elements, promoting PDGFRA overexpression and mutation. At the single-cell level, G34R/V tumors harbor dual neuronal/astroglial identity and lack oligodendroglial programs, actively repressed by GSX2/DLX-mediated cell fate specification. G34R/V may become dispensable for tumor maintenance, whereas mutant-PDGFRA is potently oncogenic. Collectively, our results open novel research avenues in deadly tumors. G34R/V gliomas are neuronal malignancies where interneuron progenitors are stalled in differentiation by G34R/V mutations and malignant gliogenesis is promoted by co-option of a potentially targetable pathway, PDGFRA signaling.

Comprehensive Genomic Analysis of Cemento-Ossifying Fibroma.

Cemento-ossifying fibroma (COF) of the jaws is currently classified as a benign mesenchymal odontogenic tumor, and only targeted approaches have been used to assess its genetic alterations. A minimal proportion of COFs harbor CDC73 somatic mutations, and copy number alterations (CNAs) involving chromosomes 7 and 12 have recently been reported in a small proportion of cases. However, the genetic background of COFs remains obscure. We used a combination of whole-exome sequencing and RNA sequencing to assess somatic mutations, fusion transcripts, and CNAs in a cohort of 12 freshly collected COFs. No recurrent fusions have been identified among the 5 cases successfully analyzed by RNA sequencing, with in-frame fusions being detected in 2 cases (MARS1::GOLT1B and PARG::BMS1 in one case and NCLN::FZR1 and NFIC::SAMD1 in the other case) and no candidate fusions identified for the remaining 3 cases. No recurrent pathogenic mutations were detected in the 11 cases that had undergone whole-exome sequencing. A KRAS p.L19F missense variant was detected in one case, and 2 CDC73 deletions were detected in another case. The other variants were of uncertain significance and included variants in PC, ACTB, DOK6, HACE1, and COL1A2 and previously unreported variants in PTPN14, ATP5F1C, APOBEC1, HDAC5, ATF7IP, PARP2, and ACTR3B. The affected genes do not clearly converge on any signaling pathway. CNAs were detected in 5/11 cases (45%), with copy gains involving chromosome 12 occurring in 3/11 cases (27%). In conclusion, no recurrent fusions or pathogenic variants have been detected in the present COF cohort, with copy gains involving chromosome 12 occurring in 27% of cases.

De novo TRPV4 Leu619Pro variant causes a new channelopathy characterised by giant cell lesions of the jaws and skull, skeletal abnormalities and polyneuropathy.

BACKGROUND: Pathogenic germline variants in Transient Receptor Potential Vanilloid 4 Cation Channel (TRPV4) lead to channelopathies, which are phenotypically diverse and heterogeneous disorders grossly divided in neuromuscular disorders and skeletal dysplasia. We recently reported in sporadic giant cell lesions of the jaws (GCLJs) novel, somatic, heterozygous, gain-of-function mutations in TRPV4, at Met713. METHODS: Here we report two unrelated women with a de novo germline p.Leu619Pro TRPV4 variant and an overlapping systemic disorder affecting all organs individually described in TRPV4 channelopathies. RESULTS: From an early age, both patients had several lesions of the nervous system including progressive polyneuropathy, and multiple aggressive giant cell-rich lesions of the jaws and craniofacial/skull bones, and other skeletal lesions. One patient had a relatively milder disease phenotype possibly due to postzygotic somatic mosaicism. Indeed, the TRPV4 p.Leu619Pro variant was present at a lower frequency (variant allele frequency (VAF)=21.6%) than expected for a heterozygous variant as seen in the other proband, and showed variable regional frequency in the GCLJ (VAF ranging from 42% to 10%). In silico structural analysis suggests that the gain-of-function p.Leu619Pro alters the ion channel activity leading to constitutive ion leakage. CONCLUSION: Our findings define a novel polysystemic syndrome due to germline TRPV4 p.Leu619Pro and further extend the spectrum of TRPV4 channelopathies. They further highlight the convergence of TRPV4 mutations on different organ systems leading to complex phenotypes which are further mitigated by possible post-zygotic mosaicism. Treatment of this disorder is challenging, and surgical intervention of the GCLJ worsens the lesions, suggesting the future use of MEK inhibitors and TRPV4 antagonists as therapeutic modalities for unmet clinical needs.

Inhibition of nuclear export restores nuclear localization and residual tumor suppressor function of truncated SMARCB1/INI1 protein in a molecular subset of atypical teratoid/rhabdoid tumors.

Loss of nuclear SMARCB1 (INI1/hSNF5/BAF47) protein expression due to biallelic mutations of the SMARCB1 tumor suppressor gene is a hallmark of atypical teratoid/rhabdoid tumors (ATRT), but the presence of cytoplasmic SMARCB1 protein in these tumors has not yet been described. In a series of 102 primary ATRT, distinct cytoplasmic SMARCB1 staining on immunohistochemistry was encountered in 19 cases (19%) and was highly over-represented in cases showing pathogenic sequence variants leading to truncation or mutation of the C-terminal part of SMARCB1 (15/19 vs. 4/83; Chi-square: 56.04, p = 1.0E-10) and, related to this, in tumors of the molecular subgroup ATRT-TYR (16/36 vs. 3/66; Chi-square: 24.47, p = 7.6E-7). Previous reports have indicated that while SMARCB1 lacks a bona fide nuclear localization signal, it harbors a masked nuclear export signal (NES) and that truncation of the C-terminal region results in unmasking of this NES leading to cytoplasmic localization. To determine if cytoplasmic localization found in ATRT is due to unmasking of NES, we generated GFP fusions of one of the SMARCB1 truncating mutations (p.Q318X) found in the tumors along with a p.L266A mutation, which was shown to disrupt the interaction of SMARCB1-NES with exportin-1. We found that while the GFP-SMARCB1(Q318X) mutant localized to the cytoplasm, the double mutant GFP-SMARCB1(Q318X;L266A) localized to the nucleus, confirming NES requirement for cytoplasmic localization. Furthermore, cytoplasmic SMARCB1(Q318X) was unable to cause senescence as determined by morphological observations and by senescence-associated ß-galactosidase assay, while nuclear SMARCB1(Q318X;L266A) mutant regained this function. Selinexor, a selective exportin-1 inhibitor, was effective in inhibiting the nuclear export of SMARCB1(Q318X) and caused rapid cell death in rhabdoid tumor cells. In conclusion, inhibition of nuclear export restores nuclear localization and residual tumor suppressor function of truncated SMARCB1. Therapies aimed at preventing nuclear export of mutant SMARCB1 protein may represent a promising targeted therapy in ATRT harboring truncating C-terminal SMARCB1 mutations.

Whole-exome sequencing reveals novel vacuolar ATPase genes' variants and variants in genes involved in lysosomal biology and autophagosomal formation in oral granular cell tumors.

BACKGROUND: Granular cell tumors (GCTs) are rare neuroectodermal soft tissue neoplasms that mainly affect the skin of the upper limbs and trunks and the oral cavity. GCTs are derived from Schwann cells and, ultrastructurally, their intracytoplasmic granules are considered autophagosomes or autophagolysosomes and are consistent with myelin accumulation. METHODS: In this study, a convenience set of 22 formalin-fixed, paraffin-embedded samples of oral GCTs, all but one sample located at the tongue, was screened for mutations by whole-exome (WES) or targeted sequencing. RESULTS: WES revealed two novel variants in genes of the vacuolar ATPase (V-ATPase) complex: ATP6AP1 frameshift c.746_749del, leading to p.P249Hfs*4, and ATP6V1A non-synonymous c.G868A, leading to p.D290N. Each of these mutations occurred in one case. With regard to the samples that were wild type for these V-ATPase variants, at least two samples presented variants in genes that are part of endosomal/lysosomal/autophagosomal networks including ABCA8, ABCC6, AGAP3, ATG9A, CTSB, DNAJC13, GALC, NPC1, SLC15A3, SLC31A2, and TMEM104. CONCLUSION: Although the mechanisms involved in oral GCT initiation and progression remain unclear, our results suggest that oral GCTs have V-ATPase variants similarly to GCTs from other tissues/organs, and additionally show variants in lysosomes/endosomes/autophagosomal genes.

Longitudinal mutational analysis of a cerebellar pilocytic astrocytoma recurring as a ganglioglioma.

A cerebellar pilocytic astrocytoma (PA) in a child recurred first with a PA histology and then with features of a ganglioglioma (GG). Molecular genetic analyses of the tumors confirmed a BRAF V600E mutation in all. They also all harbored a T202M mutation in ERK1, a kinase downstream of BRAF that is implicated in glial versus neuronal differentiation. The GG sample contained several variants that were not present in the PA samples; in particular, it had a truncating mutation in MAP2. These findings not only underscore the role of BRAF as oncogenic driver but also suggest that other genes may influence tumor morphology.

Prognostic value of medulloblastoma extent of resection after accounting for molecular subgroup: a retrospective integrated clinical and molecular analysis.

BACKGROUND: Patients with incomplete surgical resection of medulloblastoma are controversially regarded as having a marker of high-risk disease, which leads to patients undergoing aggressive surgical resections, so-called second-look surgeries, and intensified chemoradiotherapy. All previous studies assessing the clinical importance of extent of resection have not accounted for molecular subgroup. We analysed the prognostic value of extent of resection in a subgroup-specific manner. METHODS: We retrospectively identified patients who had a histological diagnosis of medulloblastoma and complete data about extent of resection and survival from centres participating in the Medulloblastoma Advanced Genomics International Consortium. We collected from resections done between April, 1997, and February, 2013, at 35 international institutions. We established medulloblastoma subgroup affiliation by gene expression profiling on frozen or formalin-fixed paraffin-embedded tissues. We classified extent of resection on the basis of postoperative imaging as gross total resection (no residual tumour), near-total resection (<1·5 cm(2) tumour remaining), or sub-total resection (≥1·5 cm(2) tumour remaining). We did multivariable analyses of overall survival and progression-free survival using the variables molecular subgroup (WNT, SHH, group 4, and group 3), age (<3 vs ≥3 years old), metastatic status (metastases vs no metastases), geographical location of therapy (North America/Australia vs rest of the world), receipt of chemotherapy (yes vs no) and receipt of craniospinal irradiation (<30 Gy or >30 Gy vs no craniospinal irradiation). The primary analysis outcome was the effect of extent of resection by molecular subgroup and the effects of other clinical variables on overall and progression-free survival. FINDINGS: We included 787 patients with medulloblastoma (86 with WNT tumours, 242 with SHH tumours, 163 with group 3 tumours, and 296 with group 4 tumours) in our multivariable Cox models of progression-free and overall survival. We found that the prognostic benefit of increased extent of resection for patients with medulloblastoma is attenuated after molecular subgroup affiliation is taken into account. We identified a progression-free survival benefit for gross total resection over sub-total resection (hazard ratio [HR] 1·45, 95% CI 1·07-1·96, p=0·16) but no overall survival benefit (HR 1·23, 0·87-1·72, p=0·24). We saw no progression-free survival or overall survival benefit for gross total resection compared with near-total resection (HR 1·05, 0·71-1·53, p=0·8158 for progression-free survival and HR 1·14, 0·75-1·72, p=0·55 for overall survival). No significant survival benefit existed for greater extent of resection for patients with WNT, SHH, or group 3 tumours (HR 1·03, 0·67-1·58, p=0·89 for sub-total resection vs gross total resection). For patients with group 4 tumours, gross total resection conferred a benefit to progression-free survival compared with sub-total resection (HR 1·97, 1·22-3·17, p=0·0056), especially for those with metastatic disease (HR 2·22, 1·00-4·93, p=0·050). However, gross total resection had no effect on overall survival compared with sub-total resection in patients with group 4 tumours (HR 1·67, 0·93-2·99, p=0·084). INTERPRETATION: The prognostic benefit of increased extent of resection for patients with medulloblastoma is attenuated after molecular subgroup affiliation is taken into account. Although maximum safe surgical resection should remain the standard of care, surgical removal of small residual portions of medulloblastoma is not recommended when the likelihood of neurological morbidity is high because there is no definitive benefit to gross total resection compared with near-total resection. FUNDING: Canadian Cancer Society Research Institute, Terry Fox Research Institute, Canadian Institutes of Health Research, National Institutes of Health, Pediatric Brain Tumor Foundation, and the Garron Family Chair in Childhood Cancer Research.

A recurrent PDGFRB mutation causes familial infantile myofibromatosis.

Infantile myofibromatosis (IM) is the most common benign fibrous tumor of soft tissues affecting young children. By using whole-exome sequencing, RNA sequencing, and targeted sequencing, we investigated germline and tumor DNA in individuals from four distinct families with the familial form of IM and in five simplex IM cases with no previous family history of this disease. We identified a germline mutation c.1681C>T (p.Arg561Cys) in platelet-derived growth factor receptor ß (PDGFRB) in all 11 affected individuals with familial IM, although none of the five individuals with nonfamilial IM had mutations in this gene. We further identified a second heterozygous mutation in PDGFRB in two myofibromas from one of the affected familial cases, indicative of a potential second hit in this gene in the tumor. PDGFR-ß promotes growth of mesenchymal cells, including blood vessels and smooth muscles, which are affected in IM. Our findings indicate p.Arg561Cys substitution in PDGFR-ß as a cause of the dominant form of this disease. They provide a rationale for further investigations of this specific mutation and gene to assess the benefits of targeted therapies against PDGFR-ß in aggressive life-threatening familial forms of the disease.

Murine diet/tissue and human brain tumorigenesis alter Mthfr/MTHFR 5'-end methylation.

Polymorphisms and decreased activity of methylenetetrahydrofolate reductase (MTHFR) are linked to disease, including cancer. However, epigenetic regulation has not been thoroughly studied. Our goal was to generate DNA methylation profiles of murine/human MTHFR gene regions and examine methylation in brain and liver tumors. Pyrosequencing in four murine tissues revealed minimal DNA methylation in the CpG island. Higher methylation was seen in liver or intestine in the CpG island shore 5' to the upstream translational start site or in another region 3' to the downstream start site. In the latter region, there was negative correlation between expression and methylation. Three orthologous regions were investigated in human MTHFR, as well as a fourth region between the two translation start sites. We found significantly increased methylation in three regions (not the CpG island) in pediatric astrocytomas compared with control brain, with decreased expression in tumors. Methylation in hepatic carcinomas was also increased in the three regions compared with normal liver, but the difference was significant for only one CpG. This work, the first overview of the Mthfr/MTHFR epigenetic landscape, suggests regulation through methylation in some regions, demonstrates increased methylation/decreased expression in pediatric astrocytomas, and should serve as a resource for future epigenetic studies.

Germline and somatic FGFR1 abnormalities in dysembryoplastic neuroepithelial tumors.

Dysembryoplastic neuroepithelial tumor (DNET) is a benign brain tumor associated with intractable drug-resistant epilepsy. In order to identify underlying genetic alterations and molecular mechanisms, we examined three family members affected by multinodular DNETs as well as 100 sporadic tumors from 96 patients, which had been referred to us as DNETs. We performed whole-exome sequencing on 46 tumors and targeted sequencing for hotspot FGFR1 mutations and BRAF p.V600E was used on the remaining samples. FISH, copy number variation assays and Sanger sequencing were used to validate the findings. By whole-exome sequencing of the familial cases, we identified a novel germline FGFR1 mutation, p.R661P. Somatic activating FGFR1 mutations (p.N546K or p.K656E) were observed in the tumor samples and further evidence for functional relevance was obtained by in silico modeling. The FGFR1 p.K656E mutation was confirmed to be in cis with the germline p.R661P variant. In 43 sporadic cases, in which the diagnosis of DNET could be confirmed on central blinded neuropathology review, FGFR1 alterations were also frequent and mainly comprised intragenic tyrosine kinase FGFR1 duplication and multiple mutants in cis (25/43; 58.1 %) while BRAF p.V600E alterations were absent (0/43). In contrast, in 53 cases, in which the diagnosis of DNET was not confirmed, FGFR1 alterations were less common (10/53; 19 %; p < 0.0001) and hotspot BRAF p.V600E (12/53; 22.6 %) (p < 0.001) prevailed. We observed overexpression of phospho-ERK in FGFR1 p.R661P and p.N546K mutant expressing HEK293 cells as well as FGFR1 mutated tumor samples, supporting enhanced MAP kinase pathway activation under these conditions. In conclusion, constitutional and somatic FGFR1 alterations and MAP kinase pathway activation are key events in the pathogenesis of DNET. These findings point the way towards existing targeted therapies.

Epigenetic dysregulation: a novel pathway of oncogenesis in pediatric brain tumors.

A remarkably large number of "epigenetic regulators" have been recently identified to be altered in cancers and a rapidly expanding body of literature points to "epigenetic addiction" (an aberrant epigenetic state to which a tumor is addicted) as a new previously unsuspected mechanism of oncogenesis. Although mutations are also found in canonical signaling pathway genes, we and others identified chromatin-associated proteins to be more commonly altered by somatic alterations than any other class of oncoprotein in several subgroups of childhood high-grade brain tumors. Furthermore, as these childhood malignancies carry fewer non-synonymous somatic mutations per case in contrast to most adult cancers, these mutations are likely drivers in these tumors. Herein, we will use as examples of this novel hallmark of oncogenesis high-grade astrocytomas, including glioblastoma, and a subgroup of embryonal tumors, embryonal tumor with multilayered rosettes (ETMR) to describe the novel molecular defects uncovered in these deadly tumors. We will further discuss evidence for their profound effects on the epigenome. The relative genetic simplicity of these tumors promises general insights into how mutations in the chromatin machinery modify downstream epigenetic signatures to drive transformation, and how to target this plastic genetic/epigenetic interface.

The Himalayas: barrier and conduit for gene flow.

The Himalayan mountain range is strategically located at the crossroads of the major cultural centers in Asia, the Middle East and Europe. Although previous Y-chromosome studies indicate that the Himalayas served as a natural barrier for gene flow from the south to the Tibetan plateau, this region is believed to have played an important role as a corridor for human migrations between East and West Eurasia along the ancient Silk Road. To evaluate the effects of the Himalayan mountain range in shaping the maternal lineages of populations residing on either side of the cordillera, we analyzed mitochondrial DNA variation in 344 samples from three Nepalese collections (Newar, Kathmandu and Tamang) and a general population of Tibet. Our results revealed a predominantly East Asian-specific component in Tibet and Tamang, whereas Newar and Kathmandu are both characterized by a combination of East and South Central Asian lineages. Interestingly, Newar and Kathmandu harbor several deep-rooted Indian lineages, including M2, R5, and U2, whose coalescent times from this study (U2, >40 kya) and previous reports (M2 and R5, >50 kya) suggest that Nepal was inhabited during the initial peopling of South Central Asia. Comparisons with our previous Y-chromosome data indicate sex-biased migrations in Tamang and a founder effect and/or genetic drift in Tamang and Newar. Altogether, our results confirm that while the Himalayas acted as a geographic barrier for human movement from the Indian subcontinent to the Tibetan highland, it also served as a conduit for gene flow between Central and East Asia.

Association of novel mutation in TRPV4 with familial nonsyndromic craniosynostosis with complete penetrance and variable expressivity.

OBJECTIVE: The aim of this study was to characterize a novel pathogenic variant in the transient receptor potential vanilloid 4 (TRPV4) gene, causing familial nonsyndromic craniosynostosis (CS) with complete penetrance and variable expressivity. METHODS: Whole-exome sequencing was performed on germline DNA of a family with nonsyndromic CS to a mean depth coverage of 300× per sample, with greater than 98% of the targeted region covered at least 25×. In this study, the authors detected a novel variant, c.496C>A in TRPV4, exclusively in the four affected family members. The variant was modeled using the structure of the TRPV4 protein from Xenopus tropicalis. In vitro assays in HEK293 cells overexpressing wild-type TRPV4 or TRPV4 p.Leu166Met were used to assess the effect of the mutation on channel activity and downstream MAPK signaling. RESULTS: The authors identified a novel, highly penetrant heterozygous variant in TRPV4 (NM_021625.4:c.496C>A) causing nonsyndromic CS in a mother and all three of her children. This variant results in an amino acid change (p.Leu166Met) in the intracellular ankyrin repeat domain distant from the Ca2+-dependent membrane channel domain. In contrast to other TRPV4 mutations in channelopathies, this variant does not interfere with channel activity as identified by in silico modeling and in vitro overexpression assays in HEK293 cells. CONCLUSIONS: Based on these findings, the authors hypothesized that this novel variant causes CS by modulating the binding of allosteric regulatory factors to TRPV4 rather than directly modifying its channel activity. Overall, this study expands the genetic and functional spectrum of TRPV4 channelopathies and is particularly relevant for the genetic counseling of CS patients.

Drosophila Tet Is Required for Maintaining Glial Homeostasis in Developing and Adult Fly Brains.

Ten-eleven translocation (TET) proteins are crucial epigenetic regulators highly conserved in multicellular organisms. TETs' enzymatic function in demethylating 5-methyl cytosine in DNA is required for proper development and TETs are frequently mutated in cancer. Recently, Drosophila melanogaster Tet (dTet) was shown to be highly expressed in developing fly brains and discovered to play an important role in brain and muscle development as well as fly behavior. Furthermore, dTet was shown to have different substrate specificity compared with mammals. However, the exact role dTet plays in glial cells and how ectopic TET expression in glial cells contributes to tumorigenesis and glioma is still not clear. Here, we report a novel role for dTet specifically in glial cell organization and number. We show that loss of dTet affects the organization of a specific glia population in the optic lobe, the "optic chiasm" glia. Additionally, we find irregularities in axon patterns in the ventral nerve cord (VNC) both, in the midline and longitudinal axons. These morphologic glia and axonal defects were accompanied by locomotor defects in developing larvae escalating to immobility in adult flies. Furthermore, glia homeostasis was disturbed in dTet-deficient brains manifesting in gain of glial cell numbers and increased proliferation. Finally, we establish a Drosophila model to understand the impact of human TET3 in glia and find that ectopic expression of hTET3 in dTet-expressing cells causes glia expansion in larval brains and affects sleep/rest behavior and the circadian clock in adult flies.

The Soliga, an isolated tribe from Southern India: genetic diversity and phylogenetic affinities.

India's role in the dispersal of modern humans can be explored by investigating its oldest inhabitants: the tribal people. The Soliga people of the Biligiri Rangana Hills, a tribal community in Southern India, could be among the country's first settlers. This forest-bound, Dravidian speaking group, lives isolated, practicing subsistence-level agriculture under primitive conditions. The aim of this study is to examine the phylogenetic relationships of the Soligas in relation to 29 worldwide, geographically targeted, reference populations. For this purpose, we employed a battery of 15 hypervariable autosomal short tandem repeat loci as markers. The Soliga tribe was found to be remarkably different from other Indian populations including other southern Dravidian-speaking tribes. In contrast, the Soliga people exhibited genetic affinity to two Australian aboriginal populations. This genetic similarity could be attributed to the 'Out of Africa' migratory wave(s) along the southern coast of India that eventually reached Australia. Alternatively, the observed genetic affinity may be explained by more recent migrations from the Indian subcontinent into Australia.

Y-STR diversity in the Himalayas.

Linguistic and ethnic diversity throughout the Himalayas suggests that this mountain range played an important role in shaping the genetic landscapes of the region. Previous Y-chromosome work revealed that the Himalayas acted as a biased bidirectional barrier to gene flow across the cordillera. In the present study, 17 Y-chromosomal short tandem repeat (Y-STR) loci included in the AmpFlSTR® Yfiler kit were analyzed in 344 unrelated males from three Nepalese populations (Tamang, Newar, and Kathmandu) and a general collection from Tibet. The latter displays the highest haplotype diversity (0.9990) followed by Kathmandu (0.9977), Newar (0.9570), and Tamang (0.9545). The overall haplotype diversity for the Himalayan populations at 17 Y-STR loci was 0.9973, and the corresponding values for the extended (11 loci) and minimal (nine loci) haplotypes were 0.9955 and 0.9942, respectively. No Y-STR profiles are shared across the four Himalayan collections at the 17-, 11-, and nine-locus resolutions considered, indicating a lack of recent gene flow among them. Phylogenetic analyses support our previous findings that Kathmandu, and to some extent Newar, received significant genetic influence from India while Tamang and Tibet exhibit limited or no gene flow from the subcontinent. A median-joining network of haplogroup O3a3c-M134 based on 15 Y-STR loci from our four Himalayan populations suggests either a male founder effect in Tamang, possibly from Tibet, or a recent bottleneck following their arrival south of the Himalayas from Tibet leading to their highly reduced Y single-nucleotide polymorphism and Y-STR diversity. The genetic uniqueness of the four Himalayan populations examined in this study merits the creation of separate databases for individual identification, parentage analysis, and population genetic studies.