bio-ucsc-api 0.2.0 → 0.2.1

data/lib/bio-ucsc/susscr2/mrna.rb

@@ -0,0 +1,81 @@
+# Copyright::
+#  Copyright (C) 2011 MISHIMA, Hiroyuki <missy at be.to / hmishima at nagasaki-u.ac.jp> 
+# License::     The Ruby licence (Ryby's / GPLv2 dual)
+#
+# In the hg18 database, this table is actually separated
+# into "chr1_*", "chr2_*", etc. This class dynamically
+# define *::Chr1_*, *::Chr2_*, etc. The
+# Rmsk.find_by_interval calls an appropreate class automatically.
+
+module Bio
+  module Ucsc
+    module SusScr2
+
+      class Mrna
+        KLASS = "Mrna"
+        KLASS_S = "mrna"
+
+        Bio::Ucsc::SusScr2::CHROMS.each do |chr|
+          class_eval %!
+            class #{chr[0..0].upcase + chr[1..-1]}_#{KLASS} < DBConnection
+              set_table_name "#{chr[0..0].downcase + chr[1..-1]}_#{KLASS_S}"
+              set_primary_key nil
+              set_inheritance_column nil
+
+              def self.find_by_interval(interval, opt = {:partial => true})
+                find_first_or_all_by_interval(interval, :first, opt)
+              end
+        
+              def self.find_all_by_interval(interval, opt = {:partial => true})
+                find_first_or_all_by_interval(interval, :all, opt)
+              end
+
+              def self.find_first_or_all_by_interval(interval, first_all, opt)
+                zstart = interval.zero_start
+                zend   = interval.zero_end
+                if opt[:partial] == true
+                  where = <<-SQL
+      tName = :chrom
+AND   bin in (:bins)
+AND ((tStart BETWEEN :zstart AND :zend)
+ OR  (tEnd BETWEEN :zstart AND :zend)
+ OR  (tStart <= :zstart AND tEnd >= :zend))
+                  SQL
+                else
+                  where = <<-SQL
+      tName = :chrom
+AND   bin in (:bins)
+AND ((tStart BETWEEN :zstart AND :zend)
+AND  (tEnd BETWEEN :zstart AND :zend))
+                  SQL
+                end
+                cond = {
+                  :chrom  => interval.chrom,
+                  :bins   => Bio::Ucsc::UcscBin.bin_all(zstart, zend),
+                  :zstart => zstart,
+                  :zend   => zend,
+                }
+                self.find(first_all,
+                          { :select => "*",
+                            :conditions => [where, cond], })
+              end
+            end
+          !
+        end # each chromosome
+
+        def self.find_by_interval(interval, opt = {:partial => true})
+          chrom = interval.chrom[0..0].upcase + interval.chrom[1..-1]
+          chr_klass = self.const_get("#{chrom}_#{KLASS}")
+          chr_klass.__send__(:find_by_interval, interval, opt)
+        end
+
+        def self.find_all_by_interval(interval, opt = {:partial => true})
+          chrom = interval.chrom[0..0].upcase + interval.chrom[1..-1]
+          chr_klass = self.const_get("#{chrom}_#{KLASS}")
+          chr_klass.__send__(:find_all_by_interval, interval, opt)
+        end
+      end # class
+
+    end # module Hg18 
+  end # module Ucsc
+end # module Bio

data/lib/bio-ucsc/table_class_detector.rb

@@ -115,9 +115,8 @@ AND  (tEnd BETWEEN :zstart AND :zend))
                 :zend => zend,
               }
                 self.find(first_all,
-                          :select => "*",
-                          :conditions => [where, cond],
-                          )
+                          { :select => "*",
+                            :conditions => [where, cond], })
               end
             end
           !
@@ -158,8 +157,8 @@ AND  (tEnd BETWEEN :zstart AND :zend))
                   :zend => zend,
                 }
                 self.find(first_all,
-                          :select => "*",
-                          :conditions => [where, cond],)
+                          { :select => "*",
+                            :conditions => [where, cond], })
               end
             end
           !
@@ -211,9 +210,8 @@ AND  (chromEnd BETWEEN :zstart AND :zend))
                 :zend => zend,
               }          
                 self.find(first_all,
-                          :select => "*",
-                          :conditions => [where, cond],
-                          )
+                          { :select => "*",
+                            :conditions => [where, cond], })
               end
             end
           !
@@ -254,8 +252,8 @@ AND  (chromEnd BETWEEN :zstart AND :zend))
                   :zend => zend,
                 }
                 self.find(first_all,
-                          :select => "*",
-                          :conditions => [where, cond],)
+                          { :select => "*",
+                            :conditions => [where, cond], })
               end
             end
           !
@@ -306,9 +304,8 @@ AND  (txEnd BETWEEN :zstart AND :zend))
                 :zend   => zend,
               }
                 self.find(first_all,
-                          :select => "*",
-                          :conditions => [where, cond],
-                          )
+                          { :select => "*",
+                            :conditions => [where, cond], })
               end
             end
           !
@@ -349,8 +346,8 @@ AND ((txStart BETWEEN :zstart AND :zend)
                   :zend => zend,
                 }
                 self.find(first_all,
-                          :select => "*",
-                          :conditions => [where, cond],)
+                          { :select => "*",
+                            :conditions => [where, cond], })
               end
             end
           !
@@ -401,9 +398,8 @@ AND  (genoEnd BETWEEN :zstart AND :zend))
                 :zend   => zend,
               }
                 self.find(first_all,
-                          :select => "*",
-                          :conditions => [where, cond],
-                          )
+                          { :select => "*",
+                            :conditions => [where, cond], })
               end
             end
           !
@@ -444,8 +440,8 @@ AND ((genoStart BETWEEN :zstart AND :zend)
                   :zend => zend,
                 }
                 self.find(first_all,
-                          :select => "*",
-                          :conditions => [where, cond],)
+                          { :select => "*",
+                            :conditions => [where, cond], })
               end
             end
           !
@@ -465,14 +461,14 @@ AND ((genoStart BETWEEN :zstart AND :zend)
       end
 
       def uphead(sym)
-        (sym.to_s[0..0].upcase + sym[1..-1])
+        (sym.to_s[0..0].upcase + sym.to_s[1..-1])
       end
 
       def downhead(sym)
         if UPPERCASED_TABLE_PREFIX.any?{|x|sym.to_s.start_with? x}
           sym.to_s
         else
-          (sym.to_s[0..0].downcase + sym[1..-1])
+          (sym.to_s[0..0].downcase + sym.to_s[1..-1])
         end
       end
 

data/lib/bio-ucsc/taegut1.rb

@@ -0,0 +1,33 @@
+#
+# = AUTOMATIC Table Definition of the Zebra finch Jul. 2008 (WUGSC 3.2.4/taeGut1) assembly
+# Copyright::   Copyright (C) 2011
+#               MISHIMA, Hiroyuki
+#               <missy at be.to / hmishima at nagasaki-u.ac.jp> 
+# License::     The Ruby licence (Ryby's / GPLv2 dual)
+#
+
+require "#{File.dirname(__FILE__)}/taegut1/db_connection"
+require "#{File.dirname(__FILE__)}/table_class_detector"
+
+module Bio
+  module Ucsc
+    module TaeGut1
+      CHROMS = 
+%w(chr1 chr1A chr1B chr2 chr3 chr4 chr4A chr5 chr6 chr7 chr8 chr9 chr10
+chr11 chr12 chr13 chr14 chr15 chr16 chr17 chr18 chr19 chr20
+chr21 chr22 chr23 chr24 chr25 chr26 chr27 chr28 chrZ chrM)
+
+      extend TableClassDetector
+
+      base = "#{File.dirname(__FILE__)}/taegut1"
+      autoload :ChainGalGal3, "#{base}/chaingalgal3"
+      autoload :ChainGalGal3Link, "#{base}/chaingalgal3link"
+      autoload :Est, "#{base}/est"
+      autoload :Gap, "#{base}/gap"
+      autoload :Gold, "#{base}/gold"
+      autoload :IntronEst, "#{base}/intronest"
+      autoload :Mrna, "#{base}/mrna"
+      autoload :Rmsk, "#{base}/rmsk"
+    end
+  end
+end

data/lib/bio-ucsc/taegut1/chaingalgal3.rb

@@ -0,0 +1,81 @@
+# Copyright::
+#  Copyright (C) 2011 MISHIMA, Hiroyuki <missy at be.to / hmishima at nagasaki-u.ac.jp> 
+# License::     The Ruby licence (Ryby's / GPLv2 dual)
+#
+# In the hg18 database, this table is actually separated
+# into "chr1_*", "chr2_*", etc. This class dynamically
+# define *::Chr1_*, *::Chr2_*, etc. The
+# Rmsk.find_by_interval calls an appropreate class automatically.
+
+module Bio
+  module Ucsc
+    module TaeGut1
+
+      class ChainGalGal3
+        KLASS = "ChainGalGal3"
+        KLASS_S = KLASS[0..0].downcase + KLASS[1..-1]
+
+        Bio::Ucsc::GalGal3::CHROMS.each do |chr|
+          class_eval %!
+            class #{chr[0..0].upcase + chr[1..-1]}_#{KLASS} < DBConnection
+              set_table_name "#{chr[0..0].downcase + chr[1..-1]}_#{KLASS_S}"
+              set_primary_key nil
+              set_inheritance_column nil
+
+              def self.find_by_interval(interval, opt = {:partial => true})
+                find_first_or_all_by_interval(interval, :first, opt)
+              end
+        
+              def self.find_all_by_interval(interval, opt = {:partial => true})
+                find_first_or_all_by_interval(interval, :all, opt)
+              end
+
+              def self.find_first_or_all_by_interval(interval, first_all, opt)
+                zstart = interval.zero_start
+                zend   = interval.zero_end
+                if opt[:partial] == true
+                  where = <<-SQL
+      tName = :chrom
+AND   bin in (:bins)
+AND ((tStart BETWEEN :zstart AND :zend)
+ OR  (tEnd BETWEEN :zstart AND :zend)
+ OR  (tStart <= :zstart AND tEnd >= :zend))
+                  SQL
+                else
+                  where = <<-SQL
+      tName = :chrom
+AND   bin in (:bins)
+AND ((tStart BETWEEN :zstart AND :zend)
+AND  (tEnd BETWEEN :zstart AND :zend))
+                  SQL
+                end
+                cond = {
+                  :chrom  => interval.chrom,
+                  :bins   => Bio::Ucsc::UcscBin.bin_all(zstart, zend),
+                  :zstart => zstart,
+                  :zend   => zend,
+                }
+                self.find(first_all,
+                          { :select => "*",
+                            :conditions => [where, cond], })
+              end
+            end
+          !
+        end # each chromosome
+
+        def self.find_by_interval(interval, opt = {:partial => true})
+          chrom = interval.chrom[0..0].upcase + interval.chrom[1..-1]
+          chr_klass = self.const_get("#{chrom}_#{KLASS}")
+          chr_klass.__send__(:find_by_interval, interval, opt)
+        end
+
+        def self.find_all_by_interval(interval, opt = {:partial => true})
+          chrom = interval.chrom[0..0].upcase + interval.chrom[1..-1]
+          chr_klass = self.const_get("#{chrom}_#{KLASS}")
+          chr_klass.__send__(:find_all_by_interval, interval, opt)
+        end
+      end # class
+
+    end # module Hg18 
+  end # module Ucsc
+end # module Bio

data/lib/bio-ucsc/taegut1/chaingalgal3link.rb

@@ -0,0 +1,81 @@
+# Copyright::
+#  Copyright (C) 2011 MISHIMA, Hiroyuki <missy at be.to / hmishima at nagasaki-u.ac.jp> 
+# License::     The Ruby licence (Ryby's / GPLv2 dual)
+#
+# In the hg18 database, this table is actually separated
+# into "chr1_*", "chr2_*", etc. This class dynamically
+# define *::Chr1_*, *::Chr2_*, etc. The
+# Rmsk.find_by_interval calls an appropreate class automatically.
+
+module Bio
+  module Ucsc
+    module TaeGut1
+
+      class ChainGalGal3Link
+        KLASS = "ChainGalGal3Link"
+        KLASS_S = KLASS[0..0].downcase + KLASS[1..-1]
+
+        Bio::Ucsc::GalGal3::CHROMS.each do |chr|
+          class_eval %!
+            class #{chr[0..0].upcase + chr[1..-1]}_#{KLASS} < DBConnection
+              set_table_name "#{chr[0..0].downcase + chr[1..-1]}_#{KLASS_S}"
+              set_primary_key nil
+              set_inheritance_column nil
+
+              def self.find_by_interval(interval, opt = {:partial => true})
+                find_first_or_all_by_interval(interval, :first, opt)
+              end
+        
+              def self.find_all_by_interval(interval, opt = {:partial => true})
+                find_first_or_all_by_interval(interval, :all, opt)
+              end
+
+              def self.find_first_or_all_by_interval(interval, first_all, opt)
+                zstart = interval.zero_start
+                zend   = interval.zero_end
+                if opt[:partial] == true
+                  where = <<-SQL
+      tName = :chrom
+AND   bin in (:bins)
+AND ((tStart BETWEEN :zstart AND :zend)
+ OR  (tEnd BETWEEN :zstart AND :zend)
+ OR  (tStart <= :zstart AND tEnd >= :zend))
+                  SQL
+                else
+                  where = <<-SQL
+      tName = :chrom
+AND   bin in (:bins)
+AND ((tStart BETWEEN :zstart AND :zend)
+AND  (tEnd BETWEEN :zstart AND :zend))
+                  SQL
+                end
+                cond = {
+                  :chrom  => interval.chrom,
+                  :bins   => Bio::Ucsc::UcscBin.bin_all(zstart, zend),
+                  :zstart => zstart,
+                  :zend   => zend,
+                }
+                self.find(first_all,
+                          { :select => "*",
+                            :conditions => [where, cond], })
+              end
+            end
+          !
+        end # each chromosome
+
+        def self.find_by_interval(interval, opt = {:partial => true})
+          chrom = interval.chrom[0..0].upcase + interval.chrom[1..-1]
+          chr_klass = self.const_get("#{chrom}_#{KLASS}")
+          chr_klass.__send__(:find_by_interval, interval, opt)
+        end
+
+        def self.find_all_by_interval(interval, opt = {:partial => true})
+chrom = interval.chrom[0..0].upcase + interval.chrom[1..-1]
+          chr_klass = self.const_get("#{chrom}_#{KLASS}")
+          chr_klass.__send__(:find_all_by_interval, interval, opt)
+        end
+      end # class
+
+    end # module Hg18 
+  end # module Ucsc
+end # module Bio

data/lib/bio-ucsc/taegut1/db_connection.rb

@@ -0,0 +1,57 @@
+#
+# = db_connection.rb - UCSC DB connection
+#
+# Copyright::   Cioyrught (C) 2011
+#               MISHIMA, Hiroyuki <missy at be.to / hmishima at nagasaki-u.ac.jp> 
+# Copyright::   Copyright (C) 2010
+#               Jan Aerts <jan.aerts@gmail.com>
+# License::     Ruby licence (Ryby's / GPLv2 dual)
+
+require 'rubygems'
+require 'active_record'
+
+module Bio
+  module Ucsc
+
+    module TaeGut1
+      DATABASE_NAME = "taeGut1"
+
+      # = DESCRIPTION
+      # Actual connection established to the UCSC mysql server.
+      class DBConnection < ActiveRecord::Base
+        # = DESCRIPTION
+        # The #connect method make the connection
+        #
+        # = USAGE
+        #  Bin::Ucsc::DB_NAME::DBConnection.connect
+        #
+        # ---
+        # *Arguments*: none
+        @@db_adapter  ||= 'mysql'
+        @@db_host     ||= 'genome-mysql.cse.ucsc.edu'
+        @@db_username ||= 'genome'
+        @@db_password ||= ''
+
+        cattr_accessor :db_adapter, :db_host, :db_username, :db_password
+
+        self.abstract_class = true
+
+        def self.default
+          @@db_adapter  = 'mysql'
+          @@db_host     = 'genome-mysql.cse.ucsc.edu'
+          @@db_username = 'genome'
+          @@db_password = ''
+        end
+
+        def self.connect
+          establish_connection({ :adapter =>  @@db_adapter,
+                                 :host =>     @@db_host,
+                                 :database => DATABASE_NAME,
+                                 :username => @@db_username,
+                                 :password => @@db_password, })
+        end
+      end
+    end # module
+    
+  end # module Ucsc
+end # module Bio

data/lib/bio-ucsc/taegut1/est.rb

@@ -0,0 +1,81 @@
+# Copyright::
+#  Copyright (C) 2011 MISHIMA, Hiroyuki <missy at be.to / hmishima at nagasaki-u.ac.jp> 
+# License::     The Ruby licence (Ryby's / GPLv2 dual)
+#
+# In the hg18 database, this table is actually separated
+# into "chr1_*", "chr2_*", etc. This class dynamically
+# define *::Chr1_*, *::Chr2_*, etc. The
+# Rmsk.find_by_interval calls an appropreate class automatically.
+
+module Bio
+  module Ucsc
+    module TaeGut1
+
+      class Est
+        KLASS = "Est"
+        KLASS_S = "est"
+
+        Bio::Ucsc::TaeGut1::CHROMS.each do |chr|
+          class_eval %!
+            class #{chr[0..0].upcase + chr[1..-1]}_#{KLASS} < DBConnection
+              set_table_name "#{chr[0..0].downcase + chr[1..-1]}_#{KLASS_S}"
+              set_primary_key nil
+              set_inheritance_column nil
+
+              def self.find_by_interval(interval, opt = {:partial => true})
+                find_first_or_all_by_interval(interval, :first, opt)
+              end
+        
+              def self.find_all_by_interval(interval, opt = {:partial => true})
+                find_first_or_all_by_interval(interval, :all, opt)
+              end
+
+              def self.find_first_or_all_by_interval(interval, first_all, opt)
+                zstart = interval.zero_start
+                zend   = interval.zero_end
+                if opt[:partial] == true
+                  where = <<-SQL
+      tName = :chrom
+AND   bin in (:bins)
+AND ((tStart BETWEEN :zstart AND :zend)
+ OR  (tEnd BETWEEN :zstart AND :zend)
+ OR  (tStart <= :zstart AND tEnd >= :zend))
+                  SQL
+                else
+                  where = <<-SQL
+      tName = :chrom
+AND   bin in (:bins)
+AND ((tStart BETWEEN :zstart AND :zend)
+AND  (tEnd BETWEEN :zstart AND :zend))
+                  SQL
+                end
+                cond = {
+                  :chrom  => interval.chrom,
+                  :bins   => Bio::Ucsc::UcscBin.bin_all(zstart, zend),
+                  :zstart => zstart,
+                  :zend   => zend,
+                }
+                self.find(first_all,
+                          { :select => "*",
+                            :conditions => [where, cond], })
+              end
+            end
+          !
+        end # each chromosome
+
+        def self.find_by_interval(interval, opt = {:partial => true})
+          chrom = interval.chrom[0..0].upcase + interval.chrom[1..-1]
+          chr_klass = self.const_get("#{chrom}_#{KLASS}")
+          chr_klass.__send__(:find_by_interval, interval, opt)
+        end
+
+        def self.find_all_by_interval(interval, opt = {:partial => true})
+          chrom = interval.chrom[0..0].upcase + interval.chrom[1..-1]
+          chr_klass = self.const_get("#{chrom}_#{KLASS}")
+          chr_klass.__send__(:find_all_by_interval, interval, opt)
+        end
+      end # class
+
+    end # module Hg18 
+  end # module Ucsc
+end # module Bio