Pyrethroid genetic resistance in the dengue vector (Aedes aegypti) in Posadas, Argentina.

Pyrethroids are extensively used to control adult populations of the arboviral vector Aedes aegypti, raising concerns regarding the increasing frequency and distribution of insecticide resistance mutations (kdr: knock-down resistance) in the voltage-gated sodium channel gene (Nav). The widespread use of pyrethroids imposes a threat to the success of mosquito control and the environment. In this study, we investigated the presence of two kdr mutations (V1016I and F1534C) in the Nav gene and their distribution across four neighborhoods in Posadas, Argentina, with different Ae. aegypti abundance and contrasting socioeconomic status (SES). Alleles at each locus were interrogated using TaqMan SNP genotyping assays in DNA extracted from adult females collected in a longitudinal study. We report the presence of both pyrethroid resistance alleles (kdr 1016I = 29.08%; kdr 1534C = 70.70%) among adult females. The frequency of combined kdr genotypes reveals that approximately 70% of local adult females have enhanced resistance to pyrethroids. Both, the proportion of resistant adult females (with at least one kdr allele in each locus) and Ae. aegypti abundance showed an uneven distribution between neighborhoods with different SES (p < 0.001). In high-SES neighborhoods, we found more mosquitoes and a higher frequency of pyrethroid resistance, possibly as a consequence of different public health interventions, social habits, and insecticide use. This is the first report of kdr mutations in Ae. Aegypti in the northeast region of Argentina. Our results focus on the need for within-population (city) distribution analyses of kdr mutations and highlight the relevance of incorporating insecticide resistance monitoring within the Integrated Vector Management initiative.

Determinants of Aedes mosquito density as an indicator of arbovirus transmission risk in three sites affected by co-circulation of globally spreading arboviruses in Colombia, Ecuador and Argentina.

BACKGROUND: The global impact of Zika virus in Latin America has drawn renewed attention to circulating mosquito-borne viruses in this region, such as dengue and chikungunya. Our objective was to assess socio-ecological factors associated with Aedes mosquito vector density as a measure of arbovirus transmission risk in three cities of potentially recent Zika virus introduction: Ibagué, Colombia; Manta, Ecuador; and Posadas, Argentina, in order to inform disease mitigation strategies. METHODS: We sampled Aedes mosquito populations in a total of 1086 households, using indoor and peridomestic mosquito collection methods, including light traps, resting traps, traps equipped with chemical attractant and aspirators. For each sampled household, we collected socio-economic data using structured questionnaires and data on microenvironmental conditions using iButton data loggers. RESULTS: A total of 3230 female Aedes mosquitoes were collected, of which 99.8% were Aedes aegypti and 0.2% were Aedes albopictus. Mean female Aedes mosquito density per household was 1.71 (standard deviation: 2.84). We used mixed-effects generalized linear Poisson regression analyses to identify predictors of Aedes density, using month, neighborhood and country as random-effects variables. Across study sites, the number of household occupants [incidence rate ratio (IRR): 1.08, 95% confidence interval (CI): 1.01-1.14], presence of entry points for mosquitoes into the household (IRR: 1.51, 95% CI: 1.30-1.76) and presence of decorative vegetation (IRR: 1.52, 95% CI: 1.22-1.88) were associated with higher Aedes density; while being in the highest wealth tertile of household wealth (IRR: 0.78, 95% CI: 0.66-0.92), knowledge of how arboviruses are transmitted (IRR: 0.94, 95% CI: 0.89-1.00) and regular emptying of water containers by occupants (IRR: 0.79, 95% CI: 0.67-0.92) were associated with lower Aedes density. CONCLUSIONS: Our study addresses the complexities of arbovirus vectors of global significance at the interface between human and mosquito populations. Our results point to several predictors of Aedes mosquito vector density in countries with co-circulation of multiple Aedes-borne viruses, and point to modifiable risk factors that may be useful for disease prevention and control.

CRISPR-based platform for carbapenemases and emerging viruses detection using Cas12a (Cpf1) effector nuclease.

CRISPR-Cas12a (also called Cpf1) has been commonly used for genomic editing, based on its ability to generate precise double-stranded DNA (dsDNA) breaks. Recently, it was demonstrated that Cas12a exhibits unspecific ssDNAse activity upon target recognition. This feature allows CRISPR-Cas to be coupled with a ssDNA reporter and generate a fast, accurate and ultrasensitive molecular detection method. Here, we demonstrate that Cas12a was able to detect DNA target sequences corresponding to carbapenemases resistance genes such as KPC, NDM and OXA. Also, with the addition of a reverse-transcription step, we were able to detect viral RNA sequences from DENV, ZIKV and HANTV genomes. In all cases, assay run time was less than two hours. Additionally, we report attomolar levels of detection. This methodology was validated using clinical samples from patients infected with Dengue virus. Reactions were visualized by detection of a fluorescent signal, as well as by the use of a simple lateral flow strip. These results indicate that Cas12a is able to detect both DNA and RNA targets, making it an appropriate and convenient tool to detect all types of pathogens.

Yerba mate (Ilex paraguariensis, A. St.-Hil.) de novo transcriptome assembly based on tissue specific genomic expression profiles.

BACKGROUND: The most common infusion in southern Latin-American countries is prepared with dried leaves of Ilex paraguariensis A. St.-Hil., an aboriginal ancestral beverage known for its high polyphenols concentration currently consumed in > 90% of homes in Argentina, in Paraguay and Uruguay. The economy of entire provinces heavily relies on the production, collection and manufacture of Ilex paraguariensis, the fifth plant species with highest antioxidant activity. Polyphenols are associated to relevant health benefits including strong antioxidant properties. Despite its regional relevance and potential biotechnological applications, little is known about functional genomics and genetics underlying phenotypic variation of relevant traits. By generating tissue specific transcriptomic profiles, we aimed to comprehensively annotate genes in the Ilex paraguariensis phenylpropanoid pathway and to evaluate differential expression profiles. RESULTS: In this study we generated a reliable transcriptome assembly based on a collection of 15 RNA-Seq libraries from different tissues of Ilex paraguariensis. A total of 554 million RNA-Seq reads were assembled into 193,897 transcripts, where 24,612 annotated full-length transcripts had complete ORF. We assessed the transcriptome assembly quality, completeness and accuracy using BUSCO and TransRate; consistency was also evaluated by experimentally validating 11 predicted genes by PCR and sequencing. Functional annotation against KEGG Pathway database identified 1395 unigenes involved in biosynthesis of secondary metabolites, 531 annotated transcripts corresponded to the phenylpropanoid pathway. The top 30 differentially expressed genes among tissue revealed genes involved in photosynthesis and stress response. These significant differences were then validated by qRT-PCR. CONCLUSIONS: Our study is the first to provide data from whole genome gene expression profiles in different Ilex paraguariensis tissues, experimentally validating in-silico predicted genes key to the phenylpropanoid (antioxidant) pathway. Our results provide essential genomic data of potential use in breeding programs for polyphenol content. Further studies are necessary to assess if the observed expression variation in the phenylpropanoid pathway annotated genes is related to variations in leaves' polyphenol content at the population scale. These results set the current reference for Ilex paraguariensis genomic studies and provide a substantial contribution to research and biotechnological applications of phenylpropanoid secondary metabolites.

Genome-wide scan for commons SNPs affecting bovine leukemia virus infection level in dairy cattle.

BACKGROUND: Bovine leukemia virus (BLV) infection is omnipresent in dairy herds causing direct economic losses due to trade restrictions and lymphosarcoma-related deaths. Milk production drops and increase in the culling rate are also relevant and usually neglected. The BLV provirus persists throughout a lifetime and an inter-individual variation is observed in the level of infection (LI) in vivo. High LI is strongly correlated with disease progression and BLV transmission among herd mates. In a context of high prevalence, classical control strategies are economically prohibitive. Alternatively, host genomics studies aiming to dissect loci associated with LI are potentially useful tools for genetic selection programs tending to abrogate the viral spreading. The LI was measured through the proviral load (PVL) set-point and white blood cells (WBC) counts. The goals of this work were to gain insight into the contribution of SNPs (bovine 50KSNP panel) on LI variability and to identify genomics regions underlying this trait. RESULTS: We quantified anti-p24 response and total leukocytes count in peripheral blood from 1800 cows and used these to select 800 individuals with extreme phenotypes in WBCs and PVL. Two case-control genomic association studies using linear mixed models (LMMs) considering population stratification were performed. The proportion of the variance captured by all QC-passed SNPs represented 0.63 (SE ± 0.14) of the phenotypic variance for PVL and 0.56 (SE ± 0.15) for WBCs. Overall, significant associations (Bonferroni's corrected -log10p > 5.94) were shared for both phenotypes by 24 SNPs within the Bovine MHC. Founder haplotypes were used to measure the linkage disequilibrium (LD) extent (r2 = 0.22 ± 0.27 at inter-SNP distance of 25-50 kb). The SNPs and LD blocks indicated genes potentially associated with LI in infected cows: i.e. relevant immune response related genes (DQA1, DRB3, BOLA-A, LTA, LTB, TNF, IER3, GRP111, CRISP1), several genes involved in cell cytoskeletal reorganization (CD2AP, PKHD1, FLOT1, TUBB5) and modelling of the extracellular matrix (TRAM2, TNXB). Host transcription factors (TFs) were also highlighted (TFAP2D; ABT1, GCM1, PRRC2A). CONCLUSIONS: Data obtained represent a step forward to understand the biology of BLV-bovine interaction, and provide genetic information potentially applicable to selective breeding programs.

Molecular characterization of yerba mate chlorosis-associated virus, a putative cytorhabdovirus infecting yerba mate (Ilex paraguariensis).

We present the molecular characterization of a new virus infecting yerba mate (Ilex paraguariensis St. Hil.) in Argentina. Deep sequencing of diseased yerba mate plants showing chlorotic linear patterns, chlorotic rings, and vein yellowing resulted in the identification of a new virus resembling plant rhabdoviruses in sequence and genome structure. We have determined the complete genome sequence of this virus, which is 12,876 nt long. Seven open reading frames (ORFs) were identified in the antigenomic orientation of the negative-sense, single-stranded viral RNA, in the order 3'-N-P-P3-P4-M-G-L-5'. Phylogenetic analysis suggested that the described virus is a new member of the genus Cytorhabdovirus, which was supported by the observation of rhabdovirus-like particles within the cytoplasm of infected yerba mate cells. The virus has been tentatively named "yerba mate chlorosis-associated virus" (YmCaV). The availability of the YmCaV genome sequence will contribute to assessing the genetic variability of this virus and determining its role in this yerba mate disease.

A unifying study of phenotypic and molecular genetic variability in natural populations of Anadenanthera colubrina var. cebil from Yungas and Paranaense biogeographic provinces in Argentina.

Anadenanthera colubrina var. cebil is a discontinuously distributed native tree species in South American subtropical forests. Thirteen quantitative traits and eight nuclear microsatellite loci were examined in individuals from two biogeographic provinces of Argentina to determine the number and composition of genetically distinguishable groups of individuals and explore possible spatial patterns of the phenotypic and genetic variability. Means of reproductive traits were higher in the Yungas than in the Paranaense biogeographic province, whereas five out of eight nonreproductive quantitative traits showed higher mean values in the latter. Variance coefficients were moderate, and there were significant differences between and within provinces. Three clusters were defined based on spatial model for cluster membership for quantitative traits. One cluster grouped the individuals from the Paranaense biogeographic province whereas the individuals from the Yungas biogeographic province grouped regarding its population of origin. Parameters of molecular genetic variability showed higher values in the Yungas than in the Paranaense biogeographic province. Observed heterozygosity was lower than expected heterozygosity in both biogeographic provinces, indicating an excess of homozygosity. The homozygosity test by Watterson and the exact test by Slatkin suggested diversifying selection for locus Ac41.1. Bayesian clustering spatial model for microsatellites loci data were performed for both all loci and for all loci excluding locus Ac41.1. In both analyses two clusters were inferred. Analysis of molecular variance revealed similar results for all genotypes and for all genotypes defined excluding locus Ac41.1. Most of the total variance is attributable to genetic variation within clusters. The presence of homogeneous clusters was detected for both the phenotypic and molecular genetic variability. Two Bayesian clustering analyses were performed according to molecular genetic data, and two clusters were inferred. Individuals were assigned to their provinces of origin. Genetic molecular variation was higher in the populations of the Yungas biogeographic province which translates into highly qualified populations for conservation. Populations from the Paranaense biogeographic province showed the highest mean value of number of seeds per fruit making them valuable as well with regard to the exploitation of management strategies as a means to recover the impacted areas where these populations are located.

Human platelet antigen typing of neonatal alloimmune thrombocytopenia patients using whole genome amplified DNA and a 5'-nuclease assay.

BACKGROUND: A serious constraint in the investigation of the human platelet antigen (HPA) status of potential neonatal alloimmune thrombocytopenia (NAIT) cases is the limited amount of DNA available from the neonates. Whole genome amplification (WGA) of these DNA samples could overcome this problem, but requires validation to ensure that it is sufficiently sensitive and accurate before its application in a clinical diagnostic setting. STUDY DESIGN AND METHODS: This study has validated the use of WGA DNA for HPA-1, -2, -3, -4, -5, and -15 genotyping with a panel of six controls and 13 previously HPA-typed samples from neonates together with parental DNA, using a 5'-nuclease (TaqMan) assay. WGA was performed using titrated amounts of genomic and WGA DNA template. HPA typing was performed on genomic and amplified DNA using a 5'-nuclease assay or polymerase chain reaction with sequence-specific primers (PCR-SSP). RESULTS: WGA DNA yields were in the suggested range of 400x to 800x, as assessed by spectrophotometry and gel analysis, and did not require further purification. HPA genotyping showed 100 percent concordance when using down to 5 ng of genomic or WGA template. CONCLUSION: This study demonstrates that WGA can be used for HPA typing using PCR-SSP or plate-based 5'-nuclease assays. The use of WGA for HPA typing in clinical samples from NAIT patients was validated with 100 percent concordance, and it is suggested that this technology can be used for other analyses where DNA amounts are limited.

Association of smoking behavior with an odorant receptor allele telomeric to the human major histocompatibility complex.

Smoking behavior has been associated in two independent European cohorts with the most common Caucasian human leukocyte antigen (HLA) haplotype (A1-B8-DR3). We aimed to test whether polymorphic members of the two odorant receptor (OR) clusters within the extended HLA complex might be responsible for the observed association, by genotyping a cohort of Hungarian women in which the mentioned association had been found. One hundred and eighty HLA haplotypes from Centre d'Etude du Polymorphisme Humain families were analyzed in silico to identify single-nucleotide polymorphisms (SNPs) within OR genes that are in linkage disequilibrium with the A1-B8-DR3 haplotype, as well as with two other haplotypes indirectly linked to smoking behavior. A nonsynonymous SNP within the OR12D3 gene (rs3749971(T)) was found to be linked to the A1-B8-DR3 haplotype. This polymorphism leads to a (97)Thr --> Ile exchange that affects a putative ligand binding region of the OR12D3 protein. Smoking was found to be associated in the Hungarian cohort with the rs3749971(T) allele (p = 1.05 x 10(-2)), with higher significance than with A1-B8-DR3 (p = 2.38 x 10(-2)). Our results link smoking to a distinct OR allele, and demonstrate that the rs3749971(T) polymorphism is associated with the HLA haplotype-dependent differential recognition of cigarette smoke components, at least among Caucasian women.

A Bayesian deconvolution strategy for immunoprecipitation-based DNA methylome analysis.

DNA methylation is an indispensible epigenetic modification required for regulating the expression of mammalian genomes. Immunoprecipitation-based methods for DNA methylome analysis are rapidly shifting the bottleneck in this field from data generation to data analysis, necessitating the development of better analytical tools. In particular, an inability to estimate absolute methylation levels remains a major analytical difficulty associated with immunoprecipitation-based DNA methylation profiling. To address this issue, we developed a cross-platform algorithm-Bayesian tool for methylation analysis (Batman)-for analyzing methylated DNA immunoprecipitation (MeDIP) profiles generated using oligonucleotide arrays (MeDIP-chip) or next-generation sequencing (MeDIP-seq). We developed the latter approach to provide a high-resolution whole-genome DNA methylation profile (DNA methylome) of a mammalian genome. Strong correlation of our data, obtained using mature human spermatozoa, with those obtained using bisulfite sequencing suggest that combining MeDIP-seq or MeDIP-chip with Batman provides a robust, quantitative and cost-effective functional genomic strategy for elucidating the function of DNA methylation.

Generation of a genomic tiling array of the human major histocompatibility complex (MHC) and its application for DNA methylation analysis.

BACKGROUND: The major histocompatibility complex (MHC) is essential for human immunity and is highly associated with common diseases, including cancer. While the genetics of the MHC has been studied intensively for many decades, very little is known about the epigenetics of this most polymorphic and disease-associated region of the genome. METHODS: To facilitate comprehensive epigenetic analyses of this region, we have generated a genomic tiling array of 2 Kb resolution covering the entire 4 Mb MHC region. The array has been designed to be compatible with chromatin immunoprecipitation (ChIP), methylated DNA immunoprecipitation (MeDIP), array comparative genomic hybridization (aCGH) and expression profiling, including of non-coding RNAs. The array comprises 7832 features, consisting of two replicates of both forward and reverse strands of MHC amplicons and appropriate controls. RESULTS: Using MeDIP, we demonstrate the application of the MHC array for DNA methylation profiling and the identification of tissue-specific differentially methylated regions (tDMRs). Based on the analysis of two tissues and two cell types, we identified 90 tDMRs within the MHC and describe their characterisation. CONCLUSION: A tiling array covering the MHC region was developed and validated. Its successful application for DNA methylation profiling indicates that this array represents a useful tool for molecular analyses of the MHC in the context of medical genomics.

An integrated resource for genome-wide identification and analysis of human tissue-specific differentially methylated regions (tDMRs).

We report a novel resource (methylation profiles of DNA, or mPod) for human genome-wide tissue-specific DNA methylation profiles. mPod consists of three fully integrated parts, genome-wide DNA methylation reference profiles of 13 normal somatic tissues, placenta, sperm, and an immortalized cell line, a visualization tool that has been integrated with the Ensembl genome browser and a new algorithm for the analysis of immunoprecipitation-based DNA methylation profiles. We demonstrate the utility of our resource by identifying the first comprehensive genome-wide set of tissue-specific differentially methylated regions (tDMRs) that may play a role in cellular identity and the regulation of tissue-specific genome function. We also discuss the implications of our findings with respect to the regulatory potential of regions with varied CpG density, gene expression, transcription factor motifs, gene ontology, and correlation with other epigenetic marks such as histone modifications.

Variation analysis and gene annotation of eight MHC haplotypes: the MHC Haplotype Project.

The human major histocompatibility complex (MHC) is contained within about 4 Mb on the short arm of chromosome 6 and is recognised as the most variable region in the human genome. The primary aim of the MHC Haplotype Project was to provide a comprehensively annotated reference sequence of a single, human leukocyte antigen-homozygous MHC haplotype and to use it as a basis against which variations could be assessed from seven other similarly homozygous cell lines, representative of the most common MHC haplotypes in the European population. Comparison of the haplotype sequences, including four haplotypes not previously analysed, resulted in the identification of >44,000 variations, both substitutions and indels (insertions and deletions), which have been submitted to the dbSNP database. The gene annotation uncovered haplotype-specific differences and confirmed the presence of more than 300 loci, including over 160 protein-coding genes. Combined analysis of the variation and annotation datasets revealed 122 gene loci with coding substitutions of which 97 were non-synonymous. The haplotype (A3-B7-DR15; PGF cell line) designated as the new MHC reference sequence, has been incorporated into the human genome assembly (NCBI35 and subsequent builds), and constitutes the largest single-haplotype sequence of the human genome to date. The extensive variation and annotation data derived from the analysis of seven further haplotypes have been made publicly available and provide a framework and resource for future association studies of all MHC-associated diseases and transplant medicine.

Germline rates of de novo meiotic deletions and duplications causing several genomic disorders.

Meiotic recombination between highly similar duplicated sequences (nonallelic homologous recombination, NAHR) generates deletions, duplications, inversions and translocations, and it is responsible for genetic diseases known as 'genomic disorders', most of which are caused by altered copy number of dosage-sensitive genes. NAHR hot spots have been identified within some duplicated sequences. We have developed sperm-based assays to measure the de novo rate of reciprocal deletions and duplications at four NAHR hot spots. We used these assays to dissect the relative rates of NAHR between different pairs of duplicated sequences. We show that (i) these NAHR hot spots are specific to meiosis, (ii) deletions are generated at a higher rate than their reciprocal duplications in the male germline and (iii) some of these genomic disorders are likely to have been underascertained clinically, most notably that resulting from the duplication of 7q11, the reciprocal of the deletion causing Williams-Beuren syndrome.

Structure and genetic relationships between Brazilian naturalized and exotic purebred goat domestic goat (Capra hircus) breeds based on microsatellites

The genetic relationships and structure of fourteen goat (Capra hircus) populations were estimated based on genotyping data from 14 goat populations (n = 410 goats) at 13 microsatellite loci. We used analysis of molecular variance (AMOVA), principal component analysis (PCA) and F statistics (F IS, F IT and F ST) to evaluate the genetic diversity (Ho, He and ad) of the goats. Genetic distances between the 14 goat populations were calculated from allelic frequency data for the 13 microsatellite markers. Moderate differentiation was observed for the populations of the undefined breeds (including the Anglo-Nubian-M breed), the naturalized Brazilian breeds (Moxotó, Canindé), the exotic purebred breeds (Alpine, Saanen, Toggenbourg and Anglo-Nubian) and the naturalized Brazilian Graúna group. Our AMOVA showed that a major portion (88.51 percent) of the total genetic variation resulted from differences between individual goats within populations, while between-populations variation accounted for the remaining 11.49 percent of genetic variation. We used a Reynolds genetic distance matrix and PCA to produce a phenogram based on the 14 goat populations and found three clusters, or groups, consisting of the goats belonging to the undefined breed, the naturalized breeds and the exotic purebred breeds. The closer proximity of the Canindé breed from the Brazilian state of Paraíba to the Graúna breed from the same state than to the genetically conserved Canindé breed from the Brazilian state of Ceará, as well as the heterozygosity values and significant deviations from Hardy-Weinberg equilibrium suggests that there was a high number of homozygotes in the populations studied, and indicates the importance of the State for the conservation of the local breeds. Cataloguing the genetic profile of Brazilian goat populations provides essential information for conservation and genetic improvements programs.

A high-resolution HLA and SNP haplotype map for disease association studies in the extended human MHC.

The proteins encoded by the classical HLA class I and class II genes in the major histocompatibility complex (MHC) are highly polymorphic and are essential in self versus non-self immune recognition. HLA variation is a crucial determinant of transplant rejection and susceptibility to a large number of infectious and autoimmune diseases. Yet identification of causal variants is problematic owing to linkage disequilibrium that extends across multiple HLA and non-HLA genes in the MHC. We therefore set out to characterize the linkage disequilibrium patterns between the highly polymorphic HLA genes and background variation by typing the classical HLA genes and >7,500 common SNPs and deletion-insertion polymorphisms across four population samples. The analysis provides informative tag SNPs that capture much of the common variation in the MHC region and that could be used in disease association studies, and it provides new insight into the evolutionary dynamics and ancestral origins of the HLA loci and their haplotypes.

Immunogenomics: molecular hide and seek.

Similar to other classical science disciplines, immunology has been embracing novel technologies and approaches giving rise to specialised sub-disciplines such as immunogenetics and, more recently, immunogenomics, which, in many ways, is the genome-wide application of immunogenetic approaches. Here, recent progress in the understanding of the immune sub-genome will be reviewed, and the ways in which immunogenomic datasets consisting of genetic and epigenetic variation, linkage disequilibrium and recombination can be harnessed for disease association and evolutionary studies will be discussed. The discussion will focus on data available for the major histocompatibility complex and the leukocyte receptor complex, the two most polymorphic regions of the human immune sub-genome.

Genetic analysis of completely sequenced disease-associated MHC haplotypes identifies shuffling of segments in recent human history.

The major histocompatibility complex (MHC) is recognised as one of the most important genetic regions in relation to common human disease. Advancement in identification of MHC genes that confer susceptibility to disease requires greater knowledge of sequence variation across the complex. Highly duplicated and polymorphic regions of the human genome such as the MHC are, however, somewhat refractory to some whole-genome analysis methods. To address this issue, we are employing a bacterial artificial chromosome (BAC) cloning strategy to sequence entire MHC haplotypes from consanguineous cell lines as part of the MHC Haplotype Project. Here we present 4.25 Mb of the human haplotype QBL (HLA-A26-B18-Cw5-DR3-DQ2) and compare it with the MHC reference haplotype and with a second haplotype, COX (HLA-A1-B8-Cw7-DR3-DQ2), that shares the same HLA-DRB1, -DQA1, and -DQB1 alleles. We have defined the complete gene, splice variant, and sequence variation contents of all three haplotypes, comprising over 259 annotated loci and over 20,000 single nucleotide polymorphisms (SNPs). Certain coding sequences vary significantly between different haplotypes, making them candidates for functional and disease-association studies. Analysis of the two DR3 haplotypes allowed delineation of the shared sequence between two HLA class II-related haplotypes differing in disease associations and the identification of at least one of the sites that mediated the original recombination event. The levels of variation across the MHC were similar to those seen for other HLA-disparate haplotypes, except for a 158-kb segment that contained the HLA-DRB1, -DQA1, and -DQB1 genes and showed very limited polymorphism compatible with identity-by-descent and relatively recent common ancestry (<3,400 generations). These results indicate that the differential disease associations of these two DR3 haplotypes are due to sequence variation outside this central 158-kb segment, and that shuffling of ancestral blocks via recombination is a potential mechanism whereby certain DR-DQ allelic combinations, which presumably have favoured immunological functions, can spread across haplotypes and populations.

A comparison of tagging methods and their tagging space.

Single-nucleotide polymorphism (SNP) tagging is widely used as a way of saving genotyping costs in association studies. A number of different tagging methods have been developed to reduce the number of markers to be genotyped while maintaining power for detecting effects on non-assayed SNPs. How the different methods perform in different settings, the degree to which they overlap and share common tags and how they differ are important questions. We investigated these questions by comparing three widely used tagging methods/algorithms--one haplotype r2-based method, one pair-wise r2-based method and one method which was based on haplotype diversity but focused on major haplotypes. Tagging efficiency was defined as the number of genotyped markers divided by the number of tagging SNPs. Tagging effectiveness was defined as the proportion of un-genotyped or 'hidden' SNPs being detected (having a pair-wise or haplotype r2 with a set of tagging SNPs over a threshold, e.g. haplotype r2> or =0.80). The ENCODE regions genotyped on the HapMap CEPH individuals were examined in this study. Tagging effectiveness was generally poor for rare SNPs than for common SNPs, for all three tagging methods. Inclusion of rare SNPs into initial HapMap scheme could enhance the performance of tags on rare hidden SNPs at the expense of increased genotyping cost. At a moderate tagging efficiency, more than 90% of hidden SNPs detected by tagging SNPs selected by one method were also detected by tagging SNPs selected by another method, and this figure could be increased to 100% if tagging efficiency was allowed to drop. These results indicate that the tagging space is highly concordant between different tagging methods, despite the fact that they often involve different sets of tagging SNPs.

A high-resolution linkage-disequilibrium map of the human major histocompatibility complex and first generation of tag single-nucleotide polymorphisms.

Autoimmune, inflammatory, and infectious diseases present a major burden to human health and are frequently associated with loci in the human major histocompatibility complex (MHC). Here, we report a high-resolution (1.9 kb) linkage-disequilibrium (LD) map of a 4.46-Mb fragment containing the MHC in U.S. pedigrees with northern and western European ancestry collected by the Centre d'Etude du Polymorphisme Humain (CEPH) and the first generation of haplotype tag single-nucleotide polymorphisms (tagSNPs) that provide up to a fivefold increase in genotyping efficiency for all future MHC-linked disease-association studies. The data confirm previously identified recombination hotspots in the class II region and allow the prediction of numerous novel hotspots in the class I and class III regions. The region of longest LD maps outside the classic MHC to the extended class I region spanning the MHC-linked olfactory-receptor gene cluster. The extended haplotype homozygosity analysis for recent positive selection shows that all 14 outlying haplotype variants map to a single extended haplotype, which most commonly bears HLA-DRB1*1501. The SNP data, haplotype blocks, and tagSNPs analysis reported here have been entered into a multidimensional Web-based database (GLOVAR), where they can be accessed and viewed in the context of relevant genome annotation. This LD map allowed us to give coordinates for the extremely variable LD structure underlying the MHC.