WO2015073707A1 - Methods of treating heart failure with agonists of hypocretin receptor 2 - Google Patents

Abstract

Methods of treating individuals for heart failure are disclosed. In particular, the invention relates to methods of treating individuals for heart failure with agonists of hypocretin receptor 2. The invention further relates to the use of genetic analysis of the single nucleotide polymorphism rs7767652, which resides immediately upstream of the HCRTR2 gene, to identify individuals in need of treatment for heart failure who are predicted to be more responsive to medical intervention.

Description

METHODS OF TREATING HEART FAILURE WITH AGONISTS OF

HYPOCRETIN RECEPTOR 2

TECHNICAL FIELD

The present invention pertains to methods of treating heart failure. In particular, the invention relates to methods of treating individuals for heart failure with agonists of hypocretin receptor 2. The invention further relates to the use of genetic analysis of the single nucleotide polymorphism rs7767652, which resides immediately upstream of the HCRTR2 gene, to identify individuals in need of treatment for heart failure who are predicted to be more responsive to medical intervention. BACKGROUND

Heart failure, a syndrome characterized by impaired cardiac function, affects over five million people in the United States alone. As the most common admitting diagnosis for patients over 65, heart failure portends a worse prognosis than most cancers (Lloyd- Jones et al. (2010) Circulation 121 :e46-e215). With the aging population, it is expected to affect 3.5% of the US population over the next 20 years (Heidenreich et al. (2012) Circulation 123:933-944).

Heart failure is the result of a wide variety of underlying conditions, most commonly coronary artery disease and hypertension (Ho et al. (1993) J. Am. Coll. Cardiol. 22:6A-13A), and can be associated with preserved or reduced left ventricular ejection fraction. Some patients in the low ejection fraction group exhibit a dramatic improvement with optimal medical therapy, which includes treatment with β- adrenergic receptor antagonists and inhibitors of the renin-angiotensin-aldosterone system. Why some patients have a dramatic response to pharmaceutical therapy whereas others decompensate and require transplantation is not well understood.

Candidate gene studies with targeted genotyping of common variants in

ADRB1 and GRK5 members of the β-adrenergic receptor signaling pathway have shown associations with long-term heart-failure survival (Cresci et al. (2009) J. Am. Coll. Cardiol. 54:432-444; Liggett et al. (2008) Nature Medicine 14:510-517). Studies with more comprehensive genotyping of cardiovascular gene candidates identified associations between common variants in HSPB7 and FRMD4B and dilated cardiomyopathy (Stark et al. (2010) PLoS Genetics 6 :e 1001167) or advanced heart failure (Cappola et al. (2010) Circ. Cardiovasc. Genet. 3: 147-154). In the limited number of genome-wide studies of heart failure performed (Larson et al. (2007) BMC Medical Genetics 8 Suppl. 1 :S5; Smith et al. (2010) Circ. Cardiovasc. Genet. 3:256- 266; Villard et al. (2011) Eur. Heart J. 32: 1065-1076), the common genomic variants that have reached genome-wide significance have not been replicated in large cohorts. Genome-wide studies of improvement in left ventricular function in response to medical intervention have yet to be performed.

There remains a need for an improved therapy for patients suffering from heart failure and diagnostic biomarkers that can be used to identify patients likely to have a favorable response to medical intervention. SUMMARY OF THE INVENTION

The present invention relates to methods of treating individuals for heart failure with agonists of hypocretin receptor 2. The invention further relates to the use of genetic analysis of the single nucleotide polymorphism rs7767652, which resides immediately upstream of the HCRTR2 gene, to identify individuals in need of treatment for heart failure who are predicted to be more responsive to medical intervention.

In one aspect, the invention includes a method for treating an individual for heart failure, the method comprising administering a therapeutically effective amount of a composition comprising at least one agonist of hypocretin receptor 2 to the individual. In one embodiment, the individual is human. The individual may also be diabetic or hypertensive.

Agonists of hypocretin receptor 2 that can be used in the practice of the invention include, but are not limited to orexin-A, orexin-B, [Alal l,D-Leul5]-orexin B, and SB-668875. In one embodiment the agonist of hypocretin receptor 2 is an orexin mimic. The agonist of hypocretin receptor 2 can be recombinantly or synthetically produced. Compositions comprising one or more agonists of hypocretin receptor 2 may be administered by any suitable method, including, but not limited to, intravenously, intra-arterially, subcutaneously, mtracerebroventricularly, intrathecally, or intranasally. In one embodiment, compositions are administered into the brain or spinal cord of the individual. In one embodiment, compositions are administered locally into the cisterna magna or the rostral ventrolateral medulla.

A therapeutically effective amount of an agonist of hypocretin receptor 2 can be administered to an individual in one or more administrations, applications or dosages. By "therapeutically effective dose or amount" of an agonist of hypocretin receptor 2 is intended an amount that, when administered, as described herein, brings about a positive therapeutic response, such as improved left ventricular function or survival of an individual treated for heart failure.

In certain embodiments, the agonist of hypocretin receptor 2 is administered according to a daily dosing regimen, such as a twice-a-week or three-times-a-week dosing regimen. Multiple cycles of the method of treatment may be administered to an individual for a time period sufficient to effect at least a partial recovery of heart function, such as a 10% or greater increase or improvement in left ventricular function (e.g., as measured by ejection fraction). Arterial blood pressure or heart rate may also be increased in the individual after administering the agonist of hypocretin receptor 2.

In certain embodiments, the method further comprises treatment with one or more other agents for treating heart failure such as, but not limited to an angiotensin converting enzyme (ACE) inhibitor, an angiotensin II receptor blocker (ARB), a beta- blocker, digoxin, a diuretic, a blood vessel dilator, potassium, magnesium, an aldactone inhibitor, a calcium channel blocker, and an inotrope.

In another aspect, the invention includes a method for enhancing

cardiovascular function of an individual with heart failure, the method comprising administering an effective amount of an agonist of hypocretin receptor 2 to the individual. Enhancing cardiovascular function may include, for example, a 10% or greater increase or improvement in left ventricular function (e.g., as measured by ejection fraction), and/or increased arterial blood pressure, and/or increased heart rate.

In another aspect, the invention includes a composition comprising an agonist of hypocretin receptor 2 for use in the treatment of heart failure. In certain

embodiments, the agonist of hypocretin receptor 2 is selected from the group consisting of orexin-A, orexin-B, [Alal l,D-Leul5]-orexin B, and SB-668875.

In another aspect, the invention includes a composition comprising an agonist of hypocretin receptor 2 for use in enhancing cardiovascular function. In certain embodiments, the agonist of hypocretin receptor 2 is selected from the group consisting of orexin-A, orexin-B, [Alal l,D-Leul5]-orexin B, and SB-668875.

In another aspect, the invention includes the use of an agonist of hypocretin receptor 2 in the manufacture of a medicament for treating an individual for heart failure. In certain embodiments, the agonist of hypocretin receptor 2 is selected from the group consisting of orexin-A, orexin-B, [Alal l ,D-Leul5]-orexin B, and SB- 668875.

In another aspect the invention includes a diagnostic method for predicting therapeutic responsiveness to treatment with a beta blocker or an ACE inhibitor for an individual with heart failure, the method comprising determining which allele is present at single nucleotide polymorphism rs7767652 in the individual, wherein the presence of the major allele at the single nucleotide polymorphism rs7767652 indicates that the individual will be more responsive to treatment for heart failure with a beta blocker or an ACE inhibitor than an individual with the minor allele at the single nucleotide polymorphism rs7767652.

These and other embodiments of the subject invention will readily occur to those of skill in the art in view of the disclosure herein.

BRIEF DESCRIPTION OF THE DRAWINGS FIG. 1 shows the primary data sources and bioinformatic analyses used for

SNV selection for the custom array. Approximately one third of the 1536 SNVs were derived from the gene expression-weighted analysis of an extreme response genome- wide association study (GWAS), one third from the most statistically significant intergenic SNPs from the extreme response GWAS, and one third from the GWAS- informed analyses of gene-coexpression networks, semantic literature networks or a curated set of SNPs manually selected from the literature.

FIGS. 2A and 2B show a Manhattan plot showing the significance of association with a 10% improvement in ejection fraction for SNPs in the customized heart failure array analysis of Phase II individuals and a regional LD plot of the statistically significant locus. FIG 2A shows SNPs plotted on the x axis according to their position on each chromosome against association with greater than 10% improvement in ejection fraction on the y axis (shown as -logio-P). The gray horizontal line represents the threshold for Bonferroni-adjusted significance (p < 3.6x l0^"5). The locus near HCRTR2 reached an association significance of P = 3.29^x 10^~5 in the Phase II analysis. FIG 2B shows an LD map (r²) of the region surrounding rs7767652 (boxed) drawn according to Phase II HapMap (release 26) CEU individuals. All SNPs are shaded according to the pairwise r value.

FIGS. 3A-3H show the functional roles of the rs7767652 allele and HCRTR2 assessed with allele-specific gene expression, measurement of human HCRTR2 gene expression in diseased tissue, measurement of ventricular function, exercise capacity and trichrome staining in HCRTR2 deficient mice and infusion of orexin A in chemically stressed mice. FIG. 3A shows a diagram of the luciferase reporter with the rs7767652 locus, which contains a TCF-4 binding site. This SNP with its flanking sequence was inserted into the multiple cloning site (MCS) of a pLuc-MCS vector, which is driven by a minimal promoter upstream from the luciferase gene. FIG. 3B shows allele-specific luciferase reporter assays at rs7767652 in C2C12 cells. The reporter construct carrying the minor allele (rs776 T-Luc) had slightly lower activity at baseline compared to the construct with the major allele (rs776 C-Luc). The major allele was highly responsive to -catenin/TCF4 overexpression, whereas the minor allele disrupted -catenin/TCF4 mediated transactivation. The TCF4 transactivation mutant (TCF4A30) blocked β-catenin effects. FIG. 3C shows human HCRTR2 gene expression measurements from microarray experiments in myocardial tissue of patients with dilated cardiomyopathy (DCM, n = 162), ischemic cardiomyopathy (ICM, n =143) and controls (n = 45) obtained from the Gene Expression Omnibus (GEO) repository. Gene expression values were graphed as mean fluorescence intensity ratios with 95% confidence intervals normalized across experiments.

HCRTR2 expression is higher in DCM tissue compared to controls (* p = 0.00017) and is higher in ICM tissue compared to controls (** p = 0.000044). FIG. 3D shows diastolic function in C57BL/6 wild-type (n =11) and HCRTR2 transcription deficient (TD) mice (n =10). The HCRTR2 TD mice had worse diastolic function as measured by a higher E/E' ratio compared to the wild-type mice (p = 0.05). FIG. 3E shows systolic function in C57BL/6 wild-type (n =11) and HCRTR2 transcription deficient (TD) mice (n = 10). There were no significant differences in ejection fraction between the two sets of mice (p = 0.94). FIG. 3F shows changes in exercise capacity before and after treatment with angiotensin II + isoproterenol chemical stress in C57BL/6 wild-type (n = 8) compared to HCRTR2 TD (n = 8) mice. The HCRTR2 TD mice had a smaller improvement in peak V0₂ (13.0 ml/kg/min +/- 3.0 vs. 4.2 ml/kg/min +/- 2.2, p =0.032) after 2 weeks of chemical stress. FIG. 3G shows representative examples of trichrome-stained myocardial sections from five C57BL/6 wild type (left) and five HCRTR2 TD (right) mice after 2 weeks of angiotensin II + isoproterenol chemical stress. The HCRTR2 TD mouse hearts had greater degree of trichrome stain (blue/red ratio 0.82 +/- 0.1 vs. 0.85 +/-0.01, p = 0.022). FIG. 3H shows systolic function in WT mice infused with either saline (n = 6) or orexin A (n = 7) for 4 weeks and stressed with angiotensin II + isoproterenol infusion during the latter 2 weeks. Mice that underwent infusion with orexin A had better systolic function after chemical stress at week 4 compared to those infused with saline (Ejection Fraction 56.6% +/- 4.5% vs. 43.6% +/- 3.3%, * p = 0.045).

FIG. 4 shows the coexpression module most significantly associated with differential gene expression in heart failure tissue. SLC25A30 is the hub of the 18- member module most significantly associated with differential gene expression in heart failure tissue compared to normal myocardial tissue.

FIG. 5 shows a quantile-quantile plot for Phase II analysis. Observed P-values of SNPs against the theoretical distribution of expected p-values were plotted from the Phase II analysis. A total of 1366 SNPs from the Phase II customized array were used to generate a plot.

FIG. 6 shows a Manhattan plot showing the significance of association with a

20% improvement in ejection fraction for SNPs in Illumina 550K assay of Phase I individuals. SNPs are plotted on the x axis according to their position on each chromosome against association with greater than 20% improvement in ejection fraction on the y axis (shown as -logio-P).

FIG. 7 shows a quantile-quantile plot for Phase I analysis. Observed P-values of SNPs against the theoretical distribution of expected /^-values were plotted from the Phase I analysis. A total of 525,220 SNPs from the Phase I Illumina 550K platform were used to generate a plot.

DETAILED DESCRIPTION OF THE INVENTION

The practice of the present invention will employ, unless otherwise indicated, conventional methods of genetics, chemistry, biochemistry, molecular biology and recombinant DNA techniques, within the skill of the art. Such techniques are explained fully in the literature. See, e.g., Single Nucleotide Polymorphisms: Methods and Protocols (Methods in Molecular Biology, A.A. Komar ed., Humana Press; 2^nd edition, 2009); Genetic Variation: Methods and Protocols (Methods in Molecular Biology, M.R. Barnes and G. Breen eds., Humana Press, 2010); Handbook of

Experimental Immunology, Vols. I-IV (D.M. Weir and C.C. Blackwell eds.,

Blackwell Scientific Publications); A.L. Lehninger, Biochemistry (Worth Publishers, Inc., current addition); Sambrook, et al, Molecular Cloning: A Laboratory Manual (2nd Edition, 1989); Methods In Enzymology (S. Colo wick and N. Kaplan eds., Academic Press, Inc.).

All publications, patents and patent applications cited herein, whether supra or infra, are hereby incorporated by reference in their entireties.

I. DEFINITIONS

In describing the present invention, the following terms will be employed, and are intended to be defined as indicated below.

It must be noted that, as used in this specification and the appended claims, the singular forms "a", "an" and "the" include plural referents unless the content clearly dictates otherwise. Thus, for example, reference to "an agonist" includes a mixture of two or more such agonists, and the like.

By "therapeutically effective dose or amount" of an agonist of hypocretin receptor 2 is intended an amount that, when administered as described herein, brings about a positive therapeutic response, such as improved left ventricular function or survival of a subject treated for heart failure.

The term "partial response" (PR) as used herein means a 10% or greater increase or improvement in left ventricular function (e.g., as measured by ejection fraction).

The term "survival" as used herein means the time from the first dose of an agonist of hypocretin receptor 2 to the time of death.

The terms "hybridize" and "hybridization" refer to the formation of complexes between nucleotide sequences which are sufficiently complementary to form complexes via Watson Crick base pairing. The terms are intended to refer to the formation of a specific hybrid between a probe and a target region. The term "derived from" is used herein to identify the original source of a molecule but is not meant to limit the method by which the molecule is made which can be, for example, by chemical synthesis or recombinant means.

"Recombinant" as used herein to describe a nucleic acid molecule means a polynucleotide of genomic, cDNA, viral, semisynthetic, or synthetic origin which, by virtue of its origin or manipulation is not associated with all or a portion of the polynucleotide with which it is associated in nature. The term "recombinant" as used with respect to a protein or polypeptide means a polypeptide produced by expression of a recombinant polynucleotide. In general, the gene of interest is cloned and then expressed in transformed organisms, as described further below. The host organism expresses the foreign gene to produce the protein under expression conditions.

"Substantially purified" generally refers to isolation of a substance

(compound, polynucleotide, oligonucleotide, protein, or polypeptide) such that the substance comprises the majority percent of the sample in which it resides. Typically in a sample, a substantially purified component comprises 50%, preferably 80%-85%, more preferably 90-95% of the sample. Techniques for purifying polynucleotides oligonucleotides and polypeptides of interest are well-known in the art and include, for example, ion-exchange chromatography, affinity chromatography and

sedimentation according to density.

By "isolated" is meant, when referring to a polypeptide, that the indicated molecule is separate and discrete from the whole organism with which the molecule is found in nature or is present in the substantial absence of other biological macro molecules of the same type. The term "isolated" with respect to a polynucleotide or oligonucleotide is a nucleic acid molecule devoid, in whole or part, of sequences normally associated with it in nature; or a sequence, as it exists in nature, but having heterologous sequences in association therewith; or a molecule disassociated from the chromosome.

The terms "polynucleotide," "oligonucleotide," "nucleic acid" and "nucleic acid molecule" are used herein to include a polymeric form of nucleotides of any length, either ribonucleotides or deoxyribonucleotides. This term refers only to the primary structure of the molecule. Thus, the term includes triple-, double- and single-stranded DNA, as well as triple-, double- and single-stranded RNA. It also includes modifications, such as by methylation and/or by capping, and unmodified forms of the polynucleotide. More particularly, the terms "polynucleotide,"

"oligonucleotide," "nucleic acid" and "nucleic acid molecule" include

polydeoxyribonucleotides (containing 2-deoxy-D-ribose), polyribonucleotides (containing D-ribose), any other type of polynucleotide which is an N- or C-glycoside of a purine or pyrimidine base, and other polymers containing nonnucleotidic backbones, for example, polyamide (e.g., peptide nucleic acids (PNAs)) and polymorpholino (commercially available from the Anti-Virals, Inc., Corvallis, Oregon, as Neugene) polymers, and other synthetic sequence-specific nucleic acid polymers providing that the polymers contain nucleobases in a configuration which allows for base pairing and base stacking, such as is found in DNA and RNA. There is no intended distinction in length between the terms "polynucleotide,"

"oligonucleotide," "nucleic acid" and "nucleic acid molecule," and these terms will be used interchangeably. Thus, these terms include, for example, 3'-deoxy-2',5'-DNA, oligodeoxyribonucleotide N3' P5' phosphoramidates, 2'-0-alkyl-substituted RNA, double- and single-stranded DNA, as well as double- and single-stranded RNA,

DNA: RNA hybrids, and hybrids between PNAs and DNA or RNA, and also include known types of modifications, for example, labels which are known in the art, methylation, "caps," substitution of one or more of the naturally occurring nucleotides with an analog, internucleotide modifications such as, for example, those with uncharged linkages (e.g., methyl phosphonates, phosphotriesters, phosphoramidates, carbamates, etc.), with negatively charged linkages (e.g., phosphorothioates, phosphorodithioates, etc.), and with positively charged linkages (e.g.,

aminoalklyphosphoramidates, aminoalkylphosphotriesters), those containing pendant moieties, such as, for example, proteins (including nucleases, toxins, antibodies, signal peptides, poly-L-lysine, etc.), those with intercalators (e.g., acridine, psoralen, etc.), those containing chelators (e.g., metals, radioactive metals, boron, oxidative metals, etc.), those containing alkylators, those with modified linkages (e.g., alpha anomeric nucleic acids, etc.), as well as unmodified forms of the polynucleotide or oligonucleotide. The term also includes locked nucleic acids (e.g., comprising a ribonucleotide that has a methylene bridge between the 2'-oxygen atom and the 4'- carbon atom). See, for example, Kurreck et al. (2002) Nucleic Acids Res. 30: 1911- 1918; Elayadi et al. (2001) Curr. Opinion Invest. Drugs 2: 558-561; Orum et al. (2001) Curr. Opinion Mol. Ther. 3: 239-243; Koshkin et al. (1998) Tetrahedron 54: 3607-3630; Obika et al. (1998) Tetrahedron Lett. 39: 5401-5404.

The terms "polymorphism," "polymorphic nucleotide," "polymorphic site" or "polymorphic nucleotide position" refer to a position in a nucleic acid that possesses the quality or character of occurring in several different forms. A nucleic acid may be naturally or non-naturally polymorphic, e.g., having one or more sequence differences (e.g., additions, deletions and/or substitutions) as compared to a reference sequence. A reference sequence may be based on publicly available information (e.g., the U.C. Santa Cruz Human Genome Browser Gateway (genome.ucsc.edu/cgi-bin/hgGateway) or the NCBI website (ncbi.nlm.nih.gov)) or may be determined by a practitioner of the present invention using methods well known in the art (e.g., by sequencing a reference nucleic acid). A nucleic acid polymorphism is characterized by two or more "alleles," or versions of the nucleic acid sequence. Typically, an allele of a polymorphism that is identical to a reference sequence is referred to as a "reference allele" and an allele of a polymorphism that is different from a reference sequence is referred to as an "alternate allele," or sometimes a "variant allele". As used herein, the term "major allele" refers to the more frequently occurring allele at a given polymorphic site, and "minor allele" refers to the less frequently occurring allele, as present in the general or study population.

The term "single nucleotide polymorphism" or "SNP" refers to a polymorphic site occupied by a single nucleotide, which is the site of variation between allelic sequences. The site is usually preceded by and followed by highly conserved sequences of the allele (e.g., sequences that vary in less than 1/100 or 1/1000 members of the populations). A single nucleotide polymorphism usually arises due to substitution of one nucleotide for another at the polymorphic site. Single nucleotide polymorphisms can also arise from a deletion of a nucleotide or an insertion of a nucleotide relative to a reference allele.

The HUGO Gene Nomenclature Committee (HGNC) has assigned unique gene symbols and names to more than 31,000 human loci, of which over 19,000 are protein coding. Genenames.org is a curated online repository of HGNC-approved gene nomenclature, and those gene names are used herein.

SNPs generally are described as having a minor allele frequency, which can vary between populations, but generally refers to the sequence variation (A,T,G, or C) that is less common than the major allele. The frequency can be obtained from dbSNP or other sources, or may be determined for a certain population using Hardy- Weinberg equilibrium (See for details see Eberle MA, Rieder MJ, Kruglyak L, Nickerson DA (2006) Allele Frequency Matching Between SNPs Reveals an Excess of Linkage Disequilibrium in Genie Regions of the Human Genome. PLoS Genet 2(9): el42. doi: 10.1371/journal.pgen.0020142; herein incorporated by reference).

SNPs described here are further detailed at the NCBI dbSNP, found online at ncbi.nlm.nih.gov/projects/SNP, from which the following brief descriptions are taken. Population data on these SNPs may be found at this site. The SNP representations are according to NCBI dbSNP convention, with the two alleles between brackets.

The term "rs7767652," as used herein, refers to the SNP at position 27 (indicated by the [C/T]) of the following sequence upstream of the HCRTR2 gene: GTTGATAAAATATTCCACATTATAAT[C/T]AAAGAAAGTACAAAAAAGTTA AATG (SEQ ID NO: l). This SNP has global minor allele frequency

(MAF)/MinorAlleleCount of T = 0.1837/400, as reported in the NCBI dbSNP.

As used herein, the term "probe" or "oligonucleotide probe" refers to a polynucleotide, as defined above, that contains a nucleic acid sequence

complementary to a nucleic acid sequence present in the target nucleic acid analyte (e.g., at SNP location). The polynucleotide regions of probes may be composed of DNA, and/or RNA, and/or synthetic nucleotide analogs. Probes may be labeled in order to detect the target sequence. Such a label may be present at the 5' end, at the 3' end, at both the 5' and 3' ends, and/or internally.

An "allele-specific probe" hybridizes to only one of the possible alleles of a SNP under suitably stringent hybridization conditions.

The term "primer" or "oligonucleotide primer" as used herein, refers to an oligonucleotide that hybridizes to the template strand of a nucleic acid and initiates synthesis of a nucleic acid strand complementary to the template strand when placed under conditions in which synthesis of a primer extension product is induced, i.e., in the presence of nucleotides and a polymerization-inducing agent such as a DNA or RNA polymerase and at suitable temperature, pH, metal concentration, and salt concentration. The primer is preferably single-stranded for maximum efficiency in amplification, but may alternatively be double-stranded. If double-stranded, the primer can first be treated to separate its strands before being used to prepare extension products. This denaturation step is typically effected by heat, but may alternatively be carried out using alkali, followed by neutralization. Thus, a "primer" is complementary to a template, and complexes by hydrogen bonding or hybridization with the template to give a primer/template complex for initiation of synthesis by a polymerase, which is extended by the addition of covalently bonded bases linked at its 3' end complementary to the template in the process of DNA or R A synthesis.

Typically, nucleic acids are amplified using at least one set of oligonucleotide primers comprising at least one forward primer and at least one reverse primer capable of hybridizing to regions of a nucleic acid flanking the portion of the nucleic acid to be amplified.

An "allele-specific primer" matches the sequence exactly of only one of the possible alleles of a SNP, hybridizes at the SNP location, and amplifies only one specific allele if it is present in a nucleic acid amplification reaction.

The term "amplicon" refers to the amplified nucleic acid product of a PCR reaction or other nucleic acid amplification process (e.g., ligase chain reaction (LGR), nucleic acid sequence based amplification (NASBA), transcription-mediated amplification (TMA), Q-beta amplification, strand displacement amplification, or target mediated amplification). Amplicons may comprise RNA or DNA depending on the technique used for amplification. For example, DNA amplicons may be generated by RT-PCR, whereas RNA amplicons may be generated by TMA/NASBA.

As used herein, the terms "label" and "detectable label" refer to a molecule capable of detection, including, but not limited to, radioactive isotopes, fluorescers, chemiluminescers, chromophores, enzymes, enzyme substrates, enzyme co factors, enzyme inhibitors, semiconductor nanoparticles, dyes, metal ions, metal sols, ligands (e.g., biotin, streptavidin or haptens) and the like. The term "fluorescer" refers to a substance or a portion thereof which is capable of exhibiting fluorescence in the detectable range. Particular examples of labels which may be used in the practice of the invention include, but are not limited to, horseradish peroxidase (HRP), SYBR® green, SYBR® gold, fluorescein, carboxyfluorescein (FAM), Alexa Fluor dyes, Cy3, Cy5, Cy7, 7-amino-4-methylcoumarin-3-acetic acid (AMCA), 5-(and-6)-carboxy-X- rhodamine, lissamine rhodamine, fluorescein-5-isothiocyanate (FITC), 7- diethylaminocoumarin-3-carboxylic acid, tetramethylrhodamine-5-(and-6)- isothiocyanate, 5-(and-6)-carboxytetramethylrhodamine, 7-hydroxycoumarin-3- carboxylic acid, 6- [fluorescein 5-(and-6)-carboxamido]hexanoic acid, N-(4,4- difluoro-5,7-dimethyl-4-bora-3a,4a diaza-3-indacenepropionic acid, eosin-5- isothiocyanate; erythrosine-5-isothiocyanate, 5-(and-6)-carboxyrhodamine 6G, CASCADE blue aectylazide, CAL Fluor Orange 560, CAL Fluor Red 610, Quasar Blue 670, tetramethyl rhodamine (TAMRA), 2', 4', 5', T- tetrachloro-4-7- dichloro fluorescein (TET), rhodamine, dansyl, umbelliferone, dimethyl acridinium ester (DMAE), Texas red, Pacific Blue, Pacific Orange, quantum dots, luminol, NADPH, and α-β-galactosidase.

The term "biological sample," as used herein, includes any cell, tissue, or bodily fluid comprising nucleic acids (e.g., genomic DNA or RNA) such as, but not limited to, blood, saliva, cells from buccal swabbing, skin, hair, biopsies of organs, amniotic fluid, various other tissues, and the like.

The terms "subject," "individual," and "patient," are used interchangeably herein and refer to any mammalian subject, particularly humans. Other subjects may include cattle, dogs, cats, guinea pigs, rabbits, rats, mice, horses, and so on. In some cases, the methods of the invention find use in experimental animals, in veterinary application, and in the development of animal models, including, but not limited to, rodents including mice, rats, and hamsters; and primates. II. Modes of Carrying Out the Invention

Before describing the present invention in detail, it is to be understood that this invention is not limited to particular formulations or process parameters as such may, of course, vary. It is also to be understood that the terminology used herein is for the purpose of describing particular embodiments of the invention only, and is not intended to be limiting.

Although a number of methods and materials similar or equivalent to those described herein can be used in the practice of the present invention, the preferred materials and methods are described herein.

The present invention is based on the discovery of a novel therapeutic methodology for effectively treating heart failure using an agonist of hypocretin receptor 2 (see Example 1). Additionally, the inventors have shown that patients having the major allele at the single nucleotide polymorphism rs7767652, which is located immediately upstream of the hypocretin receptor 2 (HCRTR2) gene, show significantly greater improvement in left ventricular function (e.g., at least a 10% improvement in ejection fraction) with treatment with beta blockers or ACE inhibitors than patients having the minor allele at rs7767652 (e.g., less than a 10% improvement in ejection fraction).

In order to further an understanding of the invention, a more detailed discussion is provided below regarding methods of treating subjects for heart failure.

A. Administration of an Agonist of Hypocretin Receptor 2

One or more agonists of hypocretin receptor 2 may be administered to an individual in need of treatment for heart failure. Exemplary agonists that can be used in the practice of the invention include orexin-A, orexin-B, [Alal l,D-Leul5]-orexin B, and SB-668875. In one embodiment, the agonist is an orexin mimic. The agonist may be recombinantly or synthetically produced.

At least one therapeutically effective cycle of treatment with an agonist of hypocretin receptor 2 will be administered to a subject for treatment of heart failure. By "therapeutically effective cycle of treatment" is intended a cycle of treatment that when administered, brings about a positive therapeutic response with respect to treatment of an individual for heart failure. Of particular interest is a cycle of treatment with an agonist of hypocretin receptor 2 that improves heart function or survival of a subject treated for heart failure. By "positive therapeutic response" is intended that the individual undergoing treatment according to the invention exhibits an improvement in one or more symptoms of heart failure, such as improved left ventricular function, increased blood pressure, or increased exercise capacity.

In certain embodiments, multiple therapeutically effective doses of compositions comprising one or more agonists of hypocretin receptor 2 and/or one or more other therapeutic agents, such as other drugs for treating heart failure will be administered. The compositions are typically, although not necessarily, administered orally, via injection (subcutaneously, intravenously, or intramuscularly), by infusion, or locally. Additional modes of administration are also contemplated, such as intracerebro ventricular, intrathecal, intracerebral, intraspinal, intranasal,

intraparenchymatous, pulmonary, rectal, transdermal, transmucosal, pericardial, intraarterial, intraocular, intraperitoneal, and so forth. The preparations according to the invention are also suitable for local treatment. For example, compositions may be administered directly by stereotactic injection into the brain. The particular preparation and appropriate method of administration are chosen to target an agonist to the hypo cretin receptor 2 receptors in the brain. In particular embodiments, compositions are administered into the brain or spinal cord of a subject. In one embodiment, a composition is administered into the cisterna magna or the rostral ventrolateral medulla of a subject.

The pharmaceutical preparation can be in the form of a liquid solution or suspension immediately prior to administration, but may also take another form such as a syrup, cream, ointment, tablet, capsule, powder, gel, matrix, suppository, or the like. The pharmaceutical compositions comprising one or more agonists of hypo cretin receptor 2 and other agents may be administered using the same or different routes of administration in accordance with any medically acceptable method known in the art.

In another embodiment of the invention, the pharmaceutical compositions comprising one or more agonists of hypo cretin receptor 2 and/or other agents are in a sustained-release formulation, or a formulation that is administered using a sustained- release device. Such devices are well known in the art, and include, for example, transdermal patches, and miniature implantable pumps that can provide for drug delivery over time in a continuous, steady-state fashion at a variety of doses to achieve a sustained-release effect with a non-sustained-release pharmaceutical composition.

The invention also provides a method for administering a conjugate comprising an agonist of hypocretin receptor 2 to a patient suffering from a condition that is responsive to treatment with an agonist of hypocretin receptor 2 contained in the conjugate or composition. The method comprises administering, via any of the herein described modes, a therapeutically effective amount of the conjugate or drug delivery system, preferably provided as part of a pharmaceutical composition. The method of administering may be used to treat any condition that is responsive to treatment with an agonist of hypocretin receptor 2. More specifically, the

compositions herein are effective in treating heart failure.

Those of ordinary skill in the art will appreciate which conditions a specific agonist of hypocretin receptor 2 can effectively treat. The actual dose to be administered will vary depending upon the age, weight, and general condition of the subject as well as the severity of the condition being treated, the judgment of the health care professional, and conjugate being administered. Therapeutically effective amounts can be determined by those skilled in the art, and will be adjusted to the particular requirements of each particular case.

Generally, a therapeutically effective amount will range from about 0.50 mg to 5 grams of an agonist of hypocretin receptor 2 inhibitor daily, more preferably from about 5 mg to 2 grams daily, even more preferably from about 7 mg to 1.5 grams daily. Preferably, such doses are in the range of 10-600 mg four times a day (QID), 200-500 mg QID, 25 - 600 mg three times a day (TID), 25-50 mg TID, 50-100 mg TID, 50-200 mg TID, 300-600 mg TID, 200-400 mg TID, 200-600 mg TID, 100 to 700 mg twice daily (BID), 100-600 mg BID, 200-500 mg BID, or 200-300 mg BID. The amount of compound administered will depend on the potency of the specific agonist of hypocretin receptor 2 and the magnitude or effect on heart function desired and the route of administration.

A purified agonist of hypocretin receptor 2 (again, preferably provided as part of a pharmaceutical preparation) can be administered alone or in combination with one or more other therapeutic agents, such as an angiotensin converting enzyme (ACE) inhibitor, an angiotensin II receptor blocker (ARB), a beta-blocker, digoxin, a diuretic, a blood vessel dilator, potassium, magnesium, an aldactone inhibitor, a calcium channel blocker, or an inotrope, or other medications used to treat a particular condition or disease according to a variety of dosing schedules depending on the judgment of the clinician, needs of the patient, and so forth. The specific dosing schedule will be known by those of ordinary skill in the art or can be determined experimentally using routine methods. Exemplary dosing schedules include, without limitation, administration five times a day, four times a day, three times a day, twice daily, once daily, three times weekly, twice weekly, once weekly, twice monthly, once monthly, and any combination thereof. Preferred compositions are those requiring dosing no more than once a day.

An agonist of hypocretin receptor 2 can be administered prior to, concurrent with, or subsequent to other agents. If provided at the same time as other agents, one or more agonists of hypocretin receptor 2 can be provided in the same or in a different composition. Thus, one or more agonists of hypocretin receptor 2 and other agents can be presented to the individual by way of concurrent therapy. By "concurrent therapy" is intended administration to a subject such that the therapeutic effect of the combination of the substances is caused in the subject undergoing therapy. For example, concurrent therapy may be achieved by administering a dose of a pharmaceutical composition comprising an agonist of hypocretin receptor 2 and a dose of a pharmaceutical composition comprising at least one other agent, such as another agonist of hypocretin receptor 2 or drug for treating heart failure, which in combination comprise a therapeutically effective dose, according to a particular dosing regimen. Similarly, one or more agonists of hypocretin receptor 2 and one or more other therapeutic agents can be administered in at least one therapeutic dose. Administration of the separate pharmaceutical compositions can be performed simultaneously or at different times (i.e., sequentially, in either order, on the same day, or on different days), as long as the therapeutic effect of the combination of these substances is caused in the subject undergoing therapy.

B. Detecting and Genotyping Gene Polymorphisms

The SNP rs7767652, which is located upstream of the HCRTR2 gene, can be used as a genetic marker for determining whether an individual with heart failure is likely to respond more favorably to medical intervention, such as treatment with a beta blocker or an ACE inhibitor. The major allele of rs7767652 has a C and the minor allele of rs7767652 has a T at position 27 (indicated by the [C/T]) of the following sequence:

GTTGATAAAATATTCCACATTATAAT[C/T]AAAGAAAGTACAAAAAAGTTA AATG (SEQ ID NO: l).

For genetic testing, a biological sample containing nucleic acids is collected from an individual in need of treatment for heart failure. The biological sample is typically blood, saliva, or cells from buccal swabbing, but can be any sample from bodily fluids, tissue or cells that contains genomic DNA or RNA of the individual. In certain embodiments, nucleic acids from the biological sample are isolated, purified, and/or amplified prior to analysis using methods well-known in the art. See, e.g.,

Green and Sambrook Molecular Cloning: A Laboratory Manual (Cold Spring Harbor Laboratory Press; 4^th edition, 2012); and Current Protocols in Molecular Biology (Ausubel ed., John Wiley & Sons, 1995); herein incorporated by reference in their entireties.

It is understood that SNPs can be detected in a sample by any suitable method known in the art. Detection of a SNP can be direct or indirect. For example, the SNP itself can be detected directly. Alternatively, the SNP can be detected indirectly from cDNAs, amplified RNAs or DNAs, or proteins expressed by the SNP. Any method that detects a single base change in a nucleic acid sample can be used. For example, allele-specific probes that specifically hybridize to a nucleic acid containing the polymorphic sequence can be used to detect SNPs. A variety of nucleic acid hybridization formats are known to those skilled in the art. For example, common formats include sandwich assays and competition or displacement assays.

Hybridization techniques are generally described in Hames, and Higgins "Nucleic Acid Hybridization, A Practical Approach," IRL Press (1985); Gall and Pardue, Proc. Natl. Acad. Sci. U.S.A., 63:378-383 (1969); and John et al Nature, 223:582-587 (1969).

Sandwich assays are commercially useful hybridization assays for detecting or isolating nucleic acids. Such assays utilize a "capture" nucleic acid covalently immobilized to a solid support and a labeled "signal" nucleic acid in solution. The clinical sample will provide the target nucleic acid. The "capture" nucleic acid and "signal" nucleic acid probe hybridize with the target nucleic acid to form a

"sandwich" hybridization complex.

In one embodiment, the allele-specific probe is a molecular beacon.

Molecular beacons are hairpin shaped oligonucleotides with an internally quenched fluorophore. Molecular beacons typically comprise four parts: a loop of about 18-30 nucleotides, which is complementary to the target nucleic acid sequence; a stem formed by two oligonucleotide regions that are complementary to each other, each about 5 to 7 nucleotide residues in length, on either side of the loop; a fluorophore covalently attached to the 5' end of the molecular beacon, and a quencher covalently attached to the 3' end of the molecular beacon. When the beacon is in its closed hairpin conformation, the quencher resides in proximity to the fluorophore, which results in quenching of the fluorescent emission from the fluorophore. In the presence of a target nucleic acid having a region that is complementary to the strand in the molecular beacon loop, hybridization occurs resulting in the formation of a duplex between the target nucleic acid and the molecular beacon. Hybridization disrupts intramolecular interactions in the stem of the molecular beacon and causes the fluorophore and the quencher of the molecular beacon to separate resulting in a fluorescent signal from the fluorophore that indicates the presence of the target nucleic acid sequence.

For SNP detection, the molecular beacon is designed to only emit fluorescence when bound to a specific allele of a SNP. When the molecular beacon probe encounters a target sequence with as little as one non-complementary nucleotide, the molecular beacon preferentially stay in its natural hairpin state and no fluorescence is observed because the fluorophore remains quenched. See, e.g., Nguyen et al. (2011) Chemistry 17(46): 13052-13058; Sato et al. (2011) Chemistry 17(41): 11650-11656; Li et al. (2011) Biosens Bioelectron. 26(5):2317-2322; Guo et al. (2012) Anal. Bioanal. Chem. 402(10):3115-3125; Wang et al. (2009) Angew. Chem. Int. Ed. Engl.

48(5):856-870; and Li et al. (2008) Biochem. Biophys. Res. Commun. 373(4):457- 461; herein incorporated by reference in their entireties.

Probes can readily be synthesized by standard techniques, e.g., solid phase synthesis via phosphoramidite chemistry, as disclosed in U.S. Patent Nos. 4,458,066 and 4,415,732, incorporated herein by reference; Beaucage et al, Tetrahedron (1992) 48:2223-2311; and Applied Biosystems User Bulletin No. 13 (1 April 1987). Other chemical synthesis methods include, for example, the phosphotriester method described by Narang et al, Meth. Enzymol. (1979) 68:90 and the phosphodiester method disclosed by Brown et al, Meth. Enzymol. (1979) 68: 109. Poly(A) or poly(C), or other non-complementary nucleotide extensions may be incorporated into polynucleotides using these same methods. Hexaethylene oxide extensions may be coupled to the polynucleotides by methods known in the art. Cload et al, J. Am. Chem. Soc. (1991) 113:6324-6326; U.S. Patent No. 4,914,210 to Levenson et al; Durand et al, Nucleic Acids Res. (1990) 18:6353-6359; and Horn et al, Tet. Lett. (1986) 27:4705-4708.

Alternatively, probes can be produced by amplification of a target nucleic acid using, e.g., polymerase chain reaction (PCR), nucleic acid sequence based

amplification (NASBA), ligase chain reaction (LCR), self-sustained sequence replication (3SR), Q-beta amplification; strand displacement amplification, or any other nucleic acid amplification method to produce a probe capable of hybridizing to the desired target sequence.

The probes may be coupled to labels for detection. There are several means known for derivatizing polynucleotides with reactive functionalities which permit the addition of a label. For example, several approaches are available for biotinylating probes so that radioactive, fluorescent, chemiluminescent, enzymatic, or electron dense labels can be attached via avidin. See, e.g., Broken et al, Nucl. Acids Res. (1978) 5:363-384 which discloses the use of ferritin-avidin-biotin labels; and Chollet et al, Nucl. Acids Res. (1985) 13: 1529-1541 which discloses biotinylation of the 5' termini of polynucleotides via an aminoalkylphosphoramide linker arm. Several methods are also available for synthesizing amino-derivatized oligonucleotides which are readily labeled by fluorescent or other types of compounds derivatized by amino- reactive groups, such as isothiocyanate, N-hydroxysuccinimide, or the like, see, e.g., Connolly, Nucl. Acids Res. (1987) 15:3131-3139, Gibson et al. Nucl. Acids Res. (1987) 15:6455-6467 and U.S. Patent No. 4,605,735 to Miyoshi et al. Methods are also available for synthesizing sulfhydryl-derivatized polynucleotides, which can be reacted with thiol- specific labels, see, e.g., U.S. Patent No. 4,757,141 to Fung et al, Connolly et al, Nucl. Acids Res. (1985) 13:4485-4502 and Spoat et al. Nucl. Acids Res. (1987) 15:4837-4848. A comprehensive review of methodologies for labeling DNA fragments is provided in Matthews et al, Anal. Biochem. (1988) 169: 1-25.

For example, probes may be fluorescently labeled. Guidance for selecting appropriate fluorescent labels can be found in Smith et al, Meth. Enzymol. (1987) 155:260-301; Karger et al, Nucl. Acids Res. (1991) 19:4955-4962; Guo et al. (2012) Anal. Bioanal. Chem. 402(10):3115-3125; and Molecular Probes Handbook, A Guide to Fluorescent Probes and Labeling Technologies, 11th edition, Johnson and Spence eds., 2010 (Molecular Probes/Life Technologies); herein incorporated by reference. Fluorescent labels include fluorescein and derivatives thereof, such as disclosed in U.S. Patent No. 4,318,846 and Lee et al, Cytometry (1989) 10: 151-164. Dyes for use in the present invention include 3-phenyl-7-isocyanatocoumarin, methyl coumarin-3- acetic acid (AMCA), acridines, such as 9-isothiocyanatoacridine and acridine orange, pyrenes, benzoxadiazoles, and stilbenes, such as disclosed in U.S. Patent

No.4, 174,384. Additional dyes include SYBR green, SYBR gold, Yakima Yellow, Texas Red, 3-(8-carboxypentyl)-3'-ethyl-5,5'-dimethyloxa-carbocyanine (CYA); 6- carboxy fluorescein (FAM); CAL Fluor Orange 560, CAL Fluor Red 610, Quasar Blue 670; 5,6-carboxyrhodamine-l lO (R110); 6-carboxyrhodamine-6G (R6G);

N',N',N',N'-tetramethyl-6-carboxyrhodamine (TAMRA); 6-carboxy-X-rhodamine (ROX); 2^*, 4^*, 5^*, 7^*, - tetrachloro-4-7-dichlorofluorescein (TET); 2^*, Ύ- dimethoxy-4^*, 5'-6 carboxyrhodamine (JOE); 6-carboxy-2',4,4',5',7,7'-hexachlorofiuorescein (HEX); Dragonfly orange; ATTO-Tec; Bodipy; ALEXA; VIC, Cy3, Cy5, and Cy7. These dyes are commercially available from various suppliers such as Life Technologies (Carlsbad, CA), Biosearch Technologies (Novato, CA), and Integrated DNA

Technolgies (Coralville, Iowa). Fluorescent labels include fluorescein and derivatives thereof, such as disclosed in U.S. Patent No. 4,318,846 and Lee et al, Cytometry

(1989) 10: 151-164, and 6-FAM, JOE, TAMRA, ROX, HEX-1, HEX-2, ZOE, TET-1 or NAN-2, and the like.

Fluorophores may be covalently attached to a particular nucleotide, for example, and the labeled nucleotide incorporated into the probe using standard techniques such as nick translation, random priming, and PCR labeling.

Alternatively, a fiuorophore may be covalently attached via a linker to a

deoxycytidine nucleotide that has been transaminated. Methods for labeling probes are described in U.S. Pat. No. 5,491,224 and Molecular Cytogenetics: Protocols and Applications (2002), Y.-S. Fan, Ed., Chapter 2, "Labeling Fluorescence In Situ Hybridization Probes for Genomic Targets," L. Morrison et al, p. 21-40, Humana Press; which are herein incorporated by reference.

One of skill in the art will recognize that other luminescent agents or dyes may be used in lieu of fluorophores as label containing moieties. Other luminescent agents, which may be used, include, for example, radio luminescent,

chemiluminescent, bioluminescent, and phosphorescent label containing moieties, as well as quantum dots. Alternatively, in situ hybridization of chromosomal probes may be employed with the use of detection moieties visualized by indirect means. Probes may be labeled with biotin or digoxygenin using routine methods known in the art, and then further processed for detection. Visualization of a biotin-containing probe may be achieved via subsequent binding of avidin conjugated to a detectable marker. Chromosomal probes hybridized to target regions may alternatively be visualized by enzymatic reactions of label moieties with suitable substrates for the production of insoluble color products. Each probe may be discriminated from other probes within the set by choice of a distinct label. A biotin-containing probe within a set may be detected via subsequent incubation with avidin conjugated to alkaline phosphatase (AP) or horseradish peroxidase (HRP) and a suitable substrate, e.g., 5- bromo-4-chloro-3-indolylphosphate and nitro blue tetrazolium (NBT) serve as substrates for alkaline phosphatase, whereas diaminobenzidine serves as a substrate for HRP.

In another embodiment, detection of the SNP sequence is performed using allele-specific amplification. In the case of PCR, amplification primers can be designed to bind to a portion of one of the disclosed genes, and the terminal base at the 3' end is used to discriminate between the major and minor alleles or mutant and wild-type forms of the genes. If the terminal base matches the major or minor allele, polymerase-dependent three prime extension can proceed. Amplification products can be detected with specific probes. This method for detecting point mutations or polymorphisms is described in detail by Sommer et al. in Mayo Clin. Proc. 64: 1361- 1372 (1989).

Tetra-primer ARMS-PCR uses two pairs of primers that can amplify two alleles of a SNP in one PCR reaction. Allele-specific primers are used that hybridize at the SNP location, but each matches perfectly to only one of the possible alleles. If a given allele is present in the PCR reaction, the primer pair specific to that allele will amplify that allele, but not the other allele of the SNP. The two primer pairs for the different alleles may be designed such that their PCR products are of significantly different length, which allows them to be distinguished readily by gel electrophoresis. See, e.g., Munoz et al. (2009) J. Microbiol. Methods. 78(2):245-246 and Chiapparino et al. (2004) Genome. 47(2):414-420; herein incorporated by reference.

SNPs may also be detected by ligase chain reaction (LCR) or ligase detection reaction (LDR). The specificity of the ligation reaction is used to discriminate between the major and minor alleles of the SNP. Two probes are hybridized at the SNP polymorphic site of a nucleic acid of interest, whereby ligation can only occur if the probes are identical to the target sequence. See e.g., Psifidi et al. (2011) PLoS One 6(l):el4560; Asari et al. (2010) Mol. Cell. Probes. 24(6):381-386; Lowe et al. (2010) Anal Chem. 82(13):5810-5814; herein incorporated by reference.

SNPs can also be detected in a biological sample by sequencing and SNP typing. In the former method, one simply carries out whole genome sequencing of a patient sample, and uses the results to detect the present sequences. Whole genome analysis is used in the field of "personal genomics," and genetic testing services exist, which provide full genome sequencing using massively parallel sequencing.

Massively parallel sequencing is described e.g. in US 5,695,934, entitled "Massively parallel sequencing of sorted polynucleotides," and US 2010/0113283 Al, entitled "Massively multiplexed sequencing." Massively parallel sequencing typically involves obtaining DNA representing an entire genome, fragmenting it, and obtaining millions of random short sequences, which are assembled by mapping them to a reference genome sequence.

Commercial services also exist which will genotype approximately 1 million

SNPs for a fixed fee. SNP analysis can be carried out by a variety of methods that do not involve massively parallel random sequencing. As described below, a

commercially available MassARRAY system can be used. This system uses matrix- assisted laser desorption ionization time-of- flight mass spectrometry (MALDI-TOF MS) coupled with single-base extension PCR for high-throughput multiplex SNP detection. Another commercial SNP system is made by Illumina. The Illumina Golden Gate assay generates SNP specific PCR products that are subsequently hybridized to beads either on a solid matrix or in solution. Three oligonucleotides are synthesized for each SNP: two allele specific oligos (ASOs) that distinguish the SNP, and a locus specific sequence (LSO) just downstream of the SNP. The ASO and LSO sequences also contain target sequences for a set of universal primers (PI through P3 in the adjacent figure), while each LSO also contains a particular address sequences (the "illumicode") complementary to sequences attached to beads.

As another example, Affymetrix SNP arrays (also used in the present examples) use multiple sets of short oligonucleotide probes for each known SNP. The design of an SNP array such as manufactured by Affymetrix and Illumina is described further in LaFamboise, "Single nucleotide polymorphism arrays: a decade of biological, computational and technological advances," Nuc. Acids Res. 37(13):4181- 4193 (2009), which provides additional description of methods for detecting SNPs.

Another technology useful in SNP analysis is PCR-dynamic allele specific hybridization (DASH), which involves dynamic heating and coincident monitoring of DNA denaturation, as disclosed by Howell et al. (Nat. Biotech. 17:87-88, 1999). A target sequence is amplified (e.g., by PCR) using one biotinylated primer. The biotinylated product strand is bound to a streptavidin-coated microtiter plate well (or other suitable surface), and the non-biotinylated strand is rinsed away with alkali wash solution. An oligonucleotide probe, specific for one allele (e.g., the wild-type allele), is hybridized to the target at low temperature. This probe forms a duplex DNA region that interacts with a double strand-specific intercalating dye. When subsequently excited, the dye emits fluorescence proportional to the amount of double-stranded DNA (probetarget duplex) present. The sample is then steadily heated while fluorescence is continually monitored. A rapid fall in fluorescence indicates the denaturing temperature of the probe-target duplex. Using this technique, a single-base mismatch between the probe and target results in a significant lowering of melting temperature (Tm) that can be readily detected.

A variety of other techniques can be used to detect polymorphisms, including but not limited to, the Invader assay with Flap endonuclease (FEN), the Serial Invasive Signal Amplification Reaction (SISAR), the oligonucleotide ligase assay, restriction fragment length polymorphism (RFLP), single-strand conformation polymorphism, temperature gradient gel electrophoresis (TGGE), and denaturing high performance liquid chromatography (DHPLC). See, for example, Molecular Analysis and Genome Discovery (R. Rap ley and S. Harbron eds., Wiley 1^st edition, 2004); Jones et al. (2009) New Phytol. 183(4):935-966; Kwok et al. (2003) Curr Issues Mol. Biol. 5(2):43-60; Munoz et al. (2009) J. Microbiol. Methods 78(2):245-246;

Chiapparino et al. (2004) Genome. 47(2):414-420; Olivier (2005) Mutat Res. 573(1- 2): 103-110; Hsu et al. (2001) Clin. Chem. 47(8): 1373-1377; Hall et al. (2000) Proc. Natl. Acad. Sci. U.S.A. 97(15):8272-8277; Li et al. (2011) J. Nanosci. Nanotechnol. 11(2):994-1003; Tang et al. (2009) Hum. Mutat. 30(10): 1460-1468; Chuang et al. (2008) Anticancer Res. 28(4A):2001-2007; Chang et al. (2006) BMC Genomics 7:30; Galeano et al. (2009) BMC Genomics 10:629; Larsen et al. (2001)

Pharmacogenomics 2(4):387-399; Yu et al. (2006) Curr. Protoc. Hum. Genet. Chapter 7: Unit 7.10; Lilleberg (2003) Curr. Opin. Drug Discov. Devel. 6(2):237-252; and U.S. Pat. Nos. 4,666,828; 4,801,531; 5,110,920; 5,268,267; 5,387,506; 5,691,153; 5,698,339; 5,736,330; 5,834,200; 5,922,542; and 5,998,137 for a description of such methods; herein incorporated by reference in their entireties.

If the polymorphism is located in the coding region of a gene of interest, the SNP can be identified indirectly by detection of the variant protein produced by the SNP. Variant proteins (i.e., containing an amino acid substitution encoded by the SNP) can be detected using antibodies specific for the variant protein. For example, immunoassays that can be used to detect variant proteins produced by SNPs include, but are not limited to, immunohistochemistry (IHC), western blotting, enzyme-linked immunosorbent assay (ELISA), radioimmunoassays (RIA), "sandwich"

immunoassays, fluorescent immunoassays, and immunoprecipitation assays, the procedures of which are well known in the art (see, e.g., Schwarz et al. (2010) Clin. Chem. Lab. Med. 48(12): 1745-1749; The Immunoassay Handbook (D.G. Wild ed., Elsevier Science; 3^rd edition, 2005); Ausubel et al, eds, 1994, Current Protocols in Molecular Biology, Vol. 1 (John Wiley & Sons, Inc., New York); Coligan Current Protocols in Immunology (1991); Harlow & Lane, Antibodies: A Laboratory Manual (1988); Handbook of Experimental Immunology, Vols. I-IV (D.M. Weir and C.C. Blackwell eds., Blackwell Scientific Publications); herein incorporated by reference herein in their entireties).

Agents for detecting which allele is present at SNP rs7767652 in a subject can be provided in the form of a kit. The kit may include one or more agents for detection of the SNP, such as allele-specific hybridization probes, PCR primers, or a SNP microarray for determining which allele is present. The kit may further comprise a container for holding a biological sample isolated from a human subject for genetic testing and printed instructions for reacting agents with the biological sample or a portion of the biological sample to detect the presence of at least one SNP in the biological sample. The agents may be packaged in separate containers. The kit may further comprise one or more control reference samples or other reagents for detecting or genotyping SNPs.

In certain embodiments, the kit further comprises reagents for performing dynamic allele-specific hybridization (DASH), Tetra-primer ARMS-PCR, a TaqMan 5 '-nuclease assay; an Invader assay with Flap endonuclease (FEN), a Serial Invasive Signal Amplification Reaction (SISAR), an oligonucleotide ligase assay, restriction fragment length polymorphism (RFLP), single-strand conformation polymorphism, temperature gradient gel electrophoresis (TGGE), denaturing high performance liquid chromatography (DHPLC), sequencing, or an immunoassay.

The kit can comprise one or more containers for compositions contained in the kit. Compositions can be in liquid form or can be lyophilized. Suitable containers for the compositions include, for example, bottles, vials, syringes, and test tubes.

Containers can be formed from a variety of materials, including glass or plastic. The kit can also comprise a package insert containing written instructions for methods of detecting SNPs and predicting future skin aging.

III. Experimental

Below are examples of specific embodiments for carrying out the present invention. The examples are offered for illustrative purposes only, and are not intended to limit the scope of the present invention in any way.

Efforts have been made to ensure accuracy with respect to numbers used (e.g., amounts, temperatures, etc.), but some experimental error and deviation should, of course, be allowed for.

Example 1

The Neuromodulator Hypocretin/Orexin Receptor 2 Regulates Cardiac Function

Heart failure is the most common admitting diagnosis for patients in the US and has a worse prognosis than most cancers (Lloyd- Jones et al. (2010) Circulation 121 :e46-e215). Although clearly heritable (Lee et al. (2006) N. Engl. J. Med.

355: 138-147), vanishingly few genetic associations have been identified. We studied 866 patients with heart failure and compared the extremes of response to medical therapy using a genome wide approach that informed the systems-based design of a customized single nucleotide variant array. Polymorphism rs7767652, upstream of the hypocretin (orexin) receptor 2 (HCRTR2) gene, is associated with improvement in heart function (odds ratio (OR) 0.40, p = 3.29^χ 10^~5) and replicated in an independent heart failure cohort. Allele-specific luciferase reporter assays demonstrated higher relative enhancer activity with the rs7767652 major allele while expression of HCRTR2 was higher in human diseased tissue compared to normal tissue. In a heart failure model, HCRTR2 deficient mice exhibited poorer cardiac function, treadmill exercise capacity and greater myocardial scarring compared to wild-type mice.

Finally, orexin, an HCRTR2 agonist was successful in rescuing heart function in a mouse model of heart failure. Here, a systems approach identifies a novel genetic contribution to human heart failure and a promising therapeutic agent efficacious in a mouse model.

Heart failure, a syndrome characterized by impaired cardiac function, affects over five million people in the US and is expected to affect 3.5% of the population over the next 20 years (Heidenreich et al. (2012) Circulation 123:933-944). Heart failure can result from a wide range of conditions, commonly coronary artery disease and hypertension (Ho et al. (1993) J. Am. Coll. Cardiol. 22:6A-13A). Few studies have explored the genetic basis of heart failure. Candidate gene studies identified associations between common variants in HSPB7 and FRMD4B and dilated cardiomyopathy (Stark et al. (2010) PLoS Genetics 6 :e 1001167) or advanced heart failure (Cappola et al. (2010) Circ. Cardiovasc. Genet. 3: 147-154). Targeted genotyping of common variants in ADRB1 and GRK5, members of the β-adrenergic receptor signaling pathway, demonstrated associations with survival (Cresci et al. (2009) J. Am. Coll. Cardiol. 54:432-444; Liggett et al. (2008) Nature Medicine

14:510-517). In the limited genome-wide studies of heart failure (Larson et al. (2007) BMC Medical Genetics 8 Suppl. 1 :S5; Smith et al. (2010) Circ. Cardiovasc. Genet. 3:256-266; Villard et al. (2011) Eur. Heart J. 32: 1065-1076; Meder et al. (2014) Eur. Heart J. 35: 1069-1077), only one common variant associated with dilated

cardiomyopathy at genome-wide significance has replicated (Meder et al, supra).

With optimal medical therapy, some exhibit a dramatic improvement in heart function while others deteriorate. Most heart failure genomic studies to date have been cross-sectional. We explored the genetic basis of extreme response to medical intervention, identifying genetic markers of dramatic improvement versus clinical deterioration towards heart transplant.

We first designed a customized single nucleotide variation (SNV) array, generating posterior probabilities by compiling publicly available myocardial gene expression data and combining this with a genome -wide association (GWA) analysis of extreme response to heart failure therapy and other resources (FIG. 1). We evaluated 866 patients with heart failure at Stanford University. Patients at the extreme of the distribution of dynamic change were genotyped using a 561,464 SNV genotype chip that was used to inform the design of a customized genotyping array. All publicly available microarray gene expression data from human cardiac tissue (Dewey et al. (2011) Circ. Cardiovasc. Genet. 4:26-35) was collected from the Gene Expression Omnibus (GEO) database (ncbi.nlm.nih.gov/geo/) and normalized using a median-absolute-deviation algorithm (Smyth et al. (2003) Methods 31 :265- 273). Differences in expression levels between failing and normal human

myocardium were measured using significance analysis of microarrays (Storey et al. (. 2003) Methods Mol. Biol. 224: 149-157). The significance value (<i-score) was used to weight (Genovese et al. (2006) Biometrika 93:509-524; Roeder et al. (2007) Genetic Epidemiology 31 :741-747) /^-values from the informative GWA analysis. The top-ranking genie SNPs were included in the customized array (Table 3). The top intergenic SNPs from the informative GWA analysis, ranked by unweighted p- values, were also included in the customized array (Table 4).

Next, gene coexpression networks with scale-free topology were created from a collated file of all human myocardial gene expression data (n = 340 expression arrays) (Dewey et al, supra). Genes were clustered using average linkage hierarchical clustering and adaptively assigned to modules using a dynamic tree cutting algorithm (Langfelder et al. (2008) Bioinformatics 24:719-720). Modules with topology shared between normal and failing myocardium were ranked according to differential gene expression (quantified by average <i-score for modular gene expression differences between normal and failing myocardium) to identify sets of genes associated with adaptation to heart failure. Common variants from the hub genes, defined as the genes in each module with maximum intra-modular

connectivity, from the highest-ranking modules were included in the customized array (Table 5). The top-ranking module with SLC25A30 as its hub is depicted in FIG. 4. Next, Pubmed abstracts were retrieved and semantic mining was used to create a gene-gene interaction network (Ashley et al. (2006) Circulation 114:2644-2654; King et al. (2005) Physiological Genomics 23: 103-118). The top common variants of the hub genes, ranked by average <i-scores (Storey et al. (2003) Methods Mol. Biol.

224: 149-157), were included in the customized array (Table 6). Finally, a curated set of SNVs derived from multiple sources including significant SNPs from prior cardiovascular GWA studies (Table 7), microRNA gene polymorphisms, and coding SNVs in genes known to be targets of pharmacologic agents used in heart failure were included. The selected SNVs were placed on a customized Illumina Golden Gate platform.

Baseline characteristics of the 591 heart failure patients, 23% of whom had an absolute improvement of 10% in their ejection fraction, genotyped with the custom array are in Table 1. An additive logistic regression model, adjusted for age, sex and race, was used to test the association between each SNV on the customized array and improved ejection fraction (Manhattan plot, FIG. 2A; quantile-quantile plot, FIG. 5). Patients with the minor allele for rs7767652 were less likely to have a response of greater than 10% in ejection fraction (OR 0.394, p = 3.29^χ 10^~5), with statistical significance below the Bonferroni cutoff of 0.05/1402 = 3.56^x 10^"5. This variant was chosen for custom array as one of the intergenic SNVs. In the informative GWA, patients with the minor allele for rs7767652 were less likely to have a response of greater than 20% in ejection fraction compared to transplanted patients (OR 0.279, p = 0.0040). The meta-analysis p-value was 9.04x l0^"6. rs7767652 is approximately 2,700 base pairs upstream from HCRTR2, within a haplotype block that encompasses several alternative splice sites (Chen et al. (2010) Biochemical J. 427:377-390) and the 5' untranslated region (UTR) ofHCRTR2 (FIG. 2B). Top results from each of the custom array selection categories are presented in Table 2 and Tables 3-8. We sought replication in 798 heart failure patients from the cross-sectional Penn Heart Failure Study (Cappola et al. (2010) Circ. Cardiovasc. Genet. 3: 147-154). After multivariate adjustment, heart failure patients with the minor allele for rs7767652 were less likely to have an EF greater than 35%> (OR 0.769, p = 0.021). To assess the potential regulatory role of the lead variant rs7767652 upstream of HCRTR2, we first mapped putative transcription factor binding sites in silico using various

bioinformatics tools, such as TRANSFAC, JASPAR, and Matlnspector (Table 8).

The minor (T) allele predicted disruption of a TCF4 binding site containing the motif ATCAAAG. We then measured allele-specific gene expression using luciferase reporters containing the predicted binding site in trans fected C2C12 myoblast cells (FIG. 3A). At baseline, transfection with the minor allele (T) construct of rs7767652 resulted in lower luciferase activity compared to the major allele (C) construct (p < 0.05) (FIG. 3B). We also observed that the minor allele significantly disrupted β- catenin/TCF4-mediated transactivation in cells overexpressing human β-catenin and TCF4 cofactors (p < 0.01). These effects were abolished in the presence of a TCF4 transactivation mutant (p < 0.005).

To assess a myocardial role for this neuropeptide receptor, we compared HCRTR2 expression in diseased hearts (dilated n = 162, ischemic n = 143) with control tissue (n = 45). Expression of HCRTR2 was greater in dilated cardiomyopathy (p = 0.00017) and ischemic cardiomyopathy samples (p = 0.000044) compared to controls (FIG. 3C).

To investigate a causal role for HCRTR2 in heart failure pathogenesis, we performed heart function testing via ultrasound in HCRTR2 transcription-disrupted (Mochizuki et al. (2011) Proc. Natl. Acad. Sci. U.S.A.108:4471-4476) (TD) mice

(n=10) compared to wild-type (n = 11). There was greater diastolic dysfunction in the HCRTR2 TD mice compared to the wild-type mice (p = 0.05, FIG. 3D), but no significant difference in systolic function (FIG. 3E). Similar studies in a small number ofHCRTR2 knockout mice (Chemelli et al. (1999) Cell 98:437-451) (n=5) compared to controls (n = 5) also demonstrated a trend towards greater diastolic dysfunction. HCRTR2 TD mice that underwent chemical stress with two weeks of angiotensin II and isoproterenol infusion had a smaller increase in treadmill exercise capacity (FIG. 3F) and greater evidence of myocardial fibrosis (FIG. 3G) compared to wild type mice.

Finally, to assess the potential of HCRTR2 as a novel therapeutic target for heart failure, wild-type mice were infused with saline (n = 6) or the HCRTR2 agonist orexin A (n = 7) for four weeks. At week two, all mice underwent infusion of angiotensin II and isoproterenol, neurohormones that mimic human heart failure. Echocardiography was performed blinded to drug infusion status at baseline, 2, 3 and 4 weeks (FIG. 3H). Mice with orexin A had better systolic function compared to controls (p = 0.045).

The HCRTR2 gene encodes a G-coupled receptor that binds hypocretin (orexin) A and B, neuropeptides involved in appetite (Lubkin et al. (1998) Biochem. Biophys. Res. Commun. 253:241-245; Sakurai et al. (1998) Cell 92: 1 page following 696) and sleep (Lin et al. (1999) Cell 98:365-376) regulation. Although highly expressed in the hypothalamus (Sakurai et al, supra), its expression has been documented in several tissues including gut (Kirchgessner et al. (1999) Neuron 24:941-951) and adrenal glands (Lopez et al. (1999) Endocrinology 140:5991-5994). Alternative splice variants of HCRTR2 are expressed in heart tissue (Chen et al.

(2010) Biochem. J. 427:377-390). Orexins provoke increases in blood pressure and heart rate via a centrally-mediated response (Shirasaka et al. (1999) Am. J. Physiol. 277:R1780-1785; Samson et al. (1999) Brain Research 831 :248-253; Matsumura et al. (2001) Hypertension 37: 1382-1387). Their effects on myocardial function have not been studied.

We used a customized gene array to identify a genomic variant in the regulatory region of HCRTR2 associated with improved left ventricular function. Allele-specific reporter assays at this locus suggested a role in regulating nearby gene transcription. Our finding that mice infused with an HCRTR2 agonist were protected from chemical stress-induced ventricular dysfunction suggests a promising role for HCRTR2 in regulation of ventricular function. Unlike the neurohormonal axis of the renin- angiotensin- aldosterone system, the effects of hypocretins on myocardial function have not yet been characterized, but they may play a role in regulating myocardial fibrosis. Additional studies to characterize central versus peripheral roles these neuropeptides play in regulation of myocardial function are beyond the scope of this work but may inform future therapeutic development.

METHODS

Study Populations

Stanford Heart Failure Genomic Analyses

Patients were recruited from Stanford University Medical Center (Stanford, California) and the Palo Alto Veterans Hospital (Palo Alto, California). 866 patients with clinically diagnosed heart failure referred for subspecialty care between 2005 and 2009 and who had an echocardiogram performed with an ejection fraction less than 55% were included. Patients who could not be contacted by telephone (n=465), or with insufficient clinical data (n = 3) were excluded. Those with congenital heart disease or with cardiomyopathy due to infiltrative disease, a peri-partum state, infection or chemotherapy were excluded. Patients whose ejection fraction improved after surgical or percutaneous revascularization, resynchronization therapy with biventricular pacing, or cardioversion were excluded. Those with a myocardial biopsy suggestive of viral cardiomyopathy who responded to viral medication therapy within 30 days of administration were excluded. Patients who suffered an acute myocardial infarction and showed subsequent improvement in their ejection fraction within 3 months of the infarct, or patients with a history of substance abuse with illicit drugs or alcohol within 6 months prior to the study were also excluded. Ejection fraction was measured using routine transthoracic echocardiography obtained by trained echo car diographers using two-dimensional scanning in the parasternal long axis, parasternal short axis and apical views. Change in ejection fraction was measured as maximal difference between lowest recorded ejection fraction and the highest subsequent recorded ejection fraction. Written informed consent was obtained from study participants in accordance with the Stanford University Internal Review Board policy.

For the informative GWAS, case subjects (n = 29) were Caucasian patients whose ejection fraction had improved by 20% or greater while on medical therapy. Control patients (n = 37) were selected from those followed in the Stanford University transplant clinic who had demonstrated lack of clinical improvement prior to transplant despite medical therapy. The case and control patients were matched by age, sex, race, medical therapy received, baseline ejection fraction on

echocardiography, type of cardiomyopathy and duration of heart failure. The 730 heart failure patients remaining underwent genotyping with the custom genotyping array. After exclusion of patients with insufficient clinical data and failure to pass quality control measures, there were a total of 591 successfully genotyped patients. Patients' charts were retrospectively reviewed and their baseline demographics, clinical characteristics and medication use were recorded. Serial measurements from clinically available echocardiograms were also recorded. Outcome data was obtained through detailed phone interviews.

UPenn Replication

For the independent replication study, patients with heart failure were recruited from the University of Pennsylvania as previously described (Cappola et al. (2010) Circ. Cardiovasc. Genet. 3: 147-154). Briefly, Caucasian heart failure patients (n = 798) were recruited from the Penn Heart Failure Study from an ongoing, NIH- sponsored (HL077101) prospective observational study of patients with advanced heart failure referred for subspecialty care at the University of Pennsylvania Health System (Philadelphia, Pennsylvania). The primary inclusion criterion was a clinical diagnosis of heart failure with abnormal left ventricular function. Extensive clinical data were collected at enrollment. Written informed consent was obtained from study participants in accordance with the Stanford University Internal Review Board policy.

Genotyping

Informative GWAS

Genomic DNA was isolated from whole blood using a commercial DNA extraction kit (Gentra Purgene Kit). The samples were genotyped using the Illumina 550K SNP platform at the Hudson Alpha Institute for Biotechnology (Huntsville, Alabama). A total of 36,244 SNPs with call rates < 90%, that had a minor allele frequency (MAF) < 1%, or deviated significantly from Hardy- Weinberg equilibrium (P < 0.001) were excluded. 525,220 common SNPs remained for association analysis. At the sample-level, 6 cases with call rates < 90% were excluded. There were no cases identified with genetic relatedness using pairwise identity by state-based analysis in PLINK (Purcell et al. (2007) Am. J. Hum. Genet. 81 :559-575).

Customized SNV Array

Genomic DNA was isolated from either whole blood or buccal swab samples using a commercial DNA extraction kit (Gentra Purgene Kit). The samples were genotyped using a customized 1536 SNP GoldenGate Assay on an Illumina BeadLab system at the HudsonAlpha Institute for Biotechnology (Huntsville, Alabama).

Genotype calls were made using Bead Studio Genotyping Module software package, version 3.1 (Illumina, San Diego,CA). A total of 134 SNPs with call rates < 50%, that had a MAF < 1%, or that deviated significantly from Hardy- Weinberg

Equilibrium (p < 0.001) were excluded. At the sample-level, 118 cases with call rates < 80%) were excluded.

UPenn Replication

Genomic DNA was isolated from peripheral blood mononuclear cells and

DNA quality was assessed utilizing optical absorbance and minigels. Genotyping analysis of the SNP selected for validation, rs7767652, was performed using the MassArray system from Sequenom. Locus-specific PCR and detection primers were designed using the MassArray Assay Design 3.0 software (Sequenom). The SNP did not deviate from Hardy- Weinberg Equilibrium (p < 0.05), the minor allele frequency was 23% (similar to HapMap CEU MAF) and 797 out of 798 samples were successfully called.

Analytical Methods

Informative GWAS

Genome-wide association analysis was performed in PLINK (Purcell et al. (2007) Am. J. Hum. Genet. 81 :559-575) using an additive logistic regression model to measure the association between each SNP and an improvement in ejection fraction of greater than 20% compared to control patients who had undergone transplant. Cases and controls were matched by age, sex and race, and no further adjustment was made for these covariates. A Manhattan plot of -log 1 OP was generated using Hap lo view (FIG. 6). Values for a quantile-quantile plot were generated using PLINK to evaluate the potential impact of population stratification (FIG. 7). Odds ratio (OR) values were measured as OR per minor allele of each SNP. The threshold for statistical (genome-wide) significance was estimated at 5 x 10 ^s . Differential Gene Expression

Microarray gene expression data was collected from every publicly available human myocardial dataset (n = 340) in the Gene Expression Omnibus database (ncbi.nlm.nih.gov/geo/). Experimental conditions and identifiable phenotypes from each dataset were also extracted. Data normalization was performed using a median- absolute-deviation algorithm (Smyth et al. (2003) Methods 31 :265-273) to allow comparison of gene expression between different platforms and experimental conditions. Gene expression data were combined and normalized using the median- absolute-deviation algorithm as described. To assess for differences in gene expression between myopathic (n = 127) tissue and control (n = 213) tissue, we ran a significance analysis of microarrays (SAM) analysis (Storey et al. (2003) Methods Mol. Biol. 224: 149-157) using the samr module in R. The statistical metric (d-score) from this analysis was subsequently used to further inform the analyses below and to rank genes and SNPs of interest for final selection in the custom genotyping array. Genie SNP Weighting

Genie SNPs from the Illumina 550K platform were defined as SNPs between 2,000 base pairs upstream of the gene start site and 2,000 base pairs downstream from the stop site. The /^-values of each genie SNP from the informative Genome- Wide Association Study were then weighted by the <i-scores of their corresponding genes from the differential gene expression analysis using the method described by

Genovesce et al. (Genovese et al. (2006) Biometrika 93:509-524). The weights are assigned as follows:

where Wj is the weight assigned to the SNPs of to the i gene, r is the strength of the weighting desired (value of r = 2 was selected), with higher values resulting in a more influential effect on the GWAS p-values, U, is the significance statistic and U_m is the average value of the test statistic. The final weighted p- value was calculated as:

The "weighted FDR" (Roeder et al. (2007) Genetic Epidemiology 31 :741 - 747) module in R was used to perform for this analysis.

Coexpression Network

Pair- wise gene-gene adjacencies were calculated by raising the pair- wise pearson correlation to an exponent of nine, where corr is the pearson correlation and j is the adjacency between gene; and gene_j:

Aj = corr(genei, gene_j)⁹ ( 1 )

This scaling exponent was empirically found to produce adjacencies that best satisfied approximate scale free topology (Jeong et al. (2000) Nature 407:651-654). Genes were clustered using average linkage hierarchical clustering and adaptively assigned to modules using the dynamic tree cut algorithm developed by Horvath, et al. (Langfelder et al. (2008) Bioinformatics 24:719-720). Mean intra-modular connectivities were calculated for each gene in each module from the intramodular adjacency, where Ki_n(i₎ is the intra-modular connectivity for gene; and A_j is the adjacency calculated as above, and N is the number of genes in the module:

K_in(l) = (∑_j≠i Ai_j )/(N-l) (2)

The gene in each module with the maximum intra-modular connectivity was subsequently chosen as the module hub gene. We then identified gene modules that exhibited significant differential expression between normal and failing human myocardium. To this end, gene modules were ranked according to the average significance of differential expression between normal and failing myocardium, or d- score (see above), of the module genes. Hub genes for the top 65 modules ranked in this manner were chosen for single-nucleotide polymorphism selection.

Semantic Network

The building of a gene-gene interaction network was automated using semantic mining of publicly available literature as previously described (Ashley et al. (2006) Circulation 114:2644-2654; King et al. (2005) Physiological Genomics

23 : 103 - 118). Publicly available Pubmed abstracts were retrieved and a search was performed for sentences containing 19,000 gene/protein terms from the gene expression microarray platforms. When multiple gene terms appear in the same sentence with an intervening verb suggestive of a relationship, such as "upregulates", "inhibits", or "binds to", an interaction is established. 9,404 genes/proteins were found to have at least one semantic interaction and a total of over 54,000 gene-gene interactions were established. The 3,514 subnetworks with fewer than 6 members were excluded from further analysis. The differential gene expression analysis described above was used to rank the semantic subnetworks by averaging the absolute (i-scores of the members in each subnet. Hub genes from the top 200 modules ranked in this manner were considered for customized gene array SNP selection.

Candidate Gene Selection

Candidate genes were selected for consideration for inclusion in the custom array from a number of sources. Genes known to be direct targets of the

pharmacologic agents used in heart failure treatment, their downstream targets, and associated downstream pathways were included. Additionally, genes known to be involved in alteration of pharmacodynamics of therapeutics for heart failure were included. Genes involved in the renin-angiotensin-aldosterone and adrenergic stimulation pathways were specifically included. A set of genes representing genes known to be responsible for monogenic, familial forms of dilated cardiomyopathy and associated members of the membrane cytoskeletal linkage and Z-disc complexes were also included. Furthermore, genes known to be involved in alteration of calcium handling and energetics were included.

Selection of SNPs from Candidate and Network Genes

A list of genes for consideration to include on the customized array was generated as described above. In order to select single nucleotide polymorphisms for these genes, we queried the Genome Variation Server (gvs.gs.washington.edu/GVS/) for tagged SNPs. We specified a Minor Allele Frequency cutoff of 0.2 from the CEU HapMap population to improve relevance of selected SNPs to our customized array genotyped population. Each candidate, hub, and nexus gene, including 5000 base pairs up and downstream from the 5 ' and 3 ' untranslated regions, was queried for tagged SNPs meeting the MAF cutoff. SNPs were selected with the tagSNP algorithm with an r² threshold of 0.8, minimum data coverage of 0.85, and minimum data coverage for clustering of 0.70. SNPs selected were then ranked by the weighted p- value from the informative GWA study. SNPs not represented on the Illumina

550K chip were excluded for further evaluation. In case of multiple SNPs within the same bin, the highest p-value followed by the highest MAF was used to rank within- bin SNPs, and only the highest-ranked SNP was selected. After selection within individual bins, tagged SNPs from each gene were ranked first by weighted /?-value followed by MAF. For each gene, a maximum of five SNPs were selected, with a maximum informative GWA p-value of 0.5 for hub and nexus gene SNPs and 0.05 for candidate gene SNPs.

Previously Identified SNPs

Literature reviews were performed evaluating evidence for and against associations between beta-blocker, ACE-inhibitor, aldosterone antagonist, antiarrhythmics, and known candidate gene polymorphisms (Wheeler et al. (2008) J. Cardiovasc. Transl. Res. 1 :25-36; Perez et al. (2008) J. Cardiovasc. Transl. Res.

1 : 155-165). Highly significant polymorphisms identified in previous genome wide association studies of coronary artery disease, diabetes mellitus, nicotine dependence, and atrial fibrillation were similarly catalogued. Additionally, SNPs identified in the Framingham Heart Study GWA associated with echo car dio graphic findings of LV size and wall thickness were catalogued (Larson et al. (2007) BMC Medical Genetics 8 Suppl. 1 :S5). SNPs identified via candidate gene sequencing (Cardiogenomics) and found to be associated with LV size and ejection fraction with p < 0.0001 in at least one comparison were further catalogued for consideration (NHLBI Program for Genomic Applications HMS. Genomics of cardiovascular development, adaptation, and remodeling. 2010).

MicroRNA Gene Polymorphisms

MicroRNAs are a unique and recently described class of transcribed RNAs that are involved in post-transcriptional regulation of mRNA via targeted degradation. Fully thirty percent of genes are thought to be targets for micro RNA-mediated regulation. Recent studies suggest microRNAs may have a significant role in the stress response to hypertrophy. We included microRNAs in our candidate gene search due to this potentially significant role in altering gene expression. Human microRNAs were catalogued from the most recently available microRNA database, miRbase (microrna.sanger.ac.uk/). Tagged SNPs were identified for genes encoding microRNAs using the chromosome locations as found on miRbase and 1000 base pairs up and downstream.

Final SNP Selection for Custom Array

A target of 1 ,536 SNPs was chosen for the creation of a customized gene array. The overall strategy was to assign approximately one third of the array to intergenic SNPs, one third to genie SNPs and one third to SNPs ascertained from the network analyses and curated lists (FIG. IB). Unless otherwise noted above, SNPs from each category were excluded if the SNP quality scores, representing likelihood of successful genotyping on the GoldenGate platform, were < 0.6 or if there was a representative SNP in linkage disequilibrium, with r² > 0.8, already selected for inclusion. The top 561 SNPs from the weighted informative GWA genie SNP analysis that met these criteria were selected (Group A). The top 540 intergenic SNPs, defined as SNPs lying further than 5000 base pairs up or downstream from the untranslated 3' and 5' regions, were also selected for custom array inclusion (Group D). SNPs from groups A and D that had CEU HapMap minor allele frequencies < 20% were excluded. The top 264 SNPs from the coexpression and semantic network analyses were selected (Group B). The final 171 SNPs were selected from the curated candidate genes and candidate SNPs as well as microRNA polymorphisms (Group C). Furthermore, SNPs selected from group C that had an informative GWA /?-value of greater than 0.5 were excluded.

Custom SNV Array Analysis

Custom array analysis was performed in PLINK using a multivariate additive logistic regression model to measure the association between each SNP and an improvement in ejection fraction of greater than 10% compared to those without documentation of such an improvement. This association was adjusted for age, sex and race. A Manhattan plot of -loglOP was generated using Haploview (FIG. 2A). Values for a quantile-quantile plot were generated using PLINK to evaluate the potential impact of population stratification (FIG. 5). OR values were measured as OR per minor allele of each SNP. The threshold for statistical significance was estimated at the Bonferroni-corrected value of 0.05/1402 = 3.56 x 10^"5.

P-values from the informative GWA and custom gene array were meta- analyzed with a weighted Z-method (Whitlock (2005) J. Evol. Biol. 18: 1368-1373) using Stouffer's Z trend, which considers sample size and effect directions, in MetaP (compute 1. lsrc . duke . edu/ so ftwares/MetaP/metap .php) .

UPenn Association Analysis

Analysis of the replication genotyping in the UPenn cohort was performed in SAS (v9.1) using a multivariate additive logistic regression model to measure association between the rs7767652 genotype and a baseline EF > 35% compared to a baseline EF < 35%. The association was adjusted for age, sex, hypertension, diabetes, renal function (glomerular filtration rate) and body mass index. The odds ratio (OR) was measured as OR per minor allele. The threshold for statistical significance in this single replication association was 0.05.

Functional Assays

Allele-specific gene expression

Oligonucleotides containing the putative regulatory elements for rs7767652-C and rs7767652-T were generated by annealing the following single-stranded sequences (rs7767652-C forward: 5'- AGCT (TTATAATCAAAGAAAG)x4 -3' (SEQ ID NO:2) and rs7767652-C reverse: 5'- GATC (CTTTCTTTGATTATAA)x4 - 3' (SEQ ID NO:3); rs7767652-T forward : 5'-AGCT (TTATAATTAAAGAAA)x4 - 3' (SEQ ID NO:4) and rs7767652-T reverse: 5'-GATC TTTCTTTAATTATAA)x4 - 3' SEQ ID NO:5) at 95 degrees for 10 minutes in annealing buffer and allowed to cool to room temperature. Double-stranded DNA fragments were then directionally subcloned into the multiple cloning site (MCS) of pLuc-MCS vector driven by a minimal promoter (Agilent). Constructs were validated by Sanger sequencing. Empty vector (pLuc-MCS), rs7767652-C or rs7767652-T and Renilla luciferase constructs, as well as empty expression construct (pcDNA3.1-CMV), human TCF4 and human β- catenin, or human TCF4A30 and β-catenin expression constructs were transfected into either C2C12, H9c2 or HEK cell lines using Lipofectamine 2000 (Life Technologies). Cells were maintained in complete Dulbecco's Modified Eagle Medium (DMEM) containing low glucose, sodium pyruvate and L-glutamine (Life Technologies) and supplemented with 10% fetal bovine serum. Media was replaced after 6 hours, and dual- luciferase activity was measured after 24 hours using a SpectraMax L

luminometer (Molecular Devices). Relative luciferase activity (firefly /Renilla luciferase ratio) is expressed as the fold change of the empty vector control (pLuc- MCS) with the empty expression construct (Empty).

HCRTR2 Gene Expression

The HCRTR2 gene expression values were collected from human left ventricular myocardial gene expression datasets in the Gene Expression Omnibus (GEO) database that were performed on the Affymetrix UA133A or UA133 Plus 2.0 arrays. Samples were identified as being control, non-ischemic (dilated)

cardiomyopathy, or ischemic cardiomyopathy. For paired samples from trials with an intervention (such as LVAD placement), only data from the pre-intervention sample was used. Individual datasets were accessed and collated by probe ID, and the combined data set was quantile-normalized using the limma package in R (function normalizeBetweenArrays). Samples without valid calls for HCRTR2 were removed. We then compared the mean gene expression values, reported as ratios of probe fluorescence to background fluorescence, between normal tissue and dilated cardiomyopathic tissue or ischemic cardiomyopathic tissue. The association between expression values and disease status was tested using Student's t-test in R.

Transcription-disrupted and Knockout Mice

Wild-type C57BL/6 14-week-old mice were obtained from Charles River International (Wilmington, MA). HCRTR2 transcription-disrupted (TD) mice were produced by inserting a loxP-flanked gene cassette into intron 1, 137 bp upstream of exon 2 which disrupted production of HCRTR2 in mice with a C57BL/6 background (Mochizuki et al, supra). These mice were provided courtesy of Dr. Thommas Scammell (Harvard Medical School) who previously demonstrated a lack of normal HCRTR2 expression in the TD mice (Mochizuki et al, supra). HCRTR2 knockout mice were created using a targeting vector to replace exon 1 with lac-Z and neo cassettes by homologous recombination in a C57BL/6J background (Chemelli et al, supra). These mice were provided courtesy of Dr. Priyattam Shiromani (Medical University of South Carolina).

Minipump Experiments

Experimental mice were anesthetized with 1.5% isoflurane by a facemask. The Alzet osmotic minipumps (Model 2002 or 2004) from the Alzet Corporation (Cupertino, CA, USA) were implanted subcutaneously in the back of mice, posterior to the scapulae. 14-week-old C57BL/6 mice were treated with either saline or orexin A (50 μg/kg/day) infusions using an Alzet (Cupertino, CA) osmotic minipump model 2004 (200 μΐ reservoir, 28-day release at 0.25 μΐ/hr) for a total of 4 weeks. At week 2, all mice were treated with angiotensin II (1.4 mg/kg/day) plus isoproterenol (15 mg/kg/day) infusion using an additional Alzet osmotic minipump, model 2002 (200 μΐ reservoir, 14-day release at 0.5 μΐ/hr) for two weeks, starting 2 weeks after initiation of the saline or orexin A infusion. Mice had echocardiography performed at baseline and weeks 2, 3 and 4.

Murine Echocardiography

For noninvasive echocardiography, mice were anesthetized with isoflurane via nose cone (maintenance: 1.5% isoflurane mixed with lL/minute 100% medical oxygen to obtain heart rate of 450 ± 50 beats per minute). Mouse body temperature was monitored with a rectal thermometer. The chest was shaved using a Nair lotion hair remover (Church & Dwight Canada, Mississauga, ON, Canada). Electrode gel was put to the four paws and tape them to the ECG electrodes. Two-dimensional, pulse Doppler and tissue Doppler flow echocardiography were performed with Vevo 2100 High-Resolution In Vivo Imaging System (VisualSonics, Toronto, ON, Canada) which has an MS550D microscan transducer with an operating frequency band of 22- 55 MHz. To minimize bias, echo cardiographs were performed blinded to drug infusion status.

Two-dimensional imaging (B-mode) was used to obtain a view of left and right ventricular movement during diastole and systole along the parasternal short axis. M-mode echocardiography was used to obtain measurements of cardiac dimensions including left ventricular internal diastolic and systolic dimensions, which were used to estimate ejection fraction, as well as septal and left ventricular posterior wall thickness. Pulsed-wave Doppler and tissue Doppler echocardiography analyses were used to evaluate diastolic function. Pulsed-wave Doppler images were used to record the trans-mitral Doppler blood flow velocities. Early (E) wave was measured as the peak velocity during filling of the ventricle, and atrial (A) wave was measured as the peak velocity during atrial contraction, at the end of diastole. Tissue Doppler imaging was used to record the peak mitral annular relaxation velocity (Ε') during diastole. All data and images were saved and analyzed using a Vevo 2100 advanced cardiovascular package software with a semi-automated quantification of cardiac function (V 1.4.0 VisualSonics, Toronto, ON, Canada).

Murine Exercise Capacity

Maximum exercise capacity was measured by using a 4-channel rodent treadmill equipped with an electrical stimulus and metabolic analyzers to measure 0₂ and C0₂ gas fractions. (Accuscan Instruments, now Omnitech Electronics, Columbus, OH, USA). Mice were familiarized with running on the treadmill before exercise testing. The exercise protocol consisted of a progressive increase in treadmill speed from 5-22.5 m/minute and incline from 0-15 degrees at 3 minute intervals until the mice reached exhaustion. The point of exhaustion was determined by the observation that mice could not keep up with the pace of the treadmill and had no response to the electrical stimulus. The reference and cage flow rate was maintained at 0.5 1/m. The flow rate, 0₂% and C0₂% of the reference and individual cage were monitored for 30 seconds at 2 minute intervals and recorded with Fusion software (version 3.9).

Exercise capacity was assessed by measuring peak V0₂. Histologic Analysis

Freshly dissected heart tissues were fixed with 4% paraformaldehyde, dehydrated, embedded in paraffin and sectioned at 4 μιη (Stanford Comparative Medicine Histology Lab). After deparaffinization and rehydration, myocardial sections were stained with Masson's trichrome stain. Digital images were acquired at 20x magnification using a Nikon 90i microscope with camera. For degree of fibrosis, blue/red ratios were estimated using RGB histograms from the Masson trichrome- stained sections analyzed in Photoshop CS6 (version 13.0).

Table 1. Baseline demographic and clinical characteristics of Phase I and Phase II heart failure cohorts.

Values are presented as total number (percentage of subjects) or mean ± standard deviation. EF = ejection fraction, NYHA Class III-IV = New York Heart Association heart failure severity class III or IV at baseline, CABG = coronary artery bypass grafting, ACE-I = angiotensin converting enzyme inhibitor, ARB = angiotensin receptor blocker, CCB = calcium channel blocker. Δ EF >20% and Δ EF >10% represent patients whose EF improved by an absolute value greater than 20 or 10 percentage points, respectively. Transplant patients are those with decompensation in their clinical status despite medical intervention and required heart transplantation. *P- values represent significance from the t-test for continuous variables and chi- square test for categorical variables between those with an improvement in EF by 20 percentage points and those who underwent transplant. ** P- value for differences between those whose EF improved by 10 percentage points and those without such an improvement. O

Table 2. Summaiy of association results for ejection fraction response in the Phase I and Phase II heart failure cohorts.

Table 2 (continued)

Top Overall SNPs are the common variants from the customized heart failure gene array with the lowest p-values. The Top Curated SNPs are the variants with the lowest p-values from a subset of variants on the customized array that had been manually pre-selected. The Phase I cohort was age, sex and race-matched. Phase II analyses were adjusted for age, sex and race. Δ EF >20% and Δ EF >10%

represent patients whose EF improved by an absolute value greater than 20 or 10 percentage points, respectively. Transplant patients are those with decompensation in their clinical status despite medical intervention and required heart transplantation. The Penn Heart Failure replication analysis was adjusted for age, sex, hypertension, diabetes, renal function (glomerular filtration rate) and body mass index. MAF = minor allele frequency, OR = odds ratio per minor allele, EF = ejection fraction, Source = group from which the variant was selected for placement on the customized array: D = Top Intergenic SNPs from phase I GWAS, A = Top weighted SNPs from phase I GWAS, B-H b ^:::: Coexpression network analysis gene hub, C-can ^::: candidate SNPs from prior studies, C-Path: variant in gene responsible for important disease pathway, C-GWA ^::: SNPs associated with, relevant phenotypes in prior studies.

Table 3. Weighted (Group A) SNPs: Summary of association results for ejection fraction response in the Phase I cohort weighted by ene expression analysis statistical values (i!-seores), and unweighted Phase II association results.

Minor Phase I Phase I Adjusted Phase II Phase II

SNP Candidate Gene Allele OR p-value rf-score Weight Phase I p OR p-value

Top SNPs Ranked by Phase I Weighted p-values

rsl 1780576 PSD3 C 0.22 0.000461 5.16 3.94 0.000117 1.08 0.6085 rs6444661 C3orf59 T 4.36 0.000436 3.21 3.5 0.000125 1.28 0.1 171 rs6061052 PDRG1 G 4.38 0.000496 4.5 3.92 0.000126 0.86 0.3514 rs6 17833 C6orfl92 C 4.1 1 0.000616 2.42 2.61 0.000236 1.23 0.178 rs6444668 C3orf59 T 3.82 0.000858 3.21 3.5 0.000245 0.91 0.5232 rsl959152 BAZ1A T 3.67 0.001062 9.87 3.95 0.000269 1.07 0.6641 rs588361 RPUSD4 G 0.26 0.001254 4.36 3.91 0.000321 1.1 0.5994 rsl0513287 ASTN2 G 0.25 0.001287 7.09 3.95 0.000326 1.08 0.6182 rs2074267 1LVBL C 0.29 0.001284 4.21 3.89 0.00033 0.8 0.1853 rsl 1731086 KIAA0746 G 6.37 0.000123 0.65 0.35 0.000349 1.12 0.4816 rs602735 RPUSD4 C 0.25 0.001371 4.36 3.91 0.000351 0.9 0.566 rs2900367 DERA C 0.25 0.001396 5.81 3.95 0.000354 1.21 0.2143 rsl049523 CYB5D2 G 5.04 0.001302 3.49 3.67 0.000354 1.23 0.2094 rs6990300 TNKS G 0.23 0.001297 3.29 3.56 0.000365 0.94 0.6812 rs35000 RASGRF2 A 3.48 0.001491 3.39 3.62 0.000412 1.06 0.6942 rs7142040 BAZ1A A 3.24 0.001666 9.87 3.95 0.000422 1.04 0.8038 rs9876221 C3orf59 G 3.41 0.001494 3.21 3.5 0.000427 1.08 0.6291 Minor Phase I Phase I Adjusted Phase II Phase II

SNP Candidate Gene Allele OR p-value rf-score Weight Phase I p OR p-value rs7773338 C6orfl 92 C 0.27 0.001206 2.42 2.61 0.000462 0.83 0.2141 rs4838267 ASTN2 G 0.25 0.001858 7.09 3.95 0.000471 0.99 0.9371

Top SNPs Ranked by Phase II unweighted p-values

rs2527366 GTF21 G 2.35 0.0164 3.33 3.58 0.004576 1.73 0.000648 rsl 546120 MAD2L1 T 0.32 0.01299 8.65 3.95 0.003293 0.61 0.007993 rs7126330 KCNQ1 T 0.36 0.006945 1.5 1.21 0.005723 1.45 0.01217 rsl 023181 PCSK5 G 2.56 0.02158 5.9 3.95 0.00547 1.45 0.01393 rs2598414 SRP68 T 2.36 0.01924 9.45 3.95 0.004877 1.42 0.01472 rsl 439977 ZAK A 2.9 0.003726 9.46 3.95 0.000944 1.42 0.01775 rs4948383 ANK3 G 3.25 0.01075 4.33 3.91 0.002752 1.45 0.02178 rsl 141371 PTPRS C 2.88 0.02253 6.3 3.95 0.00571 1 1.42 0.02187 rs8021281 SAMD4A C 0.33 0.006963 2.09 2.1 1 0.003302 1.46 0.03121 rs2272935 NPAL3 T 3.25 0.01908 9.84 3.95 0.004836 1.4 0.03137 rs709529 BBX C 0.42 0.02494 3.65 3.75 0.006649 0.71 0.03343 rs501014 AGL C 3.25 0.01 153 2.51 2.74 0.004216 0.7 0.0369 rs740058 PTPRS G 0.33 0.007253 6.3 3.95 0.001838 0.68 0.03841 rsl 77421 PCQAP A 3.96 0.003005 1.78 1.62 0.001852 0.73 0.04074 rs51 1 145 RAB22A C 0.39 0.01451 2.43 2.62 0.005529 1.33 0.05265 rsl 3417370 KBTBD10 G 0.41 0.01682 2.59 2.85 0.005893 0.74 0.05795 rs6560488 PCSK5 A 2.55 0.02434 5.9 3.95 0.00617 1.31 0.06574 rsl 858801 ATBF1 A 2.33 0.01974 3.93 3.84 0.005142 0.77 0.07262 rsl 564943 M1CAL2 C 0.4 0.01645 4.12 3.88 0.004241 1.3 0.0766 rs41 121 FAM81 B C 0.32 0.02558 4.4 3.91 0.006538 0.74 0.07832

The Phase I cohort was age, sex and race-matched. Phase I analyses tested associations for EF improvement by an absolute value greater than 20 percentage points using an additive logistic regression model. Phase II analyses tested associations for EF improvement by an absolute value greater than 10 percentage points using an additive logistic regression model adjusted for age, sex and race, d- score was the statistical value from the analysis comparing gene expression in heart failure tissue versus normal heart tissue. Weight is the weighted value based on the d- score that was used to adjust the Phase I /?-value. OR = odds ratio per minor allele, EF = ejection fraction.

Table 4. Intergenic (Group D) SNPs: Summary of association results for ejection fraction response in the Phase I and Phase II heart failure cohorts.

Phase Phase

SNP Candidate Minor Δ EF Transplant OR P Δ EF Δ EF OR P

Gene Allele >20% (n=37) >10% <10%

(n=29) MAF (n=137) (n=454)

MAF MAF MAF

Top SNPs

Ranked by

Phase I p- values

rs421500 RPS12 G 0.57 0.15 7.56 0.000069 0.44 0.41 1.13 0.4396 rs6075459 SLC24A3 A 0.45 0.12 0.000181 1 0.26 0.25 1.08 0.6731 Phase Phase

I II

SNP Candidate Minor Δ EF Transplant OR Δ EF Δ EF OR

Gene Allele >20% (n=37) >10% <10%

(n=29) MAF (n=137) (n=454)

MAF MAF MAF

5.87

rs525153 RPS12 T 0.4 0.12 4.75 0.0001895 0.32 0.26 1.32 0.08576 rs7422405 ATP6V1C2 A 0.64 0.24 5.48 0.0002242 0.46 0.45 1.06 0.7016 rs989547 DAPK2 T 0.21 0.58 0.19 0.0003092 0.41 0.46 0.84 0.2296 rsl 2337993 BNC2 T 0.07 0.35 0.14 0.0003258 0.15 0.21 0.68 0.05336 rs9829987 EPHA3 T 0.57 0.24 4.12 0.0003757 0.44 0.39 1.26 0.1 18 rs6449652 HTR1A C 0.59 0.23 4.75 0.0003876 0.38 0.35 1.15 0.3272 rsl 319025 ANXA2P3 C 0.66 0.34 3.72 0.0004225 0.5 0.44 1.28 0.1085 rs878421 GPR20 C 0.07 0.35 0.14 0.0004342 0.21 0.23 0.96 0.8267 rsl3228415 NXPH1 G 0.43 0.15 4.34 0.0004359 0.18 0.21 0.79 0.2259 rsl2091 19 C14orfl0 A 0.36 0.1 1 4.68 0.0004431 0.23 0.24 0.96 0.7981 rs2243684 C21orf34 A 0.12 0.43 0.18 0.000586 0.33 0.34 0.97 0.8221 rsl0515283 RGMB C 0.36 0.08 6.43 0.0005937 0.18 0.16 1.08 0.6903 rs6948739 CARD 1 1 G 0.22 0.57 0.22 0.00061 1 0.33 0.34 0.91 0.5702 rs2069126 SLC24A3 T 0.19 0.46 0.28 0.00061 13 0.36 0.38 0.97 0.8279 rsl 927745 FLJ10154 A 0.09 0.36 0.16 0.0006245 0.3 0.26 1.25 0.1507 rs7041706 BNC2 G 0.09 0.36 0.16 0.0006245 0.21 0.27 0.72 0.06813 rs6081741 SLC24A3 G 0.09 0.39 0.15 0.000638 0.28 0.3 0.97 0.8727 rs9358552 HDGFL1 G 0.14 0.42 0.22 0.000653 0.43 0.4 1.1 0.5138 Top SNPs

Ranked by

Phase II p- values

rs7767652 HCRTR2 T 0.14 0.36 0.28 0.004039 0.1 1 0.24 0.39 0.0000329 rs4901426 DDHD1 A 0.31 0.63 0.27 0.002487 0.52 0.43 1.54 0.003833 rsl 936602 HTR7 A 0.6 0.3 3.6 0.00203 0.38 0.49 0.67 0.006501 rs7923056 1TGA8 C 0.43 0.2 2.98 0.005915 0.46 0.38 1.45 0.01 143 rs637629 ADRBK2 G 0.53 0.23 3.85 0.003312 0.48 0.4 1.49 0.01315 rs7861 151 ADAMTSLl G 0.26 0.51 0.33 0.005119 0.5 0.41 1.45 0.01375 rsl 927496 ABCC4 G 0.57 0.3 3.12 0.001695 0.4 0.47 0.69 0.01559 rs4682910 ABHD5 T 0.4 0.19 2.82 0.005401 0.22 0.29 0.65 0.01575 rs9836395 CCDC54 G 0.62 0.31 3.72 0.001339 0.36 0.27 1.52 0.01739 rsl 875645 FAF1 A 0.28 0.53 0.34 0.004857 0.37 0.44 0.7 0.01754 rsl35893 LDOC1L C 0.5 0.24 3.1 1 0.002132 0.35 0.44 0.7 0.01832 rs2976271 PPP2R2A G 0.52 0.2 4.21 0.0009619 0.23 0.31 0.68 0.02532 rsl3252558 DOK2 T 0.14 0.38 0.26 0.003558 0.35 0.27 1.44 0.02587 rs6696250 TNRC4 c 0.63 0.28 0.001218 0.52 0.44 1.39 0.02721 Phase Phase

I II

SNP Candidate Minor Δ EF Transplant OR Δ EF Δ EF

Gene Allele >20% (n=37) >10% <10%

(n=29) MAF (n=137) (n=454)

MAF MAF MAF

4.21

rs7759442 C6orfl l 8 G 0.14 0.38 0.26 0.001991 0.35 0.27 1.44 0.02911 rsl323348 ZDHHC21 C 0.36 0.15 3.25 0.005434 0.3 0.38 0.71 0.03222 rs4978396 PTPN3 G 0.38 0.15 3.5 0.003991 0.25 0.31 0.69 0.03273 rsl 1676317 LRRTM4 A 0.33 0.61 0.31 0.005821 0.42 0.48 0.72 0.03437 rsl325376 L1G4 G 0.45 0.2 3.2 0.00391 1 0.33 0.25 1.41 0.03462 rs973330 TYRP1 T 0.6 0.35 2.81 0.004988 0.5 0.44 1.36 0.0398 rs4262135 LOC51334 T 0.57 0.28 3.37 0.003428 0.38 0.45 0.73 0.04131 rs8008076 FOXG1B G 0.45 0.19 3.48 0.003772 0.22 0.27 0.7 0.04182 rs946481 KIAA1026 G 0.22 0.55 0.23 0.0007618 0.48 0.44 1.38 0.04205 rs2197078 C2orG C 0.26 0.53 0.31 0.004852 0.53 0.47 1.36 0.0432 rs4682242 PVRL3 A 0.66 0.38 3.12 0.002765 0.52 0.45 1.38 0.04392 rs8123223 RNF24 T 0.57 0.3 3.12 0.002376 0.4 0.45 0.73 0.0464 rsl2589158 SEL1L T 0.34 0.1 1 4.21 0.002176 0.15 0.2 0.67 0.04685

The Phase I cohort was age, sex and race-matched. Phase II analyses were adjusted for age, sex and race. Δ EF >20% and Δ EF >10% represent patients whose EF improved by an absolute value greater than 20 or 10 percentage points, respectively. Transplant patients are those with decompensation in their clinical status despite medical intervention and required heart transplantation. MAF = minor

frequency, OR = odds ratio per minor allele, EF = ejection fraction.

Table 5. Hub (Group B-hub) SNPs: Summary of the hubs from the coexpression modules most significantly associated with heart failure and the association results of their corresponding SNPs with ejection fraction response in the Phase I and Phase II heart failure cohorts.

Module Average Minor Phase I Phase I Phase II Phase II

Hub Gene Size </-Score SNP Allele OR p-value OR p-value

Top Hubs Ranked by Average Module </-score

SLC25A30 18 8.521 _rs4433707 T 1.18 0.4177 1.15 0.3671

SLC25A30 18 8.521 rs9526070 c 1.55 0.1842 0.97 0.8584

UTP15 24 8.101 rsl3170849 A 1.49 0.3559 0.98 0.8934

UTP15 24 8.101 rs343120 A 0.65 0.3555 0.86 0.3587

SLC17A1 29 7.608 rsl l65196 C 0.58 0.1096 0.8 0.1588

SLC17A1 29 7.608 rsl l65210 A 2.67 0.453 0.93 0.8006

SLC17A1 29 7.608 rs3923 A 0.62 0.2992 1.06 0.6799

NUMBL 11 7.394 rs2250994 A 0.78 0.4583 0.99 0.9707

NUMBL 11 7.394 rs2561537 C 1.19 0.4598 1.1 0.5212 Module Average Minor Phase I Phase I Phase II Phase II

Hub Gene Size </-Score SNP Allele OR p-value OR p-value

LOCI 2475 47 7.386 rs2430949 A 0.79 0.3497 0.8 0.2634

DUS4L 102 7.147 rs6957510 A 3.4 0.2504 1.38 0.249

OR8D1 15 6.88 rs4935881 A 0.65 0.1795 1.07 0.6843

C1D 150 6.588 rs4671871 A 1.59 0.4773 1.17 0.4281

C1D 150 6.588 rs7582426 G 2.09 0.07674 1.34 0.07232

UGT2B7 10 6.516 rs7375178 C 1.83 0.2071 0.84 0.2596

GNG3 17 6.443 rs549312 A 0.57 0.2816 1.24 0.2148

SLC2A9 83 6.439 rs6820230 T 2.07 0.09169 0.75 0.1009

SLC2A9 83 6.439 rs6849736 A 0.42 0.448 1.24 0.3589

IL1RAP 12 6.267 rs2173911 T 1.8 0.1349 0.94 0.7164

IL1RAP 12 6.267 rs9682599 C 3.23 0.02459 1.39 0.1016

TRIM29 27 6.196 rs4936509 c 0.49 0.03817 1 0.9965

CMTM2 12 6.168 rs2290182 G 1.52 0.349 0.87 0.5073

IQCF1 122 6.161 _rs4234645 C 1.85 0.1913 1.03 0.8563

SLC37A3 13 6.147 rsl2703799 C 1.43 0.4312 0.87 0.3291

Top Hubs Ranked by Phase II /^-values.

MRO 65 5.492 rs3813089 T 1.36 0.4649 1.6 0.00748

MGC61571 32 5.787 rs7611318 G 1.53 0.354 1.41 0.02368

SLC16A9 20 5.582 rs7094971 G 1.81 0.3362 1.51 0.03891

FLJ32028 10 5.127 _rs6814401 A 2.08 0.04475 1.37 0.0438

DTX1 14 5.512 rs2701623 G 2.04 0.05177 0.74 0.06716

FLJ32028 10 5.127 rsl0034000 T 3.14 0.00502 1.33 0.06948

C1D 150 6.588 rs7582426 G 2.09 0.07674 1.34 0.07232

FBXL21 18 6.018 rs31551 T 0.42 0.03242 1.31 0.0794

MRO 65 5.492 rs2255672 c 0.66 0.4313 1.32 0.09185

SLC2A9 83 6.439 rs6820230 T 2.07 0.09169 0.75 0.1009

IL1RAP 12 6.267 rs9682599 c 3.23 0.02459 1.39 0.1016

MGC61571 32 5.787 rsl3085075 T 0.68 0.2824 1.28 0.1056

HNMT 11 5.478 rs6430676 c 1.76 0.3206 1.37 0.1283

C14orf29 26 4.977 rsl953879 T 1.89 0.4307 1.38 0.1353

SLC37A3 13 6.147 rs6947516 A 0.17 0.08895 1.38 0.1525

DTX1 14 5.512 rsl732803 T 0.67 0.2484 1.24 0.1544

SLC17A1 29 7.608 rsl l65196 c 0.58 0.1096 0.8 0.1588

EME1 32 5.351 rsl7714829 c 0.39 0.2352 0.66 0.167

ZNF238 5.165 rsl 1590416 G 0.53 0.1082 0.76 0.1869

SOCS5 5.802 rs4953410 G 1.64 0.2349 1.25 0.1959

The (i-score represents the statistical value from the analysis comparing gene expression in heart failure tissue versus normal heart tissue. The average <i-score is the absolute statistical value averaged across all members of the network module that contains the noted gene as the hub. Phase I analyses tested associations for EF improvement by an absolute value greater than 20 percentage points using an additive logistic regression model. The Phase I cohort was age, sex and race-matched. Phase II analyses tested associations for EF improvement by an absolute value greater than 10 percentage points using an additive logistic regression model adjusted for age, sex and race. OR = odds ratio per minor allele, EF = ejection fraction. Table 6. Literature Network (Group B-Lit) SNPs: Summary of the nexi from the literature networks most significantly associated with heart failure and the association results of their corresponding SNPs with ejection fraction response in the Phase I and Phase II heart failure cohorts.

Module Average Minor Phase I Phase I Phase II Phase II

Hub Gene Size </-Score SNP Allele OR p-value OR p-value

Top Hubs Ranked by Average Module </-score

A2M 8 11.094 rs226389 G 0.63 0.2145 0.78 0.1041

ITIH5 6 8.534 rsl0508332 G 0.57 0.1378 0.97 0.8269

ITIH5 6 8.534 rsl931901 A 2.19 0.08623 0.8 0.1551

ITIH5 6 8.534 rsl931902 A 0.54 0.1586 0.99 0.9338

ITIH5 6 8.534 rs7097970 C 0.6 0.1177 1.06 0.7448

ITIH1 7 8.111 rs2710323 G 0.79 0.4141 1.06 0.6725

GSPT1 7 8.047 rs33660 G 1.46 0.3279 1.04 0.7995

SLC1A7 8 7.319 rsll206093 C 0.71 0.3918 0.97 0.8894

SLC1A7 8 7.319 rsl320594 A 0.58 0.1362 0.89 0.4288

SLC1A7 8 7.319 rsl679933 G 2.53 0.01498 1.06 0.7089

SLC1A7 8 7.319 rsl769299 T 0.64 0.3578 1.03 0.8684

KCNIP2 8 7.086 rsl0883689 A 1.52 0.3043 1.13 0.4596

NODAL 10 7.085 rsl0762381 C 0.57 0.2012 0.7 0.01575

ATP7B 8 6.939 rsl801249 T 1.47 0.342 1.24 0.1611

ATP7B 8 6.939 rs754610 A 0.57 0.146 0.8 0.1206

COMMD6 6 6.827 rsl2869081 A 1.32 0.4337 0.89 0.5771

COMMD6 6 6.827 rs9543956 T 1.74 0.1207 1.21 0.2243

KCNK9 9 6.825 rslll66920 T 2.47 0.03302 1.04 0.8211

KCNK9 9 6.825 rs759656 c 2.14 0.1119 0.94 0.6957

KCNK9 9 6.825 rs882555 c 1.57 0.1725 1.09 0.6032

KCNK9 9 6.825 rs885724 c 1.92 0.1484 0.98 0.9219

KCNK9 9 6.825 rs888346 A 1.89 0.1401 0.81 0.163

PRKAG2 7 6.545 rs2727528 C 1.71 0.03734 1.07 0.6765

PRKAG2 7 6.545 rs2727541 C 1.87 0.05437 0.85 0.2853

PRKAG2 7 6.545 rs7801616 T 2.2 0.0572 1.03 0.8472

Top Hubs Ranked by Phase II p- ^■values.

RASSF6 1 5.712 rsl2507775 G 1.4 0.3837 1.53 0.01209

NODAL 10 7.085 rsl0762381 C 0.57 0.2012 0.7 0.01575

CHRM3 7 3.61 rs665159 C 1.81 0.1042 1.41 0.01675

BTAF1 9 4.106 rsl2763600 G 0.72 0.2892 1.45 0.01837

LIN7A 19 5.264 rs7969498 T 0.62 0.2028 1.42 0.02552

RAB22A 6 5.274 rs471661 c 0.53 0.1018 1.4 0.02745

MAK 8 4.585 rs2009118 T 2.15 0.1506 1.39 0.04902

PRPF4B 8 4.649 rsl2198921 c 0.62 0.1276 0.74 0.05348

BTC 11 4.873 rs4352548 c 1.77 0.18 1.36 0.06907

ANTXR1 6 6.262 rs7584948 G 1.36 0.4811 1.35 0.083

A2M 8 11.094 rs226389 G 0.63 0.2145 0.78 0.1041

NOBOX 8 5.939 rsl208208 C 1.83 0.06221 1.27 0.1071

ATP7B 8 6.939 rs754610 A 0.57 0.146 0.8 0.1206

KIR3DL2 11 6.033 rsl654644 G 1.27 0.409 1.27 0.1233

COG2 7 5.031 rs2296799 C 0.19 0.01348 1.32 0.1277

RRAGA 6 6.451 rsl0125646 G 0.76 0.4111 0.8 0.1368

EPB41L1 6 4.533 rs2746105 A 0.73 0.4998 0.77 0.1387

CYP11B1 9 6.302 rs6410 A 0.68 0.276 1.25 0.1438

ITIH5 6 8.534 rsl931901 A 2.19 0.08623 0.8 0.1551 Module Average Minor Phase I Phase I Phase II Phase II

Hub Gene Size </-Score SNP Allele OR p-value OR p-value

ANTXR1 6 6.262 rs4854550 T 0.7 0.4568 0.81 0.1572

The (i-score represents the statistical value from the analysis comparing gene expression in heart failure tissue versus normal heart tissue. The average <i-score is the absolute statistical value averaged across all members of the network module that contains the noted gene as the nexus. Phase I analyses tested associations for EF improvement by an absolute value greater than 20 percentage points using an additive logistic regression model. The Phase I cohort was age, sex and race-matched. Phase II analyses tested associations for EF improvement by an absolute value greater than 10 percentage points using an additive logistic regression model adjusted for age, sex and race. OR = odds ratio per minor allele, EF = ejection fraction.

Table 7. Replication of SNPs from other GWAS: Summary of association results for ejection fraction response in the Phase I and Phase II heart failure cohorts.

Gene Minor Phase I Phase I Phase II Phase II

SNP Symbol Phenotype Allele OR p-value OR p-value rsl2623467 NRXN1 Nicotine Dependence T 1.61 0.0267 rs7903146 TCF7L2 Type II DM T 1.43 0.3637 1.38 0.04691 rsl0511311 CD200 Atrial Fibrillation T 1.34 0.04834 rs9939609 FTO Type II DM, BMI A 1.33 0.07296 rsl3277254 CHRNB3 Nicotine Dependence G 1.34 0.08167 rs2383206 CDKN2B CAD A 1.24 0.1385 rs4506565 TCF7L2 Type II DM T 1.24 0.1606

Myocardial

rsl0757278 CDKN2B Infarction G 0.82 0.185 rs2383207 CDKN2B AAA A 0.75 0.3797 1.2 0.2178 rs662 PON1 CAD G 0.7 0.4543 1.2 0.2519 rsl 801282 PPARG Type II DM G 1.26 0.2563

LV Systolic

rsl0504543 KCNB2 Dimension A 1.91 0.1258 0.84 0.2757 rs9465871 CDKAL1 Type II DM C 1.18 0.3274 rs8055236 CDH13 CAD T 1.93 0.1739 1.18 0.3462 rs383830 APC CAD T 1.18 0.3771 rsl333049 CDKN2B CAD G 1.12 0.4386 rs 10490162 NRXN1 Nicotine Dependence G 1.07 0.8893 1.15 0.487 r_s4142041 CTNNA3 Nicotine Dependence G 1.09 0.5405 rsl0501920 CNTN5 Atrial Fibrillation G 0.89 0.5669 rsl0515869 MAT2B Heart Failure C 1.68 0.2101 1.07 0.6636 rs2836823 FLJ45139 Nicotine Dependence T 0.94 0.7026

LV Diastolic

rsl379659 SLIT2 Dimension G 1.58 0.3641 0.95 0.8026 rs5215 KCNJ11 Type II DM C 1.34 0.3122 1.04 0.8212 rsl7672135 FMN2 CAD C 0.96 0.8661 rs2221299 NRXN3 Nicotine Dependence A 1.06 0.8737 rs3807375 KCNH2 QRS Duration A (λ97 1.01 0.9458

Phase I analyses tested associations for EF improvement by an absolute value greater than 20 percentage points using an additive logistic regression model. The Phase I cohort was age, sex and race-matched. Phase II analyses tested associations for EF improvement by an absolute value greater than 10 percentage points using an additive logistic regression model adjusted for age, sex and race. OR = odds ratio per minor allele, EF = ejection fraction, DM = diabetes mellitus, BMI = body mass index, AAA

= abdominal aortic aneurysm, CAD = coronary artery disease, LV = left ventricular.

Table 8.

SNP Allele Program Transcription factor Sequence Strand Score rs7767652 C TRANSFAC LEF-1 tCAAAG 100

PMX1 TAATCa 100 TCF-4 aaTCAAAga 97.6

Matlnspector Sex determining region Y tccacATTAtaatcaaagaaag 97.9

TCF7LEF-1 ttataatCAAAgaaagt 95.9

Con e-rod homeobox-containing cattaTAATcaaagaaa 95.7 transcription factor (CRX)

PROMO LEF-1 AATCAAAG 98.6

TCF-4 TAATCAAAGA 97.2 TCF-4E ATCAAAG 96.8 GR CAAAGAA 96.2 SRY TAATCAAAG 94.9

JASPAR SRY ATAATCAAA 84.1

Nkx2.5 ATAATCA 82.7 TFSearch SRY AATCAAA 90.9

CdxA TATAATC 90.7

Octl ATTAT AATCAAAG 86.2

TRANSFAC PMX1 TAATTa 100

PMX1 tAATTA 100

DLX5 AATTAa 100

LHX3 ATTAAa 100

SATB1 ataatTAAAGaaagta 97.0

HOXA5 attaTAATTaaagaaa 95.3 Matlnspector NK6 homeobox I ttctTTAAttataat 95.6

DLX 1, 2, and 5 binding sites tccacattatAATTaaaga 100

Sex determining region Y tccacATTAtaattaaagaaagtac 98.1

Homeobox C8/Hox-3alpha ttctttaATTAtaatgtgg 95.5

LIM homeobox transcription tccacattatAATTaaagaaagt 90.0 factor 1. alpha

PROMO GR-beta AATTA 99.2

GR-beta TAATT 99.2 JASPAR SRY TTAAAGAAA 80.1

HOXA5 CTTTAATT 95.2

Pdx l TTAATT 94.0

Nkx2.5 TT AATTA 92.0

TFSearch CdxA AATTAAA 98.6

CdxA TATAATT 90.7 While the preferred embodiments of the invention have been illustrated and described, it will be appreciated that various changes can be made therein without departing from the spirit and scope of the invention.

Claims

Claims What is claimed is:

1. A method of treating an individual for heart failure, the method comprising administering a therapeutically effective amount of an agonist of hypocretin receptor

2 (HCRTR2) to the individual.

2. The method of claim 1, wherein the agonist of HCRTR2 is selected from the group consisting of orexin-A, orexin-B, [Alal l,D-Leul5]-orexin B, and SB-668875.

3. The method of claim 1, wherein the agonist of HCRTR2 is an orexin mimic.

4. The method of claim 1, wherein the agonist of HCRTR2 is recombinantly or synthetically produced.

5. The method of claim 1, wherein the agonist of HCRTR2 is administered according to a daily dosing regimen.

6. The method of claim 1, wherein the agonist of HCRTR2 is administered according to a twice-a-week or three-times-a-week dosing regimen.

7. The method of claim 1, wherein multiple cycles of the method of treatment are administered to the individual for a time period sufficient to effect at least a partial recovery of heart function.

8. The method of claim 7, wherein the time period is at least 6 months.

9. The method of claim 8, wherein the time period is at least 12 months.

10. The method of claim 7, wherein the individual shows at least a 10% improvement in left ventricular (LV) function.

11. The method of claim 10, wherein the individual shows at least a 20% improvement in left ventricular (LV) function.

12. The method of claim 7, wherein the treatment increases the expected survival time of the individual.

13. The method of claim 1, wherein the agonist is administered intravenously, intra-arterially, subcutaneously, intracerebroventricularly, intrathecally, or intranasally.

14. The method of claim 1, wherein the agonist is administered locally into the brain or spinal cord of the individual.

15. The method of claim 14, wherein the agonist is administered locally into the cisterna magna or the rostral ventrolateral medulla.

16. The method of claim 1, wherein arterial blood pressure or heart rate is increased in the individual after administering the agonist of HCRTR2.

17. The method of claim 1, further comprising treatment with one or more other agents for treating heart failure.

18. The method of claim 17, wherein the one or more other agents for treating heart failure are selected from the group consisting of an angiotensin converting enzyme (ACE) inhibitor, an angiotensin II receptor blocker (ARB), a beta-blocker, digoxin, a diuretic, a blood vessel dilator, potassium, magnesium, an aldactone inhibitor, a calcium channel blocker, and an inotrope.

19. The method of claim 1, wherein the individual is human.

20. The method of claim 1, wherein the individual is diabetic or hypertensive.

21. A method for enhancing cardiovascular function of an individual with heart failure, the method comprising administering an effective amount of an agonist of hypocretin receptor 2 (HCRTR2) to the individual.

22. The method of claim 21, wherein the agonist of HCRTR2 is selected from the group consisting of orexin-A, orexin-B, [Alal l,D-Leul5]-orexin B, and SB-668875.

23. The method of claim 21, wherein the agonist of HCRTR2 is an orexin mimic.

24. The method of claim 21, wherein the agonist of HCRTR2 is recombinantly or synthetically produced.

25. The method of claim 21, wherein the agonist is administered

intracerebroventricularly, intrathecally, or intranasally.

26. The method of claim 21, wherein the agonist is administered locally into the brain or spinal cord of the individual.

27. The method of claim 26, wherein the agonist is administered locally into the cisterna magna or the rostral ventrolateral medulla.

28. The method of claim 21, wherein the individual is administered the agonist of HCRTR2 multiple times over a period of at least 6 months.

29. The method of claim 28, wherein the individual is administered the agonist of HCRTR2 multiple times over a period of at least 12 months.

30. The method of claim 21, wherein the individual shows at least a 10% improvement in left ventricular (LV) function.

31. The method of claim 30, wherein the individual shows at least a 20% improvement in left ventricular (LV) function.

32. The method of claim 21, wherein arterial blood pressure or heart rate is increased in the individual after administering the agonist of HCRTR2.

33. A diagnostic method for predicting therapeutic responsiveness to treatment with a beta blocker or an angiotensin converting enzyme (ACE) inhibitor for an individual with heart failure, the method comprising determining which allele is present at single nucleotide polymorphism rs7767652 in the individual, wherein the presence of the major allele at the single nucleotide polymorphism rs7767652 indicates that the individual will be more responsive to treatment for heart failure with a beta blocker or an ACE inhibitor than an individual with the minor allele at the single nucleotide polymorphism rs7767652.

34. The method of claim 33, wherein the allele that is present at the single nucleotide polymorphism rs7767652 is determined by a method selected from the group consisting of dynamic allele-specific hybridization (DASH), microarray analysis, Tetra-primer ARMS-PCR, a TaqMan 5 '-nuclease assay; an Invader assay with Flap endonuclease (FEN), a Serial Invasive Signal Amplification Reaction (SISAR), an oligonucleotide ligase assay, restriction fragment length polymorphism (RFLP), single-strand conformation polymorphism, temperature gradient gel electrophoresis (TGGE), denaturing high performance liquid chromatography (DHPLC), sequencing, and immunoassay.

35. The method of claim 33, wherein the allele that is present at the single nucleotide polymorphism rs7767652 is determined using at least one set of allele- specific primers selected from the group consisting of:

a) a set of allele-specific primers comprising a forward primer comprising the sequence of SEQ ID NO:2 and a reverse primer comprising the sequence of SEQ ID NO:3; and

b) a set of allele-specific primers comprising a forward primer comprising the sequence of SEQ ID NO:4 and a reverse primer comprising the sequence of

SEQ ID NO:5.

36. A kit comprising at least one agent for determining which allele is present at the single nucleotide polymorphism (SNP) rs7767652.

37. The kit of claim 36, wherein at least one agent is an allele-specific probe that hybridizes to a nucleic acid comprising SNP rs7767652.

38. The kit of claim 36, comprising a microarray comprising the allele-specific probe.

39. The kit of claim 36, further comprising at least one set of allele-specific primers for determining which allele is present at SNP rs7767652.

40. The kit of claim 39, wherein at least one set of allele-specific primers is selected from the group consisting of:

a) a set of allele-specific primers comprising a forward primer comprising the sequence of SEQ ID NO:2 and a reverse primer comprising the sequence of SEQ ID NO:3; and

b) a set of allele-specific primers comprising a forward primer comprising the sequence of SEQ ID NO:4 and a reverse primer comprising the sequence of SEQ ID NO:5.

41. The kit of claim 36, further comprising reagents for performing dynamic allele-specific hybridization (DASH), Tetra-primer ARMS-PCR, a TaqMan 5 '- nuclease assay; an Invader assay with Flap endonuclease (FEN), a Serial Invasive Signal Amplification Reaction (SISAR), an oligonucleotide ligase assay, restriction fragment length polymorphism (RFLP), single-strand conformation polymorphism, temperature gradient gel electrophoresis (TGGE), denaturing high performance liquid chromatography (DHPLC), sequencing, or an immunoassay.

42. A composition comprising an agonist of hypocretin receptor 2 for use in the treatment of heart failure.

43. The composition of claim 42, wherein the agonist of hypocretin receptor 2 is orexin-A, orexin-B, [Alal l,D-Leul5]-orexin B, or SB-668875.

44. A composition comprising an agonist of hypocretin receptor 2 for use in enhancing cardiovascular function.

45. The composition of claim 44, wherein the agonist of hypocretin receptor 2 is orexin-A, orexin-B, [Alal l,D-Leul5]-orexin B, or SB-668875.

46. Use of an agonist of hypocretin receptor 2 (HCRTR2) in the manufacture of a medicament for treating an individual for heart failure.