Search

Bracoviruses (BVs) are endogenized nudiviruses in parasitoid wasps of the microgastroid complex (family Braconidae). Microgastroid wasps have coopted nudivirus genes to produce replication-defective virions that females use to transfer virulence genes to parasitized hosts. The microgastroid complex further consists of six subfamilies and ∼50,000 species but current understanding of BV gene inventories and organization primarily derives from analysis of two wasp species in the subfamily Microgastrinae ( Microplitis demolitor and Cotesia congregata ) that produce M. demolitor BV (MdBV) and C. congregata BV (CcBV). Notably, several genomic features of MdBV and CcBV remain conserved since divergence of M. demolitor and C. congregata ∼53 million years ago (MYA). However, it is unknown whether these conserved traits more broadly reflect BV evolution, because no complete genomes exist for any microgastroid wasps outside the Microgastrinae. In this regard, the subfamily Cheloninae is of greatest interest because it diverged earliest from the Microgastrinae (∼85 MYA) after endogenization of the nudivirus ancestor. Here, we present the complete genome of Chelonus insularis , which is an egg-larval parasitoid in the Cheloninae that produces C. insularis BV (CinsBV). We report that the inventory of nudivirus genes in C. insularis is conserved but are dissimilarly organized compared to M. demolitor and C. congregata . Reciprocally, CinsBV proviral segments share organizational features with MdBV and CcBV but virulence gene inventories exhibit almost no overlap. Altogether, our results point to the functional importance of a conserved inventory of nudivirus genes and a dynamic set of virulence genes for the successful parasitism of hosts. Our results also suggest organizational features previously identified in MdBV and CcBV are likely not essential for BV virion formation. IMPORTANCE Bracoviruses are a remarkable example of virus endogenization, because large sets of genes from a nudivirus ancestor continue to produce virions that thousands of wasp species rely upon to parasitize hosts. Understanding how these genes interact and have been coopted by wasps for novel functions is of broad interest in the study of virus evolution. This work characterizes bracovirus genome components in the parasitoid wasp Chelonus insularis , which together with existing wasp genomes captures a large portion of the diversity among wasp species that produce bracoviruses. Results provide new information about how bracovirus genome components are organized in different wasps while also providing additional insights on key features required for function.

Reorganization of topologically associated domain (TAD) is considered to be a novel mechanism for cell fate transitions. Here, we present a protocol to manipulate TAD via abscisic acid (ABA)-dependent genome linking. We use this protocol to merge two adjacent TADs and evaluate the influence on cell fate transitions. The advantages are that the manipulation does not change the genome and is reversible by withdrawing ABA. The major challenge is how to select linking loci for efficient TAD reorganization. For complete details on the use and execution of this protocol, please refer to Wang et al. (2021)., Competing Interests: The authors declare no competing interests., (© 2021 The Author(s).)

Within-species variation in genome size has been documented in many animals and plants. Despite its importance for understanding eukaryotic genome diversity, there is only sparse knowledge about how individual-level processes mediate genome size variation in populations. Here, we study a natural population of the rotifer Brachionus asplanchnoidis whose members differ up to 1.9-fold in diploid genome size, but were still able to interbreed and produce viable offspring. We show that genome size is highly heritable and can be artificially selected up or down, but not below a certain basal diploid genome size for this species. Analyses of segregation patterns in haploid males reveal that large genomic elements (several megabases in size) provide the substrate of genome size variation. These elements, and their segregation patterns, explain the generation of new genome size variants, the short-term evolutionary potential of genome size change in populations, and some seemingly paradoxical patterns, like an increase in genome size variation among highly inbred lines. Our study suggests that a conceptual model involving only two variables, 1) a basal genome size of the population, and 2) a vector containing information on additional elements that may increase genome size in this population (size, number, and meiotic segregation behavior), can effectively address most scenarios of short-term evolutionary change of genome size in a population., (© The Author(s) 2019. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Several general principles of global 3D genome organization have recently been established, including non-random positioning of chromosomes and genes in the cell nucleus, distinct chromatin compartments, and topologically associating domains (TADs). However, the extent and nature of cell-to-cell and cell-intrinsic variability in genome architecture are still poorly characterized. Here, we systematically probe heterogeneity in genome organization. High-throughput optical mapping of several hundred intra-chromosomal interactions in individual human fibroblasts demonstrates low association frequencies, which are determined by genomic distance, higher-order chromatin architecture, and chromatin environment. The structure of TADs is variable between individual cells, and inter-TAD associations are common. Furthermore, single-cell analysis reveals independent behavior of individual alleles in single nuclei. Our observations reveal extensive variability and heterogeneity in genome organization at the level of individual alleles and demonstrate the coexistence of a broad spectrum of genome configurations in a cell population., (Published by Elsevier Inc.)

Two novel HLA class II alleles, DRB4*03:01N and DQB1*03:276N, containing large deletions were identified during routine typing. Extraction of DNA encompassing the deletions was carried out with a panel of capture oligonucleotides followed by whole genome amplification. Next generation DNA sequencing was then used to characterize the sequences. DRB4*03:01N has a 16 kilobase pair deletion stretching upstream from intron 2 toward centromeric DRB8. DQB1*03:276N has two deletions separated by 844 nucleotides. The first deletion (3.7 kilobase pairs) is upstream of intron 1 and the second deletion removes 3.3 kilobase pairs further upstream towards centromeric DQA2., (Copyright © 2018 American Society for Histocompatibility and Immunogenetics. Published by Elsevier Inc. All rights reserved.)

It has been known for over a century that chromatin is not randomly distributed within the nucleus. However, the question of how DNA is folded and the influence of such folding on nuclear processes remain topics of intensive current research. A longstanding, unanswered question is whether nuclear organization is simply a reflection of nuclear processes such as transcription and replication, or whether chromatin is folded by independent mechanisms and this per se encodes function? Evidence is emerging that both may be true. Here, using the α-globin gene cluster as an illustrative model, we provide an overview of the most recent insights into the layers of genome organization across different scales and how this relates to gene activity., (© The Author 2017. Published by Oxford University Press.)

Selenocysteine (Sec) is known as the 21st amino acid, a cysteine analogue with selenium replacing sulphur. Sec is inserted co-translationally in a small fraction of proteins called selenoproteins. In selenoprotein genes, the Sec specific tRNA (tRNASec) drives the recoding of highly specific UGA codons from stop signals to Sec. Although found in organisms from the three domains of life, Sec is not universal. Many species are completely devoid of selenoprotein genes and lack the ability to synthesize Sec. Since tRNASec is a key component in selenoprotein biosynthesis, its efficient identification in genomes is instrumental to characterize the utilization of Sec across lineages. Available tRNA prediction methods fail to accurately predict tRNASec, due to its unusual structural fold. Here, we present Secmarker, a method based on manually curated covariance models capturing the specific tRNASec structure in archaea, bacteria and eukaryotes. We exploited the non-universality of Sec to build a proper benchmark set for tRNASec predictions, which is not possible for the predictions of other tRNAs. We show that Secmarker greatly improves the accuracy of previously existing methods constituting a valuable tool to identify tRNASec genes, and to efficiently determine whether a genome contains selenoproteins. We used Secmarker to analyze a large set of fully sequenced genomes, and the results revealed new insights in the biology of tRNASec, led to the discovery of a novel bacterial selenoprotein family, and shed additional light on the phylogenetic distribution of selenoprotein containing genomes. Secmarker is freely accessible for download, or online analysis through a web server at http://secmarker.crg.cat.

The genome is organized into repeating topologically associated domains (TADs), each of which is spatially isolated from its neighbor by poorly understood boundary elements thought to be conserved across cell types. Here, we show that deletion of CTCF (CCCTC-binding factor)-binding sites at TAD and sub-TAD topological boundaries that form within the HoxA and HoxC clusters during differentiation not only disturbs local chromatin domain organization and regulatory interactions but also results in homeotic transformations typical of Hox gene misregulation. Moreover, our data suggest that CTCF-dependent boundary function can be modulated by competing forces, such as the self-assembly of polycomb domains within the nucleus. Therefore, CTCF boundaries are not merely static structural components of the genome but instead are locally dynamic regulatory structures that control gene expression during development., (© 2016 Narendra et al.; Published by Cold Spring Harbor Laboratory Press.)

Tumor necrosis factor (TNF) and its superfamily (TNFSF) members are important inflammatory cytokines. Although fish have fourteen TNFSF genes, their genomic location and existence are yet to be described and confirmed in the Japanese pufferfish (Fugu) (Takifugu rubripes). Therefore, we conducted in silico identification, synteny analysis of TNFSF genes from Fugu with that of zebrafish and human TNFSF loci and their expression analysis in various tissues. We identified ten novel TNFSF genes, viz. TNFSF5 (CD40L), TNFSF6 (FasL), three TNFSF10 (TRAIL) (-1, 2 and 3), TNFSF11 (RANKlg), TNFSF12 (TWEAK), two TNFSF13B (BAFF) (1 and 2) and TNFSF14 (LIGHT) belonging to seven TNFSFs in Fugu. Several features such as existence of TNF family signature, conservation of genes in TNF loci with human and zebrafish chromosomes and phylogenetic clustering with other teleost TNFSF orthologs confirmed their identity. Fugu TNFSF genes were constitutively expressed in all eight different tissues with most of them expressed highly in liver. Fugu TNFSF10 gene has three homologs present on chromosomes 10 (TNFSF10-1), 8 (TNFSF10-2) and 2 (TNFSF10-3). Moreover, a phylogenetic analysis containing all available vertebrate (mammals, birds, reptiles, amphibians and fish) TNFSF10 orthologs showed that Fugu TNFSF10-1 and 10-3 are present in all vertebrates, whereas TNFSF10-2 was not related to any mammalian and avian orthologs. Viral double-stranded RNA mimic poly (I:C) caused an elevated expression of three Fugu TNFSF10 genes in head kidney cells at 4h indicating probable role of these genes to induce apoptosis in virus-infected cells. In conclusion, Fugu possesses genes belonging to nine TNFSFs including the newly identified seven and previously reported two, TNFSF New (TNF-N) and TNFSF2 (TNF-α). Our findings would add up information to TNFSF evolution among vertebrates., (Copyright © 2015 Elsevier B.V. All rights reserved.)

Bacillus sp. A053, isolated from the Arctic seawater, shows strong antifungal activity against plant pathogenic fungi. Here, we report the draft sequence of the approximately 4.1-Mbp genome of this strain. To the best of our knowledge, this is the first genome sequence of Bacillus strain isolated from the Arctic seawater with antifungal activity. The genome sequence may provide fundamental molecular information on elucidating the metabolic pathway of antimicrobial compound in this strain., (Copyright © 2015 Elsevier B.V. All rights reserved.)

Clostridium sulfidigenes 113A is a strictly anaerobic, rod shaped, gram positive bacterium isolated from sub-seafloor sediments associated with methane hydrates. Here, we report the first draft genome of C. sulfidigenes strain 113A, which comprises 3,717,420 bp in 96 contigs with the G+C content of 30.1%. A total of 3148 protein coding sequences were predicted. The genome annotation revealed that 113A could play an important role in biogeochemical cycles and have potential biotechnological applications such as production of organic acids and butanol., (Copyright © 2015 Elsevier B.V. All rights reserved.)

vB_DshP-R2C, a lytic phage that infects the marine bacterium Dinoroseobacter shibae DFL12(T), one of the model organisms of the Roseobacter clade, was isolated. Here we report the overall genome architecture of R2C. Genome analysis revealed that R2C is an N4-like phage with a 74.8 kb genome that contains 85 putative gene products., (Copyright © 2015. Published by Elsevier B.V.)

The complete mitochondrial genome of the intertidal barnacle Tetraclita serrata Darwin, 1854 (Crustacea: Maxillopoda: Sessilia) is presented. The genome is a circular molecule of 15,200 bp, which encodes 13 PCGs, 2 ribosomal RNA genes, and 22 transfer RNA genes. All non-coding regions are 591 bp in length, with the longest one speculated as the control region (389 bp), which is located between srRNA and trnK. The overall A+T content of the mitochondrial genome of T. serrata is 65.4%, which is lowest among all the eight mitochondrial genomes reported from sessile barnacles. There are variations of initiation and stop codons in the reported sessile barnacle mitochondrial genomes. Large-scale gene rearrangements are found in these genomes as compared to the pancrustacean ground pattern. ML and Bayesian analyses of all 15 complete mitochondrial genomes available from Maxillopoda lead to identical phylogenies. The phylogenetic tree based on mitochondrial PCGs shows that Argulus americanus (Branchiura) cluster with Armillifer armillatus (Pentastomida), distinct from all ten species from Cirripedia. Within the order Sessilia, Amphibalanus amphitrite (Balanidae) clusters with Striatobalanus amaryllis (Archaeobalanidae), and Nobia grandis (Pyrgomatidae). However, the two Megabalanus (Balanidae) are separated from the above grouping, resulting in non-monophyly of the family Balanidae. Moreover, the two Megabalanus have large-scale rearrangements as compared to the gene order shared by former three species. Therefore, both phylogenetic analysis using PCG sequences and gene order comparison suggest that Balanidae is not a monophyletic group. Given the limited taxa and moderate support values of the internal branches, the non-monophyly of the family Balanidae requires further verification., (Copyright © 2015 Elsevier B.V. All rights reserved.)

The complete mitochondrial (mt) DNA sequence of the rhodophyte Kappaphycus striatus ("Green Sacol" variety) was determined. The mtDNA is circular, 25,242 bases long (A+T content: 69.94%), and contains 50 densely packed genes comprising 93.22% of the mitochondrial genome, with genes encoded on both strands. Through comparative analysis, the overall sequence, genome structure, and organization of K. striatus mtDNA were seen to be highly similar with other fully sequenced mitochondrial genomes of the class Florideophyceae. On the other hand, certain degrees of genome rearrangements and greater sequence dissimilarities were observed for the mtDNAs of other evolutionarily distant red algae, such as those from the class Bangiophyceae and Cyanidiophyceae, compared to that of K. striatus. Furthermore, a trend was observed wherein the red algal mtDNAs tend to encode lesser number of protein-coding genes, albeit not necessarily shorter, as the organism becomes more morphologically complex. This trend is supported by the phylogenetic tree inferred from the concatenated amino acid sequences of the deduced protein products of cytochrome c oxidase subunit genes (cox1, 2, and 3)., (Copyright © 2014 Elsevier B.V. All rights reserved.)

Brevibacillus agri strain 5-2 was isolated from the formation water of a deep oil reservoir in Changqing Oilfield, China. This bacterium was found to have a capacity for degrading tetradecane, hexadecane and alkanesulfonate. To gain insights into its efficient metabolic pathway for degrading hydrocarbon and organosulfur compounds, here, we report the high quality draft genome of this strain. Two putative alkane 1-monooxygenases, one putative alkanesulfonate monooxygenase, one putative alkanesulfonate transporter, one putative sulfate permease and five putative sulfate transporters were identified in the draft genome. The genomic data of strain 5-2 may provide insights into the mechanism of microorganisms adapt to the petroleum reservoir after chemical flooding., (Copyright © 2014 Elsevier B.V. All rights reserved.)

The common carp, Cyprinus carpio, is one of the most important cyprinid species and globally accounts for 10% of freshwater aquaculture production. Here we present a draft genome of domesticated C. carpio (strain Songpu), whose current assembly contains 52,610 protein-coding genes and approximately 92.3% coverage of its paleotetraploidized genome (2n = 100). The latest round of whole-genome duplication has been estimated to have occurred approximately 8.2 million years ago. Genome resequencing of 33 representative individuals from worldwide populations demonstrates a single origin for C. carpio in 2 subspecies (C. carpio Haematopterus and C. carpio carpio). Integrative genomic and transcriptomic analyses were used to identify loci potentially associated with traits including scaling patterns and skin color. In combination with the high-resolution genetic map, the draft genome paves the way for better molecular studies and improved genome-assisted breeding of C. carpio and other closely related species.

Pyrenophora semeniperda (anamorph Drechslera campulata) is a necrotrophic fungal seed pathogen that has a wide host range within the Poaceae. One of its hosts is cheatgrass (Bromus tectorum), a species exotic to the United States that has invaded natural ecosystems of the Intermountain West. As a natural pathogen of cheatgrass, P. semeniperda has potential as a biocontrol agent due to its effectiveness at killing seeds within the seed bank; however, few genetic resources exist for the fungus. Here, the genome of P. semeniperda isolate assembled from sequence reads of 454 pyrosequencing is presented. The total assembly is 32.5 Mb and includes 11,453 gene models encoding putative proteins larger than 24 amino acids. The models represent a variety of putative genes that are involved in pathogenic pathways typically found in necrotrophic fungi. In addition, extensive rearrangements, including inter- and intrachromosomal rearrangements, were found when the P. semeniperda genome was compared to P. tritici-repentis, a related fungal species.

Background: The SBP-box gene family is specific to plants and encodes a class of zinc finger-containing transcription factors with a broad range of functions. Although SBP-box genes have been identified in numerous plants including green algae, moss, silver birch, snapdragon, Arabidopsis, rice and maize, there is little information concerning SBP-box genes, or the corresponding miR156/157, function in grapevine., Methodology/principal Findings: Eighteen SBP-box gene family members were identified in Vitis vinifera, twelve of which bore sequences that were complementary to miRNA156/157. Phylogenetic reconstruction demonstrated that plant SBP-domain proteins could be classified into seven subgroups, with the V. vinifera SBP-domain proteins being more closely related to SBP-domain proteins from dicotyledonous angiosperms than those from monocotyledonous angiosperms. In addition, synteny analysis between grape and Arabidopsis demonstrated that homologs of several grape SBP genes were found in corresponding syntenic blocks of Arabidopsis. Expression analysis of the grape SBP-box genes in various organs and at different stages of fruit development in V. quinquangularis 'Shang-24' revealed distinct spatiotemporal patterns. While the majority of the grape SBP-box genes lacking a miR156/157 target site were expressed ubiquitously and constitutively, most genes bearing a miR156/157 target site exhibited distinct expression patterns, possibly due to the inhibitory role of the microRNA. Furthermore, microarray data mining and quantitative real-time RT-PCR analysis identified several grape SBP-box genes that are potentially involved in the defense against biotic and abiotic stresses., Conclusion: The results presented here provide a further understanding of SBP-box gene function in plants, and yields additional insights into the mechanism of stress management in grape, which may have important implications for the future success of this crop.

The scaly-sided merganser (Mergus squamatus) is an endangered bird species on the IUCN Red List with the estimated global population of less than 2,500 individuals at present. In the present study, we studied the complete mitochondrial genome (mtDNA) and the phylogenetic of M. squamatus by PCR amplification and GenBank data. The genome was 16,595 bp in length and contained 37 genes (13 protein coding genes, two rRNAs, and 22 tRNAs) and a non-coding control region (D-loop). All protein-coding genes of M. squamatus mtDNA start with a typical ATG codon, except ND1, COI, and COII uses GTG as their initial codon. TAA, T- and TAG as the terminate codon occurred very commonly in the sequence. All tRNA genes can be folded into canonical cloverleaf secondary structure except for tRNA(Ser) (AGY) and tRNA(Leu) (CUN), which lose ''DHU'' arm. The genome sequences had been deposited in GenBank under accession number HQ833701. Based on the concatenated nucleotide sequences of mtDNA genes (Cyt b and D-loop), we reconstructed phylogenetic trees and discussed the phylogenetic relationships among ten Anatidae species. The results are different from the present classification, and we support Lophodytes cucullatus and Mergullus albellus to be members of the genus Mergus.

Flanking regulatory long terminal repeats (LTRs) in Human endogenous retrovirus (HERV) is a kind of typical DNA repeat that is widespread in the human genome. Currently, many algorithms have been developed to detect the latent periodicity of a wide range of DNA repeats. However, no such attempt was made for HERV LTRs. The present study focused on the investigation of the possible sequence periodic patterns in the HERV LTRs and their regulatory mechanisms. We calculated the sequence periods of 5', 3' and combined LTRs in HERVs with our devised matrix simulation algorithm. It is interesting that 5' and 3' LTRs have the same period of 7, and combined LTRs have a period of 9. These results indicated that HERV LTRs have predominant periodic patterns. Based on the obtained sequence periodicity, we constructed periodic consensus sequences of 5', 3' and combined LTRs. As to 5' and 3' LTRs with the same period - 7, we manually scanned the nucleotide bases in the corresponding positions of their periodic consensus sequences, and found some positions have the nucleotide base unchanged, such as the 1st, 5th and 7th positions. These conservative nucleotide base positions represent critical binding sites of regulatory LTRs, and may be indicative of conserved regulatory mechanisms in LRT-participating regulatory networks.

Noroviruses (NoVs) are one of the major causal agents of acute gastroenteritis in both industrial and developing countries including China. Recent studies have revealed that NoV genome is highly prone to mutation and recombination which may lead to emergence of new strains. In the present study, three full-length genomes of human NoV from China were determined and the genomic organization and recombination were analyzed. They had similar genome organization and contained three predicted ORFs, though the 5'UTR of those three strains were 2, 4 and 8 nucleotides, respectively. Phylogenetic analysis showed that the HU/GII/SHANGHAI/SH312/2008/CHN strain may be a recombinant of GII-3 capsid and GII-4 polymerase. To confirm the finding and detect the breakpoints where the recombination event occurred, we performed recombination analysis based on the genomic sequences of HU/GII/SHANGHAI/SH312/2008/CHN as the query sequence, and AB220921/NOV/JP/GII-4 and AB365435/NOV/US/GII-3 as the background sequences, using RPD software. Results indicated that the two parental strains were AB220921/NOV/JP/GII-4 and AB365435/NOV/US/GII-3. The breakpoint for this recombination event located at position 5,107 nt of the genome (in the ORF1 and ORF2 overlap).

Erythromelagia is a condition characterized by attacks of burning pain and inflammation in the extremeties. An epidemic form of this syndrome occurs in secondary students in rural China and a virus referred to as erythromelalgia-associated poxvirus (ERPV) was reported to have been recovered from throat swabs in 1987. Studies performed at the time suggested that ERPV belongs to the orthopoxvirus genus and has similarities with ectromelia virus, the causative agent of mousepox. We have determined the complete genome sequence of ERPV and demonstrated that it has 99.8% identity to the Naval strain of ectromelia virus and a slighly lower identity to the Moscow strain. Small DNA deletions in the Naval genome that are absent from ERPV may suggest that the sequenced strain of Naval was not the immediate progenitor of ERPV.

Mycoplasma bovis pneumonia in cattle has been epidemic in China since 2008. To investigate M. bovis pathogenesis, we completed genome sequencing of strain HB0801 isolated from a lesioned bovine lung from Hubei, China. The genomic plasticity was determined by comparing HB0801 with M. bovis strain ATCC® 25523™/PG45 from cow mastitis milk, Chinese strain Hubei-1 from lesioned lung tissue, and 16 other Mycoplasmas species. Compared to PG45, the genome size of HB0801 was reduced by 11.7 kb. Furthermore, a large chromosome inversion (580 kb) was confirmed in all Chinese isolates including HB0801, HB1007, a strain from cow mastitis milk, and Hubei-1. In addition, the variable surface lipoproteins (vsp) gene cluster existed in HB0801, but contained less than half of the genes, and had poor identity to that in PG45, but they had conserved structures. Further inter-strain comparisons revealed other mechanisms of gene acquisition and loss in HB0801 that primarily involved insertion sequence (IS) elements, integrative conjugative element, restriction and modification systems, and some lipoproteins and transmembrane proteins. Subsequently, PG45 and HB0801 virulence in cattle was compared. Results indicated that both strains were pathogenic to cattle. The scores of gross pathological assessment for the control group, and the PG45- and HB0801-infected groups were 3, 13 and 9, respectively. Meanwhile the scores of lung lesion for these three groups were 36, 70, and 69, respectively. In addition, immunohistochemistry detection demonstrated that both strains were similarly distributed in lungs and lymph nodes. Although PG45 showed slightly higher virulence in calves than HB0801, there was no statistical difference between the strains (P>0.05). Compared to Hubei-1, a total of 122 SNP loci were disclosed in HB0801. In conclusion, although genomic plasticity was thought to be an evolutionary advantage, it did not apparently affect virulence of M. bovis strains in cattle.

Background: Bryopsis hypnoides Lamouroux is a siphonous green alga, and its extruded protoplasm can aggregate spontaneously in seawater and develop into mature individuals. The chloroplast of B. hypnoides is the biggest organelle in the cell and shows strong autonomy. To better understand this organelle, we sequenced and analyzed the chloroplast genome of this green alga., Principal Findings: A total of 111 functional genes, including 69 potential protein-coding genes, 5 ribosomal RNA genes, and 37 tRNA genes were identified. The genome size (153,429 bp), arrangement, and inverted-repeat (IR)-lacking structure of the B. hypnoides chloroplast DNA (cpDNA) closely resembles that of Chlorella vulgaris. Furthermore, our cytogenomic investigations using pulsed-field gel electrophoresis (PFGE) and southern blotting methods showed that the B. hypnoides cpDNA had multimeric forms, including monomer, dimer, trimer, tetramer, and even higher multimers, which is similar to the higher order organization observed previously for higher plant cpDNA. The relative amounts of the four multimeric cpDNA forms were estimated to be about 1, 1/2, 1/4, and 1/8 based on molecular hybridization analysis. Phylogenetic analyses based on a concatenated alignment of chloroplast protein sequences suggested that B. hypnoides is sister to all Chlorophyceae and this placement received moderate support., Conclusion: All of the results suggest that the autonomy of the chloroplasts of B. hypnoides has little to do with the size and gene content of the cpDNA, and the IR-lacking structure of the chloroplasts indirectly demonstrated that the multimeric molecules might result from the random cleavage and fusion of replication intermediates instead of recombinational events.

The application of new and less biased methods to study the transcriptional output from genomes, such as tiling arrays and deep sequencing, has revealed that most of the genome is transcribed and that there is substantial overlap of transcripts derived from the two strands of DNA. In protein coding regions, the map of transcripts is very complex due to small transcripts from the flanking ends of the transcription unit, the use of multiple start and stop sites for the main transcript, production of multiple functional RNA molecules from the same primary transcript, and RNA molecules made by independent transcription from within the unit. In genomic regions separating those that encode proteins or highly abundant RNA molecules with known function, transcripts are generally of low abundance and short-lived. In most of these cases, it is unclear to what extent a function is related to transcription per se or to the RNA products.

Background: Among the founding mitochondrial haplogroups involved in the peopling of the Americas, haplogroup C1d has been viewed as problematic because of its phylogeny and because of the estimates of its antiquity, apparently being much younger than other founding haplogroups. Several recent analyses, based on data from the entire mitochondrial genome, have contributed to an advance in the resolution of these problems. The aim of our analysis is to compare the conclusions drawn from the available HVR-I and HVR-II data for haplogroup C1d with the ones based on whole mitochondrial genomes., Methodology/principal Findings: HVR-I and HVR-II sequences defined as belonging to haplogroup C1d by standard criteria were gathered from the literature as well as from population studies carried out in Uruguay. Sequence phylogeny was reconstructed using median-joining networks, geographic distribution of lineages was analyzed and the age of the most recent common ancestor estimated using the ρ-statistic and two different mutation rates. The putative ancestral forms of the haplogroup were found to be more widespread than the derived lineages, and the lineages defined by np 194 were found to be widely distributed and of equivalent age., Conclusions/significance: The analysis of control region sequences is found to still harbor great potential in tracing microevolutionary phenomena, especially those found to have occurred in more recent times. Based on the geographic distributions of the alleles of np 7697 and np 194, both discussed as possible basal mutations of the C1d phylogeny, we suggest that both alleles were part of the variability of the haplogroup at the time of its entrance. Moreover, based on the mutation rates of the different sites stated to be diagnostic, it is possible that the anomalies found when analyzing the haplogroup are due to paraphyly.

Background: Mobile elements (MEs) are diverse, common and dynamic inhabitants of nearly all genomes. ME transposition generates a steady stream of polymorphic genetic markers, deleterious and adaptive mutations, and substrates for further genomic rearrangements. Research on the impacts, population dynamics, and evolution of MEs is constrained by the difficulty of ascertaining rare polymorphic ME insertions that occur against a large background of pre-existing fixed elements and then genotyping them in many individuals., Results: Here we present a novel method for identifying nearly all insertions of a ME subfamily in the whole genomes of multiple individuals and simultaneously genotyping (for presence or absence) those insertions that are variable in the population. We use ME-specific primers to construct DNA libraries that contain the junctions of all ME insertions of the subfamily, with their flanking genomic sequences, from many individuals. Individual-specific "index" sequences are designed into the oligonucleotide adapters used to construct the individual libraries. These libraries are then pooled and sequenced using a ME-specific sequencing primer. Mobile element insertion loci of the target subfamily are uniquely identified by their junction sequence, and all insertion junctions are linked to their individual libraries by the corresponding index sequence. To test this method's feasibility, we apply it to the human AluYb8 and AluYb9 subfamilies. In four individuals, we identified a total of 2,758 AluYb8 and AluYb9 insertions, including nearly all those that are present in the reference genome, as well as 487 that are not. Index counts show the sequenced products from each sample reflect the intended proportions to within 1%. At a sequencing depth of 355,000 paired reads per sample, the sensitivity and specificity of ME-Scan are both approximately 95%., Conclusions: Mobile Element Scanning (ME-Scan) is an efficient method for quickly genotyping mobile element insertions with very high sensitivity and specificity. In light of recent improvements to high-throughput sequencing technology, it should be possible to employ ME-Scan to genotype insertions of almost any mobile element family in many individuals from any species.

Debaryomyces hansenii, a yeast that participates in the elaboration of foodstuff, displays important genetic diversity. Our recent phylogenetic classification of this species led to the subdivision of the species into three distinct clades. D. hansenii harbors the highest number of nuclear mitochondrial DNA (NUMT) insertions known so far for hemiascomycetous yeasts. Here we assessed the intraspecific variability of the NUMTs in this species by testing their presence/absence first in 28 strains, with 21 loci previously detected in the completely sequenced strain CBS 767(T), and second in a larger panel of 77 strains, with 8 most informative loci. We were able for the first time to structure populations in D. hansenii, although we observed little NUMT insertion variability within the clades. We determined the chronology of the NUMT insertions, which turned out to correlate with the previously defined taxonomy and provided additional evidence that colonization of nuclear genomes by mitochondrial DNA is a dynamic process in yeast. In combination with flow cytometry experiments, the NUMT analysis revealed the existence of both haploid and diploid strains, the latter being heterozygous and resulting from at least four crosses among strains from the various clades. As in the diploid pathogen Candida albicans, to which D. hansenii is phylogenetically related, we observed a differential loss of heterozygosity in the diploid strains, which can explain some of the large genetic diversity found in D. hansenii over the years.

Genome dose of active ribosome genes (ARG), average of nucleolus argyrophil structures in lymphocyte nuclei, levels of extracellular DNA (DNA(e)) concentrations and ratio of antibodies to total DNA (AB(DNA)) and ribosomal DNA (AB(DNA-rib)), and nuclease activity were determined in peripheral blood of 8 volunteered subjects (21-26 y.o.) in the experiment with 7-d DI. Results of the investigation revealed a broad individual variability ensued from heterogeneity of the group of the test-subjects as to ARG values. There was an inverse negative relationship between ARG values and increment of the ribosome genes activity index. Part of the subjected exhibited increased DNA(e) levels on completion of the experiment, whereas the others decreased the parameter demonstrating individual character of body reaction. No correlation was established between DNA(e) content and nuclease activity in blood. Concentrations of AB(DNA) and DNA AB(DNA-rib) before and after immersion were essentially unchanged; however, they were higher as compared with the control group of blood-donors. Diversity of subjects' reactions was accounted to the broad range of ARG values. Therefore, selection of test-subjects for ground-based simulation experiments should be conducted with due consideration of the parameter.

Expressed sequence tags (ESTs) represent 500-1000-bp-long sequences corresponding to mRNAs derived from different sources (cell lines, tissues, etc.). The human EST database contains over 8,000,000 sequences, with over 4,000,000,000 total nucleotides. RNA molecules are transcribed from a genomic DNA template; therefore, all ESTs should match corresponding genomes. Nevertheless, we have found in the human EST database approximately 11,000 ESTs not matching sequences in the human genome database. The presence of "trash" ESTs (TESTs) in the EST database could result from DNA or RNA contamination of the laboratory equipment, tissues, or cell lines. TESTs could also represent sequences from unidentified human genes or from species inhabiting the human body. Here, we attempt to identify the sources of human EST database contaminations. In particular, we discuss systematic contamination of the mammalian EST databases with sequences of plants.

Despite surprisingly a small number of protein-coding gene in mammalian genomes, a large variety of different RNAs is being produced. These RNAs are amazingly different in their number, size, cell localization, and mechanism of actions. Although new classes of short RNAs (sRNAs) are being continuously discovered, it is not yet obvious how many of the sRNAs are originated. Altogether, the research in the recent few years has identified an unexpectedly rich variety of mechanisms by which noncoding RNAs act, suggesting that we have identified probably only few of the many potential functional mechanism and more investigation will be needed to comprehensively understand the complex nature and biology of mammalian RNAome. Here, we focus on various aspects of the diversity of the biological role of these nonprotein-coding RNAs (ncRNAs), with emphasis on functional mechanisms recently elucidated.

One way to modulate transcription is by partitioning the chromatin fiber within the nucleus into the active or inactive domains through the establishment of higher-order chromatin structure. Such subdivision of chromatin implies the existence of insulators and boundaries that delimit differentially regulated chromosomal loci. Recently published data on transcriptional interference from the repeated component of the genome fits the classic definition of insulator/boundary activity. This review discusses the phenomena of transcriptional interference and raises the question about functionality of genomic "junk" along with the need to stimulate a dialogue on how we would define the insulators and boundaries in the light of contemporary data. Rule 19 (a) (Boundaries)"Before the toss, the umpires shall agree the boundary of the field of play with both captains. The boundary shall, if possible, be marked along its whole length" Rules of Cricket.

We report the complete cDNA sequence of the Anopheles gambiae voltage-gated sodium channel (VGSC) alpha-subunit isolated from mature adult mosquitoes. The genomic DNA contains 35 deduced exons with a predicted translation of

Several lines of evidence suggest that reiterated sequences in the human genome are targets for nonallelic homologous recombination (NAHR), which facilitates genomic rearrangements. We have used a PCR-based approach to identify breakpoint regions of rearranged structures in the human genome. In particular, we have identified intrachromosomal identical repeats that are located in reverse orientation, which may lead to chromosomal inversions. A bioinformatic workflow pathway to select appropriate regions for analysis was developed. Three such regions overlapping with known human genes, located on chromosomes 3, 15, and 19, were analyzed. The relative proportion of wild-type to rearranged structures was determined in DNA samples from blood obtained from different, unrelated individuals. The results obtained indicate that recurrent genomic rearrangements occur at relatively high frequency in somatic cells. Interestingly, the rearrangements studied were significantly more abundant in adults than in newborn individuals, suggesting that such DNA rearrangements might start to appear during embryogenesis or fetal life and continue to accumulate after birth. The relevance of our results in regard to human genomic variation is discussed.

Research programmes for constructing a 'cell factory' have been funded in several countries. In Japan, the 'Minimum genome factory' (MGF) project was launched in 2001. In this project, several model microbes have been genetically reconstructed to obtain a cell with fewer genes on a chromosome of reduced size. A microbe with a 'minimum genome' is expected to exhibit less regulation and therefore to be an ideal platform for a cell-factory system. The goal of this project is to construct such a minimum genome microbe for a cell factory. In this project, the 4.6 Mbp genome of Escherichia coli K-12 has been successfully reduced to 3.6 Mbp. The constructed reduced-genome strain, MGF-01, shows better growth and higher threonine production compared with the wild-type strain. Furthermore functional analyses of all E. coli genes have also been performed. CGH (comparative genomic hybridization) analysis revealed that about 2600 genes were commonly conserved in the 23 E. coli strains tested. This set of conserved genes was hypothesized as a core set for E. coli species. Phenotype array analysis of a nearly complete collection of single-gene knockout mutants of E. coli provided insights into E. coli metabolic networks. The data sets from the functional genomics will be used to improve design of an E. coli MGF. The present minireview summarizes the progress of the E. coli MGF project and overviews related research.

In 1997, the complete genomic DNA sequence of Bacillus subtilis (4.2 Mbp) was determined and 4100 genes were identified [Kunst, Ogasawara, Moszer, Albertini, Alloni, Azevedo, Bertero, Bessieres, Bolotin, Borchert, S. et al. (1997) Nature 90, 249-256]. In addition, B. subtilis, which shows an excellent ability to secrete proteins (enzymes) and antibiotics in large quantities outside the cell, plays an important role in industrial and medical fields. It is necessary to clarify the genes involved in the production of compounds by understanding the network of these 4100 genes and the proceeding analysis of genes of unknown functions. In promoting such a study, it is expected that the regulatory system of B. subtilis can be simplified by the creation of a Bacillus strain with a reduced genome by discriminating genes unnecessary for the production of proteins from essential genes, and deleting as many of these unnecessary genes as possible, which may help to understand this complex network of genes. We have previously distinguished essential and non-essential genes by evaluating the growth and enzyme-producing properties of strains of B. subtilis in which about 3000 genes (except 271 essential genes) have been disrupted or deleted singly, and have successfully utilized the findings from these studies in creating the MG1M strain with an approx. 1 Mbp deletion by serially deleting 17 unnecessary regions from the genome. This strain showed slightly reduced growth in enzyme-production medium, but no marked morphological changes. Moreover, we confirmed that the MG1M strain had cellulase and protease productivity comparable with that of the B. subtilis 168 strain, thus demonstrating that genome reduction does not contribute to a negative influence on enzyme productivity.

The phenotypic effects of random mutations depend on both the architecture of the genome and the gene-trait relationships. Both levels thus play a key role in the mutational variability of the phenotype, and hence in the long-term evolutionary success of the lineage. Here, by simulating the evolution of organisms with flexible genomes, we show that the need for an appropriate phenotypic variability induces a relationship between the deleteriousness of gene mutations and the quantity of non-coding sequences maintained in the genome. The more deleterious the gene mutations, the shorter the intergenic sequences. Indeed, in a shorter genome, fewer genes are affected by rearrangements (duplications, deletions, inversions, translocations) at each replication, which compensates for the higher impact of each gene mutation. This spontaneous adjustment of genome structure allows the organisms to retain the same average fitness loss per replication, despite the higher impact of single gene mutations. These results show how evolution can generate unexpected couplings between distinct organization levels.

Evolutionary conservation of DNA sequences provides a tool for the identification of functional elements in genomes. We have created a database of evolutionary conserved regions (ECRs) in vertebrate genomes, entitled ECRbase, which is constructed from a collection of whole-genome alignments produced by the ECR Browser. ECRbase features a database of syntenic blocks that recapitulate the evolution of rearrangements in vertebrates and a comprehensive collection of promoters in all vertebrate genomes generated using multiple sources of gene annotation. The database also contains a collection of annotated transcription factor binding sites (TFBSs) in evolutionary conserved and promoter elements. ECRbase currently includes human, rhesus macaque, dog, opossum, rat, mouse, chicken, frog, zebrafish and fugu genomes. It is freely accessible at http://ecrbase.dcode.org.

Pharmacologic studies have revealed that cysteinyl leukotrienes (CYSLTs) act through two receptors, cysteinyl leukotriene receptor 1 (CYSLTR1) and CYSLTR2. CYSLTR1 antagonists are widely used to treat asthma and rhinitis. In this study, we characterized the genomic structure and transcriptional regulation of CYSLTR1 and examined associations between CYSLTR1 polymorphisms and asthma/rhinitis. The experiment of rapid amplification of cDNA end revealed that CYSLTR1 contains three exons and that the entire open reading frame is located in exon 3. Reverse transcriptase-polymerase chain reaction showed that there were multiple splice variants of CYSLTR1 and that the transcript expression patterns differed from tissues and cell types. The promoter region of CYSLTR1 is from -665 to -30 bp relative to the transcription start site. We identified four polymorphisms (c.-618-434T/C, c.-618-275C/A, c.-618-136G/A, and 927C/T), and transmission disequilibrium tests revealed that none of these polymorphisms was associated with the development of asthma/rhinitis. However, the TCG and CAA haplotypes in the promoter region caused different transcriptional activity. Our findings indicate that CYSLTR1 polymorphisms are not likely to be involved in the development of asthma/rhinitis, but it is possible that these polymorphisms could influence drug responses in individuals with atopic diseases.

Previous research reveals that the genome structures of rhizobial type strains and reference strains can reflect their phylogenetic relationships. In order to further explore the potential application of genome structure as a phylogenetic marker in rhizobial natural taxonomy, this study analyzed the genome structures of 29 unclassified nodule bacteria isolated from the root nodules of leguminous trees, Robinia sp., Dalbergia spp., and A lbizia spp. and 7 rhizobial reference strains by I-CeuI cleavage, then clustered these bacteria phylogenetically based on their genome structures and compared these clusters with those based on numerical taxonomy and 16S rDNA PCR-RFLP. Eleven phylogenetic clusters were obtained. The clusters were in large part consistent with those based on numerical taxonomy and 16S rDNA PCR-RFLP. Also there are inconsistent clusters based on the above three methods. But results are completely consistent with 16S rRNA clusters. This suggested that the genome structure clustering method can be used to fastly identify root nodule isolates and detect their phylogenetic relationships. The credibility and repeatability of the results, together with the simplicity and possibility to analyze a large number of strains in a short time of the method, indicates the broad potential application of genome structure as phylogenetic marker to categorize rhizobial isolates and should in the future facilitate biodiversity studies.

Botulinum neurotoxins (BoNTXs) produced by Clostridium botulinum are among the most poisonous substances known. Of the seven types of BoNTXs, genes for type C1 and D toxins (BoNTX/C1 and D) are carried by bacteriophages. The gene for exoenzyme C3 also resides on these phages. Here, we present the complete genome sequence of c-st, a representative of BoNTX/C1-converting phages. The genome is a linear double-stranded DNA of 185,682 bp with 404-bp terminal direct repeats, the largest known temperate phage genome. We identified 198 potential protein-coding regions, including the genes for production of BoNTX/C1 and exoenzyme C3. Very exceptionally, as a viable bacteriophage, a number of insertion sequences were found on the c-st genome. By analyzing the molecular structure of the c-st genome in lysogens, we also found that it exists as a circular plasmid prophage. These features account for the unstable lysogeny of BoNTX phages, which has historically been called "pseudolysogeny." The PCR scanning analysis of other BoNTX/C1 and D phages based on the c-st sequence further revealed that BoNTX phages comprise a divergent phage family, probably generated by exchanging genomic segments among BoNTX phages and their relatives.

Until recently the study of individual DNA sequences and of total DNA content (the C-value) sat at opposite ends of the spectrum in genome biology. For gene sequencers, the vast stretches of non-coding DNA found in eukaryotic genomes were largely considered to be an annoyance, whereas genome-size researchers attributed little relevance to specific nucleotide sequences. However, the dawn of comprehensive genome sequencing has allowed a new synergy between these fields, with sequence data providing novel insights into genome-size evolution, and with genome-size data being of both practical and theoretical significance for large-scale sequence analysis. In combination, these formerly disconnected disciplines are poised to deliver a greatly improved understanding of genome structure and evolution.

Microsporidia have been known for some time to possess among the smallest genomes of any eukaryote. There is now a completely sequenced microsporidian genome, as well as several other large-scale sequencing efforts, so the nature of these genomes is becoming apparent. This paper reviews some of the characteristics of microsporidian genomes in general, and some of the recent discoveries made through comparative genomic analyses. In general, microsporidian genomes are both reduced and compacted. Reduction takes place through gene loss, which is understandable in obligate intracellular parasites that rely on their host for many metabolites. Compaction is a more complex process, and is as yet not fully understood. It is clear from genomes surveyed thus far that the remaining genes are tightly packed and that there is little non-coding sequence, resulting in some extraordinary arrangements, including overlapping genes. Compaction also seems to affect certain aspects of genome evolution, like the frequency of rearrangements. The force behind this compaction is not known, and is especially interesting in light of the fact that surveys of genomes that are significantly different in size yield similar complements of protein-coding genes. There are some interesting exceptions, including catalase, photolyase and some mitochondrial proteins, but the rarity of these raises an interesting question as to what accounts for the significant differences seen in the genome sizes among microsporidia.

The genome sequence of the microsporidian parasite Encephalitozoon cuniculi Levaditi, Nicolau et Schoen, 1923 contains about 2,000 genes that are representative of a non-redundant potential proteome composed of 1,909 protein chains. The purpose of this review is to relate some advances in the characterisation of this proteome through bioinformatics and experimental approaches. The reduced diversity of the set of E. cuniculi proteins is perceptible in all the compilations of predicted domains, orthologs, families and superfamilies, available in several public databases. The phyletic patterns of orthologs for seven eukaryotic organisms support an extensive gene loss in the fungal clade, with additional deletions in E. cuniculi. Most microsporidial orthologs are the smallest ones among eukaryotes, justifying an interest in the use of these compacted proteins to better discriminate between essential and non-essential regions. The three components of the E. cuniculi mRNA capping apparatus have been especially well characterized and the three-dimensional structure of the cap methyltransferase has been elucidated following the crystallisation of the microsporidial enzyme Ecm1. So far, our mass spectrometry-based analyses of the E. cuniculi spore proteome has led to the identification of about 170 proteins, one-quarter of these having no clearly predicted function. Immunocytochemical studies are in progress to determine the subcellular localisation of microsporidia-specific proteins. Post-translational modifications such as phosphorylation and glycosylation are expected to be soon explored.

Plasma insulin and glucose concentrations are important quantitative phenotypes related to diabetes and the metabolic syndrome. Reports purporting to identify quantitative trait loci (QTLs) that contribute to the variation in fasting insulin and glucose concentrations are discrepant. As part of the Insulin Resistance Atherosclerosis Study (IRAS) Family Study, a genome scan was performed in African-American (n = 42) and Hispanic (n = 90) extended families to identify regions that may contain positional candidate genes for fasting insulin and fasting glucose (n = 1,604 subjects). There was significant evidence for linkage of fasting insulin to the short arm of chromosome 17 (logarithm of odds [LOD] = 3.30; 54 cM between D17S1294 and D17S1299, P = 1.0 x 10(-4)). The strongest evidence for linkage over all pedigrees for fasting glucose was also observed in this region (LOD = 1.44; 58 cM, P = 9.9 x 10(-3)). The results of this study provide impetus for future positional cloning of QTLs regulating insulin and glucose levels. Identifying genes in these regions should provide insight into the nature of genetic factors regulating plasma glucose and insulin concentrations.

The transcription factor TFIID consists of TATA-binding protein (TBP) and TBP-associated factors (TAFs). TAFs are essential for modulation of transcriptional activity but the regulation of TAFs is complex and many important aspects remain unclear. In this study, we have identified and characterized five novel truncated forms of the TFIID subunit TAF4 (TAF(II)135). Analysis of the mouse gene structure revealed that all truncations were the results of alternative splicing and resulted in the loss of domains or parts of domains implicated in TAF4 functional interactions. Results from transcriptional assays showed that several of the TAF4 isoforms exerted dominant negative effects on TAF4 activity in nuclear receptor-mediated transcriptional activation. In addition, alternative TAF4 isoforms could be detected in specific cell types. Our results indicate an additional level of complexity in TAF4-mediated regulation of transcription and suggest context-specific roles for these new TAF4 isoforms in transcriptional regulation in vivo.

Cyclic nucleotide-gated (CNG) channels encoded by tax-4 and tax-2 genes are required for chemo- and thermo-sensation in Caenorhabditis elegans. Here we report the identification and the characterization of cng-3, a new CNG channel gene, found in C. elegans. CNG-3 contains six putative transmembrane regions and a cyclic nucleotide-binding domain that show high homology with CNG channels of higher animals as well as TAX-4. The expression of cng-3 is detected from early stages in worm development and restricted in five sensory neurons of amphid including AFD neuron. While a cng-3 null mutant displays normal chemotaxis to volatile odorants, the mutant worms exhibit impaired thermal tolerance. These results indicate that CNG-3, a new member of CNG channel subunits, may play a critical role in sensation or response of thermal stress in C. elegans.

Salmonella enterica serovars often have a broad host range, and some cause both gastrointestinal and systemic disease. But the serovars Paratyphi A and Typhi are restricted to humans and cause only systemic disease. It has been estimated that Typhi arose in the last few thousand years. The sequence and microarray analysis of the Paratyphi A genome indicates that it is similar to the Typhi genome but suggests that it has a more recent evolutionary origin. Both genomes have independently accumulated many pseudogenes among their approximately 4,400 protein coding sequences: 173 in Paratyphi A and approximately 210 in Typhi. The recent convergence of these two similar genomes on a similar phenotype is subtly reflected in their genotypes: only 30 genes are degraded in both serovars. Nevertheless, these 30 genes include three known to be important in gastroenteritis, which does not occur in these serovars, and four for Salmonella-translocated effectors, which are normally secreted into host cells to subvert host functions. Loss of function also occurs by mutation in different genes in the same pathway (e.g., in chemotaxis and in the production of fimbriae).