ABSTRACT
Nasopharyngeal carcinoma (NPC) is an epithelial malignancy facilitated by Epstein-Barr Virus infection. Here we resolve the major genetic influences for NPC incidence using a genome-wide association study (GWAS), independent cohort replication, and high-resolution molecular HLA class I gene typing including 4,055 study participants from the Guangxi Zhuang Autonomous Region and Guangdong province of southern China. We detect and replicate strong association signals involving SNPs, HLA alleles, and amino acid (aa) variants across the major histocompatibility complex-HLA-A, HLA -B, and HLA -C class I genes (P(HLA-A-aa-site-62) = 7.4 × 10(-29); P (HLA-B-aa-site-116) = 6.5 × 10(-19); P (HLA-C-aa-site-156) = 6.8 × 10(-8) respectively). Over 250 NPC-HLA associated variants within HLA were analyzed in concert to resolve separate and largely independent HLA-A, -B, and -C gene influences. Multivariate logistical regression analysis collapsed significant associations in adjacent genes spanning 500 kb (OR2H1, GABBR1, HLA-F, and HCG9) as proxies for peptide binding motifs carried by HLA- A*11:01. A similar analysis resolved an independent association signal driven by HLA-B*13:01, B*38:02, and B*55:02 alleles together. NPC resistance alleles carrying the strongly associated amino acid variants implicate specific class I peptide recognition motifs in HLA-A and -B peptide binding groove as conferring strong genetic influence on the development of NPC in China.
Subject(s)
Genome-Wide Association Study , HLA-A Antigens , HLA-B Antigens , Nasopharyngeal Neoplasms , Adult , Aged , Aged, 80 and over , Alleles , Asian People , Carcinoma , China , Female , Genetic Predisposition to Disease , HLA-A Antigens/genetics , HLA-A Antigens/immunology , HLA-B Antigens/genetics , HLA-B Antigens/immunology , HLA-C Antigens/genetics , HLA-C Antigens/immunology , Haplotypes , Herpesvirus 4, Human , Humans , Male , Middle Aged , Nasopharyngeal Carcinoma , Nasopharyngeal Neoplasms/genetics , Nasopharyngeal Neoplasms/immunology , Polymorphism, Single NucleotideABSTRACT
Previous studies have examined germ-line mutations to infer the processes that generate and maintain variability in microsatellite loci. Few studies, however, have examined patterns to infer processes that act on microsatellite loci over evolutionary time. Here, we examine changes in 8 dinucleotide loci across the adaptive radiation of Hawaiian honeycreepers. The loci were found to be highly variable across the radiation, and we did not detect ascertainment bias with respect to allelic diversity or allele size ranges. In examining patterns at the sequence level, we found that changes in flanking regions, repeat motifs, or repeat interruptions were often shared between closely related species and may be phylogenetically informative. Genetic distance measures based on microsatellites were strongly correlated with those based on mitochondrial DNA (mtDNA) sequences as well as with divergence time up to 3 My. Phylogenetic inferences based on microsatellite genetic distances consistently recovered 2 of the 4 honeycreeper clades observed in a tree based on mtDNA sequences but differed from the mtDNA tree in the relationships among clades. Our results confirm that microsatellite loci may be conserved over evolutionary time, making them useful in population-level studies of species that diverged from the species in which they were characterized as long as 5 Ma. Despite this, we found that their use in phylogenetic inference was limited to closely related honeycreeper species.
Subject(s)
Adaptation, Biological/genetics , Evolution, Molecular , Microsatellite Repeats/radiation effects , Passeriformes/genetics , Adaptation, Biological/radiation effects , Animals , DNA, Mitochondrial/genetics , Genetic Variation , Hawaii , Linkage Disequilibrium , Microsatellite Repeats/genetics , PhylogenyABSTRACT
Several immunodeficiencies are associated with high susceptibility to persistent and progressive human papillomavirus (HPV) infection leading to a wide range of cutaneous and mucosal lesions. However, the HPV types most commonly associated with such clinical manifestations in these patients have not been systematically defined. Here, we used virion enrichment, rolling circle amplification, and deep sequencing to identify circular DNA viruses present in skin swabs and/or wart biopsy samples from 48 patients with rare genetic immunodeficiencies, including patients with warts, hypogammaglobulinemia, infections, myelokathexis (WHIM) syndrome, or epidermodysplasia verruciformis (EV). Their profiles were compared with the profiles of swabs from 14 healthy adults and warts from 6 immunologically normal children. Individual patients were typically infected with multiple HPV types; up to 26 different types were isolated from a single patient (multiple anatomical sites, one time point). Among these, we identified the complete genomes of 83 previously unknown HPV types and 35 incomplete genomes representing possible additional new types. HPV types in the genus Gammapapillomavirus were common in WHIM patients, whereas EV patients mainly shed HPVs from the genus Betapapillomavirus. Preliminary evidence based on three WHIM patients treated with plerixafor, a leukocyte mobilizing agent, suggest that longer-term therapy may correlate with decreased HPV diversity and increased predominance of HPV types associated with childhood skin warts.IMPORTANCE Although some members of the viral family Papillomaviridae cause benign skin warts (papillomas), many human papillomavirus (HPV) infections are not associated with visible symptoms. For example, most healthy adults chronically shed Gammapapillomavirus (Gamma) virions from apparently healthy skin surfaces. To further explore the diversity of papillomaviruses, we performed viromic surveys on immunodeficient individuals suffering from florid skin warts. Our results nearly double the number of known Gamma HPV types and suggest that WHIM syndrome patients are uniquely susceptible to Gamma HPV-associated skin warts. Preliminary results suggest that treatment with the drug plerixafor may promote resolution of the unusual Gamma HPV skin warts observed in WHIM patients.
Subject(s)
DNA, Viral/genetics , Immunologic Deficiency Syndromes/complications , Papillomaviridae/classification , Papillomaviridae/isolation & purification , Papillomavirus Infections/virology , Adolescent , Adult , Aged , Aged, 80 and over , Child , Child, Preschool , DNA, Viral/chemistry , Female , Genome, Viral , High-Throughput Nucleotide Sequencing , Humans , Male , Metagenomics , Middle Aged , Mucous Membrane/virology , Nucleic Acid Amplification Techniques , Papillomaviridae/genetics , Skin/virology , Young AdultABSTRACT
Common raven (Corvus corax) populations in Mojave Desert regions of southern California and Nevada have increased dramatically over the past five decades. This growth has been attributed to increased human development in the region, as ravens have a commensal relationship with humans and feed extensively at landfills and on road-killed wildlife. Ravens, as a partially subsidized predator, also represent a problem for native desert wildlife, in particular threatened desert tortoises (Gopherus agassizii). However, it is unclear whether the more than 15-fold population increase is due to in situ population growth or to immigration from adjacent regions where ravens have been historically common. Ravens were sampled for genetic analysis at several local sites within five major areas: the West Mojave Desert (California), East Mojave Desert (southern Nevada), southern coastal California, northern coastal California (Bay Area), and northern Nevada (Great Basin). Analyses of mtDNA control region sequences reveal an increased frequency of raven 'Holarctic clade' haplotypes from south to north inland, with 'California clade' haplotypes nearly fixed in the California populations. There was significant structuring among regions for mtDNA, with high F(ST) values among sampling regions, especially between the Nevada and California samples. Analyses of eight microsatellite loci reveal a mostly similar pattern of regional population structure, with considerably smaller, but mostly significant, values. The greater mtDNA divergences may be due to lower female dispersal relative to males, lower N(e), or effects of high mutation rates on maximal values of F(ST). Analyses indicate recent population growth in the West Mojave Desert and a bottleneck in the northern California populations. While we cannot rule out in situ population growth as a factor, patterns of movement inferred from our data suggest that the increase in raven populations in the West Mojave Desert resulted from movements from southern California and the Central Valley. Ravens in the East Mojave Desert are more similar to ones from northern Nevada, indicating movement between those regions. If this interpretation of high gene flow into the Mojave Desert is correct, then efforts to manage raven numbers by local control may not be optimally effective.
Subject(s)
Crows/genetics , Demography , Genetics, Population , Animals , Base Sequence , California , Crows/physiology , DNA, Mitochondrial/genetics , Desert Climate , Gene Flow , Haplotypes/genetics , Microsatellite Repeats/genetics , Molecular Sequence Data , Nevada , Population Dynamics , Sequence Analysis, DNAABSTRACT
The native Hawaiian honeycreepers represent a classic example of adaptive radiation and speciation, but currently face one the highest extinction rates in the world. Although multiple factors have likely influenced the fate of Hawaiian birds, the relatively recent introduction of avian malaria is thought to be a major factor limiting honeycreeper distribution and abundance. We have initiated genetic analyses of class II beta chain Mhc genes in four species of honeycreepers using methods that eliminate the possibility of sequencing mosaic variants formed by cloning heteroduplexed polymerase chain reaction products. Phylogenetic analyses group the honeycreeper Mhc sequences into two distinct clusters. Variation within one cluster is high, with dN > dS and levels of diversity similar to other studies of Mhc (B system) genes in birds. The second cluster is nearly invariant and includes sequences from honeycreepers (Fringillidae), a sparrow (Emberizidae) and a blackbird (Emberizidae). This highly conserved cluster appears reminiscent of the independently segregating Rfp-Y system of genes defined in chickens. The notion that balancing selection operates at the Mhc in the honeycreepers is supported by transpecies polymorphism and strikingly high dN/dS ratios at codons putatively involved in peptide interaction. Mitochondrial DNA control region sequences were invariant in the i'iwi, but were highly variable in the 'amakihi. By contrast, levels of variability of class II beta chain Mhc sequence codons that are hypothesized to be directly involved in peptide interactions appear comparable between i'iwi and 'amakihi. In the i'iwi, natural selection may have maintained variation within the Mhc, even in the face of what appears to a genetic bottleneck.