Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Genome alignment”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,423 records · Page 79Linked to original sources

Structural characterization of a rice actin gene.

We have isolated and sequenced a full-length cDNA clone containing information for the rice actin gene RAc1. Transcript terminus mapping and sequence alignment between the RAc1 cDNA clone and a previously isolated RAc1 genomic clone were used to determine the structure of the RAc1 gene. This allowed us to make the first complete structural characterization of a plant actin gene. The analysis revealed the presence of a 5'-noncoding exon, separated by an intron, from the first translated exon of the RAc1 gene. This is one of the few reported cases of a plant gene containing such a 5'-noncoding exon. Sequence comparison between the previously isolated plant actin genes suggests that such an exon may be a common feature of plant actin gene structure. The present study also confirms that the rice actin gene family is composed of at least eight unique members.

Actins↗

Physical mapping of the E/C and grc regions of the rat major histocompatibility complex.

Alignment of class I-hybridizing cosmids from an R21 (AlBlDlEugrc+) genomic DNA library gave two contigs: one [150 kilobases (kb)] encompassed the E/C region, or a large part thereof, and the other (110 kb) contained the grc region which has genes influencing resistance to chemical carcinogens (rcc), fertility (ft), and growth (dw-3). Amplification of gene sequences in the four cosmids in the E/C region using Eu-specific and LW2 (RT1.C)-specific primers showed that each cosmid contained both Eu-like and C-like genes. They are clearly different but closely associated, and they show some variation from the prototypic E (Eu) and C (LW2) genes, respectively. Comparison of DNA from grc+ and grc- strains of rats showed that the deletion in the grc- strains was approximately 50 kb, and that it was located on two of the three cosmids in the grc-region contig. The use of specific class I probes showed that the grc region contained tandemly duplicated RT1.O-RT1.N genes and that the RT.BM1 loci lay outside of the grc region. Neither contig reacted with probes specific for class II, TNFA, Hsp70, or RT1.M genes. The data presented here and the previous data in the literature (summarized in Gill et al. 1995) suggest that the gene order in the major histocompatibility complex (MHC) and MHC-linked region of the rat is: A-E/C-grc-M.

Animals↗

Characterisation and transferability of apple SSRs to two European pear F1 populations.

European pear (Pyrus communis L.) is among the important fruit species for which only few genetic studies have been carried out. Available evidence indicates that simple sequence repeats (SSR) are very useful as molecular markers because they are codominant, highly polymorphic, abundant and reproducible. The present paper reports more than 100 apple SSR markers in two populations of European pear; a total of 41 SSR markers were then positioned on a genetic linkage map of the cross 'Passe Crassane' x 'Harrow Sweet' and 31 in the map 'Abbe Fetel' x 'Max Red Bartlett'. Syntenic relationships between pear and apple maps have been considered for the chromosomes carrying two or more SSR markers. The alignment among the two maps supports the colinearity of the two genomes with respect both to identification and to orientation of the linkage groups.

Chromosome Mapping↗

Comparative mapping on the mouse X chromosome defines a myotubular myopathy equivalent region.

The gene for X-linked myotubular myopathy (MTM1) has been localized to a 300-kb critical region in human Xq28 between IDS and GABRA3. As part of an effort to clone this gene, we developed a YAC contig on the mouse X Chromosome (Chr) which includes loci homologous to those within the human MTM1 critical region. The murine contig consists of 18 YACs and spans 2.5-3.0 Mb. We have aligned the human and murine physical maps by isolating conserved mouse genomic fragments, including CpG islands and trapped exons. We believe that the simultaneous isolation of genes from both mouse and human and continued comparative mapping will prove helpful in the eventual identification of MTM1 and other genes in the region.

Animals↗

Cloning and expression of a novel lipase gene from Bacillus sphaericus 205y.

A Bacillus sphaericus strain (205y) that produces an organic solvent-tolerant lipase was isolated in Port Dickson, Malaysia. The gene for the lipase was recovered from a genomic library and sequenced. Phylogenetic analysis was performed based on an alignment of thirteen microbial lipase sequences obtained from the NCBI database. The analysis suggested that the B. sphaericus lipase gene is a novel gene, as it is distinct from other lipase genes in Families I.4 and I.5 reported so far. Expression in Escherichia coli under the control of the lacZ promoter resulted in an eight-fold increase in enzyme activity after a 3-h induction with 1 mM IPTG. The crude enzyme thus obtained showed a slight (10%) enhancement in activity after a 30-min incubation in 25% (v/v) n-hexane at 37 degrees C, and retained 90% of its activity after a similar period in 25% (v/v) p-xylene.

Amino Acid Sequence↗

Genetic characterization of a new mammalian reovirus, type 2 Winnipeg (T2W).

We previously described isolation of a potentially new reovirus strain from the central nervous system of an 8-week-old female infant with a history of active varicella, oral thrush, hypoalbuminemia, intermittent fevers, diarrhea and feeding intolerance [Hermann et al., Ped. Inf. Dis J. 23, 373 (2004)]. This reovirus strain was tentatively identified as a member of the serotype 2 group by virus neutralization and RNA-gel electrophoresis studies and has been named type 2 Winnipeg (T2W). For this study we determined the nucleotide sequences of the T2W S1, S2, S3 and S4 genome segments to allow molecular comparison with other reoviruses. Comparative segment alignments of T2W S1 gene sequence with other reovirus S1 sequences showed T2W belongs to reovirus serotype 2. T2W S1 is most similar to the S1 genes of reovirus strains T2/Human/Netherlands/1,984 and T2/Human/Netherlands/1,973 with nucleotide identity >93%. The T2W S2 gene showed highest identity to reovirus T1 Lang S2 (approximately 75%). The T2W S3 gene showed highest identity to the S3 gene of T3/Human/Netherlands/1,983 (approximately 74%), and the T2W S4 gene showed highest identity to the T2 Jones S4 gene (approximately 73%). Pairwise protein comparisons between T2W sigma proteins and all available reovirus sigma proteins ranged from <21% identity for the sigma1 comparisons to more than 95% identity for sigma2 comparisons. The predicted T2W sigma1, sigma2 and sigma3 protein sequences were confirmed by mass spectrometry.

Amino Acid Sequence↗

Sequence and linkage analysis of the Coxiella burnetii citrate synthase-encoding gene.

The nucleotide (nt) sequence of the Coxiella burnetii citrate synthase-encoding gene (gltA), previously cloned in Escherichia coli, was determined. The nt sequence analysis revealed an open reading frame (ORF) of 1290 bp capable of coding for a protein of 430 amino acids (aa) with a deduced Mr of 48,633. Preceding an ATG start codon, a possible transcription start point (tsp) with homology to the E. coli promoter consensus was detected. A poly-purine-rich region occurred immediately upstream from the gltA reading frame and potentially serves as a ribosome-binding site. Additionally, a G + C-rich region of dyad symmetry 3' to the translational stop codon was found that could possibly function as a Rho-independent transcriptional termination signal. A large, nearly perfect, inverted repeat was identified upstream from the gltA tsp and was shown by Southern analysis to be present in multiple copies in the C. burnetii genome. The deduced aa sequence of C. burnetii GltA was optimally aligned with enzymes from various prokaryotic sources and one eukaryotic source (pig heart). Using perfect aa identity, the C. burnetii enzyme demonstrated the greatest homology with GltA from Acinetobacter anitratum (65%). Although only 26% aa identity was seen with the pig heart enzyme, many of the residues identified in ligand binding appear to be conserved. Sequencing studies of a region centered approx. 5.6 kb upstream from gltA revealed an ORF read with opposite polarity that encodes a peptide highly homologous to the C terminus of the flavoprotein subunit of E. coli succinate dehydrogenase. This report represents the first nt sequence analysis of a gene of known function from the obligate intracellular parasite, C. burnetii.

Amino Acid Sequence↗

An improved reverse transcription-polymerase chain reaction (RT-PCR) assay for the detection of two cherry flexiviruses in Prunus spp.

A one-step reverse transcription-polymerase chain reaction (RT-PCR) procedure was developed to detect Cherry green ring mottle virus (CGRMV) and Cherry necrotic rusty mottle virus (CNRMV) in woody indicators and naturally infected Prunus spp. Viral RNA suitable for RT-PCR was obtained by a simple trapping method that did not require either extraction of double-stranded RNA (dsRNA) or total RNA, availability of virus antibodies, or purification of viral particles. Consensus primers, degenerate primers and virus-specific primers, whose designs were based on alignments of available cherry flexivirus sequences, were tested to amplify viral genomic fragments of six CGRMV isolates and one CNRMV isolate. RT-PCR allowed CGRMV detection in total RNA and viral RNA preparations equivalent to 400mug and 4mug of infected leaf tissue, respectively. CGRMV was detected in tender shoots, leaves, bark and root tips, and the strongest bands were obtained using young leaves. Detection was less consistent in summer when the temperature was elevated and plant tissues were old. A direct comparison of the RT-PCR and grafting assays indicated that the RT-PCR assay is sensitive, rapid and reliable. The method will improve the routine diagnosis of cherry flexiviruses in Prunus spp.

Capsid Proteins↗

Molecular cloning and characterization of a cDNA encoding a ryegrass (Lolium perenne) ENOD40 homologue.

The ENOD40 gene, found in several leguminous and non-leguminous plant species is expressed in the pericycle of legume roots early in the nodulation process, adjacent to the protoxylem poles, but before the cortical cells divide to form the nodule itself. The ENOD40 transcript contains only short open reading frames which give rise to short peptides with a signaling function. The current work reports the cloning and analysis of ENOD40 genes from perennial ryegrass (Lolium perenne) and barley (Hordeum vulgare). Specific polymerase chain reaction (PCR) techniques lead to the isolation of a 659 bp cDNA encoding an ENOD40 homologue, designated LpE-NOD40, from a Lolium perenne stem cDNA library. In addition, a partial ENOD40 cDNA of 384 bp was isolated from barley (Hordeum vulgare) by RT-PCR cloning. The LpENOD40 transcript encodes a putative dodecapeptide, similar to that identified in ENOD40s from leguminous plants and other dicots, and also to ENOD40s from monocots. The coding sequences of ryegrass and barley ENOD40 are represented by a short open reading frame of 12 amino acids. These show a high degree of similarity to each other and to other ENOD40 sequences from monocots. The corresponding genomic DNA from a genomic ryegrass lambda library revealed that the LpENOD40 gene contains no introns. Southern blot analysis shows that the ryegrass genome contains a single copy, possibly two copies, of the gene. Alignment of the ENOD40 cDNA sequences from ryegrass and barley revealed high (77%) nucleotide homology. The ENOD40 peptides are highly conserved, not only among monocots but also on comparison with the dicot peptides. The amino acid identity of region I from ryegrass with its counterparts in maize, barley and rice is 92, 83 and 75% respectively. Expression analysis by RT-PCR demonstrates that a high level of LpENOD40 gene transcript was found expressed in stem tissue, while a lower level was detected in leaves and only a very low expression in flowers of perennial ryegrass. No LpENOD40 transcript was detected in roots.

Amino Acid Sequence↗

Structure and organization of the C4 genes.

This 200 000 Mr serum protein is coded for by at least two separate loci, C4A and C4B, which map in the HLA Class III region on chromosome 6 in man. Both loci are highly polymorphic with more than 30 alleles, including null alleles assigned to the two loci. The complete nucleotide sequence of a full length C4A cDNA clone and a substantial part of a C4b cDNA clone has shown class differences which can be used to synthesize nucleotide probes specific for C4A and C4B. Three C4 loci of approximately 16 kilobases each spaced by 10 kilobases have been identified in DNA from one individual and aligned 30 kilobases from the factor B gene by overlapping cloned genomic fragments from a cosmid library. Characterization of these genes by restriction mapping, nucleotide sequence analysis and hybridization with C4A and C4B specific synthetic oligonucleotides show that these genes are very similar.

Alleles↗

Genomic sequence analysis of Fugu rubripes CFTR and flanking genes in a 60 kb region conserving synteny with 800 kb of human chromosome 7.

To define control elements that regulate tissue-specific expression of the cystic fibrosis transmembrane regulator (CFTR), we have sequenced 60 kb of genomic DNA from the puffer fish Fugu rubripes (Fugu) that includes the CFTR gene. This region of the Fugu genome shows conservation of synteny with 800-kb sequence of the human genome encompassing the WNT2, CFTR, Z43555, and CBP90 genes. Additionally, the genomic structure of each gene is conserved. In a multiple sequence alignment of human, mouse, and Fugu, the putative WNT2 promoter sequence is shown to contain highly conserved elements that may be transcription factor or other regulatory binding sites. We have found two putative ankyrin repeat-containing genes that flank the CFTR gene. Overall sequence analysis suggests conservation of intron/exon boundaries between Fugu and human CFTR and revealed extensive homology between functional protein domains. However, the immediate 5' regions of human and Fugu CFTR are highly divergent with few conserved sequences apart from those resembling diminished cAMP response elements (CRE) and CAAT box elements. Interestingly, the polymorphic polyT tract located upstream of exon 9 is present in human and Fugu but absent in mouse. Similarly, an intron 1 and intron 9 element common to human and Fugu is absent in mouse. The euryhaline killifish CFTR coding sequence is highly homologous to the Fugu sequence, suggesting that upregulation of CFTR in that species in response to salinity may be regulated transcriptionally.

Amino Acid Sequence↗

Evolution of gene families: the multidrug resistance transporter genes in five related yeast species.

The available genomic sequences of five closely related hemiascomycetous yeast species (Kluyveromyces lactis, Kluyveromyces waltii, Candida glabrata, Ashbya (Eremothecium) gossypii with Saccharomyces cerevisiae as a reference) were analysed to identify multidrug resistance (MDR) transport proteins belonging to the ATP-binding cassette (ABC) and major facilitator superfamilies (MFS), respectively. The phylogenetic trees clearly demonstrate that a similar set of gene (sub)families already existed in the common ancestor of all five fungal species studied. However, striking differences exist between the two superfamilies with respect to the evolution of the various subfamilies. Within the ABC superfamily all six half-size transporters with six transmembrane-spanning domains (TMs) and most full-size transporters with 12 TMs have one and only one gene per genome. An exception is the PDR family, in which gene duplications and deletions have occurred independently in individual genomes. Among the MFS transporters, the DHA2 family (TC 2.A.1.3) is more variable between species than the DHA1 family (TC 2.A.1.2). Conserved gene order relationships allow to trace the evolution of most (sub)families, for which the Kluyveromyces lactis genome can serve as an optimal scaffold. Cross-species sequence alignment of orthologous upstream gene sequences led to the identification of conserved sequence motifs ("phylogenetic footprints"). Almost half of them match known sequence motifs for the MDR regulators described in S. cerevisiae. The biological significance of those and of the novel predicted motifs awaits to be confirmed experimentally.

ATP-Binding Cassette Transporters↗

Transcription initiation sites and nucleotide sequence of a herpes simplex virus 1 gene conserved in the Epstein-Barr virus genome and reported to affect the transport of viral glycoproteins.

Earlier reports have localized mutations which affect the processing and transport of herpes simplex virus 1 glycoproteins to a region located between the genes specifying glycoprotein B and the major viral DNA-binding protein (beta 8). The nucleotide sequence of this region contains a single long open reading frame encoding a 780-amino-acid protein with a predicted molecular weight of 83,845. To confirm the existence of this protein, rabbit polyclonal antibody was made against a synthetic peptide made according to the predicted sequence of a hydrophilic domain near the carboxy terminal of the protein. This antibody reacted with an infected cell protein of an apparent molecular weight of 95,500. We designated this protein infected cell protein 18.5 (ICP18.5). S1 nuclease analysis suggested that the 5.6-kilobase mRNA encoding ICP18.5 is initiated predominantly from one site, but three weaker initiation sites also seemed to occur within a 74-base-pair stretch of DNA. This gene appears to be conserved in the Epstein-Barr virus (EBV) genome, inasmuch as 174 of the 780 amino acids of ICP18.5 align with corresponding amino acids predicted by the EBV open reading frame BALF3. The EBV gene is located adjacent to the gene specifying a homolog of the herpes simplex virus 1 glycoprotein B.

Amino Acid Sequence↗

Cloning and identification of a microRNA cluster within the latency-associated region of Kaposi's sarcoma-associated herpesvirus.

MicroRNAs (miRNAs) are small, noncoding regulatory RNA molecules that bind to 3' untranslated regions (UTRs) of mRNAs to either prevent their translation or induce their degradation. Previously identified in a variety of organisms ranging from plants to mammals, miRNAs are also now known to be produced by viruses. The human gammaherpesvirus Epstein-Barr virus has been shown to encode miRNAs, which potentially regulate both viral and cellular genes. To determine whether Kaposi's sarcoma-associated herpesvirus (KSHV) encodes miRNAs, we cloned small RNAs from KSHV-positive primary effusion lymphoma-derived cells and endothelial cells. Sequence analysis revealed 11 isolated RNAs of 19 to 23 bases in length that perfectly align with KSHV. Surprisingly, all candidate miRNAs mapped to a single genomic locale within the latency-associated region of KSHV. These data suggest that viral and host cellular gene expression may be regulated by miRNAs during both latent and lytic KSHV replication.

Cloning, Molecular↗

Site-directed mutagenesis of the cytoplasmic domains of the human beta 2-adrenergic receptor. Localization of regions involved in G protein-receptor coupling.

Numerous plasma membrane-bound receptors are coupled to various effectors via a family of guanine nucleotide regulatory proteins (G proteins). Amino acid sequences of these receptors, deduced from cDNA and genomic clones, indicate the presence of seven transmembrane-spanning domains. Alignment of the available amino acid sequences of these G protein-linked receptors reveals striking homologies in regions predicted to lie near the cytoplasmic surface of the cell membrane. As these areas are likely those which interact with G proteins, we reasoned that systematic introduction of non-native sequence into these highly conserved regions of the human beta 2-adrenergic receptor would allow resolution of loci participating directly in receptor-G protein coupling. Based on this strategy, we constructed 19 mutant receptor species comprising substitutions and deletions of native sequence in the putative cytoplasmic domains of human beta 2-adrenergic receptor. By monitoring ligand binding characteristics and receptor-mediated stimulation of adenylyl cyclase, we have determined that the C-terminal portion of the third cytoplasmic loop and the N-terminal segment of the cytoplasmic tail appear to be critical for productive receptor-coupling to G proteins. In addition, we have implicated two other areas of the receptor that possibly play supportive roles in maintaining proper orientation of the G protein binding site. These comprise the second cytoplasmic loop and a conserved cysteine residue in the cytoplasmic tail.

Amino Acid Sequence↗

Alterations in organization and transcription of the mitochondrial genome of cytoplasmic male sterile sugar beet (Beta vulgaris L.).

We have constructed a physical map of the mitochondrial DNA of a cytoplasmic male sterile (CMS) sugar beet line, TK81-MS, and compared it with that published for normal fertile sugar beet (cv. TK81-O) to clarify the differences between the CMS and normal mitochondrial genomes. The TK81-MS genome is present as a single circular molecule of 481.8 kb, or as two molecules of 184.9 and 296.9 kb. The CMS genome was found to be highly rearranged relative to the normal mitochondrial genome, with at least fifteen rearrangement and/or inversion events being required to align the two DNAs. Analysis of transcription patterns of known mitochondrial genes and rearranged regions revealed six genes, coxI, coxII, atpA, atp6, rps3, and orf324, whose expression is altered in the CMS line relative to the normal line. Of these six, only the coxI transcript pattern differs between male-sterile and fertility-restored genotypes, making it likely that the coxI locus is involved in mediating CMS in sugar beet.

Base Sequence↗

Metagenomic Study of the MESA: Detection of Gemella Morbillorum and Association With Coronary Heart Disease.

BACKGROUND: Inflammation is a feature of coronary heart disease (CHD), but the role of proinflammatory microbial infection in CHD remains understudied. METHODS AND RESULTS: CHD was defined in the MESA (Multi-Ethnic Study of Atherosclerosis) as myocardial infarction (251 participants), resuscitated arrest (2 participants), and CHD death (80 participants). We analyzed sequencing reads from 4421 MESA participants in the National Heart, Lung, and Blood Institute Trans-Omics for Precision Medicine program using the PathSeq workflow of the Genome Analysis Tool Kit and a 65-gigabase microbial reference. Paired reads aligning to 840 microbes were detected in >1% of participants. The association of the presence of microbe reads with incident CHD (follow-up, ~18&#x2009;years) was examined. First, important variables were ascertained using a single regularized Cox proportional hazard model, examining change of risk as a function of presence of microbe with age, sex, education level, Life's Simple 7, and inflammation. For variables of importance, the hazard ratio (HR) was estimated in separate (unregularized) Cox proportional hazard models including the same covariates (significance threshold Bonferroni corrected P<6&#xd7;10-5, 0.05/840). Reads from 2 microbes were significantly associated with CHD: Gemella morbillorum (HR, 3.14 [95% CI, 1.92-5.12]; P=4.86&#xd7;10-6) and Pseudomonas species NFACC19-2 (HR, 3.22 [95% CI, 2.03-5.41]; P=1.58&#xd7;10-6). CONCLUSIONS: Metagenomics of whole-genome sequence reads opens a possible frontier for detection of pathogens for chronic diseases. The association of G morbillorum and Pseudomonas species reads with CHD raises the possibilities that microbes may drive atherosclerotic inflammation and that treatments for specific pathogens may provide clinical utility for CHD reduction.

Humans↗

Correcting errors in shotgun sequences.

Sequencing errors in combination with repeated regions cause major problems in shotgun sequencing, mainly due to the failure of assembly programs to distinguish single base differences between repeat copies from erroneous base calls. In this paper, a new strategy designed to correct errors in shotgun sequence data using defined nucleotide positions, DNPs, is presented. The method distinguishes single base differences from sequencing errors by analyzing multiple alignments consisting of a read and all its overlaps with other reads. The construction of multiple alignments is performed using a novel pattern matching algorithm, which takes advantage of the symmetry between indices that can be computed for similar words of the same length. This allows for rapid construction of multiple alignments, with no previous pair-wise matching of sequence reads required. Results from a C++ implementation of this method show that up to 99% of sequencing errors can be corrected, while up to 87% of the single base differences remain and up to 80% of the corrected reads contain at most one error. The results also show that the method outperforms the error correction method used in the EULER assembler. The prototype software, MisEd, is freely available from the authors for academic use.

Algorithms↗