phykit 2.1.70__tar.gz → 2.1.71__tar.gz

{phykit-2.1.70 → phykit-2.1.71}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: phykit
-Version: 2.1.70
+Version: 2.1.71
 Home-page: https://github.com/jlsteenwyk/phykit
 Author: Jacob L. Steenwyk
 Author-email: jlsteenwyk@gmail.com

{phykit-2.1.70 → phykit-2.1.71}/phykit/cli_registry.py

@@ -55,6 +55,8 @@ ALIAS_TO_HANDLER: Dict[str, str] = {
     "sum_of_pairs_score": "sum_of_pairs_score",
     "sops": "sum_of_pairs_score",
     "sop": "sum_of_pairs_score",
+    "phylo_gwas": "phylo_gwas",
+    "pgwas": "phylo_gwas",
     "vs": "variable_sites",
     # Tree aliases
     "bss": "bipartition_support_stats",

{phykit-2.1.70 → phykit-2.1.71}/phykit/phykit.py

@@ -167,6 +167,9 @@ class Phykit:
                     - calculate sum-of-pairs score between a reference and query alignment
                 thread_dna (alias: pal2nal; p2n)
                     - thread dna sequences over a protein alignment
+                phylo_gwas (alias: pgwas)
+                    - phylogenetic genome-wide association study
+                      (Pease et al. 2016 approach)
                 variable_sites (alias: vs)
                     - calculates the number and percentage of variable sites
                       in an alignment
@@ -1878,6 +1881,115 @@ class Phykit:
         _add_json_argument(parser)
         _run_service(parser, argv, VariableSites)
 
+    @staticmethod
+    def phylo_gwas(argv):
+        parser = _new_parser(
+            description=textwrap.dedent(
+                f"""\
+                {help_header}
+
+                Phylogenetic genome-wide association study following the
+                Pease et al. (2016) approach.
+
+                Performs per-site association tests between alignment
+                columns and a phenotype (categorical or continuous),
+                applies Benjamini-Hochberg FDR correction, optionally
+                classifies significant associations as monophyletic or
+                polyphyletic using a phylogenetic tree, and produces a
+                Manhattan plot.
+
+                Categorical phenotypes use Fisher's exact test (2 groups)
+                or chi-squared test (>2 groups). Continuous phenotypes use
+                point-biserial correlation. Only biallelic sites are tested.
+
+                Aliases:
+                  phylo_gwas, pgwas
+                Command line interfaces:
+                  pk_phylo_gwas, pk_pgwas
+
+                Usage:
+                phykit phylo_gwas -a <alignment> -d <phenotype> -o <output>
+                  [-t <tree>] [-p <partition>] [--alpha 0.05]
+                  [--exclude-monophyletic] [--csv <file>] [--json]
+                  [shared plot options]
+
+                Options
+                =====================================================
+                -a/--alignment              FASTA alignment file
+
+                -d/--phenotype              two-column TSV file:
+                                            taxon<tab>phenotype
+
+                -o/--output                 output Manhattan plot path
+
+                -t/--tree                   optional Newick tree for
+                                            monophyletic/polyphyletic
+                                            classification
+
+                -p/--partition              optional RAxML-style partition
+                                            file for gene annotations
+
+                --alpha                     FDR significance threshold
+                                            (default: 0.05)
+
+                --exclude-monophyletic      exclude monophyletic
+                                            associations from results
+
+                --dot-size                  scale factor for dot size
+                                            in the Manhattan plot
+                                            (default: 1.0; use 2.0
+                                            for double, 0.5 for half)
+
+                --csv                       output per-site results as CSV
+                                            to the specified file
+
+                --fig-width                 figure width in inches
+                                            (auto-scaled if omitted)
+
+                --fig-height                figure height in inches
+                                            (auto-scaled if omitted)
+
+                --dpi                       resolution in DPI
+                                            (default: 300)
+
+                --no-title                  hide the plot title
+
+                --title                     custom title text
+
+                --legend-position           legend location (e.g.,
+                                            "upper right", "none")
+
+                --ylabel-fontsize           font size for y-axis labels;
+                                            0 to hide
+
+                --xlabel-fontsize           font size for x-axis labels;
+                                            0 to hide
+
+                --title-fontsize            font size for the title
+
+                --axis-fontsize             font size for axis labels
+
+                --colors                    comma-separated colors
+                                            (hex or named)
+
+                --json                      optional argument to output
+                                            results as JSON
+                """
+            ),
+        )
+        parser.add_argument("-a", "--alignment", type=str, required=True, help=SUPPRESS, metavar="")
+        parser.add_argument("-d", "--phenotype", type=str, required=True, help=SUPPRESS, metavar="")
+        parser.add_argument("-o", "--output", type=str, required=True, help=SUPPRESS, metavar="")
+        parser.add_argument("-t", "--tree", type=str, default=None, help=SUPPRESS, metavar="")
+        parser.add_argument("-p", "--partition", type=str, default=None, help=SUPPRESS, metavar="")
+        parser.add_argument("--alpha", type=float, default=0.05, help=SUPPRESS, metavar="")
+        parser.add_argument("--exclude-monophyletic", action="store_true", default=False, help=SUPPRESS)
+        parser.add_argument("--dot-size", type=float, default=1.0, help=SUPPRESS, metavar="")
+        parser.add_argument("--csv", type=str, default=None, help=SUPPRESS, metavar="")
+        add_plot_arguments(parser)
+        _add_json_argument(parser)
+        _run_service(parser, argv, PhyloGwas)
+
     @staticmethod
     def alignment_subsample(argv):
         parser = _new_parser(
@@ -8163,6 +8275,10 @@ def dstatistic(argv=None):
     Phykit.dstatistic(sys.argv[1:])
 
 
+def phylo_gwas(argv=None):
+    Phykit.phylo_gwas(sys.argv[1:])
+
+
 def dfoil(argv=None):
     Phykit.dfoil(sys.argv[1:])
 

{phykit-2.1.70 → phykit-2.1.71}/phykit/service_factories.py

@@ -44,6 +44,7 @@ RelativeCompositionVariability = _LazyServiceFactory("phykit.services.alignment.
 RelativeCompositionVariabilityTaxon = _LazyServiceFactory("phykit.services.alignment.rcvt", "RelativeCompositionVariabilityTaxon")
 RenameFastaEntries = _LazyServiceFactory("phykit.services.alignment.rename_fasta_entries", "RenameFastaEntries")
 SumOfPairsScore = _LazyServiceFactory("phykit.services.alignment.sum_of_pairs_score", "SumOfPairsScore")
+PhyloGwas = _LazyServiceFactory("phykit.services.alignment.phylo_gwas", "PhyloGwas")
 VariableSites = _LazyServiceFactory("phykit.services.alignment.variable_sites", "VariableSites")
 
 AncestralReconstruction = _LazyServiceFactory("phykit.services.tree.ancestral_reconstruction", "AncestralReconstruction")

{phykit-2.1.70 → phykit-2.1.71}/phykit/services/alignment/__init__.py

@@ -26,6 +26,7 @@ _EXPORTS = {
     "RelativeCompositionVariabilityTaxon": "rcvt",
     "RenameFastaEntries": "rename_fasta_entries",
     "SumOfPairsScore": "sum_of_pairs_score",
+    "PhyloGwas": "phylo_gwas",
     "VariableSites": "variable_sites",
 }