Search

Protein display systems are powerful techniques used to identify protein molecules that bind with high affinity to target proteins of interest. The initial challenge in implementing a display system is the construction of a high-diversity naïve library. Here, we describe the methods to generate a designed ankyrin repeat protein (DARPin) display library using degenerate oligonucleotides. Specifically described is the construction of a single DARPin repeat module by overlap extension PCR, concatenation of the module by restriction enzyme digestion and ligation, and incorporation of the concatenated modules into a full-length DARPin sequence in a bacterial cloning or display vector containing the hydrophilic N- and C-terminal capping domains. Protocols for PCR amplification of DARPin sequences to estimate diversity of naïve and enriched libraries via next-generation sequencing are included, as is a simple Linux-based program for analysis of naïve and enriched sequences. © 2024 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Generation of a single DARPin repeat by overlap extension PCR Basic Protocol 2: Concatenation of DARPin repeats Basic Protocol 3: Ligation of internal repeats into cloning/display vector containing N- and C-terminal capping repeats Basic Protocol 4: Estimation of library size and diversity by next-generation sequencing (NGS) Basic Protocol 5: NGS analysis of naïve and enriched libraries.

Single‐cell whole‐genome sequencing (scWGS) detects cell heterogeneity at the aspect of genomic variations, which are inheritable and play an important role in life processes such as aging and cancer progression. The recent explosive development of high‐throughput single‐cell sequencing methods has enabled high‐performance heterogeneity detection through a vast number of novel strategies. Despite the limitation on total cost, technical advances in high‐throughput single‐cell whole‐genome sequencing methods are made for higher genome coverage, parallel throughput, and level of integration. This review highlights the technical advancements in high‐throughput scWGS in the aspects of strategies design, data efficiency, parallel handling platforms, and their applications on human genome. The experimental innovations, remaining challenges, and perspectives are summarized and discussed. [ABSTRACT FROM AUTHOR]

The hoverfly, Eupeodes corollae, is a worldwide natural enemy of aphids and a plant pollinator. To provide insights into the biology of this species, we examined its population genetic structure by obtaining 1.15-GB random genomic sequences using next-generation sequencing and developing genome-wide microsatellite markers. A total of 79,138 microsatellite loci were initially isolated from the genomic sequences; after strict selection and further testing of 40 primer pairs in eight individuals, 24 polymorphic microsatellites with high amplification rates were developed. These microsatellites were used to examine the population genetic structure of 96 individuals from four field populations collected across southern to northern China. The number of alleles per locus ranged from 5 to 13 with an average of 8.75; the observed and expected heterozygosity varied from 0.235 to 0.768 and from 0.333 to 0.785, respectively. Population genetic structure analysis showed weak genetic differentiation among the four geographical populations of E. corollae, suggesting a high rate of gene flow reflecting likely widespread migration of E. corollae in China. [ABSTRACT FROM AUTHOR]

Terpenes are responsible for most or all of the odor and flavor properties of Cannabis sativa, and may also impact effects users experience either directly or indirectly. We report the diversity of terpene profiles across samples bound for the Washington dispensary market. The remarkable degree of variation in terpene profiles ultimately results from action of a family of terpene synthase genes, only some of which have been described. Using a recently available genome assembly we describe 55 terpene synthases with genomic context, and tissue specific expression. The family is quite diverse from a protein similarity perspective, and subsets of the family are expressed in all tissues in the plant, including a set of root specific monoterpene synthases that could well have agronomic importance. Ultimately understanding and breeding for specific terpene profiles will require a good understanding of the gene family that underlies it. We intend for this work to serve as a foundation for that. [ABSTRACT FROM AUTHOR]

Advances in genomics have made whole genome studies increasingly feasible across the life sciences. However, new technologies and algorithmic advances do not guarantee flawless genomic sequences or annotation. Bias, errors, and artifacts can enter at any stage of the process from library preparation to annotation. When planning an experiment that utilizes a genome sequence as the basis for the design, there are a few basic checks that, if performed, may better inform the experimental design and ideally help avoid a failed experiment or inconclusive result. [ABSTRACT FROM AUTHOR]

Chinese cabbage (Brassica rapa ssp. pekinensis) is a major crop that is widely cultivated, especially in Korea, Japan, and China. With the advent of next generation sequencing technology, the cost and time required for sequencing have decreased and the development of genome research accelerated. Genome sequencing of Chinese cabbage was completed in 2011 using the variety Chiifu-401-42, and since then the genome has been continuously updated. In the present study, we conducted whole-genome sequencing of Chinese cabbage inbred line CT001, a line widely used in traditional or molecular breeding, to improve the accuracy of genetic polymorphism analysis. The constructed CT001 pseudomolecule represented 85.4% (219.8 Mb) of the Chiifu reference genome, and a total of 38,567 gene models were annotated using RNA-Seq analysis. In addition, the spontaneous mutation rate of CT001 was estimated by resequencing DNA obtained from individual plants after sexual propagation for six generations to estimate the naturally occurring variations. The CT001 pseudomolecule constructed in this study will provide valuable resources for genomic studies on Chinese cabbage. [ABSTRACT FROM AUTHOR]

Plant mitochondrial genomes are usually assembled and displayed as circular maps based on the widely-held view across the broad community of life scientists that circular genome-sized molecules are the primary form of plant mitochondrial DNA, despite the understanding by plant mitochondrial researchers that this is an inaccurate and outdated concept. Many plant mitochondrial genomes have one or more pairs of large repeats that can act as sites for inter- or intramolecular recombination, leading to multiple alternative arrangements (isoforms). Most mitochondrial genomes have been assembled using methods unable to capture the complete spectrum of isoforms within a species, leading to an incomplete inference of their structure and recombinational activity. To document and investigate underlying reasons for structural diversity in plant mitochondrial DNA, we used long-read (PacBio) and short-read (Illumina) sequencing data to assemble and compare mitochondrial genomes of domesticated (Lactuca sativa) and wild (L. saligna and L. serriola) lettuce species. We characterized a comprehensive, complex set of isoforms within each species and compared genome structures between species. Physical analysis of L. sativa mtDNA molecules by fluorescence microscopy revealed a variety of linear, branched, and circular structures. The mitochondrial genomes for L. sativa and L. serriola were identical in sequence and arrangement and differed substantially from L. saligna, indicating that the mitochondrial genome structure did not change during domestication. From the isoforms in our data, we infer that recombination occurs at repeats of all sizes at variable frequencies. The differences in genome structure between L. saligna and the two other Lactuca species can be largely explained by rare recombination events that rearranged the structure. Our data demonstrate that representations of plant mitochondrial genomes as simple, circular molecules are not accurate descriptions of their true nature and that in reality plant mitochondrial DNA is a complex, dynamic mixture of forms. [ABSTRACT FROM AUTHOR]

Aeromonas veronii is a gram-negative species abundant in aquatic environments that causes disease in humans as well as terrestrial and aquatic animals. In the current study, 41 publicly available A. veronii genomes were compared to investigate distribution of putative virulence genes, global dissemination of pathotypes, and potential mechanisms of virulence. The complete genome of A. veronii strain ML09-123 from an outbreak of motile aeromonas septicemia in farm-raised catfish in the southeastern United States was included. Dissemination of A. veronii strain types was discovered in dispersed geographical locations. Isolate ML09-123 is highly similar to Chinese isolate TH0426, suggesting the two strains have a common origin and may represent a pathotype impacting aquaculture in both countries. Virulence of strain ML09-123 in catfish in a dose-dependent manner was confirmed experimentally. Subsystem category disposition showed the majority of genomes exhibit similar distribution of genomic elements. The type I secretion system (T1SS), type II secretion system (T2SS), type 4 pilus (T4P), and flagellum core elements are conserved in all A. veronii genomes, whereas the type III secretion system (T3SS), type V secretion system (T5SS), type VI secretion system (T6SS), and tight adherence (TAD) system demonstrate variable dispersal. Distribution of mobile elements is dependent on host and geographic origin, suggesting this species has undergone considerable genetic exchange. The data presented here lends insight into the genomic variation of A. veronii and identifies a pathotype impacting aquaculture globally. [ABSTRACT FROM AUTHOR]

Cassava improvement using traditional breeding strategies is slow due to the species’ long breeding cycle. However, the use of genomic selection can lead to a shorter breeding cycle. This study aimed to estimate genetic parameters for productive traits based on pedigree (pedigree and phenotypic information) and genomic (markers and phenotypic information) analyses using biparental crosses at different stages of selection. A total of 290 clones were genotyped and phenotyped for fresh root yield (FRY), dry matter content (DMC), dry yield (DY), fresh shoot yield (FSY) and harvest index (HI). The clones were evaluated in clonal evaluation trials (CET), preliminary yield trials (PYT), advanced yield trials (AYT) and uniform yield trials (UYT), from 2013 to 2018 in ten locations. The breeding stages were analyzed as follows: one stage (CET), two stages (CET and PYT), three stages (CET, PYT and AYT) and four stages (CET, PYT, AYT and UYT). The genomic predictions were analyzed via k-fold cross-validation based on the genomic best linear unbiased prediction (GBLUP) considering a model with genetic additive effects and genotype × location interactions. Genomic and pedigree accuracies were moderate to high (0.56–0.72 and 0.62–0.78, respectively) for important starch-related traits such as DY and FRY; when considering one breeding stage (CET) with the aim of early selection, the genomic accuracies ranged from 0.60 (DMC) to 0.71 (HI). Moreover, the correlations between the genomic estimation breeding values of one-stage genomic analysis and the estimated breeding values of the four-stage (full data set) pedigree analysis were high for all traits as well as for a selection index including all traits. The results indicate great possibilities for genomic selection in cassava, especially for selection early in the breeding cycle (saving time and effort). [ABSTRACT FROM AUTHOR]

Motivation: Recently, copy number variation (CNV) has gained considerable interest as a type of genomic variation that plays an important role in complex phenotypes and disease susceptibility. Since a number of CNV detection methods have recently been developed, it is necessary to help investigators choose suitable methods for CNV detection depending on their objectives. For this reason, this study compared ten commonly used CNV detection applications, including CNVnator, ReadDepth, RDXplorer, LUMPY and Control-FREEC, benchmarking the applications by sensitivity, specificity and computational demands. Taking the DGV gold standard variants as a standard dataset, we evaluated the ten applications with real sequencing data at sequencing depths from 5X to 50X. Among the ten methods benchmarked, LUMPY performs the best for both high sensitivity and specificity at each sequencing depth. For the purpose of high specificity, Canvas is also a good choice. If high sensitivity is preferred, CNVnator and RDXplorer are better choices. Additionally, CNVnator and GROM-RD perform well for low-depth sequencing data. Our results provide a comprehensive performance evaluation for these selected CNV detection methods and facilitate future development and improvement in CNV prediction methods. [ABSTRACT FROM AUTHOR]

The frequency of epidemics caused by Dengue viruses 1–4, Zika virus and Chikungunya viruses have been on an upward trend in recent years driven primarily by uncontrolled urbanization, mobility of human populations and geographical spread of their shared vectors, Aedes aegypti and Aedes albopictus. Infections by these viruses present with similar clinical manifestations making them challenging to diagnose; this is especially difficult in regions of the world hyperendemic for these viruses. In this study, we present a targeted-enrichment methodology to simultaneously sequence the complete viral genomes for each of these viruses directly from clinical samples. Additionally, we have also developed a customized computational tool (BaitMaker) to design these enrichment baits. This methodology is robust in its ability to capture diverse sequences and is amenable to large-scale epidemiological studies. We have applied this methodology to two large cohorts: a febrile study based in Colombo, Sri Lanka taken during the 2009–2015 dengue epidemic (n = 170) and another taken during the 2016 outbreak of Zika virus in Singapore (n = 162). Results from these studies indicate that we were able to cover an average of 97.04% ± 0.67% of the full viral genome from samples in these cohorts. We also show detection of one DENV3/ZIKV co-infected patient where we recovered full genomes for both viruses. [ABSTRACT FROM AUTHOR]

Mounting evidence supports that LINE-1 (L1) retrotransposition can occur postzygotically in healthy and diseased human tissues, contributing to genomic mosaicism in the brain and other somatic tissues of an individual. However, the genomic distribution of somatic human-specific LINE-1 (L1Hs) insertions and their potential impact on carrier cells remain unclear. Here, using a PCR-based targeted bulk sequencing approach, we profiled 9,181 somatic insertions from 20 postmortem tissues from five Rett patients and their matched healthy controls. We identified and validated somatic L1Hs insertions in both cortical neurons and non-brain tissues. In Rett patients, somatic insertions were significantly depleted in exons—mainly contributed by long genes—than healthy controls, implying that cells carrying MECP2 mutations might be defenseless against a second exonic L1Hs insertion. We observed a significant increase of somatic L1Hs insertions in the brain compared with non-brain tissues from the same individual. Compared to germline insertions, somatic insertions were less sense-depleted to transcripts, indicating that they underwent weaker selective pressure on the orientation of insertion. Our observations demonstrate that somatic L1Hs insertions contribute to genomic diversity and MeCP2 dysfunction alters their genomic patterns in Rett patients. [ABSTRACT FROM AUTHOR]

Mycobacterium cosmeticum is a nontuberculous Mycobacterium recovered from different water sources including household potable water and water collected at nail salon. Individual cases of this bacterium have been reported to be associated with gastrointestinal tract infections. Here we present the first whole-genome study and comparative analysis of two new clinically-derived Mycobacterium sp. UM_RHS (referred as UM_RHS after this) and Mycobacterium sp. UM_NYF (referred as UM_NYF after this) isolated from patients in Indonesia and Malaysia respectively to have a better understanding of the biological characteristic of these isolates. Both strains are likely Mycobacterium cosmeticum as supported by the evidence from molecular phylogenetic, comparative genomic and Average Nucleotide Identity (ANI) analyses. We found the presence of a considerably large number of putative virulence genes in the genomes of UM_RHS and UM_NYF. Interestingly, we also found a horizontally transferred genomic island carrying a putative dsz operon proposing that they may have potential to perform biodesulfization of dibenzothiophene (DBT) that may be effective in cost reduction and air pollution during fuel combustion. This comparative study may provide new insights into M. cosmeticum and serve as an important reference for future functional studies of this bacterial species. [ABSTRACT FROM AUTHOR]

Illnesses caused by Shiga toxin-producing Escherichia coli (STECs) can be life threatening, such as hemolytic uremic syndrome (HUS). The STECs most frequently identified by USDA’s Microbiological Data Program (MDP) carried toxin gene subtypes stx1a and/or stx2a. Here we described the genome sequences of 331 STECs isolated from foods regulated by the FDA 2010–2017, and determined their genomic identity, serotype, sequence type, virulence potential, and prevalence of antimicrobial resistance. Isolates were selected from the MDP archive, routine food testing by FDA field labs (ORA), and food testing by a contract company. Only 276 (83%) strains were confirmed as STECs by in silico analysis. Foods from which STECs were recovered included cilantro (6%), spinach (25%), lettuce (11%), and flour (9%). Phylogenetic analysis using core genome MLST revealed these STEC genomes were highly variable, with some clustering associated with ST types and serotypes. We detected 95 different sequence types (ST); several ST were previously associated with HUS: ST21 and ST29 (O26:H11), ST11 (O157:H7), ST33 (O91:H14), ST17 (O103:H2), and ST16 (O111:H-). in silico virulome analyses showed ~ 51% of these strains were potentially pathogenic [besides stx gene they also carried eae (25%) or 26% saa (26%)]. Virulence gene prevalence was also determined: stx1 only (19%); stx2 only (66%); and stx1/sxt2 (15%). Our data form a new WGS dataset that can be used to support food safety investigations and monitor the recurrence/emergence of E. coli in foods. [ABSTRACT FROM AUTHOR]

Melaphidina aphids (Rhus-gall aphids; Eriosomatinae: Fordini) comprise five genera from eastern Asia and one monotypic genus from eastern North America. Melaphidina are unique in feeding on plant species of Rhus subgenus Rhus (Anacardiaceae), on which they form galls during the summer. The phylogenetic relationships among some species of Melaphidina aphids remain controversial. In this study, we sought to resolve the backbone phylogeny of Melaphidina aphids by sampling 15 accessions representing all six genera, all species, and all subspecies except Meitanaphis microgallis using 20 gene regions: five nuclear genes as well as 13 protein-coding genes and two rRNA genes of the mitochondrial genome. Phylogenetic analyses included Bayesian and maximum likelihood methods. Independent analyses of nuclear and mitochondrial genes returned congruent topologies, and analyses of all gene regions combined showed well-supported relationships among Melaphidina species. In particular, these were: (1) Nurudea (excluding N. ibofushi) is sister to a clade composed of the five remaining genera; (2) the monotypic North American genus Melaphis is sister to a clade comprising the four remaining genera; and (3) (Schlechtendalia + N. ibofushi) is sister to the clade (Floraphis (Meitanaphis + Kaburagia). Our results support the transfer of Meitanaphis flavogallis to Kaburagia as an additional subspecies or species, and the recognition of Floraphis as a distinct genus. This study provides important molecular resources for subsequent evolutionary studies using more nuclear genes on the Melaphidina aphids and their close allies. [ABSTRACT FROM AUTHOR]

Background: Leptospirosis is the most common zoonotic disease worldwide. The diagnostic performance of a serological test for human leptospirosis is mainly influenced by the antigen used in the test assay. An ideal serological test should cover all serovars of pathogenic leptospires with high sensitivity and specificity and use reagents that are relatively inexpensive to produce and can be used in tropical climates. Peptide-based tests fulfil at least the latter two requirements, and ORFeome phage display has been successfully used to identify immunogenic peptides from other pathogens. Methodology/Principal findings: Two ORFeome phage display libraries of the entire Leptospira spp. genomes from five local strains isolated in Malaysia or seven WHO reference strains were constructed. Subsequently, 18 unique Leptospira peptides were identified in a screen using a pool of sera from patients with acute leptospirosis. Five of these were validated by titration ELISA using different pools of patient or control sera. The diagnostic performance of these five peptides was then assessed against 16 individual sera from patients with acute leptospirosis and 16 healthy donors and was compared to that of two recombinant reference proteins from L. interrogans. This analysis revealed two peptides (SIR16-D1 and SIR16-H1) from the local isolates with good accuracy for the detection of acute leptospirosis (area under the ROC curve: 0.86 and 0.78, respectively; sensitivity: 0.88 and 0.94; specificity: 0.81 and 0.69), which was close to that of the reference proteins LipL32 and Loa22 (area under the ROC curve: 0.91 and 0.80; sensitivity: 0.94 and 0.81; specificity: 0.75 and 0.75). Conclusions/Significance: This analysis lends further support for using ORFeome phage display to identify pathogen-associated immunogenic peptides, and it suggests that this technique holds promise for the development of peptide-based diagnostics for leptospirosis and, possibly, of vaccines against this pathogen. [ABSTRACT FROM AUTHOR]

The biocide triclosan is in many consumer products and is a frequent contaminant of wastewater (WW) such that there is concern that triclosan promotes resistance to important antibiotics. This study identified functional mechanisms of triclosan resistance (TCSR) in WW metagenomes, and assessed the frequency of TCSR in WW-derived and clinical isolates of Escherichia coli and Enterococcus spp. Metagenomic DNA extracted from WW was used to profile the microbiome and construct large-insert cosmid libraries, which were screened for TCSR. Resistant cosmids were sequenced and the TCSR determinant identified by transposon mutagenesis. Wastewater Enterococcus spp. (N = 94) and E. coli (N = 99) and clinical Enterococcus spp. (N = 146) and vancomycin-resistant E. faecium (VRE; N = 149) were collected and tested for resistance to triclosan and a comprehensive drug panel. Functional metagenomic screening revealed diverse FabV homologs as major WW TCSR determinants. Resistant clones harboured sequences likely originating from Aeromonas spp., a common WW microbe. The triclosan MIC90s for E. coli, E. faecalis, and E. faecium isolates were 0.125, 32, and 32 mg/L, respectively. For E. coli, there was no correlation between the triclosan MIC and any drug tested. Negative correlations were detected between the triclosan MIC and levofloxacin resistance for E. faecalis, and between triclosan and vancomycin, teicoplanin, and ampicillin resistance for E. faecium. Thus, FabV homologs were the major contributor to the WW triclosan resistome and high-level TCSR was not observed in WW or clinical isolates. Elevated triclosan MICs were not positively correlated with antimicrobial resistance to any drug tested. [ABSTRACT FROM AUTHOR]

Genetic analyses provide a powerful tool with which to identify the biological components of historical objects. Te Tiriti o Waitangi | The Treaty of Waitangi is New Zealand’s founding document, intended to be a partnership between the indigenous Māori and the British Crown. Here we focus on an archived piece of blank parchment that has been proposed to be the missing portion of the lower parchment of the Waitangi Sheet of the Treaty. However, its physical dimensions and characteristics are not consistent with this hypothesis. We perform genetic analyses on the parchment membranes of the Treaty, plus the blank piece of parchment. We find that all three parchments were made from ewes and that the blank parchment is highly likely to be a portion cut from the lower membrane of the Waitangi Sheet because they share identical whole mitochondrial genomes, including an unusual heteroplasmic site. We suggest that the differences in size and characteristics between the two pieces of parchment may have resulted from the Treaty’s exposure to water in the early 20th century and the subsequent repair work, light exposure during exhibition or the later conservation treatments in the 1970s and 80s. The blank piece of parchment will be valuable for comparison tests to study the effects of earlier treatments and to monitor the effects of long-term display on the Treaty. [ABSTRACT FROM AUTHOR]

Genotyping-by-sequencing (GBS) provides high SNP coverage and has recently emerged as a popular technology for genetic and breeding applications in bread wheat (Triticum aestivum L.) and many other plant species. Although GBS can discover millions of SNPs, a high rate of missing data is a major concern for many applications. Accurate imputation of those missing data can significantly improve the utility of GBS data. This study compared imputation accuracies among four genome references including three wheat references (Chinese Spring survey sequence, W7984, and IWGSC RefSeq v1.0) and one barley reference genome by comparing imputed data derived from low-depth sequencing to actual data from high-depth sequencing. After imputation, the average number of imputed data points was the highest in the B genome (~48.99%). The D genome had the lowest imputed data points (~15.02%) but the highest imputation accuracy. Among the four reference genomes, IWGSC RefSeq v1.0 reference provided the most imputed data points, but the lowest imputation accuracy for the SNPs with < 10% minor allele frequency (MAF). The W7984 reference, however, provided the highest imputation accuracy for the SNPs with < 10% MAF. [ABSTRACT FROM AUTHOR]

Moraxella bovoculi is frequently cultured from the ocular secretions and conjunctiva of cattle with Infectious Bovine Keratoconjunctivitis (IBK). Previous work has shown that single nucleotide polymorphism (SNP) diversity in this species is quite high with 81,284 SNPs identified in eight genomes representing two distinct genotypes isolated from IBK affected eyes (genotype 1) and the nasopharynx of cattle without clinical IBK signs (genotype 2), respectively. The goals of this study were to identify SNPs from a collection of geographically diverse and epidemiologically unlinked M. bovoculi strains from the eyes of IBK positive cattle (n = 183) and another from the eyes of cattle (most from a single population at a single time-point) without signs of IBK (n = 63) and to characterize the genetic diversity. Strains of both genotypes were identified from the eyes of cattle without IBK signs. Only genotype 1 strains were identified from IBK affected eyes, however, these strains were isolated before the discovery of genotype 2, and the protocol for their isolation would have preferentially selected genotype 1 M. bovoculi. The core genome comprised ~74% of the whole and contained >127,000 filtered SNPs. More than 80% of these characterize diversity within genotype 1 while 23,611 SNPs (~18%) delimit the two major genotypes. Genotype 2 strains lacked a repeats-in-toxin (RTX) putative pathogenesis factor and any of ten putative antibiotic resistance genes carried within a genomic island. Within genotype 1, prevalence of these elements was 0.85 and 0.12 respectively in strains from eyes that were IBK positive. Recombination appears to be an important source of genetic diversity for genotype 1 and undermines the utility of ribosomal-locus-based species identification. The extremely high genetic diversity in genotype 1 presents a challenge to the development of an efficacious vaccine directed against them, however, several low-diversity pilin-like genes were identified. Finally, the genotype-defining SNPs described in this study are a resource that can facilitate the development of more accurate M. bovoculi diagnostic tests. [ABSTRACT FROM AUTHOR]

Dragon’s blood collected from the genus Dracaena is used as a renowned traditional medicine in various cultures worldwide. However, the genetics of the genus Dracaena and the formation mechanism of dragon’s blood remain poorly understood. Here, we generate the first draft genome reference assembly of an elite Chinese Dracaena species, Dracaena cambodiana, from next-generation sequencing data with 89.46× coverage. The reads were assembled into 2,640,704 contigs with an N50 length of 1.87 kb, and a 1.05 Gb assembly was finally assembled with 2,379,659 scaffolds. Furthermore, 97.75% of the 267,243 simple sequence repeats identified from these scaffolds were mononucleotide, dinucleotide, and trinucleotide repeats. Among all 53,700 predicted genes, 158 genes involved in cell wall and plant hormone synthesis and reactive oxygen species scavenging showed altered regulation during the formation of dragon’s blood. This study provides a genomic characterization of D. cambodiana and improves understanding of the molecular mechanism of dragon’s blood formation. This report represents the first genome-wide characterization of a Dracaena species in the Asparagaceae. [ABSTRACT FROM AUTHOR]

West Nile Virus (WNV) has been detected annually in Maricopa County, Arizona, since 2003. With this in mind, we sought to determine if contemporary strains are endemic to the county or are annually imported. As part of this effort, we developed a new protocol for tiled amplicon sequencing of WNV to efficiently attain greater than 99% coverage of 14 WNV genomes collected directly from positive mosquito pools distributed throughout Maricopa County between 2014 and 2017. Bayesian phylogenetic analyses revealed that contemporary genomes fall within two major lineages; NA/WN02 and SW/WN03. We found that all of the Arizona strains possessed an amino acid substitution known to be under positive selection, which has arisen independently at least four times in Arizona. The SW/WN03 strains exhibited transient behavior, with at least 10 separate introductions into Arizona when considering both historical and contemporary strains. However, NA/WN02 strains are geographically differentiated and appear to be endemic in Arizona, with two clades that have been circulating for four and seven years. This establishment in Maricopa County provides the first evidence of local overwintering by a WNV strain over the course of several years in Arizona. Within a national context, the placement of eleven contemporary Arizona strains in the NA/WN02 lineage indicates while WNV first entered the northeastern United States in 1999, the most ancestral extant strains of WNV are now circulating in the American southwest. [ABSTRACT FROM AUTHOR]

The availability of biological information in public databases has increased exponentially. To ensure the accuracy of this information, researchers have adopted several methods and refinements to avoid the dissemination of incorrect information; for example, several automated tools are available for annotation processes. However, manual curation ensures and enriches biological information. Additionally, the genomic finishing process is complex, resulting in increased deposition of drafts genomes. This introduces bias in other omics analyses because incomplete genomic content is used. This is also observed for complete genomes. For example, genomes generated by reference assembly may not include new products in the new sequence or errors or bias can occur during the assembly process. Thus, we developed ImproveAssembly, a tool capable of identifying new products missing from genomic sequences, which can be used for complete and draft genomes. The identified products can improve the annotation of complete genomes and drafts while significantly reducing the bias when the information is used in other omics analyses. [ABSTRACT FROM AUTHOR]

Proteins, metabolites, and 16S rRNA measurements were used to examine the community structure and functional relationships within a cellulose degrading anaerobic bioreactor. The bioreactor was seeded with bovine rumen fluid and operated with a 4 day hydraulic retention time on cellulose (avicel) as sole carbon and energy source. The reactor performance and microbial community structure was monitored during the establishment of the cellulose-degrading community. After stable operation was established in the bioreactor, the mixing intensity was increased in order to investigate the effect of a physical disruption of the microbial community structure. Finally, the original conditions were re-established to understand the stability of the microbial community after a perturbation. All factors measured were found to be inter-correlated during these three distinct phases of operation (establishment, perturbation and re-establishment). In particular, the return of community structure and function to pre-perturbed conditions suggests that propionate fermentation and acetate utilization were the explanatory factors for community establishment and re-establishment. [ABSTRACT FROM AUTHOR]

Anthocyanins have significant functions in stress tolerance in pepper (Capsicum annuum L.) and also benefit human health. Nevertheless, the key structural genes and regulatory genes involved in anthocyanin accumulation in pepper fruits are still not well understood and fine mapped. For the present study, 383 F2 plants from a cross between the green-fruited C. annuum line Z5 and the purple-fruited line Z6 was developed. Two separate bulked DNA pools were constructed with DNAs extracted from either 37 plants with high anthocyanin content or from 18 plants with no anthocyanin. A combination of specific-locus amplified fragment sequencing (SLAF-seq) and bulked segregant analysis (BSA) was used to identify candidates for regions associated with anthocyanin accumulation. We identified a total of 127,004 high-quality single nucleotide polymorphism (SNP) markers, and detected 1674 high-quality SNP markers associated with anthocyanin accumulation. Three candidate anthocyanin-associated regions including the intervals from 12.48 to 20.00 Mb, from 54.67 to 56.59 Mb, and from 192.17 to 196.82 Mb were identified within a 14.10-Mb interval on chromosome 10 containing 126 candidate genes. Based on their annotations, we identified 12 candidate genes that are predicted to be associated with anthocyanin expression. The present results provide an efficient strategy for genetic mapping of and valuable insights into the genetic mechanisms of anthocyanin accumulation in pepper fruit, and allow us to clone and functionally analyze the genes that influence anthocyanin accumulation in this species. [ABSTRACT FROM AUTHOR]

Somatic mutations promote the transformation of normal cells to cancer. Accurate identification of such mutations facilitates cancer diagnosis and treatment, but biological and technological noises, including intra-tumor heterogeneity, sample contamination, uncertainties in base sequencing and read alignment, pose a big challenge to somatic mutation discovery. A number of callers have been developed to predict them from paired tumor/normal or unpaired tumor sequencing data. However, the small size of currently available experimentally validated somatic sites limits evaluation and then improvement of callers. Fortunately, NIST reference material NA12878 genome has been well-characterized with publicly available high-confidence genotype calls, and biological and technological noises can be computationally generalized to the number of sub-clones, the VAFs, the sequencing and mapping qualities. We used BAMSurgeon to create simulated tumors by introducing somatic small variants (SNVs and small indels) into homozygous reference or wildtype sites of NA12878. We generated 135 simulated tumors from 5 pre-tumors/normals. These simulated tumors vary in sequencing and subsequent mapping error profiles, read length, the number of sub-clones, the VAF, the mutation frequency across the genome and the genomic context. Furthermore, these pure tumor/normal pairs can be mixed at desired ratios within each pair to simulate sample contamination. This database (a total size of 15 terabytes) will be of great use to benchmark somatic small variant callers and guide their improvement. [ABSTRACT FROM AUTHOR]

Background: Eastern Equine Encephalitis (EEE) (Togaviridae, Alphavirus) is a highly pathogenic mosquito-borne arbovirus that circulates in an enzootic cycle involving Culiseta melanura mosquitoes and wild Passeriformes birds in freshwater swamp habitats. Recently, the northeastern United States has experienced an intensification of virus activity with increased human involvement and northward expansion into new regions. In addition to its principal role in enzootic transmission of EEE virus among avian hosts, recent studies on the blood-feeding behavior of Cs. melanura throughout its geographic range suggest that this mosquito may also be involved in epizootic / epidemic transmission to equines and humans in certain locales. Variations in blood feeding behavior may be a function of host availability, environmental factors, and/or underlying genetic differences among regional populations. Despite the importance of Cs. melanura in transmission and maintenance of EEE virus, the genetics of this species remains largely unexplored. Methodology and principle findings: To investigate the occurrence of genetic variation in Cs. melanura, the genome of this mosquito vector was sequenced resulting in a draft genome assembly of 1.28 gigabases with a contig N50 of 93.36 kilobases. Populations of Cs. melanura from 10 EEE virus foci in the eastern North America were genotyped with double-digest RAD-seq. Following alignment of reads to the reference genome, variant calling, and filtering, 40,384 SNPs were retained for downstream analyses. Subsequent analyses revealed genetic differentiation between northern and southern populations of this mosquito species. Moreover, limited fine-scale population structure was detected throughout northeastern North America, suggesting local differentiation of populations but also a history of ancestral polymorphism or contemporary gene flow. Additionally, a genetically distinct cluster was identified predominantly at two northern sites. Conclusion and significance: This study elucidates the first evidence of fine-scale population structure in Cs. melanura throughout its eastern range and detects evidence of gene flow between populations in northeastern North America. This investigation provides the groundwork for examining the consequences of genetic variations in the populations of this mosquito species that could influence vector-host interactions and the risk of human and equine infection with EEE virus. [ABSTRACT FROM AUTHOR]

The Qinghai-Tibet plateau is a natural plague focus and is the largest such focus in China. In this area, while Marmota himalayana is the primary host, a total of 18 human plague outbreaks associated with Tibetan sheep (78 cases with 47 deaths) have been reported on the Qinghai-Tibet plateau since 1956. All of the index infectious cases had an exposure history of slaughtering or skinning diseased or dead Tibetan sheep. In this study, we sequenced and compared 38 strains of Yersinia pestis isolated from different hosts, including humans, Tibetan sheep, and M. himalayana. Phylogenetic relationships were reconstructed based on genome-wide single-nucleotide polymorphisms identified from our isolates and reference strains. The phylogenetic relationships illustrated in our study, together with the finding that the Tibetan sheep plague clearly lagged behind the M. himalayana plague, and a previous study that identified the Tibetan sheep as a plague reservoir with high susceptibility and moderate sensitivity, indicated that the human plague was transmitted from Tibetan sheep, while the Tibetan sheep plague originated from marmots. Tibetan sheep may encounter this infection by contact with dead rodents or through being bitten by fleas originating from M. himalayana during local epizootics. [ABSTRACT FROM AUTHOR]

Onion (Allium cepa L.) is an economically important vegetable crop around the world. Genetic and genomic research into various onion accessions will provide insights into the onion genome to enhance breeding strategies and improve crops. However, the onion’s large genome size means that studies of molecular markers are limited in onion. This study aimed to discover high quality single nucleotide polymorphisms (SNPs) from 192 onion inbred lines relating to short-day cultivation in Korea. Paired-end (PE) double digested restriction site-associated DNA sequencing (ddRAD-seq) was used to discover SNPs in onion. A total of 538,973,706 reads (25.9 GB), with an average of 2,658,491 high-quality reads, were generated using ddRAD-seq. With stringent filtering, 1904 SNPs were discovered based on onion reference scaffolds. Further, population structure and genetic relationship studies suggested that two well-differentiated sub-populations exist in onion lines. SNP-associated flanking sequences were also compared with a public non-redundant database for gene ontology and pathway analysis. To our knowledge, this is the first report to identify high-quality SNPs in onion based on reference sequences using the ddRAD-seq platform. The SNP markers identified will be useful for breeders and the research community to deepen their understanding, enhance breeding programs, and support the management of onion genomic resources. [ABSTRACT FROM AUTHOR]

MicroRNAs (miRNAs) are small RNAs that regulate diverse biological processes including multiple aspects of the host-pathogen interface. Consequently, miRNAs are commonly encoded by viruses that undergo long-term persistent infection. Papillomaviruses (PVs) are capable of undergoing persistent infection, but as yet, no widely-accepted PV-encoded miRNAs have been described. The incomplete understanding of PV-encoded miRNAs is due in part to lack of tractable laboratory models for most PV types. To overcome this, we have developed RNA iscovery by forced enome xpression (miDGE), a new wet bench approach to miRNA identification that screens numerous pathogen genomes in parallel. Using miDGE, we screened over 73 different PV genomes for the ability to code for miRNAs. Our results show that most PVs are unlikely to code for miRNAs and we conclusively demonstrate a lack of PV miRNA expression in cancers associated with infections of several high risk HPVs. However, we identified five different high-confidence or highly probable miRNAs encoded by four different PVs (Human PVs 17, 37, 41 and a Fringilla coelebs PV (FcPV1)). Extensive in vitro assays confirm the validity of these miRNAs in cell culture and two FcPV1 miRNAs are further confirmed to be expressed in vivo in a natural host. We show that miRNAs from two PVs (HPV41 & FcPV1) are able to regulate viral transcripts corresponding to the early region of the PV genome. Combined, these findings identify the first canonical PV miRNAs and support that miRNAs of either host or viral origin are important regulators of the PV life cycle. [ABSTRACT FROM AUTHOR]

Paulinella micropora is a rhizarian thecate amoeba, belonging to a photosynthetic Paulinella species group that has a unique organelle termed chromatophore, whose cyanobacterial origin is distinct from that of plant and algal chloroplasts. Because acquisition of the chromatophore was quite a recent event compared with that of the chloroplast ancestor, the Paulinella species are thought to be model organisms for studying the early process of primary endosymbiosis. To obtain insight into how endosymbiotically transferred genes acquire expression competence in the host nucleus, here we analyzed the 5′ end sequences of the mRNAs of P. micropora MYN1 strain with the aid of a cap-trapper cDNA library. As a result, we found that mRNAs of 27 genes, including endosymbiotically transferred genes, possessed the common 5′ end sequence of 28–33 bases that were posttranscriptionally added by spliced leader (SL) trans-splicing. We also found two subtypes of SL RNA genes encoded by the P. micropora MYN1 genome. Differing from the other SL trans-splicing organisms that usually possess poly(A)-less SL RNAs, this amoeba has polyadenylated SL RNAs. In this study, we characterize the SL trans-splicing of this unique organism and discuss the putative merits of SL trans-splicing in functional gene transfer and genome evolution. [ABSTRACT FROM AUTHOR]

To analyze a specific genome region using next-generation sequencing technologies, the enrichment of DNA libraries with targeted capture methods has been standardized. For enrichment of mitochondrial genome, a previous study developed an original targeted capture method that use baits constructed from long-range polymerase chain reaction (PCR) amplicons, common laboratory reagents, and equipment. In this study, a new targeted capture method is presented, that of bacterial artificial chromosome (BAC) double capture (BDC), modifying the previous method, but using BAC libraries as baits for sequencing a relatively large gene. We applied the BDC approach for the 214 kb autosomal region, ring finger protein 213, which is the susceptibility gene of moyamoya disease (MMD). To evaluate the reliability of BDC, cost and data quality were compared with those of a commercial kit. While the ratio of duplicate reads was higher, the cost was less than that of the commercial kit. The data quality was sufficiently the same as that of the kit. Thus, BDC can be an easy, low-cost, and useful method for analyzing individual genome regions with substantial length. [ABSTRACT FROM AUTHOR]

Transcription, replication, and repair involve interactions of specific genomic loci with many different proteins. How these interactions are orchestrated at any given location and under changing cellular conditions is largely unknown because systematically measuring protein–DNA interactions at a specific locus in the genome is challenging. To address this problem, we developed Epi-Decoder, a Tag-chromatin immunoprecipitation-Barcode-Sequencing (TAG-ChIP-Barcode-Seq) technology in budding yeast. Epi-Decoder is orthogonal to proteomics approaches because it does not rely on mass spectrometry (MS) but instead takes advantage of DNA sequencing. Analysis of the proteome of a transcribed locus proximal to an origin of replication revealed more than 400 interacting proteins. Moreover, replication stress induced changes in local chromatin proteome composition prior to local origin firing, affecting replication proteins as well as transcription proteins. Finally, we show that native genomic loci can be decoded by efficient construction of barcode libraries assisted by clustered regularly interspaced short palindromic repeats and CRISPR-associated protein 9 (CRISPR/Cas9). Thus, Epi-Decoder is an effective strategy to identify and quantify in an unbiased and systematic manner the proteome of an individual genomic locus by DNA sequencing. [ABSTRACT FROM AUTHOR]

Simultaneously detecting CRISPR-based perturbations and induced transcriptional changes in the same cell is a powerful approach to unraveling genome function. Several lentiviral approaches have been developed, some of which rely on the detection of distally located genetic barcodes as an indirect proxy of sgRNA identity. Since barcodes are often several kilobases from their corresponding sgRNAs, viral recombination-mediated swapping of barcodes and sgRNAs is feasible. Using a self-circularization-based sgRNA-barcode library preparation protocol, we estimate the recombination rate to be ~50% and we trace this phenomenon to the pooled viral packaging step. Recombination is random, and decreases the signal-to-noise ratio of the assay. Our results suggest that alternative approaches can increase the throughput and sensitivity of single-cell perturbation assays. [ABSTRACT FROM AUTHOR]

Monitoring of wild animal populations is challenging, yet reliable information about population processes is important for both management and conservation efforts. Access to molecular markers, such as SNPs, enables population monitoring through genotyping of various DNA sources. We have developed 96 high quality SNP markers for individual identification of moose (Alces alces), an economically and ecologically important top-herbivore in boreal regions. Reduced representation libraries constructed from 34 moose were high-throughput de novo sequenced, generating nearly 50 million read pairs. About 50 000 stacks of aligned reads containing one or more SNPs were discovered with the Stacks pipeline. Several quality criteria were applied on the candidate SNPs to find markers informative on the individual level and well representative for the population. An empirical validation by genotyping of sequenced individuals and additional moose, resulted in the selection of a final panel of 86 high quality autosomal SNPs. Additionally, five sex-specific SNPs and five SNPs for sympatric species diagnostics are included in the panel. The genotyping error rate was 0.002 for the total panel and probability of identities were low enough to separate individuals with high confidence. Moreover, the autosomal SNPs were highly informative also for population level analyses. The potential applications of this SNP panel are thus many including investigations of population size, sex ratios, relatedness, reproductive success and population structure. Ideally, SNP-based studies could improve today’s population monitoring and increase our knowledge about moose population dynamics. [ABSTRACT FROM AUTHOR]

Actinidia arguta is the most basal species in a phylogenetically and economically important genus in the family Actinidiaceae. To better understand the molecular basis of the Actinidia arguta chloroplast (cp), we sequenced the complete cp genome from A. arguta using Illumina and PacBio RS II sequencing technologies. The cp genome from A. arguta was 157,611 bp in length and composed of a pair of 24,232 bp inverted repeats (IRs) separated by a 20,463 bp small single copy region (SSC) and an 88,684 bp large single copy region (LSC). Overall, the cp genome contained 113 unique genes. The cp genomes from A. arguta and three other Actinidia species from GenBank were subjected to a comparative analysis. Indel mutation events and high frequencies of base substitution were identified, and the accD and ycf2 genes showed a high degree of variation within Actinidia. Forty-seven simple sequence repeats (SSRs) and 155 repetitive structures were identified, further demonstrating the rapid evolution in Actinidia. The cp genome analysis and the identification of variable loci provide vital information for understanding the evolution and function of the chloroplast and for characterizing Actinidia population genetics. [ABSTRACT FROM AUTHOR]

A bacterial artificial chromosome (BAC) library for G. mustelinum Miers ex G. Watt (AD4) was constructed. Intact nuclei from G. mustelinum (AD4) were used to isolate high molecular weight DNA, which was partially cleaved with Hind III and cloned into pSMART BAC (Hind III) vectors. The BAC library consisted of 208,182 clones arrayed in 542 384-microtiter plates, with an average insert size of 121.72 kb ranging from 100 to 150 kb. About 2% of the clones did not contain inserts. Based on an estimated genome size of 2372 Mb for G. mustelinum, the BAC library was estimated to have a total coverage of 10.50 × genome equivalents. The high capacity library of G. mustelinum will serve as a giant gene resource for map-based cloning of quantitative trait loci or genes associated with important agronomic traits or resistance to Verticillium wilt, physical mapping and comparative genome analysis. [ABSTRACT FROM AUTHOR]

Japanese eel (Anguilla japonica) constitutes one of the most important food fish in Japan; accordingly, genome sequencing and linkage mapping have been conducted for the purpose of artificial cultivation. In the next stage, integration of genomic sequences within linkage groups (LG) is required to construct higher-resolution genetic markers for quantitative trait loci mapping and selective breeding of beneficial traits in farming. In order to identify LG1-linked scaffolds from the draft genome assembly (323,776 scaffolds) reported previously, we attempted to isolate chromosomes corresponding to LG1 by flow sorting and subsequent analyses. Initially, single chromosomes were randomly collected by chromosome sorting and subjected to whole-genome amplification (WGA). A total of 60 WGA samples were screened by PCR with primers for a known LG1-linked scaffold, and five positive WGA samples were sequenced by next-generation sequencing (NGS). Following reference mapping analysis of the NGS reads, four of the five WGA samples were found to be enriched by LG1-linked sequences. These samples were cytogenetically assigned to chromosome 5 by fluorescence in situ hybridization. Using blastn searches with 82,081 contigs constructed from the NGS reads of the four WGA samples as queries, 2323 scaffolds were identified as putative LG1-linked scaffolds from the draft genome assembly. The total length of the putative LG1-linked scaffolds was 99.0 Mb, comparable to the estimated DNA amounts of chromosome 5 (91.1 Mb). These results suggest that the methodology developed herein is applicable to isolate specific chromosome DNAs and integrate unanchored scaffold sequences onto a particular LG and chromosome even in teleost fishes, in which isolation of specific chromosomes by flow sorting is generally difficult owing to similar morphologies, sizes, and GC-contents among chromosomes in the genome. The putative LG1-linked scaffolds of Japanese eel contain a total of 6833 short tandem repeats which will be available for higher-resolution linkage mapping. [ABSTRACT FROM AUTHOR]

Cerasus humilis is endemic to China and is a new fruit tree species with economic and environmental benefits, with potential developmental and utilization applications. We report the first complete chloroplast genome sequence of C. humilis. Its genome is 158,084 bp in size, and the overall GC content is 36.8%. An inverted repeats (IR) of 52,672 bp in size is separated by a large single-copy (LSC) region of 86,374 bp and a small single-copy (SSC) region of 19,038 bp. The chloroplast genome of C. humilis contains 131 genes including 90 protein-coding genes, 33 transfer RNA genes, and 8 ribosomal RNA genes. The genome has a total 510 simple sequence repeats (SSRs). Of these, 306, 149, and 55 were found in the LSC, IR, and SSC regions, respectively. In addition, a comparison of the boundaries of the LSC, SSC, and IR regions of ten other Prunus species exhibited an overall high degree of sequence similarity, with slight variations in the IR boundary region which included gene deletions, insertions, expansions, and contractions. C. humilis lost the ycf1 gene at the IRA/SSC border and it has the largest ycf1 gene at the IRB/SSC border among these Prunus species, whereas the rps19 gene was inserted at the IRB/LSC junction. Furthermore, phylogenetic reconstruction using 61 conserved coding-protein genes clustered C. humilis with Prunus tomentosa. Thus, the complete chloroplast genome sequence of C. humilis provides a rich source of genetic information for studies on Prunus taxonomy, phylogeny, and evolution, as well as lays the foundation for further development and utilization of C. humilis. [ABSTRACT FROM AUTHOR]

The dried stigmas of Crocus sativus constitute the saffron, which is considered to be the costliest spice of the world. Saffron is valuable for its constituents, which are mainly apocarotenoids. In order to enhance the production of apocarotenoids, it is imperative to understand the regulation of apocarotenoid biosynthetic pathway. In C. sativus, although the pathway has been elucidated, the information regarding the regulation of the pathwaygenes is scanty. During the present investigation, the characterization of promoters regulating the expression of two important genes i.e. CsPSY and CsUGT was performed. We successfully cloned the promoters of both the genes, which were functionally characterized in Crocus sativus and Nicotiana tabaccum. In silico analysis of the promoters demonstrated the presence of several important cis regulatory elements responding tolight, hormonesand interaction with transcription factors (TFs). Further analysis suggested the regulation of CsPSY promoter by Abscisic acid (ABA) and that of CsUGT by Gibberellic acid (GA). In addition, we also observed ABA and GA mediated modulation in the expression of significant TFs and CsPSY and CsUGT transcripts. Overall, the study addresses issues related to regulation of key genes of apocarotenoid pathway in C.sativus. [ABSTRACT FROM AUTHOR]

Discovery of copy number variations (CNVs), a major category of structural variations, have dramatically changed our understanding of differences between individuals and provide an alternate paradigm for the genetic basis of human diseases. CNVs include both copy gain and copy loss events and their detection genome-wide is now possible using high-throughput, low-cost next generation sequencing (NGS) methods. However, accurate detection of CNVs from NGS data is not straightforward due to non-uniform coverage of reads resulting from various systemic biases. We have developed an integrated platform, iCopyDAV, to handle some of these issues in CNV detection in whole genome NGS data. It has a modular framework comprising five major modules: data pre-treatment, segmentation, variant calling, annotation and visualization. An important feature of iCopyDAV is the functional annotation module that enables the user to identify and prioritize CNVs encompassing various functional elements, genomic features and disease-associations. Parallelization of the segmentation algorithms makes the iCopyDAV platform even accessible on a desktop. Here we show the effect of sequencing coverage, read length, bin size, data pre-treatment and segmentation approaches on accurate detection of the complete spectrum of CNVs. Performance of iCopyDAV is evaluated on both simulated data and real data for different sequencing depths. It is an open-source integrated pipeline available at and as Docker’s image at . [ABSTRACT FROM AUTHOR]

The growing area of European hazelnut (Corylus avellana L.) is increasing, as well as the number of producing countries, and there is a pressing need for new improved cultivars. Hazelnut conventional breeding process is slow, due to the length of juvenile phase and the high heterozygosity level. The development of genetic linkage maps and the identification of molecular markers tightly linked to QTL (quantitative trait loci) of agronomic interest are essential tools for speeding up the selection of seedlings carrying desired traits through marker-assisted selection. The objectives of this study were to enrich a previous linkage map and confirm QTL related to time of leaf budburst, using an F1 population obtained by crossing Tonda Gentile delle Langhe with Merveille de Bollwiller. Genotyping-by-Sequencing was used to identify a total of 9,999 single nucleotide polymorphism markers. Well saturated linkage maps were constructed for each parent using the double pseudo-testcross mapping strategy. A reciprocal translocation was detected in Tonda Gentile delle Langhe between two non-homologous chromosomes. Applying a bioinformatic approach, we were able to disentangle ‘pseudo-linkage’ between markers, removing markers around the translocation breakpoints and obtain a linear order of the markers for the two chromosomes arms, for each linkage group involved in the translocation. Twenty-nine QTL for time of leaf budburst were identified, including a stably expressed region on LG_02 of the Tonda Gentile delle Langhe map. The stability of these QTL and their coding sequence content indicates promise for the identification of specific chromosomal regions carrying key genes involved in leaf budburst. [ABSTRACT FROM AUTHOR]

First flower node (FFN) is an important trait for evaluating fruit earliness in pepper (Capsicum annuum L.), but the genetic mechanisms that control FFN are still poorly understood. In the present study, we developed 249 F2 plants derived from an intraspecific cross between the inbred pepper lines Z4 and Z5. Thirty plants with the highest FFN and 30 plants with the lowest FFN were chosen and their DNAs were pooled according to phenotype to construct two bulked DNA pools. Specific-locus amplified fragment sequencing (SLAF-seq) was combined with bulked segregant analysis (BSA) to identify candidate regions related to FFN. According to our genetic analysis, the FFN trait is quantitatively inherited. A total of 106,848 high-quality single nucleotide polymorphism (SNP) markers were obtained, and 393 high-quality SNP markers associated with FFN were detected. Ten candidate regions within an interval of 3.98 Mb on chromosome 12 harboring 23 candidate genes were identified as closely correlated with FFN. Five genes (CA12g15130, CA12g15160, CA12g15370, CA12g15360, and CA12g15390) are predicted based on their annotations to be associated with expression of the FFN trait. The present study demonstrates an efficient genetic mapping strategy and lays a good foundation for molecular marker-assisted breeding using SNP markers linked to FFN and for cloning and functional analysis of the key genes controlling FFN. [ABSTRACT FROM AUTHOR]

Fritillaria spp. constitute important traditional Chinese medicinal plants. Xinjiang is one of two diversity hotspots in China in which eight Fritillaria species occur, two of which are endemic to the region. Furthermore, the phylogenetic relationships of Xinjiang Fritillaria species (including F. yuminensis) within the genus are unclear. In the present study, we sequenced the chloroplast (cp) genomes of seven Fritillaria species in Xinjiang using the Illumina HiSeq platform, with the aim of assessing the global structural patterns of the seven cp genomes and identifying highly variable cp DNA sequences. These were compared to previously sequenced Fritillaria cp genomes. Phylogenetic analysis was then used to evaluate the relationships of the Xinjiang species and assess the evolution of an undivided stigma. The seven cp genomes ranged from 151,764 to 152,112 bp, presenting a traditional quadripartite structure. The gene order and gene content of the seven cp genomes were identical. A comparison of the 13 cp genomes indicated that the structure is highly conserved. Ten highly divergent regions were identified that could be valuable in phylogenetic and population genetic studies. The phylogenetic relationships of the 13 Fritillaria species inferred from the protein-coding genes, large single-copy, small single-copy, and inverted repeat regions were identical and highly resolved. The phylogenetic relationships of the species corresponded with their geographic distribution patterns, in that the north group (consisting of eight species from Xinjiang and Heilongjiang in North China) and the south group (including six species from South China) were basically divided at 40°N. Species with an undivided stigma were not monophyletic, suggesting that this trait might have evolved several times in the genus. [ABSTRACT FROM AUTHOR]

Extremely low temperatures present various challenges to life that include ice formation and effects on metabolic capacity. Psyhcrophilic microorganisms typically have an array of mechanisms to enable survival in cold temperatures. In this study, we sequenced and analysed the genome of a psychrophilic yeast isolated in the Antarctic region, Glaciozyma antarctica. The genome annotation identified 7857 protein coding sequences. From the genome sequence analysis we were able to identify genes that encoded for proteins known to be associated with cold survival, in addition to annotating genes that are unique to G. antarctica. For genes that are known to be involved in cold adaptation such as anti-freeze proteins (AFPs), our gene expression analysis revealed that they were differentially transcribed over time and in response to different temperatures. This indicated the presence of an array of adaptation systems that can respond to a changing but persistent cold environment. We were also able to validate the activity of all the AFPs annotated where the recombinant AFPs demonstrated anti-freeze capacity. This work is an important foundation for further collective exploration into psychrophilic microbiology where among other potential, the genes unique to this species may represent a pool of novel mechanisms for cold survival. [ABSTRACT FROM AUTHOR]

Flax is an important crop for oil and fiber, however, no high-density genetic maps have been reported for this species. Specific length amplified fragment sequencing (SLAF-seq) is a high-resolution strategy for large scale de novo discovery and genotyping of single nucleotide polymorphisms. In this study, SLAF-seq was employed to develop SNP markers in an F2 population to construct a high-density genetic map for flax. In total, 196.29 million paired-end reads were obtained. The average sequencing depth was 25.08 in male parent, 32.17 in the female parent, and 9.64 in each F2 progeny. In total, 389,288 polymorphic SLAFs were detected, from which 260,380 polymorphic SNPs were developed. After filtering, 4,638 SNPs were found suitable for genetic map construction. The final genetic map included 4,145 SNP markers on 15 linkage groups and was 2,632.94 cM in length, with an average distance of 0.64 cM between adjacent markers. To our knowledge, this map is the densest SNP-based genetic map for flax. The SNP markers and genetic map reported in here will serve as a foundation for the fine mapping of quantitative trait loci (QTLs), map-based gene cloning and marker assisted selection (MAS) for flax. [ABSTRACT FROM AUTHOR]

Male sexual differentiation in crustaceans is controlled by the androgenic gland (AG), a unique male endocrine organ that, in decapods, is located at the base of the 5th pereiopod. In these animals, the insulin-like androgenic gland hormone (IAG) is the major factor secreted from the AG to induce masculinization and maintain male characteristics. It has, however, recently been proposed that this hormone also plays a role in growth and ovarian development in females. In this study, we tested such a possibility by searching for the IAG gene in the marbled crayfish, a parthenogenetic animal that reproduces asexually to form an all-female genetic clone. Based on the phylogenetic relationship between the marbled crayfish and Procambarus fallax, a gonochoristic species of the same North American Cambaridae family, we searched for the IAG gene in the marbled crayfish and then fully sequenced it. The open reading frame of the gene was found to be completely identical in the two species, and their introns shared over 94% identity. It was also found that, in addition to its expression at the base of the 5th pereiopod and in the testes of male P. fallax crayfish, IAG was expressed in the muscle tissue of P. fallax males and females and even of the parthenogenetic marbled crayfish. These findings provide new insight into possible functions of IAG, in addition to its role as a masculinization-inducing factor, and also constitute the basis for a discussion of the evolutionary relationship between the above two species. [ABSTRACT FROM AUTHOR]

Apart from sharing common ancestry with chordates, sea cucumbers exhibit a unique morphology and exceptional regenerative capacity. Here we present the complete genome sequence of an economically important sea cucumber, A. japonicus, generated using Illumina and PacBio platforms, to achieve an assembly of approximately 805 Mb (contig N50 of 190 Kb and scaffold N50 of 486 Kb), with 30,350 protein-coding genes and high continuity. We used this resource to explore key genetic mechanisms behind the unique biological characters of sea cucumbers. Phylogenetic and comparative genomic analyses revealed the presence of marker genes associated with notochord and gill slits, suggesting that these chordate features were present in ancestral echinoderms. The unique shape and weak mineralization of the sea cucumber adult body were also preliminarily explained by the contraction of biomineralization genes. Genome, transcriptome, and proteome analyses of organ regrowth after induced evisceration provided insight into the molecular underpinnings of visceral regeneration, including a specific tandem-duplicated prostatic secretory protein of 94 amino acids (PSP94)-like gene family and a significantly expanded fibrinogen-related protein (FREP) gene family. This high-quality genome resource will provide a useful framework for future research into biological processes and evolution in deuterostomes, including remarkable regenerative abilities that could have medical applications. Moreover, the multiomics data will be of prime value for commercial sea cucumber breeding programs. [ABSTRACT FROM AUTHOR]

Although recent developments in DNA sequencing have allowed for great leaps in both the quality and quantity of genome assembly projects, de novo assemblies still lack the efficiency and accuracy required for studying genetic variation of individuals. Thus, efficient and accurate methods for calling and genotyping genetic variants are fundamental to studying the genomes of individuals. We study the problem of genotyping insertion variants. We assume that the location of the insertion is given, and the task is to find the insertion sequence. Insertions are the hardest structural variant to genotype, because the insertion sequence must be assembled from the reads, whereas genotyping other structural variants only requires transformations of the reference genome. The current methods for constructing insertion variants are mostly linked to variation calling methods and are only able to construct small insertions. A sub-problem in genome assembly, the gap filling problem, provides techniques that are readily applicable to insertion genotyping. Gap filling takes the context and length of a missing sequence in a genome assembly and attempts to assemble the intervening sequence. In this paper we show how tools and methods for gap filling can be used to assemble insertion variants by modeling the problem of insertion genotyping as filling gaps in the reference genome. We further give a general read filtering scheme to make the method scalable to large data sets. Our results show that gap filling methods are competitive against insertion genotyping tools. We further show that read filtering improves performance of insertion genotyping especially for long insertions. Our experiments show that on long insertions the new proposed method is the most accurate one, whereas on short insertions it has comparable performance as compared against existing tools. [ABSTRACT FROM AUTHOR]

Bacterial strains of the same species collected from different hosts frequently exhibit differences in gene content. In the ubiquitous plant pathogen Pseudomonas syringae, more than 30% of genes encoded by each strain are not conserved among strains colonizing other host species. Although they are often implicated in host specificity, the role of this large fraction of the genome in host-specific adaptation is largely unexplored. Here, we sought to relate variation in gene content between strains infecting different species to variation that persists between strains on the same host. We fully sequenced a collection of P. syringae strains collected from wild Arabidopsis thaliana populations in the Midwestern United States. We then compared patterns of variation observed in gene content within these A. thaliana-isolated strains to previously published P. syringae sequence from strains collected on a diversity of crop species. We find that strains collected from the same host, A. thaliana, differ in gene content by 21%, 2/3 the level of gene content variation observed across strains collected from different hosts. Furthermore, the frequency with which specific genes are present among strains collected within the same host and among strains collected from different hosts is highly correlated. This implies that most gene content variation is maintained irrespective of host association. At the same time, we identify specific genes whose presence is important for P. syringae’s ability to flourish within A. thaliana. Specifically, the A. thaliana strains uniquely share a genomic island encoding toxins active against plants and surrounding microbes, suggesting a role for microbe-microbe interactions in dictating the abundance within this host. Overall, our results demonstrate that while variation in the presence of specific genes can affect the success of a pathogen within its host, the majority of gene content variation is not strongly associated with patterns of host use. [ABSTRACT FROM AUTHOR]