Large-scale MHC class II genotyping of a wild lemur population by next generation sequencing

The critical role of major histocompatibility complex (MHC) genes in disease resistance, along with their putative function in sexual selection, reproduction and chemical ecology, make them an important genetic system in evolutionary ecology. Studying selective pressures acting on MHC genes in the wild nevertheless requires population-wide genotyping, which has long been challenging because of their extensive polymorphism. Here, we report on large-scale genotyping of the MHC class II loci of the grey mouse lemur (Microcebus murinus) from a wild population in western Madagascar. The second exons from MHC-DRB and -DQB of 772 and 672 individuals were sequenced, respectively, using a 454 sequencing platform, generating more than 800,000 reads. Sequence analysis, through a stepwise variant validation procedure, allowed reliable typing of more than 600 individuals. The quality of our genotyping was evaluated through three independent methods, namely genotyping the same individuals by both cloning and 454 sequencing, running duplicates, and comparing parent–offspring dyads; each displaying very high accuracy. A total of 61 (including 20 new) and 60 (including 53 new) alleles were detected at DRB and DQB genes, respectively. Both loci were non-duplicated, in tight linkage disequilibrium and in Hardy–Weinberg equilibrium, despite the fact that sequence analysis revealed clear evidence of historical selection. Our results highlight the potential of 454 sequencing technology in attempts to investigate patterns of selection shaping MHC variation in contemporary populations. The power of this approach will nevertheless be conditional upon strict quality control of the genotyping data.     

Efficacy and limits of genotyping low copy number (LCN) DNA samples by multiplex PCR of STR loci

In this study, we have evaluated the efficacy and the validity of the AmpFISTR SGM plus multiplex PCR typing system when Low Copy Number (LCN) amounts of DNA are processed. The characteristics of SGM plus profiles produced under LCN conditions were studied on the basis of heterozygote balance, between loci balance and stutter proportion based on profiles that were obtained from a variety of mock casework samples. These experiments clearly showed that LCN DNA profiles carry their own characteristic features, which must be taken into account during interpretation. Herewith, we confirmed the data of recent other studies that a comprehensive interpretation strategy is dependent upon multiple replication of the PCR using the same extract together with the proper use of extraction and amplification controls. The limitations of LCN DNA analysis were further studied in a series of single cell PCR experiments using an amplification regime of 34 PCR cycles. The allele dropout phenomenon was demonstrated to its full extent when single cells were analysed. However, the "consensus profile" which was obtained from separate single cell PCR experiments matched the actual profile of the cell donor. Single cell PCR experiments also showed that a further increase of the number of PCR cycles did not result in enhanced sensitivity and had a highly negative effect on the balance of this multiplex PCR system which hampered correct interpretation of the profile. Also, the potential of LCN typing in analysing mixtures of DNA was investigated. It was clearly shown that LCN typing had no advantages over 28 cycles amplification in the detection of the minor component of DNA-mixtures. In addition to the 34 cycles PCR amplification regime, the utility of a new approach that involved reamplification of the 28 cycle SGM plus PCR products with an extra 6 PCR cycles after the addition of fresh AmpliTaq Gold DNA Polymerase was investigated. This approach provides the scientist with an extra typing result that enhances the reliability of the consensus profile, which is commonly retrieved from two separate 34 cycle PCR results. Furthermore, the 28 + 6 cycles approach may be used to screen LCN samples for their potential to produce a 34 PCR cycle profile. Finally and as a last resort the 28 + 6 cycles approach can be used in those cases where no further extract from the crime sample is available. Finally, the potential of LCN typing was demonstrated in typing samples from non-probative and actual casework examples. From a high proportion of samples that failed to demonstrate SGM plus typing results using the standard protocol of 28 cycles, at least partial profiles could be obtained after LCN methods were used. For example, LCN typing was applied in a case where 10-year old samples from bones and teeth that were retrieved from a mass grave had to be identified. This study resulted in the positive identification of a number of victims by comparing the LCN DNA profiles with the profiles from putative relatives. The value of LCN DNA typing was further demonstrated in a strangulation case. The throat of the victim was sampled and only after 34 PCR cycles were we able to reveal that the evidential sample contained a distinct mixture of the victim's own DNA and the DNA of the defendant.     

Design of a multiplex PCR for genotyping 16 short tandem repeats in degraded DNA samples

The molecular genotyping of individuals and reconstruction of kinship through short and high polymorphic DNA markers, so-called short tandem repeats (STR), has become an important and efficient method in anthropology and forensic science. The here introduced experimental design describes a multiplex PCR capable of simultaneously amplifying 16 STRs and the sex determinant locus amelogenin in a short fragment lengths range from 84 bp to 275 bp. Thus, the design depends predominantly on the routines for DNA typing of historical samples with highly degraded ancient DNA. It is shown, that the newly designed multiplex PCR is suitable for successful typing of both forensic and historical material.     

Utilizing next generation sequencing to characterize microsatellite loci in a tropical aquatic plant species Cryptocoryne cordata var. cordata (Araceae)

Cryptocoryne cordata var. cordata (2n = 34) is an aquatic plant species distributed from the southern part of Peninsular Thailand through the Malay Peninsula. It propagates both sexually and asexually via stolons. The current study is aimed at developing nuclear microsatellite markers for the species using next generation sequencing (Roche 454 pyrosequencing) from genomic DNA. A total of 41,653 reads was generated, of which 3636 fragments contained at least one repeat motif. Seventy two primer sets in the flanking region of dinucleotide, trinucleotide and tetranucleotides repeat motifs were designed and tested for efficiency in polymerase chain reaction amplification. Using these primer sets, 11 new microsatellite marker loci were successfully amplified with unambiguous polymorphic alleles exhibited across 30 individuals examined. The number of alleles per locus ranged from 2 to 6, while observed and expected heterozygosity ranged from 0.8190 to 1.0000 and 0.5401 to 0.7548, respectively. Genotype frequencies at all loci departed significantly from Hardy–Weinberg equilibrium. Linkage disequilibrium was not detected between any pair of loci. Cross-species amplification was successful across a panel of ten Cryptocoryne species. The markers described in this study will be useful for evaluating genetic diversity within and between populations, levels of gene flow, and the population dynamics of clones. They will be of further value in the development of effective conservation programs for Cryptocoryne species.     

Microsatellite-directed selection of breeders for the next backcross generation by using a minimal number of loci

To establish a minimal number of markers for direct selection of candidate mice used for the next mating to produce congenic mice, recombination frequencies of 53 microsatellite loci on chromosomes (chr.) 1 and 19 were examined using 41 N2 mice: the donor strain was BALB/c, and recipient strain was C57BL/6J (B6J) or C57BL/6N (B6N). These markers were spaced at 0.1 to 24.2 centimorgans (cM). Among the 41 mice, B6/B6 homozygosity ranged from 18 to 24 animals (mean, 20; 2 standard deviations, 1.36) for a given locus. There was no difference in recombination frequency between chr. 1 and 19. The recombination frequency of B6J was higher than that of B6N (P < 0.05). Various densities of markers, 10 (5 markers/chr.), 8 (4 markers/chr.), and 6 (3 markers/chr.) spaced at 12.0 to 29.3, 9.0 to 45.0, and 24.5 to 53.0 cM, respectively, were selected from the 53 markers, and homozygosity was compared in each mouse. In mice with decreased homozygosity when tested using 53 markers, homozygosity differed depending on the density of the markers. The results suggested that 3 markers/chr. are sufficient for selection of the highest percentages of homozygosity but are not suitable to define mice with lower percentages of homozygosity.     

Exome sequencing followed by large-scale genotyping fails to identify single rare variants of large effect in idiopathic generalized epilepsy

Idiopathic generalized epilepsy (IGE) is a complex disease with high heritability, but little is known about its genetic architecture. Rare copy-number variants have been found to explain nearly 3% of individuals with IGE; however, it remains unclear whether variants with moderate effect size and frequencies below what are reliably detected with genome-wide association studies contribute significantly to disease risk. In this study, we compare the exome sequences of 118 individuals with IGE and 242 controls of European ancestry by using next-generation sequencing. The exome-sequenced epilepsy cases include study subjects with two forms of IGE, including juvenile myoclonic epilepsy (n = 93) and absence epilepsy (n = 25). However, our discovery strategy did not assume common genetic control between the subtypes of IGE considered. In the sequence data, as expected, no variants were significantly associated with the IGE phenotype or more specific IGE diagnoses. We then selected 3,897 candidate epilepsy-susceptibility variants from the sequence data and genotyped them in a larger set of 878 individuals with IGE and 1,830 controls. Again, no variant achieved statistical significance. However, 1,935 variants were observed exclusively in cases either as heterozygous or homozygous genotypes. It is likely that this set of variants includes real risk factors. The lack of significant association evidence of single variants with disease in this two-stage approach emphasizes the high genetic heterogeneity of epilepsy disorders, suggests that the impact of any individual single-nucleotide variant in this disease is small, and indicates that gene-based approaches might be more successful for future sequencing studies of epilepsy predisposition.     

Genome-Tagged Amplification (GTA): a PCR-based method to prepare sample-tagged amplicons from hundreds of individuals for next generation sequencing

Sampling the sequence of a relatively small fraction of the genome in large numbers of individuals is an important objective for population genetics and association genetics approaches. However, currently available ‘sequence capture’ methods either require expensive instrumentation or have problems dealing with high sample numbers and relatively small target sizes. We have developed Genome-Tagged Amplification (GTA) as a flexible PCR-based method for preparing pools of hundreds of amplicons from hundreds of samples for next generation sequencing. The method involves tagging of genomic DNA with barcode adapters at restriction sites, followed by PCR amplification from flanking DNA. It is freely scalable for both sample number and amplicon number and has no specialized equipment requirement. An optimized protocol is presented which provides a matrix of 96 × 192 combinations of samples x amplicons, corresponding to a complete 454 Titanium run. Initially, we used 454 sequencing; however, GTA could easily be adapted to Illumina sequencing platforms as read lengths have significantly increased in this system.     

Application of a multiplex PCR for the detection of protozoan pathogens of the eastern oyster Crassostrea virginica in field samples

Populations of eastern oysters Crassostrea virginica along the east coast of North America have repeatedly experienced epizootic mass mortality due to infections by protozoan parasites, and molecular diagnostic methodologies are fast becoming more widely available for the diagnosis of protozoan diseases of oysters. In this study we applied a modified version of an existing multiplex polymerase chain reaction (PCR) for detection of the eastern oyster parasites Haplosporidium nelsoni, H. costale and Perkinsus marinus from field-collected samples. We incorporated primers for DNA quality control based on the large subunit ribosomal RNA (LSU rRNA) gene of C. virginica. The multiplex PCR (MPCR) simultaneously amplified genomic DNA of C. virginica, and cloned DNA of H. nelsoni, P. marinus and H. costale. In field trial applications, we compared the performance of the MPCR to that of the conventional diagnostic techniques of histopathological tissue examination and the Ray/Mackin fluid thioglycollate medium (RMFT) assay. A total of 530 oysters were sampled from 18 sites at 12 locations along the east coast of the United States from the Gulf of Mexico to southern New England. The modified MPCR detected 21% oysters with H. nelsoni, 2% oysters with H. costale, and 40% oysters with P. marinus infections. In comparison, histopathological examination detected H. nelsoni and H. costale infections in 6 and 0.8% oysters, respectively, and the RMFT assay detected P. marinus infection in 31% oysters. The MPCR is a more sensitive diagnostic assay for detection of H. nelsoni, H. costale, and P. marinus, and incorporation of an oyster quality control product limits false negative results.     

Development of a real-time multiplex PCR assay for the detection of multiple <Emphasis Type="Italic">Salmonella</Emphasis> serotypes in chicken samples

Background  

A real-time multiplex PCR assay was developed for the detection of multiple Salmonella serotypes in chicken samples. Poultry-associated serotypes detected in the assay include Enteritidis, Gallinarum, Typhimurium, Kentucky and Dublin. The traditional cultural method according to EN ISO 6579:2002 for the detection of Salmonella in food was performed in parallel. The real-time PCR based method comprised a pre-enrichment step in Buffered Peptone Water (BPW) overnight, followed by a shortened selective enrichment in Rappaport Vasilliadis Soya Broth (RVS) for 6 hours and subsequent DNA extraction.     

FLAG assay as a novel method for real-time signal generation during PCR: application to detection and genotyping of KRAS codon 12 mutations

Real-time signal generation methods for detection and characterization of low-abundance mutations in genomic DNA are powerful tools for cancer diagnosis and prognosis. Mutations in codon 12 of the oncogene KRAS, for example, are frequently found in several types of human cancers. We have developed a novel real-time PCR technology, FLAG (FLuorescent Amplicon Generation) and adapted it for simultaneously (i) amplifying mutated codon 12 KRAS sequences, (ii) monitoring in real-time the amplification and (iii) genotyping the exact nucleotide alteration. FLAG utilizes the exceptionally thermostable endonuclease PspGI for real-time signal generation by cleavage of quenched fluorophores from the 5′-end of the PCR products and, concurrently, for selecting KRAS mutations over wild type. By including peptide-nucleic-acid probes in the reaction, simultaneous genotyping is achieved that circumvents the requirement for sequencing. FLAG enables high-throughput, closed-tube KRAS mutation detection down to ~0.1% mutant-to-wild type. The assay was validated on model systems and compared with allele-specific PCR sequencing for screening 27 cancer specimens. Diverse applications of FLAG for real-time PCR or genotyping applications in cancer, virology or infectious diseases are envisioned.     

Adaptation genomics: next generation sequencing reveals a shared haplotype for rapid early development in geographically and genetically distant populations of rainbow trout

Local adaptation occurs when a population evolves a phenotype that confers a selective advantage in its local environment, but which may not be advantageous in other habitats. Restricted gene flow and strong selection pressures are prerequisites for local adaptation. Fishes in the family Salmonidae are predicted to provide numerous examples of local adaptation because of the high fidelity of returning to spawn in their natal streams, which results in highly structured populations, and the wide diversity of environments that salmonids have colonized. These conditions are ideally suited for producing a set of specialist phenotypes, whose fitness is maximized for one specific habitat, rather than a generalist phenotype similarly viable in several environments. Understanding patterns and processes leading to local adaptations has long been a goal of evolutionary biology, but it is only recently that identifying the molecular basis for local adaptation has become feasible because of advances in genomic technologies. The study of shared adaptive phenotypes in populations that are both geographically distant and genetically distinct should reveal some of the fundamental molecular mechanisms associated with local adaptation. In this issue of Molecular Ecology, Miller et al. (2012) make a significant contribution to the development of adaptation genomics. This study suggests that salmonids use standing genetic variation to select beneficial alleles for local adaptations rather than de novo mutations in the same gene or alternative physiological pathways. Identifying the genetic basis for local adaptation has major implications for the management, conservation and potential restoration of salmonid populations.     

Evaluation of a genotyping method based on the ospA gene to detect Borrelia burgdorferi sensu lato in multiple samples of lyme borreliosis patients

In this study we have developed a new Restriction-Fragment-Length-Polymorphism (RFLP) genotyping method for rapid detection and identification of Borrelia genospecies present as unique species or as co-infection in multiple specimens obtained simultaneously from 29 individual patients affected by early or late Lyme borreliosis (LB). The target of the RFLP-genotyping was the heterogeneous plasmid located ospA gene, thus we developed a method able to detect and differentiate between six clinically relevant Borrelia genospecies circulating in Europe, B. burgdorferi sensu stricto, B. garinii, B. afzelii, B. valaisiana, B. bissettii and B. spielmanii. In this study Borrelia DNA could be detected by PCR in at least one specimen of each patient, except in one case of neuroborreliosis (NB); blood samples gave the highest sensitivity in all patient groups. The genotyping indicated that B. afzelii was present in 8 patients with skin involvement, B. garinii in 2 cases of NB and 4 cases with skin involvement, B. burgdorferi sensu stricto was detected in one patient with skin involvement and another with Lyme arthritis. Different Borrelia species in distinct specimens were identified in one patient with EM. The RFLP analysis of 11 patients revealed mixed patterns, which suggested pluri-infection with different Borrelia species.     

A de novo germline mutation of APC for inheritable colon cancer in a Chinese family using multigene next generation sequencing

Inheritable colorectal cancers (CRC) accounted for about 20% of the CRC cases, such as hereditary nonpolyposis colorectal cancer (HNPCC), Gardner syndrome and familial adenomatous polyposis (FAP). A four-generation Han Chinese family was found affected with polyposis in colons. Inferred from the pedigree structure, the disease in this family showed an autosomal dominant inheritance model. To locate the causal mutations in this family, genomic DNAs were extracted and the next generation sequencing for 5 genes relating to colon cancer performed by Ion Torrent Personal Genome Machine with a 314 chip. The reads were aligned with human reference genome hg19 to call variants in the 5 genes. After analysis, 14 variants were detected in the sequenced sample and 13 been collected in dbSNP database and assigned with a rs identification number. In these variants, 9 were synonymous, 4 missense and 1 non-sense. In them, 2 rare variants (c.694C>T in APC and c.1690A>G in MSH2) might be the putative causal mutations for familial adenomatous polyposis (FAP) since the rarity of the mutated allele in normal controls. c.694C>T was detected in only affected members and generated a premature stop codon in APC. It should be a de novo germline mutation making APC containing this stop codon as targets for nonsense-mediated mRNA decay (NMD). c.1690A>G in MSH2 was not only detected in affected members, but also in normal ones in the family. Functional prediction revealed that the amino acid affected by this variant had no effect on the function of MSH2. Here, we report a de novo germline mutation of APC as the causal variant in a Chinese family with inheritable colon cancer by the next generation sequencing.     

Validation of multiple single nucleotide variation calls by additional exome analysis with a semiconductor sequencer to supplement data of whole-genome sequencing of a human population

Background

Validation of single nucleotide variations in whole-genome sequencing is critical for studying disease-related variations in large populations. A combination of different types of next-generation sequencers for analyzing individual genomes may be an efficient means of validating multiple single nucleotide variations calls simultaneously.

Results

Here, we analyzed 12 independent Japanese genomes using two next-generation sequencing platforms: the Illumina HiSeq 2500 platform for whole-genome sequencing (average depth 32.4×), and the Ion Proton semiconductor sequencer for whole exome sequencing (average depth 109×). Single nucleotide polymorphism (SNP) calls based on the Illumina Human Omni 2.5-8 SNP chip data were used as the reference. We compared the variant calls for the 12 samples, and found that the concordance between the two next-generation sequencing platforms varied between 83% and 97%.

Conclusions

Our results show the versatility and usefulness of the combination of exome sequencing with whole-genome sequencing in studies of human population genetics and demonstrate that combining data from multiple sequencing platforms is an efficient approach to validate and supplement SNP calls.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-673) contains supplementary material, which is available to authorized users.     

Exome sequencing followed by large-scale genotyping suggests a limited role for moderately rare risk factors of strong effect in schizophrenia

Schizophrenia is a severe psychiatric disorder with strong heritability and marked heterogeneity in symptoms, course, and treatment response. There is strong interest in identifying genetic risk factors that can help to elucidate the pathophysiology and that might result in the development of improved treatments. Linkage and genome-wide association studies (GWASs) suggest that the genetic basis of schizophrenia is heterogeneous. However, it remains unclear whether the underlying genetic variants are mostly moderately rare and can be identified by the genotyping of variants observed in sequenced cases in large follow-up cohorts or whether they will typically be much rarer and therefore more effectively identified by gene-based methods that seek to combine candidate variants. Here, we consider 166 persons who have schizophrenia or schizoaffective disorder and who have had either their genomes or their exomes sequenced to high coverage. From these data, we selected 5,155 variants that were further evaluated in an independent cohort of 2,617 cases and 1,800 controls. No single variant showed a study-wide significant association in the initial or follow-up cohorts. However, we identified a number of case-specific variants, some of which might be real risk factors for schizophrenia, and these can be readily interrogated in other data sets. Our results indicate that schizophrenia risk is unlikely to be predominantly influenced by variants just outside the range detectable by GWASs. Rather, multiple rarer genetic variants must contribute substantially to the predisposition to schizophrenia, suggesting that both very large sample sizes and gene-based association tests will be required for securely identifying genetic risk factors.     

Hepatitis E virus‐based evaluation of a virion concentration method and detection of enteric viruses in environmental samples by multiplex nested RT‐PCR

Aims: The prevalence of enteric viruses in drinking and river water samples collected from Pune, India was assessed. During an outbreak of HEV in a small town near pune, water samples were screened for enteric viruses. Methods and Results: The water samples were subjected to adsorption–elution‐based virus concentration protocol followed by multiplex nested PCR. Among 64 Mutha river samples, 49 (76·56%) were positive for Hepatitis A Virus, 36 (56·25%) were positive for Rotavirus, 33 (51·56%) were positive for Enterovirus and 16 (25%) were positive for Hepatitis E Virus RNA. Only enterovirus RNA was detected in 2/662 (0·3%) drinking water samples, and the samples from the city’s water reservoir tested negative for all four viruses. HEV RNA was detected in three out of four river water samples during HEV outbreak and partial sequences from patients and water sample were identical. Conclusions: The study suggests absence of enteric viruses both in the source and in the purified water samples from Pune city, not allowing evaluation of the purification system and documents high prevalence of enteric viruses in river water, posing threat to the community. Significance and Impact of the Study: The rapid, sensitive and relatively inexpensive protocol developed for virological evaluation of water seems extremely useful and should be adapted for evaluating viral contamination of water for human consumption. This will lead to development of adequate control measures thereby reducing disease burden because of enteric viruses.     

Haplotype-specific single-locus multiplex PCR assay for molecular identification of sea-bob shrimp,Xiphopenaeus kroyeri (Heller, 1862), cryptic species from the Southwest Atlantic using a DNA pooling strategy for simultaneous identification of multiple samples

This work describes a single-locus multiplex PCR assay based on partial COI mitochondrial gene polymorphisms for identification of two Atlantic cryptic species of the sea-bob shrimp, Xiphopenaeus kroyeri (Heller, 1862) that were recently identified using molecular approaches. Reliable identification of cryptic species of Xiphopenaeus spp. has fundamental implications for management and conservation of the sea-bob shrimp fishery stocks. The assay was developed based on sequence polymorphisms of 130 specimens of both species, comprising samples from Venezuela to the southern coast of Brazil, validated by the amplification of 368 adult shrimp samples from nine different locations and confirmed by direct sequencing. The methodology has been optimized to enable the identification of equi-molar mixtures of DNA from up to 10 individuals by PCR reaction, allowing the fast and cheap identification of many specimens for large scale studies on fisheries biology and population genetics. The DNA pooling strategy enabled the identification of a new locality of occurrence of Xiphopenaeus sp. II in the Brazilian coast, Caravelas, indicating that the species distribution may be continuous on the coast, and not disjoint as observed so far.