Pooled DNA Research Articles

Genome-wide detection of copy number variants in European autochthonous and commercial pig breeds by whole-genome sequencing of DNA pools identified breed-characterising copy number states.

In this study, we identified copy number variants (CNVs) in 19 European autochthonous pig breeds and in two commercial breeds (Italian Large White and Italian Duroc) that represent important genetic resources for this species. The genome of 725 pigs was sequenced using a breed-specific DNA pooling approach (30-35 animals per pool) obtaining an average depth per pool of 42×. This approach maximised CNV discovery as well as the related copy number states characterising, on average, the analysed breeds. By mining more than 17.5billion reads, we identified a total of 9592 CNVs (~683 CNVs per breed) and 3710 CNV regions (CNVRs; 1.15% of the reference pig genome), with an average of 77 CNVRs per breed that were considered as private. A few CNVRs were analysed in more detail, together with other information derived from sequencing data. For example, the CNVR encompassing the KIT gene was associated with coat colour phenotypes in the analysed breeds, confirming the role of the multiple copies in determining breed-specific coat colours. The CNVR covering the MSRB3 gene was associated with ear size in most breeds. The CNVRs affecting the ELOVL6 and ZNF622 genes were private features observed in the Lithuanian Indigenous Wattle and in the Turopolje pig breeds respectively. Overall, the genome variability unravelled here can explain part of the genetic diversity among breeds and might contribute to explain their origin, history and adaptation to a variety of production systems.

A New Design for the Fixed Primer Regions in an Oligonucleotide Library for SELEX Aptamer Screening.

Single-stranded DNA or RNA oligonucleotides, often called aptamers, can bind to a molecular target with both high affinity and selectivity due to their distinct three-dimensional structures. A technique called systematic evolution of ligands by exponential enrichment (SELEX) is used to screen aptamers from a random DNA or RNA pool, or library. The traditionally-designed oligonucleotides in libraries contain a randomized central region along with a fixed primer region at each end for amplifying target-bound central sequences. The single-stranded forward and reverse primer sequences may interfere with target-binding to the central region, resulting in a partial or complete loss of high-affinity aptamers during the SELEX process. To address this issue, researchers have modified the traditional oligonucleotide libraries and developed new types of oligonucleotide libraries; however, these approaches come with various limitations. The author proposes a new design that uses a conformation-changeable sequence as primers, which may open a new avenue for developing an optimized aptamer sequence with both high affinity for, and selective binding to, a particular target via SELEX.

Culicidae-centric metabarcoding through targeted use of D2 ribosomal DNA primers.

A practical limitation to many metabarcoding initiatives is that sampling methods tend to collect many non-target taxa, which become “amplicon noise” that can saturate Next Generation Sequencing results and lead to both financial and resource inefficiencies. An available molecular tool that can significantly decrease these non-target amplicons and decrease the need for pre-DNA-extraction sorting of bycatch is the design of PCR primers tailored to the taxa under investigation. We assessed whether the D2 extension segment of the 28S ribosomal operon can limit this shortcoming within the context of mosquito (Culicidae) monitoring. We designed PCR primers that are fully conserved across mosquitos and exclude from amplification most other taxa likely to be collected with current sampling apparatuses. We show that, given enough sequencing depth, D2 is an effective marker for the detection of mosquito sequences within mock genomic DNA pools. As few as 3,050 quality-filtered Illumina reads were able to recover all 17 species in a bulk pool containing as little as 0.2% of constituent DNA from single taxa. We also mixed these mosquito DNA pools with high concentrations of non-Culicidae bycatch DNA and show that the component mosquito species are generally still recoverable and faithful to their original relative frequencies. Finally, we show that there is little loss of fidelity in abundance parameters when pools from degraded DNA samples were sequenced using the D2 primers.

High-Resolution Mapping of the Novel Early Leaf Senescence Gene Els2 in Common Wheat.

Early leaf senescence negatively impacts the grain yield in wheat (Triticum aestivum L.). Induced mutants provide an important resource for mapping and cloning of genes for early leaf senescence. In our previous study, Els2, a single incomplete dominance gene, that caused early leaf senescence phenotype in the wheat mutant LF2099, had been mapped on the long arm of chromosome 2B. The objective of this study was to develop molecular markers tightly linked to the Els2 gene and construct a high-resolution map surrounding the Els2 gene. Three tightly linked single-nucleotide polymorphism (SNP) markers were obtained from the Illumina Wheat 90K iSelect SNP genotyping array and converted to Kompetitive allele-specific polymerase chain reaction (KASP) markers. To saturate the Els2 region, the Axiom® Wheat 660K SNP array was used to screen bulked extreme phenotype DNA pools, and 9 KASP markers were developed. For fine mapping of the Els2 gene, these KASP markers and previously identified polymorphic markers were analyzed in a large F2 population of the LF2099 × Chinese Spring cross. The Els2 gene was located in a 0.24-cM genetic region flanked by the KASP markers AX-111643885 and AX-111128667, which corresponded to a physical interval of 1.61 Mb in the Chinese Spring chromosome 2BL containing 27 predicted genes with high confidence. The study laid a foundation for a map-based clone of the Els2 gene controlling the mutation phenotype and revealing the molecular regulatory mechanism of wheat leaf senescence.

Development of a Large Gene-Associated SSR Marker Set and in-Depth Genetic Characterization in Scarlet Sage.

Salvia splendens, scarlet or tropical sage, is a tender perennial herbaceous flowering plant popularly grown in public and private gardens all over the world. In this study, we developed a set of simple sequence repeats (SSRs) from genome-wide sequences to assess the genetic diversity and population structure among 112 cultivars. We obtained 364,379 SSRs by mining scarlet sage’s recently published whole genome sequence; 14,545 gene-associated SSR loci were identified in 2 kb gene flanking regions. Among the 768 gene-associated SSR primer sets we screened, 576 loci successfully amplified in DNA pools of 3–4 different cultivars, of which 271 remained polymorphic when tested across eight individual plants. We searched for the related gene functions attributable to these gene-associated SSRs using diverse databases, resulting in 259 Non-redundant matching sequences, 205 individual Gene Ontology (GO) terms, 236 assigned to eukaryotic orthologous groups, and 67 KEGG-annotated (Kyoto Encyclopedia of Genes and Genomes) sequences. We finally selected 41 polymorphic SSR loci to infer genetic diversity and population structure among 112 S. splendens accessions. Based on the developed gene-associated SSRs, clustering analyses consistently revealed two distinct genetic groups within the core collection of S. splendens cultivars. This work developed and characterized an exhaustive set of genome-wide gene-associated SSR markers for scarlet sage. These SSRs can provide species identification, genetic diversity and population structure information for S. splendens, and will therefore be important tools for the management and protection of S. splendens germplasm.

Accurate prediction of somatic variants using deep learning model.

e13659 Background: Efficient and accurate identification of somatic variant is important for understanding the formation, progression, and treatment of cancer. It is necessary to conduct manual review by Integrative Genomic Viewer (IGV) in traditional variant calling process. However, the traditional manual is heavy workload when evaluating tumor with a high variant burden. In this study, a new convolutional neural network (CNN) method was created to train models for somatic mutation identification, which was suitable for Panel sequencing platform with different tumor purities. Methods: A total of 1000 tumor samples from next generation sequencing (NGS)-based genetic testing by a College of American Pathologists (CAP) accredited and Clinical Laboratory Improvement Amendments (CLIA) certified laboratory. Through variant calling program, like GATK, the candidate mutation locations were identified and standardized by manual confirmation. For each candidate mutation location, reads of both tumor and control tissue were extracted. A 2-dimensional feature matrix M of size (2k+1) * 32 in each candidate base was created. The rows of 2k+1 represented the length of candidate region, and the 32 columns included the reads coverage frequency, mapping quality messages, and genome local scores of different tumor and control tissues. CNN model, which includes nine convolutional layers structured by Temporal Convolutional Networks (TCN) but with a different structure to adapt to the proposed input matrix, was used for training. The training data set including manually validated sequence data was used as benchmark test, and optimized by Stochastic Gradient Descent (SGD) optimizer with a learning rate of 0.01 was used for training. Results: The validation data set included 15 mixed samples which were composed of different proportions of known cell lines and real mixed blood samples. The pooled DNA contained 2,359 somatic variants, with expected variant allele frequencies ranged from 3% to 97% in each pool. The overall sensitivity and positive predictive value (PPV) of single nucleotide variants (SNVs) were 99.3% and 99.8%, respectively. Conclusions: A novel and sensitive computational tool for somatic variation detection in DNA Panel sequencing was developed. Our result showed that the deep learning CNN model could call variant in Panel sequencing data.

In silico validation of pooled genotyping strategies for genomic evaluation in Angus cattle.

In this study, we aimed to assess the value of genotyping DNA pools as a strategy to generate accurate and cost-effective genomic estimated breeding values (GEBV) of sires in multi-sire mating systems. In order to do that, we used phenotypic records of 2,436 Australian Angus cattle from 174 sires, including yearling weight (YWT; N = 1,589 records), coat score (COAT; N = 2,026 records), and Meat Standards Australia marbling score (MARB; N = 1,304 records). Phenotypes were adjusted for fixed effects and age at measurement and pools of 2, 5, 10, 15, 20, and 25 animals were explored. Pools were created either by phenotype or at random. When pools were created at random, 10 replicates were examined to provide a measure of sampling variation. The relative accuracy of each pooling strategy was measured by the Pearson correlation coefficient between the sire's GEBV with pooled progeny and the GEBV using individually genotyped progeny. Random pools allow the computation of sire GEBV that are, on average, moderately correlated (i.e., r > 0.5 at pool sizes [PS] ≤ 10) with those obtained without pooling. However, for pools assigned at random, the difference between the best and the worst relative accuracy obtained out of the 10 replicates was as high as 0.41 for YWT, 0.36 for COAT, and 0.61 for MARB. This uncertainty associated with the relative accuracy of GEBV makes randomly assigning animals to pools an unreliable approach. In contrast, pooling by phenotype allowed the estimation of sires' GEBV with a relative accuracy ≥ 0.9 at PS < 10 for all three phenotypes. Moreover, even with larger PS, the lowest relative accuracy obtained was 0.88 (YWT, PS = 20). In agreement with results using simulated data, we conclude that pooling by phenotype is a robust approach to implementing genomic evaluation using commercial herd data, and PS larger than 10 individuals can be considered.

De Novo Genome Assembly of the Meadow Brown Butterfly, Maniola jurtina.

Meadow brown butterflies (Maniola jurtina) on the Isles of Scilly represent an ideal model in which to dissect the links between genotype, phenotype and long-term patterns of selection in the wild - a largely unfulfilled but fundamental aim of modern biology. To meet this aim, a clear description of genotype is required. Here we present the draft genome sequence of M. jurtina to serve as a founding genetic resource for this species. Seven libraries were constructed using pooled DNA from five wild caught spotted females and sequenced using Illumina, PacBio RSII and MinION technology. A novel hybrid assembly approach was employed to generate a final assembly with an N50 of 214 kb (longest scaffold 2.9 Mb). The sequence assembly described here predicts a gene count of 36,294 and includes variants and gene duplicates from five genotypes. Core BUSCO (Benchmarking Universal Single-Copy Orthologs) gene sets of Arthropoda and Insecta recovered 90.5% and 88.7% complete and single-copy genes respectively. Comparisons with 17 other Lepidopteran species placed 86.5% of the assembled genes in orthogroups. Our results provide the first high-quality draft genome and annotation of the butterfly M. jurtina.

Technical note: Droplet digital PCR precisely and accurately quantifies sex skew in bovine semen

Quantifying the relative population of sperm cells bearing the X or Y chromosome in a sexed-semen sample has historically been limited to methods that are either low throughput and sensitive to user subjectivity (e.g., fluorescence in situ hybridization), conterminous (using the same technology to generate and confirm the sex skew), or relatively insensitive (including quantitative PCR with a change detection threshold of 2×). Customers pay a premium for sexed semen and should have access to reliable sex skew data, generated by an accurate, precise test that is orthogonal to the method used to generate sexed semen. Droplet digital PCR (ddPCR) has the capacity to provide an accurate and precise sex skew quantitation by subdividing a pool of template DNA into nanoliter-scale droplets containing either 1 or 0 copies of template DNA. Then PCR amplification is conducted in the droplets, and the number of copies of the amplicon of interest can be counted as the number of fluorescence-positive droplets based on classic quantitative PCR fluorescent reporters. We have optimized and validated a multiplexed ddPCR assay that uses this copy counting method to quantify the sex skew (ratio of X or Y chromosomes) in frozen-thawed bovine sexed semen. The assay interrogates at least 1,000 cells per sample well, quantifying X and Y chromosome copy numbers along with an autosomal gene, GAPDH, used as an internal assay control to confirm total cells counted. The ddPCR sex skew assay achieved a 0.5-percentage-point variance for %X or %Y with a broad linear detection range, from 10 to 95% X, and provided reproducible skew values across a range of 9 to 27 ng of genomic DNA input. This approach overcomes some limitations of other sex skew assays by quantifying absolute X and Y chromosome copy numbers, thus providing a rigorous, independent assessment of sex-skewed semen.

Assessing breed integrity of G\xf6ttingen Minipigs

BackgroundGöttingen Minipigs (GMP) is the smallest commercially available minipig breed under a controlled breeding scheme and is globally bred in five isolated colonies. The genetic isolation harbors the risk of stratification which might compromise the identity of the breed and its usability as an animal model for biomedical and human disease. We conducted whole genome re-sequencing of two DNA-pools per colony to assess genomic differentiation within and between colonies. We added publicly available samples from 13 various pig breeds and discovered overall about 32 M loci, ~ 16 M. thereof variable in GMPs. Individual samples were virtually pooled breed-wise. FST between virtual and DNA pools, a phylogenetic tree, principal component analysis (PCA) and evaluation of functional SNP classes were conducted. An F-test was performed to reveal significantly differentiated allele frequencies between colonies. Variation within a colony was quantified as expected heterozygosity.ResultsPhylogeny and PCA showed that the GMP is easily discriminable from all other breads, but that there is also differentiation between the GMP colonies. Dependent on the contrast between GMP colonies, 4 to 8% of all loci had significantly different allele frequencies. Functional annotation revealed that functionally non-neutral loci are less prone to differentiation. Annotation of highly differentiated loci revealed a couple of deleterious mutations in genes with putative effects in the GMPs .ConclusionDifferentiation and annotation results suggest that the underlying mechanisms are rather drift events than directed selection and limited to neutral genome regions. Animal exchange seems not yet necessary. The Relliehausen colony appears to be the genetically most unique GMP sub-population and could be a valuable resource if animal exchange is required to maintain uniformity of the GMP.

A novel metabarcoded 18S ribosomal DNA sequencing tool for the detection of Plasmodium species in malaria positive patients

Various PCR based methods have been described for the diagnosis of malaria, but most depend on the use of Plasmodium species-specific probes and primers; hence only the tested species are identified and there is limited available data on the true circulating species diversity. Sensitive diagnostic tools and platforms for their use are needed to detect Plasmodium species in both clinical cases and asymptomatic infections that contribute to disease transmission. We have recently developed for the first time a novel high throughput ‘haemoprotobiome’ metabarcoded DNA sequencing method and applied it for the quantification of haemoprotozoan parasites (Theleria and Babesia) of livestock. Here, we describe a novel, high throughput method using an Illumina MiSeq platform to demonstrate the proportions of Plasmodium species in metabarcoded DNA samples derived from human malaria patients. Plasmodium falciparum and Plasmodium vivax positive control gDNA was used to prepare mock DNA pools of parasites to evaluate the detection threshold of the assay for each of the two species. The different mock pools demonstrate the accurate detection ability and to show the proportions of each of the species being present. We then applied the assay to malaria-positive human samples to show the species composition of Plasmodium communities in the Punjab province of Pakistan and in the Afghanistan-Pakistan tribal areas. The diagnostic performance of the deep amplicon sequencing method was compared to an immunochromatographic assay that is widely used in the region. The deep amplicon sequencing showed that P. vivax was present in 69.8%, P. falciparum in 29.5% and mixed infection in 0.7% patients examined. The immunochromatographic assay showed that P. vivax was present in 65.6%, P. falciparum in 27.4%, mixed infection 0.7% patients and 6.32% malaria-positive cases were negative in immunochromatographic assay, but positive in the deep amplicon sequencing. Overall, metabarcoded DNA sequencing demonstrates better diagnostic performance, greatly increasing the estimated prevalence of Plasmodium infection. The next-generation sequencing method using metabarcoded DNA has potential applications in the diagnosis, surveillance, treatment, and control of Plasmodium infections, as well as to study the parasite biology.

Polymorphisms in JAK2 Gene are Associated with Production Traits and Mastitis Resistance in Dairy Cattle

Abstract The present study was designed to investigate the effects of single nucleotide polymorphisms (SNPs) in the JAK2 gene on the production and mastitis related traits in dairy cattle. Blood and milk samples were collected from 201 lactating dairy cattle of three breeds, i.e. Holstein Friesian (HF), Jersey (J) and Achai (A) and their crosses maintained at well-established dairy farms in Khyber Pakhtunkhwa, Pakistan. Generalized linear model was used to evaluate the association between genotypes and the studied traits. A DNA pool was made from randomly selected 30 samples which revealed three SNPs, i.e. SNP 1 in 5’ upstream region (G&gt;A, rs379754157), SNP 2 in intron 15 (A&gt;G, rs134192265), and SNP 3 in exon 20 (A&gt;G, rs110298451) that were further validated in the population under study using SNaPshot technique. Of the three SNPs, SNP 1 did not obey Hardy-Weinberg equilibrium (P&lt;0.05). SNP 2 and SNP 3 were found to be in strong linkage disequilibrium and allele G was highly prevalent compared to allele A in these SNPs. in SNP 1, the GG genotype was associated with significantly (P&lt;0.01) higher SCC, whereas SNP 2 and SNP 3 were significantly (P&lt;0.01) associated with higher lactose percentage compared to the other geno-types. The haplogroups association analysis revealed that H1H2 (GG GG AG) has significantly lower SCC than H2H2 (GG GG GG). The results infer that JAK2 could be an important candidate gene and the studied SNPs might be useful genetic markers for production and mastitis related traits.

Digital DNA microarray generation on glass substrates

In this work we show how DNA microarrays can be produced batch wise on standard microscope slides in a fast, easy, reliable and cost-efficient way. Contrary to classical microarray generation, the microarrays are generated via digital solid phase PCR. We have developed a cavity-chip system made of a PDMS/aluminum composite which allows such a solid phase PCR in a scalable and easy to handle manner. For the proof of concept, a DNA pool composed of two different DNA species was used to show that digital PCR is possible in our chips. In addition, we demonstrate that DNA microarray generation can be realized with different laboratory equipment (slide cycler, manually in water baths and with an automated cartridge system). We generated multiple microarrays and analyzed over 13,000 different monoclonal DNA spots to show that there is no significant difference between the used equipment. To show the scalability of our system we also varied the size and number of the cavities located in the array region up to more than 30,000 cavities with a volume of less than 60 pL per cavity. With this method, we present a revolutionary tool for novel DNA microarrays. Together with new established label-free measurement systems, our technology has the potential to give DNA microarray applications a new boost.

Human papillomavirus prevalence in oral potentially malignant disorders: Systematic review and meta‐analysis

We aimed to provide pooled estimates of human papillomavirus (HPV) prevalence in oral potentially malignant disorders (OPMD) and evaluate the impact of presence of epithelial dysplasia. We searched PubMed, Embase, and Cochrane Library databases for studies that examined the prevalence of HPV DNA in OPMD tested by polymerase chain reaction (PCR). Across 52 eligible studies (2,677 cases), we found an overall pooled HPV prevalence of 22.5% (95% confidence interval [CI] 16.6-29.0). Between-study heterogeneity was 93%. When stratified by subgroup, the pooled HPV prevalence in leukoplakia (1,232 cases) was 20.2% (95% CI 11.2-31.1), lichen planus (767 cases) 23.0% (95% CI 15.0-32.2), oral submucous fibrosis (238 cases) 28.6% (95% CI 23.0-34.5), proliferative verrucous leukoplakia (60 cases) 24.7% (95% CI 1.8-62.0), and OPMD unspecified (377 cases) 25.4% (95% CI 16.2-35.8). Information on presence of epithelial dysplasia was available in 19 studies, and the results did not vary substantially between non-dysplastic and dysplastic samples. HPV16 was the predominant genotype among HPV-positive OPMD cases (48.2%, 95% CI 31.4-65.2). We found a pooled HPV DNA prevalence of 22.5% in OPMD cases with great between-study heterogeneity. The HPV prevalence appeared to be comparable across subgroups and independent of epithelial dysplasia.

Role of Hepatitis B virus capsid phosphorylation in nucleocapsid disassembly and covalently closed circular DNA formation.

Hepatitis B virus (HBV) delivers a partially double-stranded, relaxed circular (RC) DNA genome in complete virions to the host cell nucleus for conversion to the covalently closed circular (CCC) DNA, which establishes and sustains viral infection. An overlength pregenomic RNA (pgRNA) is then transcribed from CCC DNA and packaged into immature nucleocapsids (NCs) by the viral core (HBc) protein. pgRNA is reverse transcribed to produce RC DNA in mature NCs, which are then enveloped and secreted as complete virions, or delivered to the nucleus to replenish the nuclear CCC DNA pool. RC DNA, whether originating from extracellular virions or intracellular mature NCs, must be released upon NC disassembly (uncoating) for CCC DNA formation. HBc is known to undergo dynamic phosphorylation and dephosphorylation at its C-terminal domain (CTD) to facilitate pgRNA packaging and reverse transcription. Here, two putative phosphorylation sites in the HBc N-terminal domain (NTD), S44 and S49, were targeted for genetic and biochemical analysis to assess their potential roles in viral replication. The NTD mutant that mimics the non-phosphorylated state (N2A) was competent in all steps of viral replication tested from capsid assembly, pgRNA packaging, reverse transcription, to virion secretion, except for a decrease in CCC DNA formation. On the other hand, the phosphor-mimetic mutant N2E showed a defect in the early step of pgRNA packaging but enhanced the late step of mature NC uncoating and consequently, increased CCC DNA formation. N2E also enhanced phosphorylation in CTD and possibly elsewhere in HBc. Furthermore, inhibition of the cyclin-dependent kinase 2 (CDK2), which is packaged into viral capsids, could block CCC DNA formation. These results prompted us to propose a model whereby rephosphorylation of HBc at both NTD and CTD by the packaged CDK2, following CTD dephosphorylation during NC maturation, facilitates uncoating and CCC DNA formation by destabilizing mature NCs.

Comparative study of population genomic approaches for mapping colony-level traits.

Social insect colonies exhibit colony-level phenotypes such as social immunity and task coordination, which are produced by the individual phenotypes. Mapping the genetic basis of such phenotypes requires associating the colony-level phenotype with the genotypes in the colony. In this paper, we examine alternative approaches to DNA extraction, library construction, and sequencing for genome wide association studies (GWAS) of colony-level traits using a population sample of Cataglyphis niger ants. We evaluate the accuracy of allele frequency estimation from sequencing a pool of individuals (pool-seq) from each colony using either whole-genome sequencing or reduced representation genomic sequencing. Based on empirical measurement of the experimental noise in sequenced DNA pools, we show that reduced representation pool-seq is drastically less accurate than whole-genome pool-seq. Surprisingly, normalized pooling of samples did not result in greater accuracy than un-normalized pooling. Subsequently, we evaluate the power of the alternative approaches for detecting quantitative trait loci (QTL) of colony-level traits by using simulations that account for an environmental effect on the phenotype. Our results can inform experimental designs and enable optimizing the power of GWAS depending on budget, availability of samples and research goals. We conclude that for a given budget, sequencing un-normalized pools of individuals from each colony provides optimal QTL detection power.

An enzymatic on/off switch-mediated assay for KRAS hotspot point mutation detection of circulating tumor DNA.

BackgroundTo detect the mutations of KRAS gene in colorectal cancer patients and other cancer patients, it is of value to develop non‐invasive, sensitive, specific, easy, and low‐cost assays.MethodsTemplates harboring hotspot mutations of the KRAS gene were constructed, and primers were designed for evaluation of the specificity, and sensitivity of detection system consisted of exonuclease polymerase‐mediated on/off switch; then, gel electrophoresis and real‐time PCR were performed for verification. The assay was verified by testing the DNA pool of normal controls and circulating DNA (ctDNA) samples from 14 tumor patients, as compared to Sanger sequencing.ResultsA specific and sensitive assay consisted of exonuclease polymerase‐mediated on/off switch, and multiplex real‐time PCR method has been established. This assay could detect <100 copies of KRAS mutation in more than 10 million copies of wild‐type KRAS gene fragments. This assay was applied to test KRAS gene mutations in three cases of fourteen ctDNA samples, and the results were consistent with Sanger sequencing. However, this PCR‐based assay was more sensitive and easier to be interpreted.ConclusionThis assay can detect the presence of KRAS hotspot mutations in clinical circulating tumor DNA samples. The assay has a potential to be used in early diagnosis of colorectal cancer as well as other types of cancer.

Genotyping-by-sequencing of pooled drone DNA for the management of living honeybee (Apis mellifera) queens in commercial beekeeping operations in New Zealand

The absence of a full pedigree can hinder selective breeding efforts. In honeybees, definitive maternity and especially paternity of queens is difficult to determine, even under managed mating schemes (e.g. using artificial insemination) due to the negative effects of single-drone mating on colony fitness. Here we genotyped 388 living queens from two beekeeping operations using Genotyping-by-Sequencing (GBS). We evaluate two methods to call single-nucleotide polymorphism (SNPs), Tassel 5 and Stacks, for their ability to supply SNPs that can recover known relationships. While Stacks discovered more SNPs (29,433), SNPs called with Tassel 5 (16,757) were found to be more accurate for the derivation of relationships. This methodology presents a low-cost genotyping approach and can be used to support commercial honeybee breeding schemes.

Identification of coronavirus sequences in carp cDNA from Wuhan, China.

Severe acute respiratory syndrome (SARS)‐like coronavirus sequences were identified in two separate complementary DNA (cDNA) pools. The first pool was from a Carassius auratus (crusian carp) cell line and the second was from Ctenopharyngodon idella (grass carp) head kidney tissue. BLAST analysis suggests that these sequences belong to SARS‐like coronaviruses, and that they are not evolutionarily conserved in other species. Investigation of the submitting laboratories revealed that two laboratories from the Institute of Hydrobiology at the Chinese Academy of Sciences in Wuhan, China performed the research and submitted the cDNA libraries to GenBank. This institution is very close in proximity to the Wuhan South China Seafood Wholesale Market where SARS‐CoV‐2 first amplified in the human population. It is possible that these sequences are an artifact of the bioinformatics pipeline that was used. It is also possible that SARS‐like coronaviruses are a common environmental pathogen in the region that may be in aquatic habitats.

Lack of evidence for CDK12 as an ovarian cancer predisposing gene.

CDK12 variants were investigated as a genetic susceptibility to ovarian cancer in a series of 416 unrelated and consecutive patients with ovarian carcinoma and who carry neither germline BRCA1 nor BRCA2 pathogenic variant. The presence of CDK12 variants was searched in germline DNA by massive parallel sequencing on pooled DNAs. The lack of detection of deleterious variants and the observed proportion of missense variants in the series of ovarian carcinoma patients as compared with all human populations strongly suggests that CDK12 is not an ovarian cancer predisposing gene.