A genome-wide survey of switchgrass genome structure and organization

The perennial grass, switchgrass (Panicum virgatum L.), is a promising bioenergy crop and the target of whole genome sequencing. We constructed two bacterial artificial chromosome (BAC) libraries from the AP13 clone of switchgrass to gain insight into the genome structure and organization, initiate functional and comparative genomic studies, and assist with genome assembly. Together representing 16 haploid genome equivalents of switchgrass, each library comprises 101,376 clones with average insert sizes of 144 (HindIII-generated) and 110 kb (BstYI-generated). A total of 330,297 high quality BAC-end sequences (BES) were generated, accounting for 263.2 Mbp (16.4%) of the switchgrass genome. Analysis of the BES identified 279,099 known repetitive elements, >50,000 SSRs, and 2,528 novel repeat elements, named switchgrass repetitive elements (SREs). Comparative mapping of 47 full-length BAC sequences and 330K BES revealed high levels of synteny with the grass genomes sorghum, rice, maize, and Brachypodium. Our data indicate that the sorghum genome has retained larger microsyntenous regions with switchgrass besides high gene order conservation with rice. The resources generated in this effort will be useful for a broad range of applications.     

A first survey of the rye (Secale cereale) genome composition through BAC end sequencing of the short arm of chromosome 1R

Background

Thellungiella halophila (also known as Thellungiella salsuginea) is a model halophyte with a small plant size, short life cycle, and small genome. It easily undergoes genetic transformation by the floral dipping method used with its close relative, Arabidopsis thaliana. Thellungiella genes exhibit high sequence identity (approximately 90% at the cDNA level) with Arabidopsis genes. Furthermore, Thellungiella not only shows tolerance to extreme salinity stress, but also to chilling, freezing, and ozone stress, supporting the use of Thellungiella as a good genomic resource in studies of abiotic stress tolerance.

Results

We constructed a full-length enriched Thellungiella (Shan Dong ecotype) cDNA library from various tissues and whole plants subjected to environmental stresses, including high salinity, chilling, freezing, and abscisic acid treatment. We randomly selected about 20 000 clones and sequenced them from both ends to obtain a total of 35 171 sequences. CAP3 software was used to assemble the sequences and cluster them into 9569 nonredundant cDNA groups. We named these cDNAs "RTFL" (RIKEN Thellungiella Full-Length) cDNAs. Information on functional domains and Gene Ontology (GO) terms for the RTFL cDNAs were obtained using InterPro. The 8289 genes assigned to InterPro IDs were classified according to the GO terms using Plant GO Slim. Categorical comparison between the whole Arabidopsis genome and Thellungiella genes showing low identity to Arabidopsis genes revealed that the population of Thellungiella transport genes is approximately 1.5 times the size of the corresponding Arabidopsis genes. This suggests that these genes regulate a unique ion transportation system in Thellungiella.

Conclusion

As the number of Thellungiella halophila (Thellungiella salsuginea) expressed sequence tags (ESTs) was 9388 in July 2008, the number of ESTs has increased to approximately four times the original value as a result of this effort. Our sequences will thus contribute to correct future annotation of the Thellungiella genome sequence. The full-length enriched cDNA clones will enable the construction of overexpressing mutant plants by introduction of the cDNAs driven by a constitutive promoter, the complementation of Thellungiella mutants, and the determination of promoter regions in the Thellungiella genome.     

A set of BAC clones spanning the human genome

Using the human bacterial artificial chromosome (BAC) fingerprint-based physical map, genome sequence assembly and BAC end sequences, we have generated a fingerprint-validated set of 32855 BAC clones spanning the human genome. The clone set provides coverage for at least 98% of the human fingerprint map, 99% of the current assembled sequence and has an effective resolving power of 79 kb. We have made the clone set publicly available, anticipating that it will generally facilitate FISH or array-CGH-based identification and characterization of chromosomal alterations relevant to disease.     

Isolating promoters of multigene family members from the polyploid sugarcane genome by PCR-based walking in BAC DNA

The availability of a wider range of promoters for regulated expression in valuable transgenic crops would benefit functional genomics studies and current biotechnology programs aimed at improved productivity. Polymerase chain reaction (PCR)-based genome walking techniques are commonly used to isolate promoters or 5' flanking genomic regions adjacent to known cDNA sequences in genomes that are not yet completely sequenced. However, these techniques are problematic when applied directly to DNA isolated from crops with highly complex and large genomes. An adaptor ligation-mediated PCR-based BAC genome walking method is described here for the efficient isolation of promoters of multigene family members, such as the putative defense and fiber biosynthesis DIRIGENT genes that are abundant in the stem of sugarcane, a species with a highly polyploid genome. The advantage of this method is the efficient and specific amplification of the target promoter using BAC genomic DNA as template for the adaptor ligation-mediated PCR walking.     

A chromosome-scale assembly of the bilberry genome identifies a complex locus controlling berry anthocyanin composition

Bilberry (Vaccinium myrtillus L.) belongs to the Vaccinium genus, which includes blueberries (Vaccinium spp.) and cranberry (V. macrocarpon). Unlike its cultivated relatives, bilberry remains largely undomesticated, with berry harvesting almost entirely from the wild. As such, it represents an ideal target for genomic analysis, providing comparisons with the domesticated Vaccinium species. Bilberry is prized for its taste and health properties and has provided essential nutrition for Northern European indigenous populations. It contains high concentrations of phytonutrients, with perhaps the most important being the purple colored anthocyanins, found in both skin and flesh. Here, we present the first bilberry genome assembly, comprising 12 pseudochromosomes assembled using Oxford Nanopore (ONT) and Hi-C Technologies. The pseudochromosomes represent 96.6% complete BUSCO genes with an assessed LAI score of 16.3, showing a high conservation of synteny against the blueberry genome. Kmer analysis showed an unusual third peak, indicating the sequenced samples may have been from two individuals. The alternate alleles were purged so that the final assembly represents only one haplotype. A total of 36,404 genes were annotated after nearly 48% of the assembly was masked to remove repeats. To illustrate the genome quality, we describe the complex MYBA locus, and identify the key regulating MYB genes that determine anthocyanin production. The new bilberry genome builds on the genomic resources and knowledge of Vaccinium species, to help understand the genetics underpinning some of the quality attributes that breeding programs aspire to improve. The high conservation of synteny between bilberry and blueberry genomes means that comparative genome mapping can be applied to transfer knowledge about marker-trait association between these two species, as the loci involved in key characters are orthologous.     

A genome-wide analysis of FRT-like sequences in the human genome

Efficient and precise genome manipulations can be achieved by the Flp/FRT system of site-specific DNA recombination. Applications of this system are limited, however, to cases when target sites for Flp recombinase, FRT sites, are pre-introduced into a genome locale of interest. To expand use of the Flp/FRT system in genome engineering, variants of Flp recombinase can be evolved to recognize pre-existing genomic sequences that resemble FRT and thus can serve as recombination sites. To understand the distribution and sequence properties of genomic FRT-like sites, we performed a genome-wide analysis of FRT-like sites in the human genome using the experimentally-derived parameters. Out of 642,151 identified FRT-like sequences, 581,157 sequences were unique and 12,452 sequences had at least one exact duplicate. Duplicated FRT-like sequences are located mostly within LINE1, but also within LTRs of endogenous retroviruses, Alu repeats and other repetitive DNA sequences. The unique FRT-like sequences were classified based on the number of matches to FRT within the first four proximal bases pairs of the Flp binding elements of FRT and the nature of mismatched base pairs in the same region. The data obtained will be useful for the emerging field of genome engineering.     

A BAC clone fingerprinting approach to the detection of human genome rearrangements

We present a method, called fingerprint profiling (FPP), that uses restriction digest fingerprints of bacterial artificial chromosome clones to detect and classify rearrangements in the human genome. The approach uses alignment of experimental fingerprint patterns to in silico digests of the sequence assembly and is capable of detecting micro-deletions (1-5 kb) and balanced rearrangements. Our method has compelling potential for use as a whole-genome method for the identification and characterization of human genome rearrangements.     

A genome-wide association study of psoriasis and psoriatic arthritis identifies new disease loci

A genome-wide association study was performed to identify genetic factors involved in susceptibility to psoriasis (PS) and psoriatic arthritis (PSA), inflammatory diseases of the skin and joints in humans. 223 PS cases (including 91 with PSA) were genotyped with 311,398 single nucleotide polymorphisms (SNPs), and results were compared with those from 519 Northern European controls. Replications were performed with an independent cohort of 577 PS cases and 737 controls from the U.S., and 576 PSA patients and 480 controls from the U.K.. Strongest associations were with the class I region of the major histocompatibility complex (MHC). The most highly associated SNP was rs10484554, which lies 34.7 kb upstream from HLA-C (P = 7.8x10(-11), GWA scan; P = 1.8x10(-30), replication; P = 1.8x10(-39), combined; U.K. PSA: P = 6.9x10(-11)). However, rs2395029 encoding the G2V polymorphism within the class I gene HCP5 (combined P = 2.13x10(-26) in U.S. cases) yielded the highest ORs with both PS and PSA (4.1 and 3.2 respectively). This variant is associated with low viral set point following HIV infection and its effect is independent of rs10484554. We replicated the previously reported association with interleukin 23 receptor and interleukin 12B (IL12B) polymorphisms in PS and PSA cohorts (IL23R: rs11209026, U.S. PS, P = 1.4x10(-4); U.K. PSA: P = 8.0x10(-4); IL12B:rs6887695, U.S. PS, P = 5x10(-5) and U.K. PSA, P = 1.3x10(-3)) and detected an independent association in the IL23R region with a SNP 4 kb upstream from IL12RB2 (P = 0.001). Novel associations replicated in the U.S. PS cohort included the region harboring lipoma HMGIC fusion partner (LHFP) and conserved oligomeric golgi complex component 6 (COG6) genes on chromosome 13q13 (combined P = 2x10(-6) for rs7993214; OR = 0.71), the late cornified envelope gene cluster (LCE) from the Epidermal Differentiation Complex (PSORS4) (combined P = 6.2x10(-5) for rs6701216; OR 1.45) and a region of LD at 15q21 (combined P = 2.9x10(-5) for rs3803369; OR = 1.43). This region is of interest because it harbors ubiquitin-specific protease-8 whose processed pseudogene lies upstream from HLA-C. This region of 15q21 also harbors the gene for SPPL2A (signal peptide peptidase like 2a) which activates tumor necrosis factor alpha by cleavage, triggering the expression of IL12 in human dendritic cells. We also identified a novel PSA (and potentially PS) locus on chromosome 4q27. This region harbors the interleukin 2 (IL2) and interleukin 21 (IL21) genes and was recently shown to be associated with four autoimmune diseases (Celiac disease, Type 1 diabetes, Grave's disease and Rheumatoid Arthritis).     

A genome-wide survey of human thioredoxin and glutaredoxin family pseudogenes

The thioredoxin/glutaredoxin family consists of small heat-stable proteins that have a highly conserved CXXC active site and that participate in the regulation of many redox reactions. We have searched the human genome sequence to find putative pseudogenes (non-functional copies of protein-coding genes) for all known members of this family. This survey has resulted in the identification of seven processed pseudogenes for human Trx1 and two more for human Grx1. No evidence for the presence of processed pseudogenes has been found for the remaining members of this family. In addition, we have been unable to detect any non-processed pseudogenes derived from any member of the family in the human genome. The seven thioredoxin pseudogenes can be divided into two groups: Trx1-psi2, -psi3, -psi4, -psi5 and -psi6 arose from the functional ancestor, whereas Trx1-psi1 and -psi7 originated from Trx1-psi2 and -psi6, respectively. In all cases, the pseudogenes originated after the human/rodent radiation as shown by phylogenetic analysis. This is also the case for Grx1-psi1 and Grx1-psi2, which are placed between rodent and human sequences in the phylogenetic tree. Our study provides a molecular record of the recent evolution of these two genes in the hominid lineage.     

A genome-wide CRISPR screen identifies interactors of the autophagy pathway as conserved coronavirus targets

Over the past 20 years, 3 highly pathogenic human coronaviruses (HCoVs) have emerged—Severe Acute Respiratory Syndrome Coronavirus (SARS-CoV), Middle East Respiratory Syndrome Coronavirus (MERS-CoV), and, most recently, Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2)—demonstrating that coronaviruses (CoVs) pose a serious threat to human health and highlighting the importance of developing effective therapies against them. Similar to other viruses, CoVs are dependent on host factors for their survival and replication. We hypothesized that evolutionarily distinct CoVs may exploit similar host factors and pathways to support their replication cycles. Herein, we conducted 2 independent genome-wide CRISPR/Cas-9 knockout (KO) screens to identify MERS-CoV and HCoV-229E host dependency factors (HDFs) required for HCoV replication in the human Huh7 cell line. Top scoring genes were further validated and assessed in the context of MERS-CoV and HCoV-229E infection as well as SARS-CoV and SARS-CoV-2 infection. Strikingly, we found that several autophagy-related genes, including TMEM41B, MINAR1, and the immunophilin FKBP8, were common host factors required for pan-CoV replication. Importantly, inhibition of the immunophilin protein family with the compounds cyclosporine A, and the nonimmunosuppressive derivative alisporivir, resulted in dose-dependent inhibition of CoV replication in primary human nasal epithelial cell cultures, which recapitulate the natural site of virus replication. Overall, we identified host factors that are crucial for CoV replication and demonstrated that these factors constitute potential targets for therapeutic intervention by clinically approved drugs.

This study identifies essential host dependency factors for human coronavirus replication, showing that these can be directly targeted by clinically approved inhibitors and that treatment leads to effective inhibition of coronavirus replication in primary human nasal epithelial cell cultures.     

A genome-wide survey of photoreceptor and circadian genes in the coral, Acropora digitifera

Corals exhibit circadian behaviors, but little is known about the molecular mechanisms underlying the regulation of these behaviors. We surveyed the recently decoded genome of the coral, Acropora digitifera, for photoreceptor and circadian genes, using molecular phylogenetic analyses. Our search for photoreceptor genes yielded seven opsin and three cryptochrome genes. Two genes from each family likely underwent tandem duplication in the coral lineage. We also found the following A. digitifera orthologs to Drosophila and mammalian circadian clock genes: four clock, one bmal/cycle, three pdp1-like, one creb/atf, one sgg/zw3, two ck2alpha, one dco (csnk1d/cnsk1e), one slim/BTRC, and one grinl. No vrille, rev-ervα/nr1d1, bhlh2, vpac2, adcyap1, or adcyaplr1 orthologs were found. Intriguingly, in spite of an extensive survey, we also failed to find homologs of period and timeless, although we did find one timeout gene. In addition, the coral genes were compared to orthologous genes in the sea anemone, Nematostella vectensis. Thus, the coral and sea anemone genomes share a similar repertoire of circadian clock genes, although A. digitifera contains more clock genes and fewer photoreceptor genes than N. vectensis. This suggests that the circadian clock system was established in a common ancestor of corals and sea anemones, and was diversified by tandem gene duplications and the loss of paralogous genes in each lineage. It will be interesting to determine how the coral circadian clock functions without period.     

A comparative, BAC end sequence enabled map of the genome of the American mink (Neovison vison)

In this report we present the results of the analysis of approximately 2.7 Mb of genomic information for the American mink (Neovison vison) derived through BAC end sequencing. Our study, which encompasses approximately 1/1000^th of the mink genome, suggests that simple sequence repeats (SSRs) are less common in the mink than in the human genome, whereas the average GC content of the mink genome is slightly higher than that of its human counterpart. The 2.7 Mb mink genomic dataset also contained 2,416 repeat elements (retroids and DNA transposons) occupying almost 31% of the sequence space. Among repeat elements, LINEs were over-represented and endogenous viruses (aka LTRs) under-represented in comparison to the human genome. Finally, we present a virtual map of the mink genome constructed with reference to the human and canine genome assemblies using a comparative genomics approach and incorporating over 200 mink BESs with unique hits to the human genome.     

A genome-wide survey does not show the genetic distinctiveness of Basques

Basques are a cultural isolate, and, according to mainly allele frequencies of classical polymorphisms, also a genetic isolate. We investigated the differentiation of Spanish Basques from the rest of Iberian populations by means of a dense, genome-wide SNP array. We found that F _ST distances between Spanish Basques and other populations were similar to those between pairs of non-Basque populations. The same result is found in a PCA of individuals, showing a general distinction between Iberians and other South Europeans independently of being Basques. Pathogen-mediated natural selection may be responsible for the high differentiation previously reported for Basques at very specific genes such as ABO, RH, and HLA. Thus, Basques cannot be considered a genetic outlier under a general genome scope and interpretations on their origin may have to be revised.     

A set of simple-sequence repeat (SSR) markers covering the Prunus genome

A set of 109 microsatellite primer pairs recently developed for peach and cherry have been studied in the almond x peach F(2) progeny previously used to construct a saturated Prunus map containing mainly restriction fragment length polymorphism markers. All but one gave amplification products, and 87 (80%) segregated in the progeny and detected 96 loci. The resulting Prunus map contains a total of 342 markers covering a total distance of 522 cM. The approximate position of nine additional simple sequence repeats (SSRs) was established by comparison with other almond and peach maps. SSRs were placed in all the eight linkage groups of this map, and their distribution was relatively even, providing a genome-wide coverage with an average density of 5.4 cM/SSR. Twenty-four single-locus SSRs, highly polymorphic in peach, and each falling within 24 evenly spaced approximately 25-cM regions covering the whole Prunus genome, are proposed as a 'genotyping set' useful as a reference for fingerprinting, pedigree and genetic analysis of this species.     

A first generation physical map of the medaka genome in BACs essential for positional cloning and clone-by-clone based genomic sequencing

In order to realize the full potential of the medaka as a model system for developmental biology and genetics, characterized genomic resources need to be established, culminating in the sequence of the medaka genome. To facilitate the map-based cloning of genes underlying induced mutations and to provide templates for clone-based genomic sequencing, we have created a first-generation physical map of the medaka genome in bacterial artificial chromosome (BAC) clones. In particular, we exploited the synteny to the closely related genome of the pufferfish, Takifugu rubripes, by marker content mapping. As a first step, we clustered 103,144 public medaka EST sequences to obtain a set of 21,121 non-redundant sequence entities. Avoiding oversampling of gene-dense regions, 11,254 of EST clusters were successfully matched against the draft sequence of the fugu genome, and 2363 genes were selected for the BAC map project. We designed 35mer oligonucleotide probes from the selected genes and hybridized them against 64,500 BAC clones of strains Cab and Hd-rR, representing 14-fold coverage of the medaka genome. Our data set is further supplemented with 437 results generated from PCR-amplified inserts of medaka cDNA clones and BAC end-fragment markers. Our current, edited, first generation medaka BAC map consists of 902 map segments that cover about 74% of the medaka genome. The map contains 2721 markers. Of these, 2534 are from expressed sequences, equivalent to a non-redundant set of 2328 loci. The 934 markers (724 different) are anchored to the medaka genetic map. Thus, genetic map assignments provide immediate access to underlying clones and contigs, simplifying molecular access to candidate gene regions and their characterization.     

Sequence-ready 1-Mb YAC, BAC and cosmid contigs covering the distal imprinted region of mouse chromosome 7.

We have constructed approximately 1-Mb contigs of yeast artificial chromosome (YAC), bacterial artificial chromosome (BAC) and cosmid clones covering the imprinted region in mouse chromosome band 7F4/F5. This region is syntenic to human chromosome 11p15.5, which is associated with Beckwith-Wiedemann syndrome (BWS) and certain childhood and adult tumors. These contigs provide the basis for genomic sequencing, identification of genes and their regulatory elements, and functional studies in transgenic and knockout mice, which should be of help to understand not only the mechanisms of imprinting but also the molecular events involved in the genesis of BWS and tumors.     

BAC library construction and BAC end sequencing of five Drosophila species: the comparative map with the D. melanogaster genome

We constructed and characterized arrayed bacterial artificial chromosome (BAC) libraries of five Drosophila species (D. melanogaster, D. simulans, D. sechellia, D. auraria, and D. ananassae), which are genetically well characterized in the studies of meiosis, evolution, population genetics, and developmental biology. The BAC libraries comprise 8,000 to 12,500 clones for each species, estimated to cover the most of the genomes. We sequenced both ends of most of these BAC clones with a success rate of 91%. Of these, 53,701 clones consisting of non-repetitive BAC end sequences (BESs) were mapped with reference of the public D. melanogaster genome sequences. The BES mapping estimated that the BAC libraries of D. auraria and D. ananassae covered 47% and 57% of the D. melanogaster genome, respectively, and those of D. melanogaster, D. sechellia, and D. simulans covered 94-97%. The low coverage by BESs of D. auraria and D. ananassae may be due to the high sequence divergence with D. melanogaster. From the comparative BES mapping, 111 possible breakpoints of chromosomal rearrangements were identified in these four species. The breakpoints of the major chromosome rearrangement between D. simulans and D. melanogaster on the third chromosome were determined within 20 kb in 84E and 30 kb in 93E/F. Corresponding breakpoints were also identified in D. sechellia. The BAC clones described here will be an important addition to the Drosophila genomic resources.