Peptidomic Analysis Reveals Seasonal Neuropeptide and Peptide Hormone Changes in the Hypothalamus and Pituitary of a Hibernating Mammal

Abstract

During the winter, hibernating mammals undergo extreme changes in physiology, which allow them to survive several months without access to food. These animals enter a state of torpor, which is characterized by decreased metabolism, near-freezing body temperatures, and a dramatically reduced heart rate. The neurochemical basis of this regulation is largely unknown. Based on prior evidence suggesting that the peptide-rich hypothalamus plays critical roles in hibernation, we hypothesized that changes in specific cell–cell signaling peptides (neuropeptides and peptide hormones) underlie physiological changes during torpor/arousal cycles. To test this hypothesis, we used a mass spectrometry-based peptidomics approach to examine seasonal changes of endogenous peptides that occur in the hypothalamus and pituitary of a model hibernating mammal, the thirteen-lined ground squirrel (Ictidomys tridecemlineatus). In the pituitary, we observed changes in several distinct peptide hormones as animals prepared for torpor in October, exited torpor in March, and progressed from spring (March) to fall (August). In the hypothalamus, we observed an overall increase in neuropeptides in October (pre-torpor), a decrease as the animal entered torpor, and an increase in a subset of neuropeptides during normothermic interbout arousals. Notable changes were observed for feeding regulatory peptides, opioid peptides, and several peptides without well-established functions. Overall, our study provides critical insight into changes in endogenous peptides in the hypothalamus and pituitary during mammalian hibernation that were not available from transcriptomic measurements. Understanding the molecular basis of the hibernation phenotype may pave the way for future efforts to employ hibernation-like strategies for organ preservation, combating obesity, and treatment of stroke.

Project and data background

We hypothesized that changes in cell–cell signaling peptides from the hypothalamus and pituitary may be responsible for regulating critical physiological changes during the hibernation season. To test this hypothesis, we applied a non-targeted liquid chromatography–mass spectrometry (LC–MS) and LC-tandem mass spectrometry (LC–MS/MS) “peptidomics” approach to identify and quantify endogenous neuropeptides and peptide hormones as a function of season and hibernation stage in the thirteen-lined ground squirrel. LC–MS and LC–MS/MS-based peptidomics allow the detection, identification, and quantification of peptides without preselection, allowing for findings not possible with targeted antibody-based methods.

Thirteen-lined ground squirrel hibernation cycle. A. Active thirteen-lined ground squirrel. B. Torpid thirteen-lined ground squirrel. C. Overview of the annual hibernation cycle and timeline for tissue sampling. The central circle represents months of the year progressing clockwise, with the general active period in white and the hibernation states of torpor and IBA in gray, as indicated in the center. Colored outer regions and text boxes indicate time points for sample collection for this study. For time points in pre-hibernation and hibernation seasons, the measured range of internal body temperatures is indicated in parentheses.

Pituitary

Label-free peptidomics

Total number of samples after LC-MS: 37

• August: 9

• October: 9

• IBA: 6

• Torpor: 5

• March: 8

Total number of feature provided: 124

Three step data processing:

• Batch correction

• Log2 transformation

• EigenMS normalization

Batch correction

Batch correction using package statTarget evaluates the missing values and a feature will be kept if it has non-zero value for at least 50% of samples (statTarget default is 80%) in any one group (remove 31 features). It then imputes missing values for the present features and QC-based signal correction. At last, it removes features CV% > 50% (default, removes 3 features).

Code

pheno <- "CNS_meta4samples_pituitary.csv"
dfile = "CNS_data4samples_pituitary.csv"
#labels <- read.csv("CNS_label4samples_pituitary.csv", row.names = "SampleName")

statTarget::shiftCor(pheno, dfile,  QCspan = 0.25, Frule = 0.5,
                     degree = 2,imputeM = "KNN", ntree=500, coCV = 50)

statTarget: Signal Correction Start... Time: Mon Nov 20 19:24:07 2023 

* Step 1: Data File Checking Start..., Time:  Mon Nov 20 19:24:07 2023 

 Data Link 
 metaFile: CNS_meta4samples_pituitary.csv 
 profileFile: CNS_data4samples_pituitary.csv 

 44 Meta Samples vs 44 Profile samples
 The Meta samples list (*NA, missing data from the Profile File)

 [1] "20220302_QC1" "20220302_QC2" "Pituitary_53" "Pituitary_31" "Pituitary_38"
 [6] "Pituitary_63" "Pituitary_55" "Pituitary_51" "Pituitary_50" "Pituitary_34"
[11] "Pituitary_59" "Pituitary_43" "20220302_QC3" "Pituitary_66" "Pituitary_26"
[16] "Pituitary_19" "Pituitary_24" "Pituitary_41" "Pituitary_57" "Pituitary_47"
[21] "Pituitary_52" "Pituitary_23" "20220303_QC1" "Pituitary_58" "Pituitary_44"
[26] "Pituitary_37" "Pituitary_39" "Pituitary_46" "Pituitary_65" "Pituitary_56"
[31] "Pituitary_48" "Pituitary_35" "Pituitary_61" "20220304_QC1" "Pituitary_49"
[36] "Pituitary_60" "Pituitary_54" "Pituitary_45" "Pituitary_42" "Pituitary_62"
[41] "Pituitary_25" "Pituitary_64" "20220305_QC1" "20220305_QC2"

 Meta-information: 

    Class No.
1  August   9
2     IBA   6
3   March   8
4 October   9
5      QC   7
6  Torpor   5
  Batch No.
1     1  44

 Metabolic profile information: 

               no.
QC and samples  44
Metabolites    128

* Step 2: Evaluation of Missing Value... 

 The number of missing value before QC based signal correction:  732
 The number of filtered variables using the modified  50 % rule :   31 

* Step 3: Imputation start... 

 The imputation method was set at 'KNN'
 The number of missing value after imputation:  0
 Imputation Finished! 
 
* Step 4: QC-based Signal Correction Start... Time:  Mon Nov 20 19:24:07 2023 

 The Signal Correction method was set at QC-RFSC 

  |                                                                            
  |                                                                      |   0%
  |                                                                            
  |=                                                                     |   1%
  |                                                                            
  |=                                                                     |   2%
  |                                                                            
  |==                                                                    |   3%
  |                                                                            
  |===                                                                   |   4%
  |                                                                            
  |====                                                                  |   5%
  |                                                                            
  |====                                                                  |   6%
  |                                                                            
  |=====                                                                 |   7%
  |                                                                            
  |======                                                                |   8%
  |                                                                            
  |=======                                                               |   9%
  |                                                                            
  |=======                                                               |  10%
  |                                                                            
  |========                                                              |  11%
  |                                                                            
  |=========                                                             |  12%
  |                                                                            
  |=========                                                             |  14%
  |                                                                            
  |==========                                                            |  15%
  |                                                                            
  |===========                                                           |  16%
  |                                                                            
  |============                                                          |  17%
  |                                                                            
  |============                                                          |  18%
  |                                                                            
  |=============                                                         |  19%
  |                                                                            
  |==============                                                        |  20%
  |                                                                            
  |===============                                                       |  21%
  |                                                                            
  |===============                                                       |  22%
  |                                                                            
  |================                                                      |  23%
  |                                                                            
  |=================                                                     |  24%
  |                                                                            
  |==================                                                    |  25%
  |                                                                            
  |==================                                                    |  26%
  |                                                                            
  |===================                                                   |  27%
  |                                                                            
  |====================                                                  |  28%
  |                                                                            
  |====================                                                  |  29%
  |                                                                            
  |=====================                                                 |  30%
  |                                                                            
  |======================                                                |  31%
  |                                                                            
  |=======================                                               |  32%
  |                                                                            
  |=======================                                               |  33%
  |                                                                            
  |========================                                              |  34%
  |                                                                            
  |=========================                                             |  35%
  |                                                                            
  |==========================                                            |  36%
  |                                                                            
  |==========================                                            |  38%
  |                                                                            
  |===========================                                           |  39%
  |                                                                            
  |============================                                          |  40%
  |                                                                            
  |============================                                          |  41%
  |                                                                            
  |=============================                                         |  42%
  |                                                                            
  |==============================                                        |  43%
  |                                                                            
  |===============================                                       |  44%
  |                                                                            
  |===============================                                       |  45%
  |                                                                            
  |================================                                      |  46%
  |                                                                            
  |=================================                                     |  47%
  |                                                                            
  |==================================                                    |  48%
  |                                                                            
  |==================================                                    |  49%
  |                                                                            
  |===================================                                   |  50%
  |                                                                            
  |====================================                                  |  51%
  |                                                                            
  |====================================                                  |  52%
  |                                                                            
  |=====================================                                 |  53%
  |                                                                            
  |======================================                                |  54%
  |                                                                            
  |=======================================                               |  55%
  |                                                                            
  |=======================================                               |  56%
  |                                                                            
  |========================================                              |  57%
  |                                                                            
  |=========================================                             |  58%
  |                                                                            
  |==========================================                            |  59%
  |                                                                            
  |==========================================                            |  60%
  |                                                                            
  |===========================================                           |  61%
  |                                                                            
  |============================================                          |  62%
  |                                                                            
  |============================================                          |  64%
  |                                                                            
  |=============================================                         |  65%
  |                                                                            
  |==============================================                        |  66%
  |                                                                            
  |===============================================                       |  67%
  |                                                                            
  |===============================================                       |  68%
  |                                                                            
  |================================================                      |  69%
  |                                                                            
  |=================================================                     |  70%
  |                                                                            
  |==================================================                    |  71%
  |                                                                            
  |==================================================                    |  72%
  |                                                                            
  |===================================================                   |  73%
  |                                                                            
  |====================================================                  |  74%
  |                                                                            
  |====================================================                  |  75%
  |                                                                            
  |=====================================================                 |  76%
  |                                                                            
  |======================================================                |  77%
  |                                                                            
  |=======================================================               |  78%
  |                                                                            
  |=======================================================               |  79%
  |                                                                            
  |========================================================              |  80%
  |                                                                            
  |=========================================================             |  81%
  |                                                                            
  |==========================================================            |  82%
  |                                                                            
  |==========================================================            |  83%
  |                                                                            
  |===========================================================           |  84%
  |                                                                            
  |============================================================          |  85%
  |                                                                            
  |=============================================================         |  86%
  |                                                                            
  |=============================================================         |  88%
  |                                                                            
  |==============================================================        |  89%
  |                                                                            
  |===============================================================       |  90%
  |                                                                            
  |===============================================================       |  91%
  |                                                                            
  |================================================================      |  92%
  |                                                                            
  |=================================================================     |  93%
  |                                                                            
  |==================================================================    |  94%
  |                                                                            
  |==================================================================    |  95%
  |                                                                            
  |===================================================================   |  96%
  |                                                                            
  |====================================================================  |  97%
  |                                                                            
  |===================================================================== |  98%
  |                                                                            
  |===================================================================== |  99%
  |                                                                            
  |======================================================================| 100%

 Calculation of CV distribution of raw peaks (QC)...

        CV<5%   CV<10%   CV<15%  CV<20%   CV<25%   CV<30%   CV<35%  CV<40%
Batch_1     0 3.092784 22.68041 37.1134 54.63918 60.82474 69.07216 74.2268
Total       0 3.092784 22.68041 37.1134 54.63918 60.82474 69.07216 74.2268
          CV<45%   CV<50%   CV<55%   CV<60%   CV<65%   CV<70% CV<75% CV<80%
Batch_1 77.31959 87.62887 91.75258 94.84536 95.87629 98.96907    100    100
Total   77.31959 87.62887 91.75258 94.84536 95.87629 98.96907    100    100
        CV<85% CV<90% CV<95% CV<100%
Batch_1    100    100    100     100
Total      100    100    100     100

 Calculation of CV distribution of corrected peaks (QC)...

           CV<5%   CV<10%   CV<15%   CV<20%   CV<25%  CV<30%   CV<35%   CV<40%
Batch_1 1.030928 20.61856 41.23711 59.79381 68.04124 74.2268 82.47423 90.72165
Total   1.030928 20.61856 41.23711 59.79381 68.04124 74.2268 82.47423 90.72165
          CV<45%   CV<50%   CV<55% CV<60% CV<65% CV<70% CV<75% CV<80% CV<85%
Batch_1 92.78351 96.90722 98.96907    100    100    100    100    100    100
Total   92.78351 96.90722 98.96907    100    100    100    100    100    100
        CV<90% CV<95% CV<100%
Batch_1    100    100     100
Total      100    100     100

* Step 5:  Removal of the features (CV% > 50%) Mon Nov 20 19:24:08 2023 

 No. of removed features: 3
 Feature name: LAGAPEPAEP KRLGVLFNPYYDPLQWKSSRFE TPEAQAAFQKVVAGVANALSHKYH 

 Output Link: C:/Users/qiuha/OneDrive/projects/MaxHub/projects/2023-07-03-CNS-groundsquirrel 

 Correction Finished! Time:  Mon Nov 20 19:24:08 2023 

 #################################### 
 # Software Version: statTarget 2.0 + # 
 # Email: luanhm@sustech.edu.cn# 
 #################################### 

Code

fh = read.csv(file = "statTarget/shiftCor/After_shiftCor/shift_sample_cor.csv",header = TRUE, row.names = 1, stringsAsFactors = FALSE, sep = ",") %>%
        t(.) %>%
        as.data.frame()
labels_d1 <- as.matrix(fh[1,])
N_feature = nrow(fh[-1,])

# change row annotation for heatmap
fh_mod = read.csv(file = "statTarget/shiftCor/After_shiftCor/shift_sample_cor_mod_pituitary.csv",header = TRUE, row.names = 1, stringsAsFactors = FALSE, sep = ",") %>%
        t(.) %>%
        as.data.frame()

feature_name = as.matrix(fh_mod[-1, "sample"]) %>%
        `rownames<-`(rownames(fh_mod)[-1]) %>%
        `colnames<-`("name") %>%
        as.data.frame() %>%
        rownames_to_column(., var = "variable")

• Features = 94

Data preview (after batch correction)

Code

#rownames(fh)[-1] == feature_name$variable
rownames(fh)[-1] = feature_name$name
# display table
fh %>% 
        #mutate(across(-variable, round, digit=2)) %>%
        knitr::kable(., digits = 3, "html") %>%
        kableExtra::kable_styling(bootstrap_options = c("striped", "hover", "condensed", "responsive")) %>%
        kableExtra::scroll_box(width = "100%", height = "500px")

	Pituitary_53	Pituitary_31	Pituitary_38	Pituitary_63	Pituitary_55	Pituitary_51	Pituitary_50	Pituitary_34	Pituitary_59	Pituitary_43	Pituitary_66	Pituitary_26	Pituitary_19	Pituitary_24	Pituitary_41	Pituitary_57	Pituitary_47	Pituitary_52	Pituitary_23	Pituitary_58	Pituitary_44	Pituitary_37	Pituitary_39	Pituitary_46	Pituitary_65	Pituitary_56	Pituitary_48	Pituitary_35	Pituitary_61	Pituitary_49	Pituitary_60	Pituitary_54	Pituitary_45	Pituitary_42	Pituitary_62	Pituitary_25	Pituitary_64
class	August	Torpor	IBA	October	August	August	August	Torpor	October	March	October	Torpor	Torpor	Torpor	March	August	March	August	IBA	October	March	IBA	IBA	March	October	August	March	IBA	October	August	October	August	March	March	October	IBA	October
AGAPEPAEPMELAKPSA	8440.891791	775.554963	894.443289	11.425360	1548.747707	16807.763410	7436.475845	24.346540	1116.529464	4772.051234	477.122095	1347.690917	4.423521	779.355400	5978.243266	14785.218490	1030.061559	1617.177966	963.767877	26.824983	4079.207738	1310.612505	1056.466008	1104.657872	109.499447	475.793453	1154.425233	1146.309064	952.829079	6014.112942	728.606820	6652.094453	40.168600	45.421386	133.549590	1696.183331	980.962340
LAGAPEPA	1637.378702	575.595032	975.630047	45.861858	618.251226	6120.353519	5723.721568	99.165909	48.814586	1434.796982	125.069634	8881.278769	163.538557	227.448079	872.068168	8555.658726	130.762483	1051.025765	15.254078	69.059286	2.970965	4940.096814	33.457682	40.529443	64.864598	2213.437898	82.727210	2029.214809	201.302627	5663.122020	55.439455	4929.551350	91.342434	73.612294	253.726328	5986.048232	218.966671
AFGPRGPGLQL	550.29629	649.65194	755.18678	345.80726	1221.06916	791.64920	881.68783	812.51054	20.84877	2777.31318	356.64258	456.75401	801.66424	688.52076	804.50951	548.14269	590.88412	333.61744	475.63361	624.53244	1184.96964	565.58643	362.17408	2421.97165	4249.72287	478.94991	1558.92600	1171.09202	2395.91454	1019.61708	2884.94196	2131.58046	1049.11926	453.65046	596.25238	180.22440	826.22788
AAMDLEL	525.06319	13.00213	874.75070	1218.98064	7439.16463	25.49921	297.86676	924.28122	689.19719	375.21130	629.74732	124.37370	784.78327	957.32957	2253.14606	671.64965	449.49625	303.24242	1481.36860	829.00897	611.15013	675.69158	736.71143	1012.66076	22.83013	952.48399	416.59759	57.17709	442.89607	552.91798	4101.33316	1482.06788	482.06852	814.28190	1113.41676	116.09223	137.89368
AELDQLLHY	1521.42903	1002.28099	319.04781	1904.33387	2971.66541	2247.60675	72.35338	1899.71042	2508.03824	4280.23291	3344.65898	858.16806	2704.74526	157.78970	796.25832	5154.79758	665.90566	331.78645	4495.39135	1127.25751	2517.78311	2150.77456	234.06884	3391.50560	1430.31875	360.19220	745.14195	244.52844	212.58989	1280.93174	706.79407	947.91799	2080.38023	1587.68654	613.10502	795.80551	718.56146
EDREEASLQDRQYASHHT	42.17504	797.29683	1235.93663	645.17823	184.07864	547.59149	489.51312	1215.11341	1156.37250	2059.00835	1077.36160	1069.12980	382.01498	980.73374	247.32133	586.13359	1602.37731	665.29413	1551.51409	172.45288	2739.57538	257.36413	1022.67087	2859.40863	697.95067	710.02483	364.97251	1410.20775	1408.38193	385.27649	767.71181	135.72669	703.49989	707.01072	632.72335	446.80660	519.89247
EDREEASLQDRQYASHHTTE	103.494268	527.248428	902.432299	530.758483	213.092060	538.638119	498.095470	977.873970	800.188114	1931.523871	531.038140	765.766642	414.387373	1069.247181	562.047648	645.154991	1606.570019	673.070093	1508.743502	235.332843	2955.934267	252.633601	1234.983948	4391.329206	679.960741	666.583467	662.772621	1163.094189	1411.795110	497.315730	1398.673698	307.393847	818.341676	207.556735	2.944484	534.580232	1040.775925
EEASLQDRQYASHHT	55.977706	577.709867	1326.664190	835.710334	55.430721	403.754716	396.416524	1575.441576	1555.499722	1816.940059	1266.284003	1069.694237	416.330001	941.406267	183.235430	222.398505	1248.631166	519.670512	1470.596814	252.783118	2282.001923	248.859735	1279.058948	2761.773044	1049.763295	208.921965	406.380606	1545.591591	1714.273409	102.133829	990.366681	71.689021	441.925256	110.825801	4.302888	505.348092	725.029992
EEASLQDRQYASHHTTE	139.71983	495.46860	1072.65147	572.76324	217.54310	527.15196	549.87897	1318.03658	940.23415	1734.41444	650.58705	787.96181	481.14419	988.25398	436.81115	623.64710	1790.77073	630.12132	1380.09498	267.59218	2384.29342	310.02315	1187.27888	3739.72553	1056.69198	580.88166	690.62213	1345.06364	2049.28193	368.99512	1596.60559	342.84747	566.92150	136.15521	25.30894	501.69637	1215.65739
EESEESEEGTTSEVT	254.28797	396.82302	1454.76710	524.26478	459.59644	223.76097	445.93341	2431.03205	803.17354	370.86905	1407.42960	397.53098	1538.47669	8100.06053	118.40360	403.21272	2990.34082	412.19670	2156.28479	1791.83580	2406.82879	323.59815	1303.54333	2270.37937	2172.79854	148.13757	283.37044	2946.79497	2340.76900	521.42110	2066.92534	424.30620	365.28014	65.94815	2086.61914	760.79018	987.37709
ELDQLLHY	1029.73311	1510.97697	1067.39391	669.81819	983.48923	628.20442	681.85894	1020.06360	1299.43517	1795.94991	1401.02161	460.30632	2154.10252	801.67790	1359.26169	856.45478	1468.28770	319.17416	1340.78192	788.31253	1803.85903	552.30906	790.46238	3595.49952	2729.17864	301.69341	2465.12908	1420.89800	1957.30849	97.86531	992.35921	417.34480	593.74461	32.96241	915.30199	194.77194	1188.81152
FLGEAYHHVPE	493.4311	322.7630	657.7509	369.3847	734.5341	629.2542	320.0933	677.0220	115.5980	2320.1347	374.6594	439.0358	949.0371	444.6169	1248.6644	1500.9906	300.2621	131.0724	1904.5890	1034.1041	670.1074	2005.1242	255.0434	1771.6900	1546.4188	127.9063	471.0338	291.5775	496.3285	612.4399	453.2683	162.8636	1206.0531	1271.5075	23.3822	1071.1865	326.5323
FLGEAYHHVPES	481.21481	654.24171	1047.65011	329.80621	842.32271	515.26592	682.86542	1326.22183	225.60144	2239.06766	721.41903	694.97339	1288.01876	504.01045	525.07038	528.91362	841.68651	1096.58151	1512.49735	1144.77095	954.72033	1152.73119	258.74070	1417.68697	1635.10268	72.67321	628.49496	643.30132	2065.45078	441.04958	787.10100	206.27630	957.96976	128.69579	213.47122	515.66261	285.22421
FLGEAYHHVPESQ	428.70360	758.26145	912.96618	406.40473	846.08793	79.05816	783.88383	1395.71003	249.02890	2313.96918	465.60644	809.88547	809.16139	749.01573	464.38649	554.87005	890.71067	299.68694	1294.83834	919.60470	841.72658	593.31631	185.36898	2804.12490	1713.15637	153.82720	1060.05182	1188.88783	2069.10972	84.42699	803.96634	495.43250	994.35460	367.69513	267.12389	564.01480	505.96548
FLGEAYHHVPESQRD	247.7085	488.8710	753.6601	482.4274	675.5019	372.5328	446.2955	883.9051	223.1465	2342.6750	325.4790	343.9048	700.3629	781.2719	1024.8392	527.7075	1813.8460	450.1383	1419.6284	851.5824	1643.1506	599.7167	362.1726	5741.6535	1112.5673	210.4716	2430.3198	963.9576	1414.9978	564.0439	1589.6277	677.4201	1246.0614	687.5427	20.1341	443.1950	530.4494
FLGEAYHHVPESQRDKA	63.83861	425.14360	708.01010	691.15817	345.96203	410.27961	421.94871	834.57931	232.34500	2402.56414	180.76880	548.61649	472.53191	820.97403	962.28249	745.17544	1806.78991	262.12347	1347.75389	681.88412	1137.75698	1544.45893	672.18398	5446.53105	645.23605	201.77124	1564.56968	609.19995	1340.20771	585.85442	2231.79484	614.44507	2304.55415	305.02820	949.50425	754.30831	762.21780
GEAYHHVPESQRD	108.5881	662.3954	1698.3930	352.3198	189.7030	106.8252	472.2040	1208.0028	446.6579	2028.1897	909.4631	782.0409	1480.9752	1081.6284	1589.6561	1685.8281	911.1899	160.9587	1582.7503	427.2400	1697.3173	1132.9635	1480.3471	3353.7126	987.9464	314.1462	894.3933	1224.2174	1607.2474	341.8750	2879.8256	530.0516	524.4536	330.8959	743.9444	876.1253	804.6718
GEAYHHVPESQRDKA	11.92236	395.62409	813.88227	501.06086	239.91843	400.09200	253.31832	734.02854	386.50858	1920.88127	174.43463	1008.40384	987.85891	684.55000	1872.89073	1242.10266	1235.12391	114.96755	1106.49652	428.55056	995.45208	2367.28765	2066.11269	4548.26675	720.33087	120.43822	805.45484	311.28189	1244.85174	432.56575	3258.74429	327.14693	1474.92964	251.63912	925.21924	457.54335	567.04523
GVLFNPYYDPLQWKSSRFE	6611.4346	875.9025	685.7050	472.3530	865.2550	1074.7683	3886.0569	331.1997	960.9155	7846.9287	543.5063	2637.2926	1005.2123	218.1399	1472.0790	6479.8435	1142.6181	203.2473	1050.6886	394.2170	1674.5674	2205.6678	960.0407	4358.9356	985.1988	1677.7208	2466.2344	881.2165	4798.2964	1318.4202	1715.2404	11978.9482	471.7387	2172.5147	1197.1231	1155.3646	1948.9287
LDQLLHY	275.9542	793.0110	895.8305	451.5764	1004.0843	279.7714	288.4017	1994.8197	932.9811	3206.3428	544.0862	219.0557	2486.7365	749.9878	1264.3156	621.4381	1052.6288	534.3379	1757.2866	1132.4615	2160.3010	1447.9382	865.8969	3079.2076	1420.9880	163.3341	1491.3648	724.8362	846.6669	262.0112	1123.4937	387.9258	1457.3647	1070.4640	1303.5664	455.8696	751.2653
LGEAYHHVPESQRD	167.21461	550.89299	761.28101	416.57572	664.14201	43.32686	192.79973	947.61855	127.53569	3021.89913	167.23181	190.43521	841.13448	878.00125	954.61050	721.53158	1938.28268	257.33678	1739.43515	878.69710	1862.44146	827.99034	184.08044	5794.52841	1276.08095	86.16034	2340.22125	730.11514	1150.07058	346.23813	1836.12807	822.88674	2173.45214	518.13479	738.44704	487.49897	39.91444
LGVLFNPYYDPLQWKSSRFE	8669.6495	627.6412	657.2180	2093.6500	4464.1618	2311.5561	10826.9857	1080.3823	1941.9735	9129.7677	1663.8509	2553.1065	417.4635	312.6697	442.5947	12800.4487	730.3036	806.9694	1341.2673	835.1249	3117.1786	2264.1320	176.6967	11557.5027	2057.3074	4560.8705	4188.6763	2091.2422	5271.7920	7873.1617	5597.6518	16651.9515	1866.7894	1832.6857	1731.6016	509.6936	1664.2611
Q(-17.03)KIAEKFSQ	394.65207	989.09485	1563.36256	174.54546	76.71573	631.76771	10.75832	2022.87482	332.67193	793.28098	239.85917	236.07875	396.79905	2330.69967	864.55441	241.44595	794.72469	581.56294	1167.23519	1036.85109	964.04812	412.62181	1028.77867	1979.95156	1382.01711	345.09604	1434.80009	298.36758	2180.02907	367.47702	2715.75645	74.14704	769.91394	125.13626	888.36737	313.74587	99.00356
Q(-17.03)KIAEKFSQR(-.98)	13.555279	412.189102	857.423054	63.616254	38.934850	530.182135	5.441377	1219.043920	366.004642	704.230341	77.965252	315.974786	418.366163	1271.987762	1415.228254	444.460457	1853.511233	230.219022	1428.640272	1358.385843	2389.541328	752.584749	1658.803085	3418.286575	1078.201212	706.032686	1213.248992	202.746358	896.524996	731.196235	3804.701044	120.987800	1791.039468	406.218462	311.066852	624.413222	103.382018
Q(-17.03)YDRVAELDQLLHY	16423.77064	220.46894	324.73862	3250.03108	3600.84988	3194.66819	13497.81109	854.97455	2641.47061	8033.12524	2085.40169	1450.06221	386.99210	271.41982	898.83112	15508.98794	849.99339	310.55343	840.38469	908.70988	2369.15980	2064.59586	46.19556	4049.32491	4933.21084	3718.76707	3083.31826	673.65307	3248.15809	5186.08688	4082.30745	33985.28329	3112.48130	4548.34755	95.67132	659.50459	567.47456
QYDRVAELDQLLHY	1779.169075	110.810896	1984.625763	3073.594982	4481.431768	1519.311376	2896.246777	1735.430427	2199.744258	3831.964478	3092.654594	548.393817	703.703883	504.263075	1197.232466	2186.854596	380.145695	372.932879	1400.544479	606.524804	4054.252621	1498.332456	383.825462	4.910719	2056.050241	1195.321366	578.218830	1101.266090	1118.642384	1406.298599	1456.922656	6738.520580	2090.537237	1501.578752	60.034898	425.761408	1175.979262
RVLTEEEEKELENL	1496.1039	1489.5749	1863.2634	425.9575	1257.1078	4316.7888	1761.4553	3408.0382	1564.9475	1699.1262	884.2519	553.1472	803.3353	3182.5054	2451.1502	1203.1286	362.0642	271.8096	5294.0775	4540.5105	1756.6278	1046.7054	1467.0762	1164.8802	1742.0443	1118.4210	1500.4058	1210.1720	3468.4995	2327.7286	10675.1755	408.1506	1036.4430	686.7855	126.4312	660.8108	225.2719
SQGGEPGAYLTPDT	933.8392	1007.3850	902.4255	354.2801	965.4500	302.4190	470.7034	875.1652	573.3923	1174.4731	320.0198	857.8360	878.6701	642.4561	959.2307	394.1289	1324.9163	692.6524	1408.6957	472.0319	937.2687	696.1563	500.3454	3518.6506	1927.9244	142.7497	391.3537	774.4224	1072.9653	410.5230	1002.5809	773.8843	627.6872	299.4090	568.1871	517.6635	1022.8996
SQGGEPGAYLTPDTREE	457.2856	436.5969	630.6995	330.7102	1144.8602	462.1013	501.0232	494.7017	579.0398	1321.0060	389.4226	565.9996	897.9201	861.9230	621.8995	617.8421	1222.0180	673.9228	1294.0884	1031.7024	769.4819	487.2670	506.1190	1807.9632	798.6428	182.9074	1293.4814	440.0699	730.5911	519.7184	907.3238	591.6321	688.3643	319.6933	894.5599	385.4963	774.6026
SSQEGNPPSEERGHTLQEPEEAKVGE	270.4985	730.6174	615.5500	436.6446	1081.0976	975.8356	1258.7788	1002.5869	569.9957	1896.9526	492.3425	1096.8643	638.1843	1216.7679	690.9403	1828.5668	1834.4817	762.8766	2071.3178	582.5196	2029.7862	853.9988	583.4784	3717.7003	2620.9343	788.7831	1875.7158	1705.9751	1000.8422	1656.7842	940.2866	1263.6121	2044.5158	190.9979	187.9178	1690.1970	929.7172
VAELDQLLHY	710.6329	1661.2074	995.6014	2024.5493	3404.2238	1685.2024	1091.5948	1560.6975	1872.9170	137.2283	1753.6776	503.8729	1332.1929	549.8334	1522.4090	1996.1246	540.6041	499.3167	2132.7675	989.0957	2696.9916	1340.7765	367.5124	5604.4396	1701.2064	644.0655	1015.1564	1223.6588	1170.9386	1159.5714	1486.2551	2652.1020	1270.3575	1139.8518	1164.6195	225.5111	539.7864
YDRVAELDQLLHY	1921.2349	779.4586	1183.3663	4642.7824	5231.8957	2200.0387	1612.0324	1169.5422	3437.6512	3720.9585	4187.5864	612.2748	571.2310	169.8196	2095.6228	8518.2837	193.6967	124.2603	2001.3383	1399.9847	4114.0251	3023.2594	116.8445	1047.7075	3320.1231	1317.3165	406.6363	125.3643	536.7377	3125.2496	1696.5310	7766.6467	3266.3952	4275.8072	1328.8320	866.0230	740.5959
Q(-17.03)EPGAPAAGM	1182.6064	576.3984	776.1936	614.1799	1189.9451	2163.6382	1409.6190	809.8199	537.4106	1427.3021	551.5436	765.6882	685.3389	691.1032	868.2769	1018.9770	1526.9570	1232.3924	1170.6636	871.0869	1397.0738	767.7241	608.4399	1653.6092	2759.6232	433.5372	148.7911	2045.6867	2017.5743	2077.3394	2125.9028	1808.6240	1003.3029	707.0681	1967.4411	649.9396	1036.2081
SLDSPAGPA	1407.60227	597.20486	899.90567	599.69070	876.35173	1004.43049	1218.87335	92.19355	1193.21654	744.44072	564.87617	814.56617	647.64301	167.62533	62.61696	1110.08840	4962.65211	2782.84682	60.08567	216.29559	1330.02735	729.76293	1297.97498	3689.57339	6127.75839	692.13749	6119.97577	5076.68455	1459.45084	451.28754	2567.89114	442.38150	1960.05130	56.60028	956.03057	2475.50683	3288.53589
SLDSPAGPAE	224.688672	776.741851	1139.862167	323.229867	785.439721	229.639017	1063.428607	110.943635	1542.458008	368.718257	434.823797	921.740524	389.621240	341.448954	39.152612	6.308664	5868.176667	2822.689795	6.518965	128.809849	779.917445	721.876138	1145.982574	3966.032611	4610.137183	168.784033	6080.818746	2725.968495	1275.506683	195.326074	1657.595108	164.686308	2047.773879	235.508453	5.756192	1898.797220	2854.147910
VKVGVNGF	2212.1589	879.3813	312.4563	1385.2188	4011.5349	2161.7012	354.8655	307.8283	2990.9895	735.4303	721.8823	1468.7655	541.6038	467.4043	2533.7810	728.5705	1485.0691	121.0391	897.2713	890.5916	1422.0956	1996.4345	787.9567	2124.6160	8452.9563	1242.5262	1293.7715	212.3951	314.7743	756.8240	293.1550	1231.0724	2685.2269	4635.4415	1491.9776	1439.5450	945.5468
Q(-17.03)HWSYGLRPG(-.98)	821.423755	887.395308	543.252364	573.822647	896.330038	156.486655	1322.969829	336.261309	271.589693	1232.051280	749.754688	55.270095	778.814142	349.856089	132.108752	89.240132	2185.141985	2174.950579	64.780758	326.457993	1046.780803	271.666059	7.128048	1345.333136	1344.761419	558.409425	2109.564757	1479.411528	2179.336941	608.336491	1081.103300	174.328387	1556.556669	994.782829	33.617020	544.336962	529.181934
SVSTVLTSKYR	511.74364	1066.83668	629.23461	630.64485	3291.97912	2792.32956	194.22572	502.86604	1242.77294	843.25636	742.65424	1058.82814	245.98463	877.25843	579.91026	329.10116	1844.81809	81.93880	1037.88314	2626.49562	815.30701	1796.57417	713.04349	1496.32049	2738.00545	385.33324	856.21166	73.93288	148.02640	1169.23286	597.59938	832.91949	1969.21699	1570.60247	1803.43038	1538.27855	1226.10399
GVANALSHKYH	84.97186	1420.81691	714.41442	1198.17922	1798.47488	473.04146	230.27930	859.93232	2164.51731	1565.97138	1204.75999	820.22776	218.86975	1253.04894	836.43356	1065.06708	4211.75240	83.51805	1722.24996	1795.43793	1781.37783	2311.58571	72.17162	2251.25115	2241.27363	1538.07529	61.26901	476.15935	116.28334	1298.83855	158.88758	911.07471	2812.79265	2195.76618	60.91228	19.43212	1171.83679
KVVAGVANALSHKYH	1295.3962	1228.5311	519.0045	855.6822	2369.9914	1028.9242	221.7337	476.4842	1179.8338	1920.9000	797.2765	925.3959	1419.1548	1188.2542	1135.3516	1072.8624	4026.0229	1482.9382	2048.9855	1439.2888	2041.0056	2698.7278	1232.3757	2661.1546	2109.3469	1183.9586	1319.4228	377.8791	1695.2464	2227.9271	1706.0109	1186.8453	3442.8773	2827.0645	1501.2990	1312.7818	968.7908
VANALSHKYH	67.77349	1510.54698	825.89552	2233.12137	901.70499	350.61214	348.92269	887.73587	4385.58244	1105.96771	2140.00307	690.91956	286.35218	1264.09210	840.82224	1769.60580	2323.03594	49.58380	894.60897	811.80865	1662.20278	3205.02220	57.44170	1938.65977	1856.28805	2007.90450	1601.26430	317.38474	53.90293	1569.80784	192.17078	1132.54046	1147.00757	1722.93609	64.90289	934.87208	1522.23353
GRTLYGFGG	210.79106	1047.61329	292.41303	1247.20614	2161.54613	2053.18357	245.41516	636.79697	1765.78190	168.56984	1112.69824	117.12027	732.53098	781.40538	1745.69030	395.69164	2406.19717	349.26655	1260.23134	1421.68460	1503.41041	1741.02772	127.54586	1569.82459	6212.41099	455.63762	262.09716	609.21820	587.59834	846.15612	400.16026	955.21457	3212.91745	1672.19179	1510.66823	96.54891	56.25776
RTLYGFGG	223.6174	978.2271	220.2125	1533.7512	2525.2230	1701.0789	174.2894	557.3945	2168.6592	301.2160	1100.3548	119.3319	641.3648	791.2343	1525.5741	331.1254	2565.5384	322.1529	1092.1647	1331.8070	1327.8681	1812.3511	191.5573	1191.3973	7186.2700	550.2879	283.0906	614.8341	626.6266	576.0633	562.7893	746.3481	3190.7879	2513.4609	1604.2847	148.7276	113.1286
TLYGFGG	1078.42158	2632.59413	1322.19361	2237.88478	4281.83058	2700.18873	3833.76464	228.52010	2696.80900	433.17355	2303.16981	36.58765	5557.08767	219.91817	1941.22982	825.04401	902.70937	3843.43252	449.70379	408.15663	505.09863	1527.71901	1161.16800	872.80191	5251.73941	509.61539	4088.48295	1983.02510	2220.11131	319.16476	991.16706	502.27567	1601.03847	3974.15696	375.33475	562.29194	1314.31080
AADQDLGPEAPPEGVLGA	2901.8157	672.7262	177.7539	358.4347	1313.7772	1262.3843	4857.9957	1052.5964	763.6487	3256.6101	755.3845	1450.1382	495.3376	116.4219	1540.4689	1655.2416	771.4656	174.7108	459.6791	491.2544	816.9146	1326.2971	518.6031	918.1944	1280.8434	1068.5097	2000.2317	1555.2968	840.4997	2130.0988	450.0622	3194.4151	1136.8468	230.8847	1079.4181	1843.1214	2223.0568
AADQDLGPEAPPEGVLGAL	20153.3480	956.6158	204.6107	6386.2001	8540.2739	3845.5771	22177.9979	2953.2585	3032.4969	15476.4061	2385.2506	2401.1125	164.0238	557.3158	2009.3125	9596.4939	827.6753	981.9257	997.2533	290.5575	3928.4518	1742.7590	1025.8917	5205.2513	8052.2519	6171.2306	3032.4938	1175.8996	2206.6297	8989.4532	939.2032	36212.0728	3388.6153	10166.8262	885.1045	2434.9761	2690.4065
LENPQPQAPARRLLPP	1190.76922	842.61542	439.98537	491.07796	930.34859	801.49427	1392.92750	188.72075	137.78095	977.42039	739.44097	977.16104	916.94428	422.38212	825.49809	780.67822	836.42948	1571.69028	262.91801	475.12640	781.69005	818.23824	206.25306	1274.53917	3190.16541	828.48681	1070.73824	1148.82049	1426.29262	2312.41875	1083.13822	781.48364	1992.09689	493.57180	34.18693	1305.72102	627.64834
VKMALQQEGFD	1433.4862	827.2746	407.0671	238.9144	2355.4541	5705.0479	2638.0575	629.5865	1136.9658	783.1829	257.5033	1726.0843	621.0280	334.7879	1295.9967	1052.0031	1172.2126	1299.0673	504.1724	721.1416	730.3526	1022.2065	388.7564	396.4622	3675.0981	1000.5473	842.4559	1938.2006	1100.3952	4473.3412	579.8407	1830.7952	1325.4213	370.3733	385.7423	880.4956	217.6216
GGFMRGL	1496.46263	1688.20867	591.83332	1498.25623	496.09793	177.63674	542.55803	1035.70028	934.90147	767.92367	1078.73139	1587.19706	5060.28470	90.51220	133.56343	253.05359	835.96457	2446.87459	656.92193	467.99198	792.96502	3460.92931	70.88582	1551.98913	4250.68241	589.34582	1106.41402	2006.17226	114.91533	648.92810	669.08974	883.41027	4533.29436	2048.69724	1950.28559	1273.68088	839.50276
Q(-17.03)LEDEAKELQ	775.99320	1479.36541	710.35774	459.31607	1647.03835	818.38893	161.52018	148.02224	306.23259	1303.45631	909.55371	1257.71787	1962.34358	213.01920	1692.68030	210.58887	2057.22799	4083.38718	177.80203	413.99111	1572.62460	1326.20146	1193.50844	1279.52295	5273.62608	767.82309	2320.95430	2982.26184	1476.29250	1064.10736	1366.63939	729.70625	1780.75972	578.13398	91.66786	1705.47619	4848.67783
SPQLEDE	276.940288	1396.146722	868.421715	1621.192868	541.077833	157.469711	1376.129368	173.778454	536.193917	599.786281	2129.672108	3.747196	2007.199859	182.917465	1630.921342	151.337529	2475.792186	3013.991442	168.921508	344.846822	977.798290	1594.029523	1164.211456	2069.255579	7458.458534	1465.719104	2979.703605	5069.150264	2789.936732	404.259081	1911.965961	407.764109	1594.143458	149.209407	50.698346	1416.144893	2577.922942
SPQLEDEAKEL	277.786221	1036.769599	654.393133	543.158332	1201.051478	532.904772	1153.984688	93.904069	227.779413	1256.263483	848.883078	11.723532	1718.711467	249.621916	232.122231	143.639516	2338.913313	3901.151115	230.170538	424.933721	1030.205023	1016.289722	9.938066	2206.949104	5638.847764	624.672814	1898.221903	1798.849490	1476.675831	906.315700	1807.089223	465.656026	2584.407912	542.809015	136.917168	1221.109312	2302.560668
SPQLEDEAKELQ	241.81852	927.04902	593.46943	412.70948	1069.11015	533.98925	801.28466	73.37816	214.17205	1128.25290	664.97569	858.80506	1447.58946	209.51772	185.81760	141.15500	1599.47307	2603.37298	284.19503	307.22308	614.62298	848.57225	844.03596	1204.66812	3106.45944	334.82910	1899.97455	1372.41877	1147.46899	452.55092	905.11587	234.48943	1240.27119	314.72025	158.90466	1060.96160	1903.28535
VGRPEWWMDYQ	1358.65238	1102.29997	889.81150	1231.63184	2599.58091	639.78954	1571.67682	247.10198	403.28590	1725.18970	2640.00610	44.69051	3200.06445	131.43542	129.04138	1423.58847	1377.74487	1899.69496	146.03149	528.95820	1174.19960	1615.34573	1051.83839	1615.29359	5583.56003	1107.01914	1930.71736	1734.28048	2425.28805	1193.97178	2159.53222	1164.79296	1945.54173	636.37644	864.14885	947.59896	1895.08407
YGGFM(+15.99)RF	308.72641	1857.00611	431.99899	1104.60674	3090.62302	1078.02620	1193.75619	144.04006	363.70816	499.25315	1458.25799	1265.20513	1295.66668	234.68305	162.85851	33.11038	2946.16832	2884.89510	343.08152	1159.78803	1007.22268	1484.29441	1485.36170	1133.40160	10249.66731	838.89131	3793.75036	4006.98441	4505.30649	534.95961	3558.76461	772.80113	2392.42032	373.12693	427.22210	1326.31802	2252.70854
YGGFM(+15.99)RGL	379.2099	2089.9072	731.6304	841.8555	2775.9645	374.2918	935.5794	83.4596	402.6122	1038.4256	1462.9972	1481.3787	1453.9900	210.7689	119.9386	1644.0987	3069.9812	3360.3432	374.7593	924.2952	817.3762	1649.4453	1548.2935	2332.4758	9752.8839	825.9417	3705.4005	4033.7214	3343.9755	635.2441	3495.1601	835.5263	3050.4199	329.4857	602.6551	2150.2488	2565.2627
YGGFMRF	497.01736	1072.45165	657.85941	390.74855	1517.08828	244.21196	1554.66620	112.81281	158.01486	1131.18459	946.43495	25.09619	2253.31242	114.42392	140.31986	66.76026	1874.88530	2053.41614	90.99779	339.80553	464.00538	570.87299	10.20850	2234.56324	5548.33404	626.14708	3144.52475	1639.39121	1949.17177	374.88962	1006.41738	379.70159	1174.01710	198.78394	20.12067	735.51331	1280.08874
YGGFMRGL	455.524465	1076.660757	731.799562	503.240530	885.183573	347.129756	1229.231896	64.953071	104.692304	976.237469	746.840484	13.584416	1732.230280	175.666854	166.674002	51.737575	1598.639631	2114.666742	142.632741	345.057006	1087.930780	1217.931336	9.304989	1577.833891	5638.117818	713.249485	2884.506078	1915.576318	1685.077620	825.545519	1258.356088	298.362835	1707.331613	515.077918	11.420737	723.209313	626.059157
EIGDEENSAKFPI(-.98)	684.13642	1410.96250	646.02682	302.45713	3245.14893	1194.65850	1627.76406	705.96206	829.37014	855.18367	660.46934	159.38188	1395.32122	824.27309	94.08887	212.89589	1681.87644	2244.52753	767.08280	420.17673	1534.93003	867.05916	72.27757	307.27640	2052.56475	814.77669	2291.23308	1853.18938	1223.98844	506.21449	744.46896	871.47618	2913.11112	347.84685	1713.70450	1711.82257	1297.38821
DDGPYRMEHFRWGSPPKD	888.0572	963.9873	778.7072	159.2630	1088.9537	3044.2984	2095.5458	682.8386	701.4690	7495.9115	412.2031	1752.1296	526.3064	357.5323	2394.6282	1438.4015	767.3897	365.3827	805.7882	745.8070	1334.7368	2395.1000	501.5079	666.1333	389.5712	1328.1925	523.2651	1088.8715	488.6027	8686.6042	227.8335	1777.8577	1591.5709	364.4288	1053.8615	1302.8425	157.8620
FKNAIIKNAH	268.18564	2262.00275	630.90863	222.67547	459.62972	538.33252	2842.81733	1521.00816	2097.98308	1872.80321	741.44975	1235.61020	431.51476	654.35859	2376.64313	2678.26587	2196.37466	19.90206	545.49302	683.86008	834.08123	1213.78134	1969.22179	1076.24002	1284.69057	2033.28923	1179.15512	2507.17002	451.98455	4786.01655	923.92415	1309.30975	1885.42008	263.12043	984.08458	1817.17792	1008.37974
FKNAIIKNAHKKGQ	1068.24973	616.22640	565.04669	360.62832	346.00555	77.54253	1326.01402	1852.36027	1833.90254	978.62750	4537.85561	261.39455	696.48980	1427.84678	943.04308	4570.06188	147.62783	671.19572	990.27794	849.31405	2653.58271	635.24536	2415.98368	1810.29062	228.90026	1771.95713	26.62369	458.20485	426.99343	1281.29942	287.09892	904.95755	90.10210	262.74748	463.27945	475.99949	47.21013
GGPEFRDDGAEPGPRE(-.98)	184.34248	409.66628	516.77914	838.62365	2540.15618	1655.25502	222.58258	519.13462	1381.73306	1554.81341	575.42936	220.52535	163.16835	447.08426	1992.40359	1841.38729	1067.70061	162.55123	1192.38880	3808.54640	3507.88639	1130.17613	72.16027	1693.86708	3408.05571	1422.56663	352.24162	255.41762	534.41627	1162.01400	350.11235	604.42017	1723.70927	184.02595	257.81121	102.00221	117.24573
GPEFRDDGAEPGPRE(-.98)	230.43919	491.78614	586.33114	1742.38708	2758.03902	1447.44269	239.64908	486.38715	1509.67262	1083.83336	1165.46144	161.95845	157.08193	404.09396	1897.64600	2055.16359	1128.97156	165.90879	1231.57062	1824.44855	4294.58862	1313.28439	92.07857	1940.39295	3481.14670	1865.25933	422.14152	306.39618	610.68631	943.73315	372.71368	598.73971	1605.97972	1744.80208	475.44030	162.96999	144.46238
KYVMGHFRWD	666.83001	919.03956	488.17867	79.26367	2276.84373	3537.79734	1259.29414	891.59384	820.76288	1853.21539	132.03477	1190.99449	216.61804	316.29118	1862.93105	847.15950	796.49801	365.55210	476.94614	1044.53757	942.44799	1692.66198	824.08453	257.70230	646.71090	908.59956	335.65493	1020.51415	171.39937	3264.04992	495.70035	1228.04311	1074.60414	287.91871	799.49838	1156.08301	545.24816
KYVMGHFRWDRF(-.98)	523.94947	610.84660	208.12384	103.28573	1553.99133	3179.39383	1145.10540	690.21057	787.75085	1489.83337	220.14246	1178.58549	222.54323	231.15933	2435.86516	1234.28548	847.97340	347.19952	568.89702	447.92820	1414.73048	1804.43796	604.92168	638.78623	195.89304	1100.33812	520.56015	443.48029	291.18103	5564.39222	135.57800	1092.37310	2722.20930	357.51755	1007.42049	1317.40521	80.29325
PEFRDDGAEPGPRE(-.98)	101.5571	596.9674	452.4494	1990.4450	2191.1090	1305.3927	318.2228	753.1191	1962.0430	1404.0559	709.4040	205.6463	160.8840	479.5994	1850.9203	1693.1350	1543.1788	225.1386	1603.9442	1974.4034	4349.6843	1117.5126	888.3964	1639.6116	2423.2787	1500.7441	666.3537	185.0070	875.8161	802.2767	332.5011	340.0470	1188.1347	1039.8707	188.8059	242.9633	173.2729
RPVKVYPNGAEDES(+79.97)AESFPLEF	11070.911012	1024.212935	598.930738	803.337809	1890.399514	5483.960323	12819.912157	1280.582473	2013.762122	2.263158	985.079174	985.226118	177.458969	149.420336	1288.865876	11322.215087	137.203384	350.236262	452.653906	1318.013793	525.399139	1313.939201	89.573796	398.426466	4241.817379	1952.205252	655.557869	1316.604545	1001.470239	4223.601885	559.535447	8715.438528	575.801552	411.634553	98.938620	314.634411	853.556459
RPVKVYPNGAEDESAESFPLEF	8413.8285	413.0344	338.3428	661.9594	3161.6411	11184.3257	3404.5756	2516.1064	1870.7808	2733.8307	998.1490	2845.4641	510.6685	106.7193	3105.8971	2646.3593	181.0622	109.3660	690.8507	1919.2375	442.4930	311.2595	120.4821	1401.7923	1986.6106	585.1426	425.6389	859.3504	318.6119	1767.5116	420.8891	7855.0619	150.8136	317.5392	130.3227	436.2797	860.9705
S(+42.01)YSMEHFRWGKPV(-.98)	1445.26051	1034.17101	594.73437	216.16841	1849.27024	5106.88113	2373.08811	762.46943	751.09516	3995.77067	222.74034	1434.68463	221.32407	273.21215	2159.03240	1318.43543	561.81899	201.60802	621.99264	1744.30736	867.82676	1931.20450	354.59992	1123.50072	1586.33718	822.79776	772.28777	1137.27354	556.57733	6273.38074	212.44451	3022.66171	2174.82932	304.06473	20.09392	1107.82212	42.29011
SQTPLVTLF	2026.3571	1338.3774	1623.4480	2445.3891	2433.8217	2375.7449	1749.4797	2302.1809	958.0414	1496.5210	1965.0195	509.4767	966.5813	992.6191	1519.5531	921.3968	731.3100	782.0178	1802.1732	1235.6303	1699.5273	241.5109	1093.2346	676.9668	3519.4247	369.1602	3336.4293	1910.7103	1857.1677	1348.9777	1192.4735	2850.4817	440.4337	347.0073	3960.2972	18.1852	1008.0168
SQTPLVTLFKNAIIKNAHKKGQ	117.48860	757.42141	280.49213	1837.08784	1090.60960	1926.86155	1054.30282	550.62580	1280.11986	26577.13347	2518.14077	162.30858	95.95411	35.02182	219.53943	1019.87965	1126.77236	786.73696	1416.79001	893.44103	2665.32701	550.16239	871.62217	4885.21030	790.21630	515.94608	1189.60516	972.49127	1543.78200	430.22349	235.00719	1618.50745	722.11728	549.73388	525.34170	826.92977	1804.46577
SYSMEHF	251.8199	454.5423	1294.8463	720.5086	1186.6389	1207.4403	477.5258	1406.0815	662.9903	1143.6044	873.6375	595.1338	544.0245	1035.8847	1347.5163	388.6288	1444.8139	1321.3567	1109.8358	1271.2487	1666.5523	596.2820	1446.0838	1652.4600	1869.7001	247.4647	752.6563	1594.6875	1185.4428	772.3699	1695.7522	1509.0425	687.0539	1371.3856	2411.2158	385.8057	1447.4501
SYSMEHFRWGKPV(-.98)	871.230465	500.403496	643.408852	162.443722	1217.779121	1200.968282	792.456632	678.602887	282.251955	4061.782860	364.162704	1137.518104	442.042424	373.295725	1345.266997	1461.102210	1093.807251	862.250912	1337.728993	685.261637	2349.972008	1047.818175	392.938806	1006.289154	482.209703	482.141417	581.026236	416.758618	566.741233	899.759385	377.526805	602.579019	621.020594	1414.679349	4.007785	846.135939	735.753967
TPLVTLFKNAIIKNAHKKGQ	39.017319	59.253908	324.591693	5997.692736	4015.516224	1727.033128	265.104454	59.435508	2515.269223	1015.082744	614.686639	1326.689093	582.734372	51.543614	2577.540941	3281.249169	145.988688	5.532145	599.201862	606.678534	2751.414666	911.412250	708.273628	2352.920690	1273.930689	3405.909709	130.065672	115.108949	123.705923	592.024646	125.156847	1279.954375	1360.806783	8028.490536	1589.838004	1710.198870	951.774009
Y(+42.01)GGFM(+15.99)TSEKSQTPLVTLFKNAIIKNAH	14162.883902	2365.057731	158.584889	150.033360	1631.640545	9011.993295	14121.987855	459.118373	2756.309415	2474.392467	375.559132	1208.715658	1036.564495	749.624185	2058.697338	2656.479389	668.447605	511.980078	519.999158	855.164940	436.039162	3514.722470	158.205674	895.194346	3636.670696	2895.218123	2333.079053	2940.203022	744.520437	5817.330372	325.354293	13247.741149	1609.137381	196.111785	8.994823	2012.047551	640.937726
Y(+42.01)GGFMTSEKSQTPLVTLFKNAIIKNA	23585.955975	382.829697	561.830452	2.812048	3284.611173	19705.363557	23676.182460	387.076220	1236.833564	5527.907107	443.781787	2161.675193	54.842816	484.070825	288.411734	6380.381734	136.686925	162.916960	133.572199	430.824586	148.733429	1542.313664	560.927937	437.555832	2188.827943	2285.938157	586.843613	1073.933933	324.560096	10578.781438	137.602105	21344.819396	1093.705647	135.836551	865.835574	1553.389076	501.464470
Y(+42.01)GGFMTSEKSQTPLVTLFKNAIIKNAH	10562.507025	1313.873053	326.827403	185.917882	1569.881824	6231.741689	11887.322857	875.327942	1576.701996	4496.915948	369.626719	1895.314932	116.868448	128.287745	1991.375605	3063.330306	588.820515	322.615369	504.834166	624.502360	931.560661	3005.680615	318.987323	1286.978004	3566.243306	3087.327742	2298.889006	3011.368614	702.229602	12637.042044	197.312580	8671.848646	1288.368399	153.314102	1.266135	1189.672749	92.427463
Y(+42.01)GGFMTSEKSQTPLVTLFKNAIIKNAHKKGQ	1665.118109	2358.463974	462.164207	323.549390	124.093425	2876.484699	1764.176098	1967.862107	737.734284	3099.549506	918.660769	536.229356	273.353854	411.823698	2090.703456	994.902695	458.077857	652.567267	1210.835641	229.356185	1213.990431	2246.736673	807.261405	565.522042	624.540004	756.459013	1449.365160	919.152149	290.393644	5464.374704	324.698935	3339.414959	716.832285	137.239539	5.673444	1017.026609	257.101993
YGGFM(+15.99)TSEKSQTPLVTLFKNAIIKNAHKKGQ	855.61804	1564.62538	679.22428	3299.93371	988.48249	419.93398	1529.47214	1928.35277	1212.12873	2146.84750	3696.69073	189.13095	111.49999	1318.82335	375.64435	1651.57970	2915.53495	2794.01381	2534.75999	887.30389	1905.95113	297.58898	1149.79025	5203.24189	555.50651	1204.74838	6925.62624	2044.63784	3871.58467	346.47124	1266.41939	1445.10787	11.67290	590.22729	86.98747	465.23928	6528.60896
YGGFMTSE	698.4797	461.6090	926.2728	579.2999	1302.5862	1772.6989	726.7832	1327.9245	470.9490	4159.3878	559.5050	505.2308	551.4393	1257.9355	1113.2142	752.3002	1126.7745	1627.2060	1244.2337	840.4334	1524.7202	771.6526	1253.4835	1401.1738	2089.1147	255.5741	2800.7054	1989.6387	1203.5297	1502.7101	1267.6343	3296.1996	546.6067	937.7522	2665.3345	332.3717	1065.0374
YGGFMTSEKSQTPLVTLFKNAIIKNA	2951.5013	1539.9506	646.0024	689.7253	3335.0018	6548.3851	3512.8243	1398.3273	1727.1455	6010.2293	826.5009	2729.5885	184.4494	250.1372	1686.4381	7816.5240	817.2334	1608.1524	1620.6477	794.3135	2752.2311	3901.2476	426.4441	1762.2194	11793.6100	2887.3495	2554.2672	2294.5725	3100.5964	5287.8632	2175.2867	6908.9752	1365.0428	1473.6719	2038.7831	2393.6143	1054.3216
YGGFMTSEKSQTPLVTLFKNAIIKNAH	873.61538	2534.52617	1003.51082	474.94346	2619.84130	3612.68479	997.01071	1014.50509	686.67347	1824.02489	219.21242	904.18906	904.63712	304.38916	2327.82683	4316.03166	86.29545	606.95202	709.87681	444.40690	1193.59729	2882.42940	574.67644	201.99292	4608.70178	1080.22942	531.03424	2120.79391	1045.03194	2485.87129	730.94920	3774.92429	530.61881	405.83311	1057.42788	2154.77732	503.54405
YGGFMTSEKSQTPLVTLFKNAIIKNAHKKGQ	393.11266	830.16951	900.89913	1823.70062	969.88503	515.22282	1392.10677	2702.47749	2313.65046	3007.59031	3636.37601	605.44693	861.46873	863.07440	729.28256	820.22952	1335.47654	998.89814	1162.14875	437.77553	2709.20596	369.74234	1043.47270	3537.03313	310.27799	734.82812	3108.06961	817.92246	1168.32549	368.95374	819.46657	998.11605	331.29404	184.76103	60.03794	342.45879	2596.88959
ATLESVFQELGKLTGPSSQKRE	3294.04044	903.90235	935.07711	1895.27679	777.60628	2856.17627	1788.73418	959.81735	1316.80474	3300.47857	2771.12015	2015.65181	2279.14060	416.85423	2890.40177	1192.75292	1197.34697	2522.08608	5283.46813	1174.28619	1925.02997	2127.06071	1720.28421	4274.33010	1432.29168	191.30353	2295.86514	96.88988	2327.35350	216.47887	7822.12955	542.79891	430.69764	707.59705	721.83937	727.12987	1677.23823
MYEENSRDNPF	305.1265	567.2840	946.3929	379.4790	776.3466	847.9241	598.9437	973.3149	1012.5519	1420.9099	547.8667	565.5895	981.6504	751.4806	837.0393	311.9473	1301.8650	708.3824	1629.0183	609.6255	1354.0874	727.9994	1024.1684	1403.1404	2867.3674	196.5511	650.2069	1118.6491	2479.7351	434.4935	3143.0369	874.4292	525.1122	656.3637	332.9595	307.8804	417.7585
TNEIVEEQYTPQSL	2314.7539	1336.9040	3433.1293	425.7653	1255.2410	3517.1165	290.5042	1436.5918	2870.4362	3490.4357	279.9286	417.9799	1378.7863	1181.1678	2541.7544	2009.4243	703.1639	42.0380	5958.2769	3529.5288	1472.9574	1687.4461	2177.8225	3152.0782	3019.2665	656.9153	3505.1688	275.2392	4599.9721	751.0174	13376.8058	1084.1515	969.0547	677.7245	2963.3904	777.1403	136.0975
FPTLGGSQDKSLHN	664.2153	370.2542	557.7496	221.4898	597.3154	606.4768	371.3213	588.3091	372.4494	1565.7994	345.6432	698.4697	672.1565	578.0652	1128.3113	509.1754	1963.2686	864.1597	720.0531	927.4844	1430.8851	873.7425	926.7753	1476.3045	2076.4818	389.9772	2021.4806	802.8369	1225.4325	954.9039	1629.4122	434.4113	1162.6729	435.7855	1296.2335	571.3601	589.7155
TLGGSQDKSLHN	9.40251	592.52172	802.17080	561.50195	115.96206	159.53973	349.82833	1054.58435	2160.66247	1168.83266	714.69657	360.78445	885.07056	850.57742	1887.46821	1041.36511	2003.38339	294.67852	547.65671	370.16717	3213.77105	889.56429	2026.84314	2339.60508	1338.58510	596.74598	1144.12295	602.37438	1464.63799	323.35547	1665.69287	228.59093	710.10797	443.24938	434.63305	51.13244	390.57529
SVNPYLQGQRL	485.96282	411.30195	733.85188	38.58792	86.95171	817.31076	22.67851	1159.42085	513.10949	496.82094	53.22789	370.91560	607.23846	808.71197	2386.96739	196.65309	1552.83320	125.40403	1046.38220	829.46815	1821.59823	529.74321	1004.27342	1352.80673	1174.56028	277.67287	1059.35390	145.33931	1015.85492	439.73768	3538.90011	79.60038	884.72503	1232.92605	1142.59679	288.76396	30.83659
SVNPYLQGQRLDNVV	3815.52145	610.10195	851.37226	2049.88684	2490.87854	6634.58109	3752.73163	751.16707	1720.10395	2761.64171	1618.63698	1535.62826	17.36857	405.98752	774.78039	3009.91874	429.25437	506.77740	1262.37252	1164.57054	2597.84915	685.81527	242.01545	2403.89090	4948.88733	2261.88994	2771.05251	1594.16713	2287.64083	4374.91678	991.23702	5924.71792	728.08695	1855.16075	1542.86184	609.48502	624.73813
SVNPYLQGQRLDNVVA	2694.0518	647.5007	673.2173	1477.8779	2094.4568	3090.5216	2525.7932	720.1189	1055.1988	3077.7936	1142.0299	1543.8883	530.1360	405.8454	1761.6106	1651.9937	1074.6270	621.1074	1130.0799	913.4880	3208.9678	978.3610	404.5577	3778.7130	3976.7559	1346.4645	2928.0116	1292.0485	2099.3698	3167.2281	1268.3325	3335.1360	1377.1888	1314.7793	1541.2214	341.3430	602.0742
SANSNPAM(+15.99)APRE	1260.112708	2407.584605	1322.967798	1138.071781	781.401032	21.081017	916.543676	30.822463	568.511184	688.646946	1801.814999	1498.654531	958.876436	45.279093	5.431026	35.691451	4102.192382	2906.255234	270.550403	162.012907	1802.270129	1201.215803	1611.176928	4154.590623	3845.167459	946.164790	1575.399704	2388.857335	6502.607984	125.496761	3124.517422	31.557790	2130.711937	297.813342	1347.675117	2533.847678	1858.003876
SANSNPAMAPRE	216.84933	1064.92141	1023.03393	632.29622	633.01556	56.89700	1342.72703	64.95270	324.50022	1817.11808	731.77210	32.48293	844.89077	48.98649	187.64042	84.60936	3282.09110	1949.56949	99.86112	430.44818	1768.89111	620.53899	15.36907	3199.86188	3633.86203	1038.03566	2087.23858	3026.89867	2434.76128	309.08295	1282.94845	86.84305	1642.58491	474.95798	126.71698	834.54089	1060.56441

Log2 transformation

Data transformation applies a mathematical transformation on individual values themselves. For mass spec data, log transformation is a good choice, as it reduces or removes the skewness of mass spec data.

Code

plot_grid(hist_log2, qq_log2, pca_log2, nrow = 1)

Code

#rownames(log2_d1) == feature_name$variable
rownames(log2_d1) = feature_name$name
# display table
log2_d1 %>% 
        mutate(across(everything(), round, digit=2)) %>%        
        knitr::kable(., digits = 3, "html") %>%
        kableExtra::kable_styling(bootstrap_options = c("striped", "hover", "condensed", "responsive")) %>%
        kableExtra::scroll_box(width = "100%", height = "500px")

	Pituitary_53	Pituitary_31	Pituitary_38	Pituitary_63	Pituitary_55	Pituitary_51	Pituitary_50	Pituitary_34	Pituitary_59	Pituitary_43	Pituitary_66	Pituitary_26	Pituitary_19	Pituitary_24	Pituitary_41	Pituitary_57	Pituitary_47	Pituitary_52	Pituitary_23	Pituitary_58	Pituitary_44	Pituitary_37	Pituitary_39	Pituitary_46	Pituitary_65	Pituitary_56	Pituitary_48	Pituitary_35	Pituitary_61	Pituitary_49	Pituitary_60	Pituitary_54	Pituitary_45	Pituitary_42	Pituitary_62	Pituitary_25	Pituitary_64
AGAPEPAEPMELAKPSA	13.04	9.60	9.80	3.51	10.60	14.04	12.86	4.61	10.12	12.22	8.90	10.40	2.15	9.61	12.55	13.85	10.01	10.66	9.91	4.75	11.99	10.36	10.05	10.11	6.77	8.89	10.17	10.16	9.90	12.55	9.51	12.70	5.33	5.51	7.06	10.73	9.94
LAGAPEPA	10.68	9.17	9.93	5.52	9.27	12.58	12.48	6.63	5.61	10.49	6.97	13.12	7.35	7.83	9.77	13.06	7.03	10.04	3.93	6.11	1.57	12.27	5.06	5.34	6.02	11.11	6.37	10.99	7.65	12.47	5.79	12.27	6.51	6.20	7.99	12.55	7.77
AFGPRGPGLQL	9.10	9.34	9.56	8.43	10.25	9.63	9.78	9.67	4.38	11.44	8.48	8.84	9.65	9.43	9.65	9.10	9.21	8.38	8.89	9.29	10.21	9.14	8.50	11.24	12.05	8.90	10.61	10.19	11.23	9.99	11.49	11.06	10.03	8.83	9.22	7.49	9.69
AAMDLEL	9.04	3.70	9.77	10.25	12.86	4.67	8.22	9.85	9.43	8.55	9.30	6.96	9.62	9.90	11.14	9.39	8.81	8.24	10.53	9.70	9.26	9.40	9.52	9.98	4.51	9.90	8.70	5.84	8.79	9.11	12.00	10.53	8.91	9.67	10.12	6.86	7.11
AELDQLLHY	10.57	9.97	8.32	10.90	11.54	11.13	6.18	10.89	11.29	12.06	11.71	9.75	11.40	7.30	9.64	12.33	9.38	8.37	12.13	10.14	11.30	11.07	7.87	11.73	10.48	8.49	9.54	7.93	7.73	10.32	9.47	9.89	11.02	10.63	9.26	9.64	9.49
EDREEASLQDRQYASHHT	5.40	9.64	10.27	9.33	7.52	9.10	8.94	10.25	10.18	11.01	10.07	10.06	8.58	9.94	7.95	9.20	10.65	9.38	10.60	7.43	11.42	8.01	10.00	11.48	9.45	9.47	8.51	10.46	10.46	8.59	9.58	7.08	9.46	9.47	9.31	8.80	9.02
EDREEASLQDRQYASHHTTE	6.69	9.04	9.82	9.05	7.74	9.07	8.96	9.93	9.64	10.92	9.05	9.58	8.69	10.06	9.13	9.33	10.65	9.39	10.56	7.88	11.53	7.98	10.27	12.10	9.41	9.38	9.37	10.18	10.46	8.96	10.45	8.26	9.68	7.70	1.56	9.06	10.02
EEASLQDRQYASHHT	5.81	9.17	10.37	9.71	5.79	8.66	8.63	10.62	10.60	10.83	10.31	10.06	8.70	9.88	7.52	7.80	10.29	9.02	10.52	7.98	11.16	7.96	10.32	11.43	10.04	7.71	8.67	10.59	10.74	6.67	9.95	6.16	8.79	6.79	2.11	8.98	9.50
EEASLQDRQYASHHTTE	7.13	8.95	10.07	9.16	7.77	9.04	9.10	10.36	9.88	10.76	9.35	9.62	8.91	9.95	8.77	9.28	10.81	9.30	10.43	8.06	11.22	8.28	10.21	11.87	10.05	9.18	9.43	10.39	11.00	8.53	10.64	8.42	9.15	7.09	4.66	8.97	10.25
EESEESEEGTTSEVT	7.99	8.63	10.51	9.03	8.84	7.81	8.80	11.25	9.65	8.53	10.46	8.63	10.59	12.98	6.89	8.66	11.55	8.69	11.07	10.81	11.23	8.34	10.35	11.15	11.09	7.21	8.15	11.52	11.19	9.03	11.01	8.73	8.51	6.04	11.03	9.57	9.95
ELDQLLHY	10.01	10.56	10.06	9.39	9.94	9.30	9.41	9.99	10.34	10.81	10.45	8.85	11.07	9.65	10.41	9.74	10.52	8.32	10.39	9.62	10.82	9.11	9.63	11.81	11.41	8.24	11.27	10.47	10.93	6.61	9.95	8.71	9.21	5.04	9.84	7.61	10.22
FLGEAYHHVPE	8.95	8.33	9.36	8.53	9.52	9.30	8.32	9.40	6.85	11.18	8.55	8.78	9.89	8.80	10.29	10.55	8.23	7.03	10.90	10.01	9.39	10.97	7.99	10.79	10.59	7.00	8.88	8.19	8.96	9.26	8.82	7.35	10.24	10.31	4.55	10.06	8.35
FLGEAYHHVPES	8.91	9.35	10.03	8.37	9.72	9.01	9.42	10.37	7.82	11.13	9.49	9.44	10.33	8.98	9.04	9.05	9.72	10.10	10.56	10.16	9.90	10.17	8.02	10.47	10.68	6.18	9.30	9.33	11.01	8.78	9.62	7.69	9.90	7.01	7.74	9.01	8.16
FLGEAYHHVPESQ	8.74	9.57	9.83	8.67	9.72	6.30	9.61	10.45	7.96	11.18	8.86	9.66	9.66	9.55	8.86	9.12	9.80	8.23	10.34	9.84	9.72	9.21	7.53	11.45	10.74	7.27	10.05	10.22	11.01	6.40	9.65	8.95	9.96	8.52	8.06	9.14	8.98
FLGEAYHHVPESQRD	7.95	8.93	9.56	8.91	9.40	8.54	8.80	9.79	7.80	11.19	8.35	8.43	9.45	9.61	10.00	9.04	10.82	8.81	10.47	9.73	10.68	9.23	8.50	12.49	10.12	7.72	11.25	9.91	10.47	9.14	10.63	9.40	10.28	9.43	4.33	8.79	9.05
FLGEAYHHVPESQRDKA	6.00	8.73	9.47	9.43	8.43	8.68	8.72	9.70	7.86	11.23	7.50	9.10	8.88	9.68	9.91	9.54	10.82	8.03	10.40	9.41	10.15	10.59	9.39	12.41	9.33	7.66	10.61	9.25	10.39	9.19	11.12	9.26	11.17	8.25	9.89	9.56	9.57
GEAYHHVPESQRD	6.76	9.37	10.73	8.46	7.57	6.74	8.88	10.24	8.80	10.99	9.83	9.61	10.53	10.08	10.63	10.72	9.83	7.33	10.63	8.74	10.73	10.15	10.53	11.71	9.95	8.30	9.80	10.26	10.65	8.42	11.49	9.05	9.03	8.37	9.54	9.77	9.65
GEAYHHVPESQRDKA	3.58	8.63	9.67	8.97	7.91	8.64	7.98	9.52	8.59	10.91	7.45	9.98	9.95	9.42	10.87	10.28	10.27	6.85	10.11	8.74	9.96	11.21	11.01	12.15	9.49	6.91	9.65	8.28	10.28	8.76	11.67	8.35	10.53	7.98	9.85	8.84	9.15
GVLFNPYYDPLQWKSSRFE	12.69	9.77	9.42	8.88	9.76	10.07	11.92	8.37	9.91	12.94	9.09	11.36	9.97	7.77	10.52	12.66	10.16	7.67	10.04	8.62	10.71	11.11	9.91	12.09	9.94	10.71	11.27	9.78	12.23	10.36	10.74	13.55	8.88	11.09	10.23	10.17	10.93
LDQLLHY	8.11	9.63	9.81	8.82	9.97	8.13	8.17	10.96	9.87	11.65	9.09	7.78	11.28	9.55	10.30	9.28	10.04	9.06	10.78	10.15	11.08	10.50	9.76	11.59	10.47	7.35	10.54	9.50	9.73	8.03	10.13	8.60	10.51	10.06	10.35	8.83	9.55
LGEAYHHVPESQRD	7.39	9.11	9.57	8.70	9.38	5.44	7.59	9.89	6.99	11.56	7.39	7.57	9.72	9.78	9.90	9.49	10.92	8.01	10.76	9.78	10.86	9.69	7.52	12.50	10.32	6.43	11.19	9.51	10.17	8.44	10.84	9.68	11.09	9.02	9.53	8.93	5.32
LGVLFNPYYDPLQWKSSRFE	13.08	9.29	9.36	11.03	12.12	11.17	13.40	10.08	10.92	13.16	10.70	11.32	8.71	8.29	8.79	13.64	9.51	9.66	10.39	9.71	11.61	11.14	7.47	13.50	11.01	12.16	12.03	11.03	12.36	12.94	12.45	14.02	10.87	10.84	10.76	8.99	10.70
Q(-17.03)KIAEKFSQ	8.62	9.95	10.61	7.45	6.26	9.30	3.43	10.98	8.38	9.63	7.91	7.88	8.63	11.19	9.76	7.92	9.63	9.18	10.19	10.02	9.91	8.69	10.01	10.95	10.43	8.43	10.49	8.22	11.09	8.52	11.41	6.21	9.59	6.97	9.80	8.29	6.63
Q(-17.03)KIAEKFSQR(-.98)	3.76	8.69	9.74	5.99	5.28	9.05	2.44	10.25	8.52	9.46	6.28	8.30	8.71	10.31	10.47	8.80	10.86	7.85	10.48	10.41	11.22	9.56	10.70	11.74	10.07	9.46	10.24	7.66	9.81	9.51	11.89	6.92	10.81	8.67	8.28	9.29	6.69
Q(-17.03)YDRVAELDQLLHY	14.00	7.78	8.34	11.67	11.81	11.64	13.72	9.74	11.37	12.97	11.03	10.50	8.60	8.08	9.81	13.92	9.73	8.28	9.71	9.83	11.21	11.01	5.53	11.98	12.27	11.86	11.59	9.40	11.67	12.34	12.00	15.05	11.60	12.15	6.58	9.37	9.15
QYDRVAELDQLLHY	10.80	6.79	10.95	11.59	12.13	10.57	11.50	10.76	11.10	11.90	11.59	9.10	9.46	8.98	10.23	11.09	8.57	8.54	10.45	9.24	11.99	10.55	8.58	2.30	11.01	10.22	9.18	10.10	10.13	10.46	10.51	12.72	11.03	10.55	5.91	8.73	10.20
RVLTEEEEKELENL	10.55	10.54	10.86	8.73	10.30	12.08	10.78	11.73	10.61	10.73	9.79	9.11	9.65	11.64	11.26	10.23	8.50	8.09	12.37	12.15	10.78	10.03	10.52	10.19	10.77	10.13	10.55	10.24	11.76	11.18	13.38	8.67	10.02	9.42	6.98	9.37	7.82
SQGGEPGAYLTPDT	9.87	9.98	9.82	8.47	9.92	8.24	8.88	9.77	9.16	10.20	8.32	9.74	9.78	9.33	9.91	8.62	10.37	9.44	10.46	8.88	9.87	9.44	8.97	11.78	10.91	7.16	8.61	9.60	10.07	8.68	9.97	9.60	9.29	8.23	9.15	9.02	10.00
SQGGEPGAYLTPDTREE	8.84	8.77	9.30	8.37	10.16	8.85	8.97	8.95	9.18	10.37	8.61	9.14	9.81	9.75	9.28	9.27	10.26	9.40	10.34	10.01	9.59	8.93	8.98	10.82	9.64	7.51	10.34	8.78	9.51	9.02	9.83	9.21	9.43	8.32	9.81	8.59	9.60
SSQEGNPPSEERGHTLQEPEEAKVGE	8.08	9.51	9.27	8.77	10.08	9.93	10.30	9.97	9.15	10.89	8.94	10.10	9.32	10.25	9.43	10.84	10.84	9.58	11.02	9.19	10.99	9.74	9.19	11.86	11.36	9.62	10.87	10.74	9.97	10.69	9.88	10.30	11.00	7.58	7.55	10.72	9.86
VAELDQLLHY	9.47	10.70	9.96	10.98	11.73	10.72	10.09	10.61	10.87	7.10	10.78	8.98	10.38	9.10	10.57	10.96	9.08	8.96	11.06	9.95	11.40	10.39	8.52	12.45	10.73	9.33	9.99	10.26	10.19	10.18	10.54	11.37	10.31	10.15	10.19	7.82	9.08
YDRVAELDQLLHY	10.91	9.61	10.21	12.18	12.35	11.10	10.65	10.19	11.75	11.86	12.03	9.26	9.16	7.41	11.03	13.06	7.60	6.96	10.97	10.45	12.01	11.56	6.87	10.03	11.70	10.36	8.67	6.97	9.07	11.61	10.73	12.92	11.67	12.06	10.38	9.76	9.53
Q(-17.03)EPGAPAAGM	10.21	9.17	9.60	9.26	10.22	11.08	10.46	9.66	9.07	10.48	9.11	9.58	9.42	9.43	9.76	9.99	10.58	10.27	10.19	9.77	10.45	9.58	9.25	10.69	11.43	8.76	7.22	11.00	10.98	11.02	11.05	10.82	9.97	9.47	10.94	9.34	10.02
SLDSPAGPA	10.46	9.22	9.81	9.23	9.78	9.97	10.25	6.53	10.22	9.54	9.14	9.67	9.34	7.39	5.97	10.12	12.28	11.44	5.91	7.76	10.38	9.51	10.34	11.85	12.58	9.43	12.58	12.31	10.51	8.82	11.33	8.79	10.94	5.82	9.90	11.27	11.68
SLDSPAGPAE	7.81	9.60	10.15	8.34	9.62	7.84	10.05	6.79	10.59	8.53	8.76	9.85	8.61	8.42	5.29	2.66	12.52	11.46	2.70	7.01	9.61	9.50	10.16	11.95	12.17	7.40	12.57	11.41	10.32	7.61	10.69	7.36	11.00	7.88	2.53	10.89	11.48
VKVGVNGF	11.11	9.78	8.29	10.44	11.97	11.08	8.47	8.27	11.55	9.52	9.50	10.52	9.08	8.87	11.31	9.51	10.54	6.92	9.81	9.80	10.47	10.96	9.62	11.05	13.05	10.28	10.34	7.73	8.30	9.56	8.20	10.27	11.39	12.18	10.54	10.49	9.89
Q(-17.03)HWSYGLRPG(-.98)	9.68	9.79	9.09	9.16	9.81	7.29	10.37	8.39	8.09	10.27	9.55	5.79	9.61	8.45	7.05	6.48	11.09	11.09	6.02	8.35	10.03	8.09	2.83	10.39	10.39	9.13	11.04	10.53	11.09	9.25	10.08	7.45	10.60	9.96	5.07	9.09	9.05
SVSTVLTSKYR	9.00	10.06	9.30	9.30	11.68	11.45	7.60	8.97	10.28	9.72	9.54	10.05	7.94	9.78	9.18	8.36	10.85	6.36	10.02	11.36	9.67	10.81	9.48	10.55	11.42	8.59	9.74	6.21	7.21	10.19	9.22	9.70	10.94	10.62	10.82	10.59	10.26
GVANALSHKYH	6.41	10.47	9.48	10.23	10.81	8.89	7.85	9.75	11.08	10.61	10.23	9.68	7.77	10.29	9.71	10.06	12.04	6.38	10.75	10.81	10.80	11.17	6.17	11.14	11.13	10.59	5.94	8.90	6.86	10.34	7.31	9.83	11.46	11.10	5.93	4.28	10.19
KVVAGVANALSHKYH	10.34	10.26	9.02	9.74	11.21	10.01	7.79	8.90	10.20	10.91	9.64	9.85	10.47	10.21	10.15	10.07	11.98	10.53	11.00	10.49	11.00	11.40	10.27	11.38	11.04	10.21	10.37	8.56	10.73	11.12	10.74	10.21	11.75	11.47	10.55	10.36	9.92
VANALSHKYH	6.08	10.56	9.69	11.12	9.82	8.45	8.45	9.79	12.10	10.11	11.06	9.43	8.16	10.30	9.72	10.79	11.18	5.63	9.81	9.66	10.70	11.65	5.84	10.92	10.86	10.97	10.64	8.31	5.75	10.62	7.59	10.15	10.16	10.75	6.02	9.87	10.57
GRTLYGFGG	7.72	10.03	8.19	10.28	11.08	11.00	7.94	9.31	10.79	7.40	10.12	6.87	9.52	9.61	10.77	8.63	11.23	8.45	10.30	10.47	10.55	10.77	6.99	10.62	12.60	8.83	8.03	9.25	9.20	9.72	8.64	9.90	11.65	10.71	10.56	6.59	5.81
RTLYGFGG	7.80	9.93	7.78	10.58	11.30	10.73	7.45	9.12	11.08	8.23	10.10	6.90	9.33	9.63	10.58	8.37	11.33	8.33	10.09	10.38	10.37	10.82	7.58	10.22	12.81	9.10	8.15	9.26	9.29	9.17	9.14	9.54	11.64	11.30	10.65	7.22	6.82
TLYGFGG	10.07	11.36	10.37	11.13	12.06	11.40	11.90	7.84	11.40	8.76	11.17	5.19	12.44	7.78	10.92	9.69	9.82	11.91	8.81	8.67	8.98	10.58	10.18	9.77	12.36	8.99	12.00	10.95	11.12	8.32	9.95	8.97	10.64	11.96	8.55	9.14	10.36
AADQDLGPEAPPEGVLGA	11.50	9.39	7.47	8.49	10.36	10.30	12.25	10.04	9.58	11.67	9.56	10.50	8.95	6.86	10.59	10.69	9.59	7.45	8.84	8.94	9.67	10.37	9.02	9.84	10.32	10.06	10.97	10.60	9.72	11.06	8.81	11.64	10.15	7.85	10.08	10.85	11.12
AADQDLGPEAPPEGVLGAL	14.30	9.90	7.68	12.64	13.06	11.91	14.44	11.53	11.57	13.92	11.22	11.23	7.36	9.12	10.97	13.23	9.69	9.94	9.96	8.18	11.94	10.77	10.00	12.35	12.98	12.59	11.57	10.20	11.11	13.13	9.88	15.14	11.73	13.31	9.79	11.25	11.39
LENPQPQAPARRLLPP	10.22	9.72	8.78	8.94	9.86	9.65	10.44	7.56	7.11	9.93	9.53	9.93	9.84	8.72	9.69	9.61	9.71	10.62	8.04	8.89	9.61	9.68	7.69	10.32	11.64	9.69	10.06	10.17	10.48	11.18	10.08	9.61	10.96	8.95	5.10	10.35	9.29
VKMALQQEGFD	10.49	9.69	8.67	7.90	11.20	12.48	11.37	9.30	10.15	9.61	8.01	10.75	9.28	8.39	10.34	10.04	10.20	10.34	8.98	9.49	9.51	10.00	8.60	8.63	11.84	9.97	9.72	10.92	10.10	12.13	9.18	10.84	10.37	8.53	8.59	9.78	7.77
GGFMRGL	10.55	10.72	9.21	10.55	8.95	7.47	9.08	10.02	9.87	9.58	10.08	10.63	12.31	6.50	7.06	7.98	9.71	11.26	9.36	8.87	9.63	11.76	6.15	10.60	12.05	9.20	10.11	10.97	6.84	9.34	9.39	9.79	12.15	11.00	10.93	10.31	9.71
Q(-17.03)LEDEAKELQ	9.60	10.53	9.47	8.84	10.69	9.68	7.34	7.21	8.26	10.35	9.83	10.30	10.94	7.73	10.73	7.72	11.01	12.00	7.47	8.69	10.62	10.37	10.22	10.32	12.36	9.58	11.18	11.54	10.53	10.06	10.42	9.51	10.80	9.18	6.52	10.74	12.24
SPQLEDE	8.11	10.45	9.76	10.66	9.08	7.30	10.43	7.44	9.07	9.23	11.06	1.91	10.97	7.52	10.67	7.24	11.27	11.56	7.40	8.43	9.93	10.64	10.19	11.01	12.86	10.52	11.54	12.31	11.45	8.66	10.90	8.67	10.64	7.22	5.66	10.47	11.33
SPQLEDEAKEL	8.12	10.02	9.35	9.09	10.23	9.06	10.17	6.55	7.83	10.29	9.73	3.55	10.75	7.96	7.86	7.17	11.19	11.93	7.85	8.73	10.01	9.99	3.31	11.11	12.46	9.29	10.89	10.81	10.53	9.82	10.82	8.86	11.34	9.08	7.10	10.25	11.17
SPQLEDEAKELQ	7.92	9.86	9.21	8.69	10.06	9.06	9.65	6.20	7.74	10.14	9.38	9.75	10.50	7.71	7.54	7.14	10.64	11.35	8.15	8.26	9.26	9.73	9.72	10.23	11.60	8.39	10.89	10.42	10.16	8.82	9.82	7.87	10.28	8.30	7.31	10.05	10.89
VGRPEWWMDYQ	10.41	10.11	9.80	10.27	11.34	9.32	10.62	7.95	8.66	10.75	11.37	5.48	11.64	7.04	7.01	10.48	10.43	10.89	7.19	9.05	10.20	10.66	10.04	10.66	12.45	10.11	10.91	10.76	11.24	10.22	11.08	10.19	10.93	9.31	9.76	9.89	10.89
YGGFM(+15.99)RF	8.27	10.86	8.75	10.11	11.59	10.07	10.22	7.17	8.51	8.96	10.51	10.31	10.34	7.87	7.35	5.05	11.52	11.49	8.42	10.18	9.98	10.54	10.54	10.15	13.32	9.71	11.89	11.97	12.14	9.06	11.80	9.59	11.22	8.54	8.74	10.37	11.14
YGGFM(+15.99)RGL	8.57	11.03	9.51	9.72	11.44	8.55	9.87	6.38	8.65	10.02	10.51	10.53	10.51	7.72	6.91	10.68	11.58	11.71	8.55	9.85	9.67	10.69	10.60	11.19	13.25	9.69	11.86	11.98	11.71	9.31	11.77	9.71	11.57	8.36	9.24	11.07	11.32
YGGFMRF	8.96	10.07	9.36	8.61	10.57	7.93	10.60	6.82	7.30	10.14	9.89	4.65	11.14	6.84	7.13	6.06	10.87	11.00	6.51	8.41	8.86	9.16	3.35	11.13	12.44	9.29	11.62	10.68	10.93	8.55	9.98	8.57	10.20	7.64	4.33	9.52	10.32
YGGFMRGL	8.83	10.07	9.52	8.98	9.79	8.44	10.26	6.02	6.71	9.93	9.54	3.76	10.76	7.46	7.38	5.69	10.64	11.05	7.16	8.43	10.09	10.25	3.22	10.62	12.46	9.48	11.49	10.90	10.72	9.69	10.30	8.22	10.74	9.01	3.51	9.50	9.29
EIGDEENSAKFPI(-.98)	9.42	10.46	9.34	8.24	11.66	10.22	10.67	9.46	9.70	9.74	9.37	7.32	10.45	9.69	6.56	7.73	10.72	11.13	9.58	8.71	10.58	9.76	6.18	8.26	11.00	9.67	11.16	10.86	10.26	8.98	9.54	9.77	11.51	8.44	10.74	10.74	10.34
DDGPYRMEHFRWGSPPKD	9.79	9.91	9.60	7.32	10.09	11.57	11.03	9.42	9.45	12.87	8.69	10.77	9.04	8.48	11.23	10.49	9.58	8.51	9.65	9.54	10.38	11.23	8.97	9.38	8.61	10.38	9.03	10.09	8.93	13.08	7.83	10.80	10.64	8.51	10.04	10.35	7.30
FKNAIIKNAH	8.07	11.14	9.30	7.80	8.84	9.07	11.47	10.57	11.03	10.87	9.53	10.27	8.75	9.35	11.21	11.39	11.10	4.31	9.09	9.42	9.70	10.25	10.94	10.07	10.33	10.99	10.20	11.29	8.82	12.22	9.85	10.35	10.88	8.04	9.94	10.83	9.98
FKNAIIKNAHKKGQ	10.06	9.27	9.14	8.49	8.43	6.28	10.37	10.86	10.84	9.93	12.15	8.03	9.44	10.48	9.88	12.16	7.21	9.39	9.95	9.73	11.37	9.31	11.24	10.82	7.84	10.79	4.73	8.84	8.74	10.32	8.17	9.82	6.49	8.04	8.86	8.89	5.56
GGPEFRDDGAEPGPRE(-.98)	7.53	8.68	9.01	9.71	11.31	10.69	7.80	9.02	10.43	10.60	9.17	7.78	7.35	8.80	10.96	10.85	10.06	7.34	10.22	11.90	11.78	10.14	6.17	10.73	11.73	10.47	8.46	8.00	9.06	10.18	8.45	9.24	10.75	7.52	8.01	6.67	6.87
GPEFRDDGAEPGPRE(-.98)	7.85	8.94	9.20	10.77	11.43	10.50	7.90	8.93	10.56	10.08	10.19	7.34	7.30	8.66	10.89	11.01	10.14	7.37	10.27	10.83	12.07	10.36	6.52	10.92	11.77	10.87	8.72	8.26	9.25	9.88	8.54	9.23	10.65	10.77	8.89	7.35	7.17
KYVMGHFRWD	9.38	9.84	8.93	6.31	11.15	11.79	10.30	9.80	9.68	10.86	7.04	10.22	7.76	8.31	10.86	9.73	9.64	8.51	8.90	10.03	9.88	10.73	9.69	8.01	9.34	9.83	8.39	10.00	7.42	11.67	8.95	10.26	10.07	8.17	9.64	10.18	9.09
KYVMGHFRWDRF(-.98)	9.03	9.25	7.70	6.69	10.60	11.63	10.16	9.43	9.62	10.54	7.78	10.20	7.80	7.85	11.25	10.27	9.73	8.44	9.15	8.81	10.47	10.82	9.24	9.32	7.61	10.10	9.02	8.79	8.19	12.44	7.08	10.09	11.41	8.48	9.98	10.36	6.33
PEFRDDGAEPGPRE(-.98)	6.67	9.22	8.82	10.96	11.10	10.35	8.31	9.56	10.94	10.46	9.47	7.68	7.33	8.91	10.85	10.73	10.59	7.81	10.65	10.95	12.09	10.13	9.80	10.68	11.24	10.55	9.38	7.53	9.77	9.65	8.38	8.41	10.21	10.02	7.56	7.92	7.44
RPVKVYPNGAEDES(+79.97)AESFPLEF	13.43	10.00	9.23	9.65	10.88	12.42	13.65	10.32	10.98	1.18	9.94	9.94	7.47	7.22	10.33	13.47	7.10	8.45	8.82	10.36	9.04	10.36	6.49	8.64	12.05	10.93	9.36	10.36	9.97	12.04	9.13	13.09	9.17	8.69	6.63	8.30	9.74
RPVKVYPNGAEDESAESFPLEF	13.04	8.69	8.40	9.37	11.63	13.45	11.73	11.30	10.87	11.42	9.96	11.47	9.00	6.74	11.60	11.37	7.50	6.77	9.43	10.91	8.79	8.28	6.91	10.45	10.96	9.19	8.73	9.75	8.32	10.79	8.72	12.94	7.24	8.31	7.03	8.77	9.75
S(+42.01)YSMEHFRWGKPV(-.98)	10.50	10.01	9.22	7.76	10.85	12.32	11.21	9.57	9.55	11.96	7.80	10.49	7.79	8.09	11.08	10.36	9.13	7.66	9.28	10.77	9.76	10.92	8.47	10.13	10.63	9.68	9.59	10.15	9.12	12.62	7.73	11.56	11.09	8.25	4.33	10.11	5.40
SQTPLVTLF	10.98	10.39	10.66	11.26	11.25	11.21	10.77	11.17	9.90	10.55	10.94	8.99	9.92	9.96	10.57	9.85	9.51	9.61	10.82	10.27	10.73	7.92	10.09	9.40	11.78	8.53	11.70	10.90	10.86	10.40	10.22	11.48	8.78	8.44	11.95	4.18	9.98
SQTPLVTLFKNAIIKNAHKKGQ	6.88	9.56	8.13	10.84	10.09	10.91	10.04	9.10	10.32	14.70	11.30	7.34	6.58	5.13	7.78	9.99	10.14	9.62	10.47	9.80	11.38	9.10	9.77	12.25	9.63	9.01	10.22	9.93	10.59	8.75	7.88	10.66	9.50	9.10	9.04	9.69	10.82
SYSMEHF	7.98	8.83	10.34	9.49	10.21	10.24	8.90	10.46	9.37	10.16	9.77	9.22	9.09	10.02	10.40	8.60	10.50	10.37	10.12	10.31	10.70	9.22	10.50	10.69	10.87	7.95	9.56	10.64	10.21	9.59	10.73	10.56	9.42	10.42	11.24	8.59	10.50
SYSMEHFRWGKPV(-.98)	9.77	8.97	9.33	7.34	10.25	10.23	9.63	9.41	8.14	11.99	8.51	10.15	8.79	8.54	10.39	10.51	10.10	9.75	10.39	9.42	11.20	10.03	8.62	9.97	8.91	8.91	9.18	8.70	9.15	9.81	8.56	9.24	9.28	10.47	2.00	9.72	9.52
TPLVTLFKNAIIKNAHKKGQ	5.29	5.89	8.34	12.55	11.97	10.75	8.05	5.89	11.30	9.99	9.26	10.37	9.19	5.69	11.33	11.68	7.19	2.47	9.23	9.24	11.43	9.83	9.47	11.20	10.32	11.73	7.02	6.85	6.95	9.21	6.97	10.32	10.41	12.97	10.63	10.74	9.89
Y(+42.01)GGFM(+15.99)TSEKSQTPLVTLFKNAIIKNAH	13.79	11.21	7.31	7.23	10.67	13.14	13.79	8.84	11.43	11.27	8.55	10.24	10.02	9.55	11.01	11.38	9.38	9.00	9.02	9.74	8.77	11.78	7.31	9.81	11.83	11.50	11.19	11.52	9.54	12.51	8.35	13.69	10.65	7.62	3.17	10.97	9.32
Y(+42.01)GGFMTSEKSQTPLVTLFKNAIIKNA	14.53	8.58	9.13	1.49	11.68	14.27	14.53	8.60	10.27	12.43	8.79	11.08	5.78	8.92	8.17	12.64	7.09	7.35	7.06	8.75	7.22	10.59	9.13	8.77	11.10	11.16	9.20	10.07	8.34	13.37	7.10	14.38	10.10	7.09	9.76	10.60	8.97
Y(+42.01)GGFMTSEKSQTPLVTLFKNAIIKNAH	13.37	10.36	8.35	7.54	10.62	12.61	13.54	9.77	10.62	12.13	8.53	10.89	6.87	7.00	10.96	11.58	9.20	8.33	8.98	9.29	9.86	11.55	8.32	10.33	11.80	11.59	11.17	11.56	9.46	13.63	7.62	13.08	10.33	7.26	0.34	10.22	6.53
Y(+42.01)GGFMTSEKSQTPLVTLFKNAIIKNAHKKGQ	10.70	11.20	8.85	8.34	6.96	11.49	10.78	10.94	9.53	11.60	9.84	9.07	8.09	8.69	11.03	9.96	8.84	9.35	10.24	7.84	10.25	11.13	9.66	9.14	9.29	9.56	10.50	9.84	8.18	12.42	8.34	11.71	9.49	7.10	2.50	9.99	8.01
YGGFM(+15.99)TSEKSQTPLVTLFKNAIIKNAHKKGQ	9.74	10.61	9.41	11.69	9.95	8.71	10.58	10.91	10.24	11.07	11.85	7.56	6.80	10.37	8.55	10.69	11.51	11.45	11.31	9.79	10.90	8.22	10.17	12.35	9.12	10.23	12.76	11.00	11.92	8.44	10.31	10.50	3.55	9.21	6.44	8.86	12.67
YGGFMTSE	9.45	8.85	9.86	9.18	10.35	10.79	9.51	10.37	8.88	12.02	9.13	8.98	9.11	10.30	10.12	9.56	10.14	10.67	10.28	9.71	10.57	9.59	10.29	10.45	11.03	8.00	11.45	10.96	10.23	10.55	10.31	11.69	9.09	9.87	11.38	8.38	10.06
YGGFMTSEKSQTPLVTLFKNAIIKNA	11.53	10.59	9.34	9.43	11.70	12.68	11.78	10.45	10.75	12.55	9.69	11.41	7.53	7.97	10.72	12.93	9.67	10.65	10.66	9.63	11.43	11.93	8.74	10.78	13.53	11.50	11.32	11.16	11.60	12.37	11.09	12.75	10.41	10.53	10.99	11.22	10.04
YGGFMTSEKSQTPLVTLFKNAIIKNAH	9.77	11.31	9.97	8.89	11.36	11.82	9.96	9.99	9.42	10.83	7.78	9.82	9.82	8.25	11.18	12.08	6.43	9.25	9.47	8.80	10.22	11.49	9.17	7.66	12.17	10.08	9.05	11.05	10.03	11.28	9.51	11.88	9.05	8.66	10.05	11.07	8.98
YGGFMTSEKSQTPLVTLFKNAIIKNAHKKGQ	8.62	9.70	9.82	10.83	9.92	9.01	10.44	11.40	11.18	11.55	11.83	9.24	9.75	9.75	9.51	9.68	10.38	9.96	10.18	8.77	11.40	8.53	10.03	11.79	8.28	9.52	11.60	9.68	10.19	8.53	9.68	9.96	8.37	7.53	5.91	8.42	11.34
ATLESVFQELGKLTGPSSQKRE	11.69	9.82	9.87	10.89	9.60	11.48	10.80	9.91	10.36	11.69	11.44	10.98	11.15	8.70	11.50	10.22	10.23	11.30	12.37	10.20	10.91	11.05	10.75	12.06	10.48	7.58	11.16	6.60	11.18	7.76	12.93	9.08	8.75	9.47	9.50	9.51	10.71
MYEENSRDNPF	8.25	9.15	9.89	8.57	9.60	9.73	9.23	9.93	9.98	10.47	9.10	9.14	9.94	9.55	9.71	8.29	10.35	9.47	10.67	9.25	10.40	9.51	10.00	10.45	11.49	7.62	9.34	10.13	11.28	8.76	11.62	9.77	9.04	9.36	8.38	8.27	8.71
TNEIVEEQYTPQSL	11.18	10.38	11.75	8.73	10.29	11.78	8.18	10.49	11.49	11.77	8.13	8.71	10.43	10.21	11.31	10.97	9.46	5.39	12.54	11.79	10.52	10.72	11.09	11.62	11.56	9.36	11.78	8.10	12.17	9.55	13.71	10.08	9.92	9.40	11.53	9.60	7.09
FPTLGGSQDKSLHN	9.38	8.53	9.12	7.79	9.22	9.24	8.54	9.20	8.54	10.61	8.43	9.45	9.39	9.18	10.14	8.99	10.94	9.76	9.49	9.86	10.48	9.77	9.86	10.53	11.02	8.61	10.98	9.65	10.26	9.90	10.67	8.76	10.18	8.77	10.34	9.16	9.20
TLGGSQDKSLHN	3.23	9.21	9.65	9.13	6.86	7.32	8.45	10.04	11.08	10.19	9.48	8.49	9.79	9.73	10.88	10.02	10.97	8.20	9.10	8.53	11.65	9.80	10.99	11.19	10.39	9.22	10.16	9.23	10.52	8.34	10.70	7.84	9.47	8.79	8.76	5.68	8.61
SVNPYLQGQRL	8.92	8.68	9.52	5.27	6.44	9.67	4.50	10.18	9.00	8.96	5.73	8.53	9.25	9.66	11.22	7.62	10.60	6.97	10.03	9.70	10.83	9.05	9.97	10.40	10.20	8.12	10.05	7.18	9.99	8.78	11.79	6.31	9.79	10.27	10.16	8.17	4.95
SVNPYLQGQRLDNVV	11.90	9.25	9.73	11.00	11.28	12.70	11.87	9.55	10.75	11.43	10.66	10.58	4.12	8.67	9.60	11.56	8.75	8.99	10.30	10.19	11.34	9.42	7.92	11.23	12.27	11.14	11.44	10.64	11.16	12.10	9.95	12.53	9.51	10.86	10.59	9.25	9.29
SVNPYLQGQRLDNVVA	11.40	9.34	9.39	10.53	11.03	11.59	11.30	9.49	10.04	11.59	10.16	10.59	9.05	8.66	10.78	10.69	10.07	9.28	10.14	9.84	11.65	9.93	8.66	11.88	11.96	10.39	11.52	10.34	11.04	11.63	10.31	11.70	10.43	10.36	10.59	8.42	9.23
SANSNPAM(+15.99)APRE	10.30	11.23	10.37	10.15	9.61	4.40	9.84	4.95	9.15	9.43	10.82	10.55	9.91	5.50	2.44	5.16	12.00	11.50	8.08	7.34	10.82	10.23	10.65	12.02	11.91	9.89	10.62	11.22	12.67	6.97	11.61	4.98	11.06	8.22	10.40	11.31	10.86
SANSNPAMAPRE	7.76	10.06	10.00	9.30	9.31	5.83	10.39	6.02	8.34	10.83	9.52	5.02	9.72	5.61	7.55	6.40	11.68	10.93	6.64	8.75	10.79	9.28	3.94	11.64	11.83	10.02	11.03	11.56	11.25	8.27	10.33	6.44	10.68	8.89	6.99	9.70	10.05

EigenMS normalization

After imputation using knn (k = 7), the data was normalized using EigenMS to account for sample-to-sample variability.

EigenMS normalization preserves the treatment group differences in the data by estimating treatment effects with an ANOVA model, then uses singular value decomposition on the model residual matrix to identify and remove the bias.

Code

plot_grid(hist_eigenms, qq_eigenms, pca_eigenms, nrow = 1)

Code

norm_d1_mod = norm_d1[-1,] %>%
        as.data.frame() %>%
        rownames_to_column(., var = "rowname") %>%
        mutate(across(-rowname, as.numeric)) %>%
        column_to_rownames(., var = "rowname")
#write.csv(norm_d1_mod, file = "after_norm_data.csv")

#rownames(norm_d1_mod) == feature_name$variable
rownames(norm_d1_mod) = feature_name$name

# display table
norm_d1_mod %>% 
        mutate(across(everything(), round, digit=2)) %>%
        knitr::kable(., digits = 3, "html") %>%
        kableExtra::kable_styling(bootstrap_options = c("striped", "hover", "condensed", "responsive")) %>%
        kableExtra::scroll_box(width = "100%", height = "500px")

	Pituitary_53	Pituitary_31	Pituitary_38	Pituitary_63	Pituitary_55	Pituitary_51	Pituitary_50	Pituitary_34	Pituitary_59	Pituitary_43	Pituitary_66	Pituitary_26	Pituitary_19	Pituitary_24	Pituitary_41	Pituitary_57	Pituitary_47	Pituitary_52	Pituitary_23	Pituitary_58	Pituitary_44	Pituitary_37	Pituitary_39	Pituitary_46	Pituitary_65	Pituitary_56	Pituitary_48	Pituitary_35	Pituitary_61	Pituitary_49	Pituitary_60	Pituitary_54	Pituitary_45	Pituitary_42	Pituitary_62	Pituitary_25	Pituitary_64
AGAPEPAEPMELAKPSA	12.03	9.92	10.59	6.66	13.59	12.75	11.30	4.23	9.24	8.51	9.25	6.88	5.87	9.45	11.35	12.35	10.60	12.74	10.60	5.48	11.92	11.52	8.28	9.18	7.20	10.74	8.94	9.67	7.42	12.04	8.01	11.65	7.99	9.38	8.49	10.34	8.72
LAGAPEPA	8.74	8.21	11.22	7.57	10.60	11.67	10.90	8.25	5.13	8.47	7.36	8.96	8.93	9.75	9.28	13.71	7.93	12.96	7.11	6.94	3.42	10.97	6.13	6.70	4.96	12.27	5.28	9.74	7.79	11.91	6.93	11.21	5.78	6.42	6.62	9.56	6.13
AFGPRGPGLQL	9.41	8.75	9.50	10.25	10.05	8.95	10.09	9.84	5.16	10.93	9.56	9.36	9.49	9.49	10.30	9.34	8.80	8.27	9.02	9.16	10.07	8.75	9.14	10.24	10.00	9.55	10.10	9.61	9.78	9.40	10.00	11.15	9.76	11.03	9.74	7.76	10.62
AAMDLEL	10.02	5.12	9.15	8.77	12.37	4.89	9.73	8.42	9.16	9.86	9.13	8.67	9.42	8.39	10.08	8.06	8.92	7.61	7.92	8.71	8.24	10.39	7.90	9.86	6.63	9.10	10.09	7.72	9.65	9.35	11.95	10.82	9.62	8.36	8.73	8.86	8.48
AELDQLLHY	11.35	10.00	8.39	10.05	10.59	10.57	6.93	10.05	10.39	12.29	11.51	10.24	11.44	7.58	8.88	11.45	10.26	10.77	10.73	9.16	10.86	9.77	9.26	11.88	9.31	8.24	10.81	8.56	8.71	9.71	10.19	9.24	10.41	9.91	10.30	10.25	10.85
EDREEASLQDRQYASHHT	5.79	9.90	10.00	9.34	7.84	9.26	8.85	9.71	10.18	10.73	9.88	10.69	8.86	9.31	8.10	8.81	10.34	8.57	9.89	7.52	10.76	8.81	9.30	10.70	9.83	9.52	8.39	10.51	9.64	8.77	8.78	7.27	10.42	10.50	10.68	9.63	9.00
EDREEASLQDRQYASHHTTE	7.59	8.86	9.49	8.62	8.32	9.47	7.86	8.84	9.01	9.64	7.86	10.58	9.47	9.56	9.92	8.74	10.30	8.79	9.52	8.35	10.05	8.83	10.35	10.32	8.80	9.58	8.96	9.39	8.63	9.23	9.22	8.22	11.57	10.32	7.69	10.29	9.36
EEASLQDRQYASHHT	6.62	8.95	10.04	9.16	6.53	9.23	7.27	9.60	9.96	9.47	8.93	11.00	9.50	9.39	8.46	7.33	9.84	8.13	9.68	8.66	9.68	8.98	10.31	9.62	9.60	7.89	8.08	9.63	8.81	7.09	8.73	6.18	10.84	9.47	8.62	10.11	8.48
EEASLQDRQYASHHTTE	7.78	8.84	9.79	9.34	8.40	9.23	8.30	9.55	9.61	9.61	8.71	10.37	9.62	9.41	9.43	8.81	10.39	8.56	9.67	8.47	10.06	9.15	9.98	10.27	9.56	9.51	8.91	9.74	9.17	8.69	9.25	8.47	10.78	9.65	9.13	10.02	9.81
EESEESEEGTTSEVT	8.68	8.78	10.01	10.04	8.81	7.66	9.38	10.79	10.46	8.51	11.16	10.10	10.38	12.05	7.42	8.46	10.80	7.12	10.16	10.65	10.45	9.16	9.59	9.78	10.67	7.47	7.85	11.57	9.65	8.94	9.16	9.24	9.19	8.05	11.59	10.87	10.84
ELDQLLHY	10.41	10.35	9.90	10.52	10.17	8.96	9.37	9.70	10.68	10.09	10.87	9.38	11.37	9.32	10.87	9.56	10.14	7.92	10.04	9.67	10.26	9.40	9.59	10.64	10.38	8.74	10.83	10.06	9.42	6.36	8.52	8.79	9.83	7.24	11.49	8.27	10.62
FLGEAYHHVPE	9.69	7.94	9.43	8.38	8.91	8.75	8.50	8.74	6.19	10.75	8.37	9.27	10.08	9.17	10.19	9.98	8.75	8.91	9.97	9.46	8.82	9.85	9.57	10.30	8.69	7.11	9.47	8.07	8.86	8.65	8.80	6.79	9.99	11.03	7.18	10.58	9.29
FLGEAYHHVPES	9.45	8.86	9.87	9.44	9.46	8.61	9.58	10.22	8.27	10.70	10.05	10.31	10.17	8.91	9.70	9.11	9.35	9.90	10.27	10.04	9.43	9.93	8.67	9.37	8.90	6.60	8.98	8.81	9.71	8.38	8.33	7.77	9.92	9.00	9.40	9.57	8.90
FLGEAYHHVPESQ	9.08	9.09	9.76	9.22	9.62	6.16	9.40	10.38	8.14	10.72	8.97	10.16	9.59	9.65	9.48	9.27	9.54	8.11	10.32	9.94	9.38	9.01	8.20	10.67	9.41	7.55	9.71	9.61	10.06	6.20	8.87	8.97	10.05	9.98	10.01	9.38	9.16
FLGEAYHHVPESQRD	8.80	8.49	9.29	9.13	9.20	8.48	8.56	9.19	7.76	10.60	8.13	9.63	9.45	9.44	10.78	8.86	10.51	8.60	9.71	9.74	9.73	9.19	9.27	11.11	8.53	7.93	11.06	9.29	9.11	8.96	9.48	9.42	10.85	11.50	8.11	9.71	9.42
FLGEAYHHVPESQRDKA	6.34	8.71	9.37	10.77	8.46	8.05	9.22	9.44	8.26	10.80	8.34	9.48	9.13	9.34	9.91	9.19	10.66	8.14	9.90	9.10	9.82	10.63	9.27	11.60	8.40	8.14	10.50	9.30	9.35	8.72	9.90	9.26	11.37	9.88	9.86	10.20	10.53
GEAYHHVPESQRD	7.08	9.58	10.64	9.32	8.00	6.33	8.99	9.68	8.78	10.27	10.13	9.75	11.24	9.58	10.45	10.06	9.78	7.45	9.93	8.61	10.18	10.61	9.99	10.87	9.74	8.71	9.71	10.37	9.60	8.18	10.43	8.95	9.89	9.97	10.42	10.54	10.09
GEAYHHVPESQRDKA	4.07	8.91	9.59	10.47	8.03	7.69	8.81	8.87	8.77	10.32	8.45	10.26	10.61	8.85	10.36	9.36	10.37	7.57	9.03	8.12	9.44	11.28	10.60	11.27	8.73	7.48	9.86	8.76	9.27	8.09	10.33	8.16	10.97	9.73	9.50	9.85	10.57
GVLFNPYYDPLQWKSSRFE	12.04	9.65	9.99	10.13	10.85	9.38	11.15	8.50	9.32	11.34	9.19	9.34	11.52	8.25	9.90	12.20	10.73	9.45	10.78	8.90	11.05	10.93	9.88	12.18	9.65	11.50	10.85	9.43	11.70	9.98	10.74	12.84	9.43	12.16	10.60	9.42	10.34
LDQLLHY	8.70	9.65	9.54	9.12	9.43	7.88	8.90	10.63	10.23	11.99	9.59	8.87	10.85	9.20	10.45	9.11	9.85	8.82	9.98	9.69	10.67	10.38	9.93	11.04	9.75	7.34	10.80	9.74	9.37	7.77	9.43	8.77	10.33	10.65	10.36	9.62	10.62
LGEAYHHVPESQRD	8.32	8.88	9.23	9.97	8.51	4.69	8.80	9.49	7.77	11.79	8.60	9.19	9.13	9.38	10.25	9.29	10.58	7.91	9.68	9.02	10.28	9.23	8.11	11.37	8.35	6.73	11.41	9.65	9.14	7.73	9.29	9.87	10.64	10.72	9.61	10.10	7.29
LGVLFNPYYDPLQWKSSRFE	13.00	8.91	9.74	11.88	12.54	10.53	12.89	9.96	10.38	11.90	10.71	10.34	9.68	8.80	8.57	13.26	9.98	11.40	10.60	9.73	11.57	10.54	8.25	13.21	9.76	12.75	11.87	10.56	11.78	12.45	12.28	13.39	11.10	12.10	12.45	8.70	10.66
Q(-17.03)KIAEKFSQ	9.45	10.29	10.13	9.61	6.08	8.40	4.90	10.38	9.52	9.59	9.64	9.36	8.59	10.01	9.78	7.29	9.02	8.25	8.81	9.29	9.17	9.28	9.04	9.43	9.53	9.02	10.45	8.79	9.33	7.88	8.94	6.62	9.99	9.49	8.52	9.97	8.73
Q(-17.03)KIAEKFSQR(-.98)	4.91	9.13	9.29	7.98	4.80	7.82	4.31	9.23	9.30	9.45	8.06	9.95	8.80	9.14	10.07	7.68	10.62	7.92	8.40	9.22	10.25	9.71	10.11	10.22	8.79	10.01	10.74	8.59	8.36	8.58	9.58	7.05	11.07	11.04	7.20	11.32	9.47
Q(-17.03)YDRVAELDQLLHY	14.27	7.09	8.88	11.42	11.77	10.99	12.91	9.23	9.90	11.57	10.22	9.55	9.65	9.18	9.47	13.26	10.80	11.62	9.53	9.64	10.93	9.45	7.72	11.88	10.12	12.22	12.00	8.77	11.75	11.68	12.73	13.90	11.55	12.87	10.65	9.01	9.12
QYDRVAELDQLLHY	11.24	6.66	11.19	9.79	11.84	10.65	10.94	9.98	9.52	11.51	10.16	8.88	9.99	9.58	9.74	10.39	9.52	10.83	9.67	9.01	11.51	9.60	10.03	2.58	10.48	9.87	10.04	10.10	11.11	10.37	11.88	11.92	11.14	9.68	9.37	8.79	9.94
RVLTEEEEKELENL	11.51	10.70	10.73	9.60	10.24	11.18	11.39	10.56	10.24	9.89	10.20	9.83	10.47	11.12	10.82	8.98	8.79	9.39	10.66	11.42	9.73	9.86	10.74	8.93	9.44	10.60	11.03	10.59	10.65	10.47	12.11	8.25	10.75	11.51	9.01	10.81	9.30
SQGGEPGAYLTPDT	10.08	9.72	9.75	9.40	9.92	7.92	8.97	9.75	9.52	9.82	8.80	10.06	9.81	9.25	10.25	8.65	10.13	9.29	10.41	8.85	9.67	9.40	9.17	11.11	9.94	7.52	8.31	9.29	9.15	8.41	9.07	9.65	9.38	9.59	9.79	9.29	10.42
SQGGEPGAYLTPDTREE	9.05	8.78	9.17	9.20	10.16	8.56	9.32	8.84	9.59	10.23	9.18	9.55	9.82	9.44	9.41	9.17	10.01	9.05	10.06	9.87	9.38	9.11	8.75	10.27	9.20	7.77	10.19	8.81	8.80	8.81	8.95	9.34	9.57	9.34	9.59	9.03	10.16
SSQEGNPPSEERGHTLQEPEEAKVGE	8.47	9.05	9.26	9.67	10.11	9.51	10.10	9.74	9.19	10.07	9.17	10.40	9.61	10.35	9.86	10.70	10.74	9.98	10.83	9.18	10.55	9.46	9.89	10.91	9.76	10.10	10.57	10.16	8.78	10.30	8.91	10.14	11.27	9.49	9.78	11.08	10.23
VAELDQLLHY	10.36	10.70	9.74	10.43	10.78	10.45	10.90	9.92	10.65	7.60	10.81	10.23	9.93	8.99	10.45	10.49	9.31	9.73	9.74	9.20	10.81	9.70	9.41	12.12	9.70	9.06	10.82	10.68	10.49	9.81	10.49	11.24	9.89	10.06	11.15	8.73	10.38
YDRVAELDQLLHY	11.52	9.63	10.61	11.04	11.43	10.25	11.27	9.25	10.15	11.78	11.56	8.91	9.69	8.15	9.67	11.82	9.14	10.99	9.58	9.29	11.77	9.62	8.80	10.73	10.38	10.13	10.33	7.75	10.74	10.76	12.28	11.74	10.85	10.67	11.50	9.97	10.87
Q(-17.03)EPGAPAAGM	10.35	9.01	9.61	10.39	10.13	10.54	10.85	9.68	9.45	10.21	9.86	9.71	9.48	9.39	9.81	9.92	10.49	10.52	10.10	9.53	10.41	9.35	9.42	10.26	10.47	9.15	7.10	10.94	10.33	10.59	10.28	10.79	9.80	10.53	10.56	9.55	10.75
SLDSPAGPA	9.68	7.84	9.89	11.29	10.25	10.26	8.98	8.21	11.77	8.80	9.83	9.36	8.65	8.09	8.08	12.12	11.00	9.11	8.60	9.09	11.12	9.69	10.91	11.11	10.85	10.19	10.47	10.19	8.47	9.00	9.94	9.46	10.57	8.20	10.78	9.89	10.32
SLDSPAGPAE	7.16	7.35	10.10	8.96	9.84	9.22	7.43	8.82	11.98	7.81	8.07	10.25	7.13	9.72	8.90	5.73	10.75	7.95	6.21	9.21	10.14	9.54	12.02	10.92	9.71	7.75	9.85	7.94	8.05	8.43	9.75	8.31	10.64	10.33	7.55	9.02	8.61
VKVGVNGF	11.05	9.49	8.70	10.16	11.21	10.46	8.92	8.39	10.97	9.74	9.65	9.87	8.99	9.77	10.63	9.36	11.44	9.41	9.86	9.16	10.99	9.31	11.18	12.00	11.89	10.20	11.14	8.00	9.66	8.95	9.45	9.60	10.05	10.81	9.64	9.86	10.66
Q(-17.03)HWSYGLRPG(-.98)	10.03	8.08	8.76	8.68	8.77	8.32	9.28	9.43	9.06	10.71	8.99	7.58	7.58	9.36	9.80	8.63	9.88	8.76	7.50	9.29	10.00	7.40	5.17	9.55	7.74	8.88	9.87	8.37	10.07	9.66	9.61	8.21	9.59	11.04	9.29	8.45	8.10
SVSTVLTSKYR	9.10	9.96	9.59	9.63	11.00	10.59	8.43	8.91	9.97	9.88	10.16	9.67	7.97	10.30	8.39	7.97	11.62	8.60	9.69	10.54	10.00	9.44	10.57	11.12	10.26	8.68	10.54	6.71	8.12	9.43	9.83	9.14	9.83	9.89	9.36	10.40	11.54
GVANALSHKYH	7.82	10.35	9.42	8.18	9.37	8.57	8.45	8.32	9.57	10.94	9.21	11.01	7.62	10.66	9.14	8.93	13.10	9.24	8.54	9.66	9.75	9.49	8.46	11.01	9.39	9.98	7.67	9.44	8.07	9.77	8.37	9.02	10.94	10.24	9.66	5.40	11.66
KVVAGVANALSHKYH	10.63	10.13	8.90	10.11	10.56	9.73	8.52	8.95	10.65	11.41	10.27	10.55	9.87	10.19	10.27	10.20	11.87	10.51	10.73	10.08	11.01	10.95	10.69	11.27	10.24	10.19	10.59	8.73	10.77	10.82	10.43	10.34	11.07	11.50	9.63	10.62	10.88
VANALSHKYH	6.82	10.16	9.94	8.94	8.89	8.41	8.12	8.94	10.35	10.00	9.57	9.61	8.30	11.24	9.30	10.12	12.34	8.64	8.78	9.12	10.18	9.93	8.28	11.31	9.43	10.45	11.89	8.28	7.16	10.30	9.31	9.19	9.66	9.51	10.06	9.95	10.80
GRTLYGFGG	9.04	9.76	7.88	9.84	8.93	10.33	9.81	8.72	10.96	8.78	10.88	9.07	8.07	9.73	10.71	8.36	11.55	9.76	8.58	9.01	10.08	8.99	8.99	10.43	10.21	8.38	9.44	9.96	10.06	8.85	8.67	9.80	9.90	10.06	10.31	7.70	8.54
RTLYGFGG	8.94	9.64	7.53	10.04	9.30	10.19	9.09	8.68	11.24	9.58	10.73	8.83	7.92	9.84	10.56	8.25	11.62	9.51	8.68	9.09	10.04	9.14	9.49	10.19	10.67	8.64	9.42	9.86	10.23	8.42	9.34	9.46	9.92	10.46	10.32	8.07	9.19
TLYGFGG	10.40	10.47	10.04	10.00	10.92	12.31	11.63	8.47	12.04	9.87	10.69	6.80	10.59	8.28	12.55	11.14	9.08	10.19	9.47	9.04	9.03	10.04	11.66	9.67	11.17	8.39	11.68	9.99	11.29	8.73	10.21	9.61	9.57	11.41	10.27	8.82	9.98
AADQDLGPEAPPEGVLGA	10.70	8.85	8.19	9.71	11.15	9.59	11.33	10.55	9.03	10.15	9.68	8.30	10.18	7.87	10.22	10.66	10.24	9.61	10.15	9.27	10.35	9.57	9.81	10.21	9.30	10.84	10.48	9.89	9.46	10.57	9.19	10.85	10.06	8.62	10.54	9.55	10.44
AADQDLGPEAPPEGVLGAL	13.86	9.63	8.33	12.24	13.66	11.60	13.35	11.42	10.15	12.63	10.30	9.34	8.68	10.09	10.28	12.73	10.73	12.66	10.46	8.38	12.18	9.91	11.02	12.94	12.56	12.86	11.73	9.82	11.66	12.90	11.12	14.14	11.95	13.05	11.87	10.32	10.48
LENPQPQAPARRLLPP	10.27	8.56	8.87	9.43	9.59	9.71	9.61	8.14	7.39	9.39	9.42	10.16	9.37	9.55	10.99	10.56	9.35	10.48	9.08	9.40	9.68	8.91	9.37	9.74	9.42	10.00	9.29	8.73	9.56	11.04	9.68	9.63	10.47	10.33	7.94	9.74	8.82
VKMALQQEGFD	10.40	8.98	8.99	9.32	11.01	11.65	11.39	9.67	10.31	8.92	8.77	10.22	9.45	9.10	10.49	10.22	10.42	11.74	9.56	9.33	9.83	8.88	9.89	8.43	9.73	10.58	9.48	10.34	9.57	11.40	8.80	10.44	9.64	9.70	8.86	9.28	8.36
GGFMRGL	10.32	9.47	9.31	9.91	7.61	7.92	8.88	11.32	10.47	10.67	10.04	11.12	10.43	7.83	8.36	9.77	9.44	10.95	10.97	9.12	10.53	10.17	8.48	11.38	10.17	8.78	9.87	9.91	7.78	9.36	10.45	10.05	9.98	9.61	10.90	8.92	9.45
Q(-17.03)LEDEAKELQ	9.31	9.27	9.48	9.51	10.50	10.10	6.36	8.35	9.13	10.13	9.88	10.56	9.99	8.54	12.50	9.33	10.18	10.57	9.26	9.57	11.04	9.99	11.50	9.88	10.57	9.83	9.92	9.85	9.49	10.22	9.90	9.95	10.14	10.38	8.35	9.73	11.29
SPQLEDE	8.76	8.99	9.16	10.09	8.10	8.52	9.50	8.16	10.24	9.82	10.48	4.40	8.86	7.88	13.58	9.23	9.74	8.30	8.34	9.37	9.48	10.58	11.90	9.68	10.65	10.19	10.34	10.31	9.98	9.25	9.86	9.71	10.14	8.76	10.24	10.48	10.51
SPQLEDEAKEL	8.69	8.05	8.94	9.11	8.55	9.78	9.72	7.77	9.33	11.18	9.86	5.96	8.08	8.97	10.85	9.59	9.84	9.58	9.27	9.33	10.07	8.77	6.20	10.15	8.80	9.08	9.83	8.67	9.46	9.91	10.02	9.76	9.55	10.31	10.44	9.72	11.10
SPQLEDEAKELQ	7.62	8.61	9.13	9.45	9.88	9.60	8.70	7.43	8.88	10.07	9.53	10.23	9.36	8.37	9.51	8.93	9.59	9.32	10.02	9.23	9.68	9.57	10.78	9.69	9.95	8.61	9.48	8.69	8.95	9.09	9.08	8.50	9.65	9.62	8.90	9.10	9.90
VGRPEWWMDYQ	10.55	8.96	9.55	10.51	10.51	9.83	10.30	8.89	9.78	11.31	11.58	6.77	10.02	7.58	8.87	12.10	9.46	9.04	8.35	9.55	10.40	10.17	11.43	10.15	10.56	10.03	10.06	9.43	10.51	10.37	10.50	10.86	9.89	10.06	11.11	9.41	10.65
YGGFM(+15.99)RF	8.04	9.29	8.63	11.12	11.02	10.56	9.45	8.68	10.01	9.15	10.97	11.15	8.69	8.73	9.70	7.24	10.29	9.26	10.54	11.10	10.52	10.01	12.12	9.54	10.91	9.95	10.38	10.00	10.84	9.21	10.91	10.36	10.05	9.97	10.15	9.30	10.43
YGGFM(+15.99)RGL	8.32	9.67	9.45	10.92	11.12	8.84	9.15	7.68	9.97	9.93	11.03	11.08	9.30	8.44	8.90	12.50	10.52	9.87	10.43	10.68	10.14	10.30	11.85	10.54	11.10	10.05	10.44	10.22	10.34	9.37	10.80	10.32	10.71	9.98	10.48	10.14	10.70
YGGFMRF	9.41	7.74	9.04	8.62	9.21	8.86	9.37	8.17	8.63	10.46	9.57	6.77	8.70	8.13	10.57	8.73	9.42	8.56	8.44	9.46	8.87	7.98	6.59	9.92	8.41	9.23	10.10	7.86	9.38	8.77	9.12	9.40	8.73	9.51	9.52	8.67	9.49
YGGFMRGL	9.58	7.73	9.10	8.68	8.00	9.33	9.35	7.20	8.01	10.57	9.23	6.40	8.01	8.73	10.82	8.29	9.25	8.78	8.67	9.22	9.95	8.79	6.82	9.39	8.10	9.25	10.29	8.25	9.41	9.81	9.52	9.05	9.00	10.65	8.87	8.92	8.91
EIGDEENSAKFPI(-.98)	9.43	9.32	9.19	8.77	10.72	10.47	10.61	10.55	10.89	10.39	9.92	8.33	8.81	10.36	8.18	9.35	9.91	9.73	10.86	9.04	11.04	8.97	7.68	8.02	8.97	9.66	10.43	9.69	9.78	8.94	9.11	10.34	10.10	8.91	11.02	10.07	10.41
DDGPYRMEHFRWGSPPKD	9.53	9.77	10.13	8.69	10.35	10.41	11.22	9.39	9.01	11.86	9.40	9.34	10.10	9.03	10.32	9.86	10.39	11.09	9.79	9.11	10.73	10.25	9.59	9.60	7.45	11.06	9.25	10.21	8.91	12.23	7.89	9.99	10.32	9.14	9.26	9.92	8.01
FKNAIIKNAH	7.75	10.97	9.91	9.01	9.45	8.00	11.18	10.41	10.23	9.39	9.83	8.49	10.25	9.97	10.29	10.60	12.02	7.16	9.28	9.19	9.92	9.40	11.57	10.19	9.27	11.75	10.30	11.21	8.62	11.46	9.99	9.39	11.01	8.97	10.37	10.33	10.20
FKNAIIKNAHKKGQ	10.84	10.34	8.90	7.47	8.75	6.20	10.81	9.28	9.79	9.73	11.47	8.53	10.47	9.46	8.73	10.47	7.72	10.16	7.74	9.14	10.24	10.10	10.12	10.37	9.24	10.55	5.73	10.03	8.89	10.37	8.08	9.50	7.94	8.03	10.14	10.49	6.15
GGPEFRDDGAEPGPRE(-.98)	8.95	8.60	8.92	9.79	10.07	9.54	9.13	7.70	9.78	10.55	9.61	9.12	7.43	8.78	10.31	9.56	10.83	10.01	7.96	10.46	10.78	8.63	7.88	9.95	9.26	10.56	9.84	8.67	9.11	9.06	8.03	8.55	10.24	8.36	9.88	8.16	9.42
GPEFRDDGAEPGPRE(-.98)	9.13	8.95	9.07	9.91	10.00	9.80	9.15	7.77	9.80	10.60	10.19	8.73	6.98	8.73	10.26	9.98	10.93	9.70	8.18	9.50	11.26	8.88	8.21	10.68	10.00	10.53	10.27	9.05	10.07	9.09	8.87	8.66	9.89	10.36	10.33	8.57	9.30
KYVMGHFRWD	9.08	9.67	9.42	7.81	11.35	10.66	10.57	9.92	9.46	9.99	7.91	8.90	8.62	8.82	10.09	9.26	10.31	10.77	9.16	9.62	10.31	9.79	10.25	8.23	8.17	10.51	8.52	10.09	7.35	10.84	8.90	9.58	9.63	8.80	8.46	9.70	9.83
KYVMGHFRWDRF(-.98)	8.84	9.29	8.16	7.90	10.64	10.47	10.73	9.35	9.27	9.93	8.64	9.00	8.64	8.25	10.17	9.57	10.55	10.94	9.03	8.16	10.84	9.84	9.70	9.71	6.72	10.64	9.49	9.25	8.49	11.58	7.25	9.37	10.90	8.63	8.30	10.09	7.36
PEFRDDGAEPGPRE(-.98)	7.93	9.52	8.60	10.12	10.22	9.85	9.27	8.16	10.12	10.64	9.24	9.02	7.45	8.55	10.20	9.45	11.20	9.50	8.38	9.86	11.00	9.43	10.63	10.10	10.28	10.29	10.69	8.36	10.09	9.11	8.31	7.95	10.30	10.15	9.55	9.46	9.13
RPVKVYPNGAEDES(+79.97)AESFPLEF	13.73	9.42	9.62	9.05	10.38	11.95	13.33	10.00	9.86	0.60	9.34	9.48	7.84	8.22	10.04	13.09	8.02	11.21	8.60	10.00	8.94	8.76	8.57	8.87	10.24	10.97	9.97	10.04	10.60	11.49	10.11	12.22	8.58	8.48	9.23	7.96	10.02
RPVKVYPNGAEDESAESFPLEF	12.98	8.52	9.15	9.72	12.26	12.35	11.06	10.57	9.09	9.49	9.44	9.46	11.04	7.59	10.32	10.01	8.97	10.97	9.07	10.57	8.64	7.06	8.16	10.59	9.70	9.85	9.28	9.68	8.46	9.96	9.46	11.47	7.70	9.06	9.64	8.42	9.80
S(+42.01)YSMEHFRWGKPV(-.98)	10.87	9.28	9.77	8.82	10.72	10.99	11.03	9.08	8.54	10.37	7.99	9.59	8.94	9.07	10.58	9.53	10.17	11.34	8.93	10.22	9.52	9.13	10.64	9.69	7.66	10.46	9.97	9.68	8.65	11.45	7.66	10.37	10.79	9.91	7.21	10.01	6.35
SQTPLVTLF	11.46	10.70	10.47	11.78	11.29	10.86	11.30	10.57	9.99	10.37	11.34	9.57	10.21	9.36	10.33	9.22	9.45	9.58	9.83	9.90	10.18	8.26	9.60	8.72	11.67	8.71	11.90	11.30	10.22	10.16	9.33	11.49	9.30	9.44	12.12	5.12	10.80
SQTPLVTLFKNAIIKNAHKKGQ	7.03	9.41	8.24	10.67	10.55	10.92	9.28	8.69	9.66	13.77	10.62	7.03	7.33	5.26	7.83	9.61	10.33	10.28	10.29	10.05	10.94	9.16	10.08	11.78	9.26	9.21	10.08	9.50	10.04	8.75	7.78	10.32	10.28	10.06	11.86	9.82	10.28
SYSMEHF	8.22	9.05	10.15	10.02	10.17	10.04	9.41	10.25	9.74	10.28	10.28	9.73	9.03	9.55	10.35	8.39	10.28	9.90	9.62	10.08	10.47	9.54	10.00	10.28	10.89	8.05	9.59	10.93	9.77	9.47	10.02	10.75	9.59	11.00	10.56	9.16	11.12
SYSMEHFRWGKPV(-.98)	10.25	8.65	9.41	6.81	10.49	10.23	8.77	8.70	7.17	10.92	7.51	10.18	9.55	8.78	10.56	9.98	10.41	10.85	9.84	9.57	10.44	9.80	9.52	9.25	7.97	9.06	9.22	8.14	8.53	9.74	8.49	8.74	10.11	11.67	6.39	10.09	9.11
TPLVTLFKNAIIKNAHKKGQ	5.49	6.33	8.92	10.62	11.12	10.03	8.78	5.03	9.27	10.32	8.51	9.26	9.87	6.54	9.21	10.26	9.20	7.14	7.93	7.96	11.62	7.83	11.03	12.91	10.26	11.18	9.15	8.23	9.89	8.53	9.68	8.95	9.35	9.79	9.96	10.52	10.96
Y(+42.01)GGFM(+15.99)TSEKSQTPLVTLFKNAIIKNAH	13.60	9.79	8.07	8.14	10.95	12.35	12.29	9.10	10.40	9.16	8.08	8.63	11.12	11.23	11.38	11.47	10.16	12.15	10.17	10.10	8.93	9.96	10.06	9.48	8.54	12.38	10.69	9.87	8.75	11.72	8.63	12.56	10.37	9.53	7.81	9.77	8.70
Y(+42.01)GGFMTSEKSQTPLVTLFKNAIIKNA	13.17	7.85	10.45	5.03	13.14	12.19	13.71	9.35	9.62	9.53	10.08	6.99	8.33	10.42	6.95	12.09	8.38	11.82	8.96	8.83	8.44	9.07	10.13	9.22	8.91	12.97	8.49	9.31	7.51	11.96	7.03	12.85	9.86	9.21	8.71	8.67	8.88
Y(+42.01)GGFMTSEKSQTPLVTLFKNAIIKNAH	13.53	9.06	9.16	8.37	11.15	11.55	11.86	9.34	8.94	9.33	7.70	9.20	8.74	8.56	10.98	10.91	10.32	12.48	9.37	9.49	9.47	9.72	11.20	9.54	8.24	12.64	10.92	9.99	8.32	12.65	7.72	11.56	10.74	9.94	6.90	9.54	6.06
Y(+42.01)GGFMTSEKSQTPLVTLFKNAIIKNAHKKGQ	10.98	10.77	9.21	8.44	7.65	11.07	9.58	10.21	8.21	9.69	8.90	8.19	9.64	9.19	10.89	9.12	9.48	11.53	9.93	8.05	9.57	10.68	10.73	8.36	7.92	10.12	10.40	9.09	7.25	12.06	8.22	10.81	10.57	8.98	7.43	10.08	7.44
YGGFM(+15.99)TSEKSQTPLVTLFKNAIIKNAHKKGQ	10.19	11.01	8.97	10.59	11.23	9.81	9.08	9.91	9.62	10.03	10.15	8.21	7.75	9.37	9.21	10.17	10.89	9.45	10.54	10.76	9.53	10.28	8.73	10.92	10.70	10.15	12.09	10.46	10.33	9.42	9.37	10.80	6.21	11.01	11.61	9.98	10.91
YGGFMTSE	9.50	8.95	9.81	10.40	10.61	10.33	9.85	10.29	9.29	11.65	9.88	8.99	9.42	9.95	10.04	9.32	9.96	10.51	10.10	9.58	10.46	9.87	9.78	9.93	10.74	8.43	11.22	11.05	9.39	10.27	9.32	11.74	9.38	11.09	10.69	8.74	10.61
YGGFMTSEKSQTPLVTLFKNAIIKNA	11.31	10.04	9.86	11.01	11.95	11.57	11.61	10.58	10.42	11.26	10.33	10.11	8.48	8.73	10.33	12.60	10.30	13.05	11.12	9.42	11.72	10.81	9.86	10.67	11.61	12.32	11.20	10.74	11.08	11.50	10.86	11.98	10.08	11.86	11.49	10.66	10.53
YGGFMTSEKSQTPLVTLFKNAIIKNAH	9.74	10.83	10.42	10.11	11.13	10.72	10.24	10.06	9.15	10.09	8.50	8.96	10.34	8.99	10.71	11.76	7.11	11.70	9.62	8.30	10.52	10.11	10.50	7.74	10.21	10.66	9.29	10.92	10.00	10.36	9.54	11.17	8.29	9.35	9.90	10.67	9.92
YGGFMTSEKSQTPLVTLFKNAIIKNAHKKGQ	9.08	9.87	9.57	9.78	10.70	9.72	9.28	10.52	10.41	10.68	10.38	9.70	10.52	9.23	9.92	9.16	10.16	9.16	9.45	9.37	10.31	9.67	9.49	10.76	9.00	9.43	11.31	9.23	9.18	9.12	9.23	9.98	10.22	8.80	10.50	9.24	10.15
ATLESVFQELGKLTGPSSQKRE	11.86	10.14	9.76	11.32	10.20	11.39	10.67	9.42	10.27	11.08	11.41	11.01	11.84	8.15	11.36	9.67	10.10	11.00	11.83	10.28	10.41	11.82	9.91	11.39	10.91	7.84	10.99	6.70	10.32	7.80	12.13	9.09	9.79	10.65	10.39	10.13	10.66
MYEENSRDNPF	8.81	9.00	9.68	9.45	9.51	9.35	9.51	9.54	10.29	10.08	9.57	9.95	10.01	9.22	10.04	8.03	10.07	9.27	10.04	9.05	9.83	9.58	10.14	9.42	10.38	7.96	9.21	9.96	10.12	8.44	10.38	9.84	9.37	11.11	9.65	9.06	9.47
TNEIVEEQYTPQSL	11.93	10.78	11.78	10.97	10.11	10.00	9.82	9.50	11.55	11.06	9.85	8.90	11.40	9.64	10.11	9.43	10.01	7.59	10.78	10.45	9.95	10.14	11.01	10.71	9.99	10.19	12.53	9.12	11.26	8.23	12.16	9.51	9.99	11.42	10.17	10.97	9.78
FPTLGGSQDKSLHN	9.39	8.34	9.08	9.42	9.28	8.70	8.94	9.43	9.31	10.32	9.50	9.59	9.37	9.01	10.40	9.13	10.58	9.38	9.68	9.76	10.54	9.80	9.70	9.94	10.11	9.13	10.55	9.47	9.23	9.51	9.48	8.93	10.06	10.24	9.41	9.32	9.89
TLGGSQDKSLHN	4.15	9.63	9.24	9.25	6.84	7.13	8.97	8.97	11.01	9.99	9.55	9.76	10.10	8.81	10.78	9.12	10.79	7.86	7.44	8.11	10.55	10.42	10.41	10.02	10.29	9.29	10.48	9.68	9.56	8.21	9.50	7.92	10.46	10.22	10.39	7.25	9.53
SVNPYLQGQRL	9.56	9.13	9.30	7.66	6.05	8.24	6.46	9.64	9.92	9.01	7.88	9.35	9.39	8.79	10.51	6.71	10.55	7.48	8.62	8.51	10.48	8.94	9.32	9.52	9.12	8.77	10.50	8.22	9.01	7.74	9.90	6.34	9.56	11.99	7.13	9.53	7.65
SVNPYLQGQRLDNVV	11.80	9.16	10.22	11.75	11.79	11.82	11.58	9.17	9.83	10.10	10.67	9.22	5.50	9.12	8.75	10.68	9.60	11.57	10.14	9.94	11.31	8.72	8.51	11.21	11.39	11.74	11.66	10.62	10.99	11.46	10.10	11.64	9.82	11.70	11.70	9.05	9.49
SVNPYLQGQRLDNVVA	11.48	9.09	9.61	11.33	11.05	10.94	11.37	9.36	9.82	10.93	10.51	10.17	9.55	8.96	10.54	10.36	10.42	10.63	10.05	9.59	11.62	9.31	9.31	11.66	10.80	10.82	11.59	10.21	10.69	11.10	10.05	11.28	10.31	11.22	11.14	8.39	9.75
SANSNPAM(+15.99)APRE	9.53	9.52	10.10	10.89	9.40	5.78	8.35	7.24	11.31	9.96	10.97	11.41	7.66	6.31	5.69	8.40	10.01	7.01	11.40	9.13	11.75	10.58	11.54	11.56	10.53	9.92	8.24	8.64	11.14	7.86	10.71	6.40	10.04	9.36	11.44	9.60	8.78
SANSNPAMAPRE	8.03	8.08	9.65	9.23	8.10	6.83	9.37	7.41	9.75	11.41	9.30	6.97	7.32	6.66	10.67	8.98	10.23	8.19	8.53	9.77	10.95	8.45	6.50	10.74	8.70	9.85	9.60	9.14	10.02	8.63	9.62	7.37	9.29	10.19	10.77	8.86	9.18

Peak intensity plots

EigenMS removes sample-to-sample variability.

Intensity after log2 transformation

Intensity after EigenMS normalization

Univariate analysis

Univariate statistic analysis will proceed for each pairs of comparison. For each pair,

Step 1, differential analysis;

• t-test

• Wilcoxon rank-sum test

For all statistical tests, the Benjamini-Hochberg (BH) procedure was applied to correct for multiple hypothesis testing.

Step 2, fold change, both linear and log2;

Step 3, regulation shown in volcano plot.

• A feature is considered unregulated when FC(log2) > 0.6 (or FC(lin) > 1.5);

• A feature is considered downregulated when FC(log2) < -0.6 (or FC(lin) < -1.5).

Code

d5 <- rbind.data.frame(labels_d1, norm_d1_mod) 
d5_mod <- t(d5) %>%
  as.data.frame(.) %>%
  rownames_to_column(., var = "rowname") %>%
  as_tibble() %>%
  dplyr::rename(., Label = class) %>%
  mutate(across(-c(Label,rowname), as.numeric)) %>%
  #rename_with(str_trim)
  column_to_rownames(., var = "rowname")

#colnames(d5_mod)[-1] == feature_name$variable
colnames(d5_mod)[-1] = feature_name$name

grps = as.factor(t(labels_d1))
trt= levels(grps)

combo_list = do.call("c", lapply(seq_along(time), function(i) combn(trt, 2, FUN = list))) # top 10 are unique combo
#combo_list = combo_list[1:10]

feature_meta = read.csv(file = "feature_meta_pituitary.csv")

folder = paste(getwd(), "/pituitary/", sep = "")

August vs IBA

Full stat

Code

contrast = combo_list[[1]]

df = d5_mod %>%
    dplyr::filter(., Label %in% contrast ) %>%
    droplevels()

temp_d <- df %>%
        mutate(across(-Label, ~ 2^(.) )) %>% # d3_mod is log-scale number, this step will de-log and change data back to linear number
        group_by(Label) %>%
        summarise(across(everything(), mean))

print(paste("Order of Fold Change:", temp_d$Label[1], "over", temp_d$Label[2], sep = " "))

[1] "Order of Fold Change: August over IBA"

Code

uni_res <- do_univariate(df)

uni_res <- uni_res %>%
        rowwise() %>%
        mutate(padj = min(c(BHT, BHW))) %>%
        # get the lowest padj
        ungroup() %>%
        mutate(`-log10padj` = -log10(padj))

uni_res_annotation = uni_res %>% left_join(., feature_meta, by = c("variable" = "Peptide"))

uni_res_annotation %>% 
        knitr::kable(., digits = 3, "html") %>%
        kableExtra::kable_styling(bootstrap_options = c("striped", "hover", "condensed", "responsive")) %>%
        kableExtra::scroll_box(width = "100%", height = "500px")

variable	pT	BHT	pW	BHW	FC(lin)	FC(log2)	padj	-log10padj	Accession	Name	GN
AADQDLGPEAPPEGVLGA	0.012	0.035	0.012	0.035	0.482	-1.053	0.035	1.454	tr|I3NGD6|I3NGD6_ICTTR	ProSAAS	PCSK1N
AADQDLGPEAPPEGVLGAL	0.000	0.001	0.000	0.003	0.119	-3.072	0.001	2.966	tr|I3NGD6|I3NGD6_ICTTR	ProSAAS	PCSK1N
AAMDLEL	0.592	0.704	0.456	0.579	0.430	-1.218	0.579	0.237	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
AELDQLLHY	0.561	0.676	0.529	0.654	0.635	-0.654	0.654	0.184	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
AFGPRGPGLQL	0.134	0.233	0.224	0.351	0.609	-0.716	0.233	0.633	tr|I3MMW2|I3MMW2_ICTTR	Chromogranin-A	CHGA
AGAPEPAEPMELAKPSA	0.005	0.019	0.002	0.008	0.269	-1.893	0.008	2.102	tr|I3MXH3|I3MXH3_ICTTR	Vasopressin-neurophysin 2-copeptin-like	AVP
ATLESVFQELGKLTGPSSQKRE	0.934	0.954	0.955	0.955	1.184	0.244	0.954	0.021	tr|I3N6V4|I3N6V4_ICTTR	Secretogranin-2	SCG2
DDGPYRMEHFRWGSPPKD	0.053	0.122	0.113	0.213	0.546	-0.872	0.122	0.915	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
EDREEASLQDRQYASHHT	0.010	0.030	0.008	0.025	2.267	1.181	0.025	1.609	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
EDREEASLQDRQYASHHTTE	0.011	0.031	0.018	0.047	1.935	0.952	0.031	1.503	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
EEASLQDRQYASHHT	0.000	0.001	0.001	0.006	4.550	2.186	0.001	3.194	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
EEASLQDRQYASHHTTE	0.000	0.002	0.002	0.008	2.043	1.030	0.002	2.626	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
EESEESEEGTTSEVT	0.002	0.011	0.002	0.008	3.628	1.859	0.008	2.102	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
EIGDEENSAKFPI(-.98)	0.331	0.472	0.456	0.579	0.799	-0.323	0.472	0.326	tr|I3MQ49|I3MQ49_ICTTR	Pro-melanin concentrating hormone	PMCH
ELDQLLHY	0.229	0.341	0.328	0.467	1.289	0.367	0.341	0.467	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
FKNAIIKNAH	0.366	0.514	0.607	0.713	1.108	0.148	0.514	0.289	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
FKNAIIKNAHKKGQ	0.795	0.839	0.328	0.467	0.790	-0.340	0.467	0.331	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
FLGEAYHHVPE	0.065	0.139	0.113	0.213	1.831	0.872	0.139	0.857	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
FLGEAYHHVPES	0.108	0.198	0.113	0.213	1.515	0.599	0.198	0.702	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
FLGEAYHHVPESQ	0.059	0.128	0.088	0.188	1.770	0.824	0.128	0.893	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
FLGEAYHHVPESQRD	0.003	0.012	0.008	0.025	1.530	0.613	0.012	1.924	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
FLGEAYHHVPESQRDKA	0.003	0.013	0.000	0.003	2.422	1.276	0.003	2.467	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
FPTLGGSQDKSLHN	0.031	0.082	0.050	0.119	1.278	0.353	0.082	1.089	tr|I3MPL6|I3MPL6_ICTTR	Secretogranin-3	SCG3
GEAYHHVPESQRD	0.000	0.003	0.002	0.008	3.477	1.798	0.003	2.580	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
GEAYHHVPESQRDKA	0.005	0.019	0.003	0.011	3.990	1.996	0.011	1.942	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
GGFMRGL	0.533	0.658	0.607	0.713	1.167	0.223	0.658	0.181	tr|I3MGR7|I3MGR7_ICTTR	Proenkephalin-A	PENK
GGPEFRDDGAEPGPRE(-.98)	0.001	0.007	0.003	0.011	0.437	-1.195	0.007	2.138	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
GPEFRDDGAEPGPRE(-.98)	0.004	0.017	0.003	0.011	0.514	-0.961	0.011	1.942	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
GRTLYGFGG	0.197	0.313	0.328	0.467	0.699	-0.516	0.313	0.504	tr|I3N7M0|I3N7M0_ICTTR	Histone H4	LOC106145262
GVANALSHKYH	0.441	0.583	0.864	0.902	0.840	-0.251	0.583	0.235	tr|A0A287CUK4|A0A287CUK4_ICTTR	Hemoglobin subunit beta-S/F-like	LOC101967841
GVLFNPYYDPLQWKSSRFE	0.070	0.146	0.145	0.252	0.412	-1.280	0.146	0.836	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
KVVAGVANALSHKYH	0.893	0.933	0.456	0.579	1.042	0.059	0.579	0.237	tr|A0A287CUK4|A0A287CUK4_ICTTR	Hemoglobin subunit beta-S/F-like	LOC101967841
KYVMGHFRWD	0.097	0.188	0.181	0.304	0.623	-0.682	0.188	0.727	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
KYVMGHFRWDRF(-.98)	0.035	0.088	0.066	0.152	0.483	-1.051	0.088	1.055	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
LAGAPEPA	0.036	0.090	0.036	0.089	0.235	-2.089	0.089	1.051	tr|I3MXH3|I3MXH3_ICTTR	Vasopressin-neurophysin 2-copeptin-like	AVP
LDQLLHY	0.000	0.003	0.000	0.003	2.350	1.233	0.003	2.580	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
LENPQPQAPARRLLPP	0.001	0.006	0.005	0.018	0.491	-1.025	0.006	2.224	tr|I3NGD6|I3NGD6_ICTTR	ProSAAS	PCSK1N
LGEAYHHVPESQRD	0.039	0.092	0.066	0.152	1.922	0.943	0.092	1.036	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
LGVLFNPYYDPLQWKSSRFE	0.000	0.003	0.003	0.011	0.156	-2.680	0.003	2.580	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
MYEENSRDNPF	0.016	0.045	0.018	0.047	1.605	0.682	0.045	1.346	tr|I3N6V4|I3N6V4_ICTTR	Secretogranin-3	SCG2
PEFRDDGAEPGPRE(-.98)	0.758	0.829	0.776	0.848	0.935	-0.097	0.829	0.082	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
Q(-17.03)EPGAPAAGM	0.141	0.236	0.181	0.304	0.737	-0.440	0.236	0.627	tr|I3MIV6|I3MIV6_ICTTR	Carboxypeptidase E	CPE
Q(-17.03)HWSYGLRPG(-.98)	0.055	0.123	0.018	0.047	0.466	-1.103	0.047	1.326	tr|I3MCA0|I3MCA0_ICTTR	Progonadoliberin-1	GNRH1
Q(-17.03)KIAEKFSQ	0.007	0.022	0.012	0.035	2.551	1.351	0.022	1.652	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
Q(-17.03)KIAEKFSQR(-.98)	0.006	0.021	0.008	0.025	3.821	1.934	0.021	1.683	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
Q(-17.03)LEDEAKELQ	0.474	0.611	0.864	0.902	1.237	0.307	0.611	0.214	tr|I3MGR7|I3MGR7_ICTTR	Proenkephalin-A	PENK
Q(-17.03)YDRVAELDQLLHY	0.000	0.000	0.000	0.003	0.066	-3.912	0.000	3.769	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
QYDRVAELDQLLHY	0.030	0.081	0.026	0.067	0.516	-0.956	0.067	1.175	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
RPVKVYPNGAEDES(+79.97)AESFPLEF	0.000	0.001	0.000	0.003	0.099	-3.337	0.001	3.194	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
RPVKVYPNGAEDESAESFPLEF	0.000	0.003	0.000	0.003	0.143	-2.801	0.003	2.580	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
RTLYGFGG	0.494	0.627	0.689	0.771	0.867	-0.206	0.627	0.202	tr|I3N7M0|I3N7M0_ICTTR	Histone H4	LOC106145262
RVLTEEEEKELENL	0.416	0.559	0.607	0.713	1.027	0.039	0.559	0.253	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
S(+42.01)YSMEHFRWGKPV(-.98)	0.007	0.022	0.012	0.035	0.486	-1.041	0.022	1.652	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
SANSNPAM(+15.99)APRE	0.005	0.019	0.005	0.018	4.009	2.003	0.018	1.744	tr|I3MP38|I3MP38_ICTTR	Somatostatin	SST
SANSNPAMAPRE	0.790	0.839	0.689	0.771	1.103	0.142	0.771	0.113	tr|I3MP38|I3MP38_ICTTR	Somatostatin	SST
SLDSPAGPA	0.938	0.954	0.955	0.955	0.832	-0.265	0.954	0.021	tr|I3M5J1|I3M5J1_ICTTR	Corticoliberin	CRH
SLDSPAGPAE	0.236	0.346	0.224	0.351	3.406	1.768	0.346	0.460	tr|I3M5J1|I3M5J1_ICTTR	Corticoliberin	CRH
SPQLEDE	0.096	0.188	0.088	0.188	2.513	1.329	0.188	0.727	tr|I3MGR7|I3MGR7_ICTTR	Proenkephalin-A	PENK
SPQLEDEAKEL	0.176	0.286	0.145	0.252	0.671	-0.575	0.252	0.599	tr|I3MGR7|I3MGR7_ICTTR	Proenkephalin-A	PENK
SPQLEDEAKELQ	0.129	0.228	0.145	0.252	1.600	0.678	0.228	0.641	tr|I3MGR7|I3MGR7_ICTTR	Proenkephalin-A	PENK
SQGGEPGAYLTPDT	0.104	0.195	0.224	0.351	1.377	0.462	0.195	0.710	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
SQGGEPGAYLTPDTREE	0.666	0.764	0.955	0.955	1.055	0.078	0.764	0.117	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
SQTPLVTLF	0.204	0.320	0.272	0.406	0.577	-0.794	0.320	0.495	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
SQTPLVTLFKNAIIKNAHKKGQ	0.944	0.954	0.689	0.771	0.853	-0.230	0.771	0.113	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
SSQEGNPPSEERGHTLQEPEEAKVGE	0.640	0.742	0.864	0.902	1.171	0.228	0.742	0.129	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
SVNPYLQGQRL	0.006	0.021	0.036	0.089	2.152	1.106	0.021	1.683	tr|I3MEQ9|I3MEQ9_ICTTR	Neuroendocrine protein 7B2	SCG5
SVNPYLQGQRLDNVV	0.002	0.009	0.000	0.003	0.281	-1.833	0.003	2.467	tr|I3MEQ9|I3MEQ9_ICTTR	Neuroendocrine protein 7B3	SCG5
SVNPYLQGQRLDNVVA	0.001	0.007	0.000	0.003	0.366	-1.448	0.003	2.467	tr|I3MEQ9|I3MEQ9_ICTTR	Neuroendocrine protein 7B4	SCG5
SVSTVLTSKYR	0.786	0.839	0.388	0.529	1.166	0.221	0.529	0.276	tr|A0A287D1R9|A0A287D1R9_ICTTR	Hemoglobin subunit alpha	LOC101956162
SYSMEHF	0.219	0.332	0.388	0.529	1.302	0.380	0.332	0.479	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
SYSMEHFRWGKPV(-.98)	0.414	0.559	0.388	0.529	0.769	-0.379	0.529	0.276	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
TLGGSQDKSLHN	0.110	0.199	0.113	0.213	2.468	1.303	0.199	0.701	tr|I3MPL6|I3MPL6_ICTTR	Secretogranin-4	SCG3
TLYGFGG	0.539	0.658	0.529	0.654	0.684	-0.548	0.654	0.184	tr|I3N7M0|I3N7M0_ICTTR	Histone H4	LOC106145262
TNEIVEEQYTPQSL	0.098	0.188	0.113	0.213	1.656	0.728	0.188	0.727	tr|I3N6V4|I3N6V4_ICTTR	Secretogranin-4	SCG2
TPLVTLFKNAIIKNAHKKGQ	0.978	0.978	0.776	0.848	0.869	-0.203	0.848	0.072	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
VAELDQLLHY	0.084	0.171	0.088	0.188	0.633	-0.659	0.171	0.766	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
VANALSHKYH	0.708	0.783	0.955	0.955	1.007	0.010	0.783	0.106	tr|A0A287CUK4|A0A287CUK4_ICTTR	Hemoglobin subunit beta-S/F-like	LOC101967841
VGRPEWWMDYQ	0.209	0.322	0.145	0.252	0.668	-0.582	0.252	0.599	tr|I3MGR7|I3MGR7_ICTTR	Proenkephalin-A	PENK
VKMALQQEGFD	0.000	0.003	0.001	0.006	0.353	-1.501	0.003	2.580	tr|I3MEG8|I3MEG8_ICTTR	Proprotein convertase subtilisin/kexin type 2	PCSK2
VKVGVNGF	0.446	0.583	0.456	0.579	0.808	-0.307	0.579	0.237	tr|A0A287D887|A0A287D887_ICTTR	Glyceraldehyde-3-phosphate dehydrogenase	GAPDH
Y(+42.01)GGFM(+15.99)TSEKSQTPLVTLFKNAIIKNAH	0.000	0.001	0.000	0.003	0.170	-2.554	0.001	2.857	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
Y(+42.01)GGFMTSEKSQTPLVTLFKNAIIKNA	0.000	0.000	0.000	0.003	0.108	-3.208	0.000	3.769	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
Y(+42.01)GGFMTSEKSQTPLVTLFKNAIIKNAH	0.000	0.002	0.002	0.008	0.209	-2.255	0.002	2.654	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
Y(+42.01)GGFMTSEKSQTPLVTLFKNAIIKNAHKKGQ	0.502	0.629	0.272	0.406	0.619	-0.692	0.406	0.391	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
YDRVAELDQLLHY	0.006	0.021	0.002	0.008	0.336	-1.574	0.008	2.102	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
YGGFM(+15.99)RF	0.329	0.472	0.456	0.579	1.650	0.723	0.472	0.326	tr|I3MGR7|I3MGR7_ICTTR	Proenkephalin-A	PENK
YGGFM(+15.99)RGL	0.411	0.559	0.272	0.406	1.058	0.081	0.406	0.391	tr|I3MGR7|I3MGR7_ICTTR	Proenkephalin-A	PENK
YGGFM(+15.99)TSEKSQTPLVTLFKNAIIKNAHKKGQ	0.612	0.720	0.864	0.902	0.871	-0.199	0.720	0.143	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
YGGFMRF	0.039	0.092	0.008	0.025	0.569	-0.815	0.025	1.609	tr|I3MGR7|I3MGR7_ICTTR	Proenkephalin-A	PENK
YGGFMRGL	0.150	0.247	0.113	0.213	0.674	-0.569	0.213	0.671	tr|I3MGR7|I3MGR7_ICTTR	Proenkephalin-A	PENK
YGGFMTSE	0.704	0.783	0.689	0.771	0.823	-0.281	0.771	0.113	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
YGGFMTSEKSQTPLVTLFKNAIIKNA	0.000	0.003	0.000	0.003	0.350	-1.513	0.003	2.580	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
YGGFMTSEKSQTPLVTLFKNAIIKNAH	0.139	0.236	0.224	0.351	0.688	-0.539	0.236	0.627	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
YGGFMTSEKSQTPLVTLFKNAIIKNAHKKGQ	0.707	0.783	0.607	0.713	0.889	-0.169	0.713	0.147	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC

Significant stat

Code

# get DE features only tibble
uni_res_filt <- uni_res %>%
        dplyr::filter(BHT < FDR & BHW < FDR) %>%
        mutate(status = case_when(`FC(log2)` < - 0.6 ~ "Down",
                                  `FC(log2)` > 0.6 ~ "Up") )

if (nrow(uni_res_filt) == 0){
        print(paste("There is no significant differentiated features between",
                    temp_d$Label[2], "and", temp_d$Label[1], sep = " "))
} else {
        uni_res_filt = uni_res_filt %>%
                left_join(., feature_meta, by = c("variable" = "Peptide"))
        
        uni_res_filt %>% 
                knitr::kable(., digits = 3, "html") %>%
                kableExtra::kable_styling(bootstrap_options = c("striped", "hover", "condensed", "responsive")) %>%
                kableExtra::scroll_box(width = "100%", height = "500px")
}

variable	pT	BHT	pW	BHW	FC(lin)	FC(log2)	padj	-log10padj	status	Accession	Name	GN
AADQDLGPEAPPEGVLGA	0.012	0.035	0.012	0.035	0.482	-1.053	0.035	1.454	Down	tr|I3NGD6|I3NGD6_ICTTR	ProSAAS	PCSK1N
AADQDLGPEAPPEGVLGAL	0.000	0.001	0.000	0.003	0.119	-3.072	0.001	2.966	Down	tr|I3NGD6|I3NGD6_ICTTR	ProSAAS	PCSK1N
AGAPEPAEPMELAKPSA	0.005	0.019	0.002	0.008	0.269	-1.893	0.008	2.102	Down	tr|I3MXH3|I3MXH3_ICTTR	Vasopressin-neurophysin 2-copeptin-like	AVP
EDREEASLQDRQYASHHT	0.010	0.030	0.008	0.025	2.267	1.181	0.025	1.609	Up	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
EDREEASLQDRQYASHHTTE	0.011	0.031	0.018	0.047	1.935	0.952	0.031	1.503	Up	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
EEASLQDRQYASHHT	0.000	0.001	0.001	0.006	4.550	2.186	0.001	3.194	Up	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
EEASLQDRQYASHHTTE	0.000	0.002	0.002	0.008	2.043	1.030	0.002	2.626	Up	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
EESEESEEGTTSEVT	0.002	0.011	0.002	0.008	3.628	1.859	0.008	2.102	Up	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
FLGEAYHHVPESQRD	0.003	0.012	0.008	0.025	1.530	0.613	0.012	1.924	Up	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
FLGEAYHHVPESQRDKA	0.003	0.013	0.000	0.003	2.422	1.276	0.003	2.467	Up	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
GEAYHHVPESQRD	0.000	0.003	0.002	0.008	3.477	1.798	0.003	2.580	Up	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
GEAYHHVPESQRDKA	0.005	0.019	0.003	0.011	3.990	1.996	0.011	1.942	Up	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
GGPEFRDDGAEPGPRE(-.98)	0.001	0.007	0.003	0.011	0.437	-1.195	0.007	2.138	Down	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
GPEFRDDGAEPGPRE(-.98)	0.004	0.017	0.003	0.011	0.514	-0.961	0.011	1.942	Down	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
LDQLLHY	0.000	0.003	0.000	0.003	2.350	1.233	0.003	2.580	Up	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
LENPQPQAPARRLLPP	0.001	0.006	0.005	0.018	0.491	-1.025	0.006	2.224	Down	tr|I3NGD6|I3NGD6_ICTTR	ProSAAS	PCSK1N
LGVLFNPYYDPLQWKSSRFE	0.000	0.003	0.003	0.011	0.156	-2.680	0.003	2.580	Down	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
MYEENSRDNPF	0.016	0.045	0.018	0.047	1.605	0.682	0.045	1.346	Up	tr|I3N6V4|I3N6V4_ICTTR	Secretogranin-3	SCG2
Q(-17.03)KIAEKFSQ	0.007	0.022	0.012	0.035	2.551	1.351	0.022	1.652	Up	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
Q(-17.03)KIAEKFSQR(-.98)	0.006	0.021	0.008	0.025	3.821	1.934	0.021	1.683	Up	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
Q(-17.03)YDRVAELDQLLHY	0.000	0.000	0.000	0.003	0.066	-3.912	0.000	3.769	Down	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
RPVKVYPNGAEDES(+79.97)AESFPLEF	0.000	0.001	0.000	0.003	0.099	-3.337	0.001	3.194	Down	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
RPVKVYPNGAEDESAESFPLEF	0.000	0.003	0.000	0.003	0.143	-2.801	0.003	2.580	Down	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
S(+42.01)YSMEHFRWGKPV(-.98)	0.007	0.022	0.012	0.035	0.486	-1.041	0.022	1.652	Down	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
SANSNPAM(+15.99)APRE	0.005	0.019	0.005	0.018	4.009	2.003	0.018	1.744	Up	tr|I3MP38|I3MP38_ICTTR	Somatostatin	SST
SVNPYLQGQRLDNVV	0.002	0.009	0.000	0.003	0.281	-1.833	0.003	2.467	Down	tr|I3MEQ9|I3MEQ9_ICTTR	Neuroendocrine protein 7B3	SCG5
SVNPYLQGQRLDNVVA	0.001	0.007	0.000	0.003	0.366	-1.448	0.003	2.467	Down	tr|I3MEQ9|I3MEQ9_ICTTR	Neuroendocrine protein 7B4	SCG5
VKMALQQEGFD	0.000	0.003	0.001	0.006	0.353	-1.501	0.003	2.580	Down	tr|I3MEG8|I3MEG8_ICTTR	Proprotein convertase subtilisin/kexin type 2	PCSK2
Y(+42.01)GGFM(+15.99)TSEKSQTPLVTLFKNAIIKNAH	0.000	0.001	0.000	0.003	0.170	-2.554	0.001	2.857	Down	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
Y(+42.01)GGFMTSEKSQTPLVTLFKNAIIKNA	0.000	0.000	0.000	0.003	0.108	-3.208	0.000	3.769	Down	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
Y(+42.01)GGFMTSEKSQTPLVTLFKNAIIKNAH	0.000	0.002	0.002	0.008	0.209	-2.255	0.002	2.654	Down	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
YDRVAELDQLLHY	0.006	0.021	0.002	0.008	0.336	-1.574	0.008	2.102	Down	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
YGGFMTSEKSQTPLVTLFKNAIIKNA	0.000	0.003	0.000	0.003	0.350	-1.513	0.003	2.580	Down	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC

Volcano plot

Code

# volcano plot
volcano_plot(uni_res_annotation, feature_col = "variable", fdr = FDR, log2fc = LOG2FC, save = FALSE)

Heatmap

Code

# heatmap
anno <- data.frame(Label = as.factor(t(d5)[, "class"])) %>%
        filter(Label %in% contrast ) %>%
        droplevels()

hm <- d5[-1, d5[1,] %in% contrast ] %>%
        rownames_to_column("variable") %>%
        mutate(across(-variable, as.numeric)) #%>%
        #left_join(., feature_meta, by = c("variable" = "Feature"))

# select color palette of current pair
selected_colors = anno_colors[levels(anno[[1]])]

heatmap(hm, feature_col = "variable", sample_anno = anno, sample_colors = selected_colors, feature_anno = feature_meta, feature_label = "Peptide", rowname_switch = TRUE, save = FALSE)

August vs March

Full stat

Code

contrast = combo_list[[2]]

df = d5_mod %>%
    dplyr::filter(., Label %in% contrast ) %>%
    droplevels()

temp_d <- df %>%
        mutate(across(-Label, ~ 2^(.) )) %>% # d3_mod is log-scale number, this step will de-log and change data back to linear number
        group_by(Label) %>%
        summarise(across(everything(), mean))

print(paste("Order of Fold Change:", temp_d$Label[1], "over", temp_d$Label[2], sep = " "))

[1] "Order of Fold Change: August over March"

Code

uni_res <- do_univariate(df)

uni_res <- uni_res %>%
        rowwise() %>%
        mutate(padj = min(c(BHT, BHW))) %>%
        # get the lowest padj
        ungroup() %>%
        mutate(`-log10padj` = -log10(padj))

uni_res_annotation = uni_res %>% left_join(., feature_meta, by = c("variable" = "Peptide"))

uni_res_annotation %>% 
        knitr::kable(., digits = 3, "html") %>%
        kableExtra::kable_styling(bootstrap_options = c("striped", "hover", "condensed", "responsive")) %>%
        kableExtra::scroll_box(width = "100%", height = "500px")

variable	pT	BHT	pW	BHW	FC(lin)	FC(log2)	padj	-log10padj	Accession	Name	GN
AADQDLGPEAPPEGVLGA	0.078	0.127	0.036	0.070	0.674	-0.570	0.070	1.152	tr|I3NGD6|I3NGD6_ICTTR	ProSAAS	PCSK1N
AADQDLGPEAPPEGVLGAL	0.021	0.048	0.036	0.070	0.482	-1.053	0.048	1.321	tr|I3NGD6|I3NGD6_ICTTR	ProSAAS	PCSK1N
AAMDLEL	0.727	0.781	0.743	0.785	0.624	-0.680	0.781	0.107	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
AELDQLLHY	0.229	0.296	0.370	0.446	1.533	0.617	0.296	0.529	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
AFGPRGPGLQL	0.136	0.196	0.093	0.141	1.409	0.495	0.141	0.852	tr|I3MMW2|I3MMW2_ICTTR	Chromogranin-A	CHGA
AGAPEPAEPMELAKPSA	0.001	0.006	0.002	0.006	0.248	-2.011	0.006	2.213	tr|I3MXH3|I3MXH3_ICTTR	Vasopressin-neurophysin 2-copeptin-like	AVP
ATLESVFQELGKLTGPSSQKRE	0.174	0.241	0.370	0.446	1.245	0.316	0.241	0.619	tr|I3N6V4|I3N6V4_ICTTR	Secretogranin-2	SCG2
DDGPYRMEHFRWGSPPKD	0.317	0.382	0.277	0.356	0.752	-0.412	0.356	0.448	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
EDREEASLQDRQYASHHT	0.007	0.021	0.021	0.047	3.147	1.654	0.021	1.677	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
EDREEASLQDRQYASHHTTE	0.001	0.004	0.001	0.003	2.874	1.523	0.003	2.545	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
EEASLQDRQYASHHT	0.000	0.002	0.001	0.003	3.938	1.978	0.002	2.680	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
EEASLQDRQYASHHTTE	0.000	0.003	0.001	0.003	2.415	1.272	0.003	2.545	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
EESEESEEGTTSEVT	0.278	0.348	0.423	0.498	1.832	0.873	0.348	0.459	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
EIGDEENSAKFPI(-.98)	0.522	0.577	0.673	0.727	0.940	-0.089	0.577	0.239	tr|I3MQ49|I3MQ49_ICTTR	Pro-melanin concentrating hormone	PMCH
ELDQLLHY	0.086	0.138	0.046	0.079	1.957	0.969	0.079	1.101	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
FKNAIIKNAH	0.364	0.417	0.606	0.670	1.116	0.158	0.417	0.380	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
FKNAIIKNAHKKGQ	0.133	0.195	0.046	0.079	0.503	-0.991	0.079	1.101	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
FLGEAYHHVPE	0.011	0.027	0.011	0.029	2.349	1.232	0.027	1.563	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
FLGEAYHHVPES	0.071	0.119	0.236	0.308	1.548	0.631	0.119	0.923	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
FLGEAYHHVPESQ	0.006	0.018	0.001	0.004	2.518	1.332	0.004	2.355	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
FLGEAYHHVPESQRD	0.000	0.000	0.000	0.001	4.087	2.031	0.000	4.650	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
FLGEAYHHVPESQRDKA	0.000	0.001	0.000	0.001	4.343	2.119	0.001	3.141	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
FPTLGGSQDKSLHN	0.000	0.000	0.000	0.001	2.250	1.170	0.000	5.140	tr|I3MPL6|I3MPL6_ICTTR	Secretogranin-3	SCG3
GEAYHHVPESQRD	0.001	0.004	0.001	0.004	3.056	1.611	0.004	2.423	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
GEAYHHVPESQRDKA	0.001	0.004	0.000	0.001	4.852	2.279	0.001	3.066	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
GGFMRGL	0.179	0.245	0.236	0.308	1.491	0.576	0.245	0.612	tr|I3MGR7|I3MGR7_ICTTR	Proenkephalin-A	PENK
GGPEFRDDGAEPGPRE(-.98)	0.101	0.154	0.093	0.141	1.558	0.640	0.141	0.852	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
GPEFRDDGAEPGPRE(-.98)	0.002	0.006	0.002	0.009	1.895	0.923	0.006	2.209	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
GRTLYGFGG	0.036	0.072	0.036	0.070	1.904	0.929	0.070	1.152	tr|I3N7M0|I3N7M0_ICTTR	Histone H4	LOC106145262
GVANALSHKYH	0.055	0.102	0.046	0.079	3.878	1.955	0.079	1.101	tr|A0A287CUK4|A0A287CUK4_ICTTR	Hemoglobin subunit beta-S/F-like	LOC101967841
GVLFNPYYDPLQWKSSRFE	0.873	0.922	0.815	0.851	0.839	-0.253	0.851	0.070	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
KVVAGVANALSHKYH	0.006	0.018	0.004	0.012	1.886	0.916	0.012	1.936	tr|A0A287CUK4|A0A287CUK4_ICTTR	Hemoglobin subunit beta-S/F-like	LOC101967841
KYVMGHFRWD	0.061	0.108	0.046	0.079	0.580	-0.787	0.079	1.101	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
KYVMGHFRWDRF(-.98)	0.489	0.553	0.541	0.606	0.794	-0.333	0.553	0.257	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
LAGAPEPA	0.000	0.001	0.000	0.001	0.043	-4.541	0.001	3.141	tr|I3MXH3|I3MXH3_ICTTR	Vasopressin-neurophysin 2-copeptin-like	AVP
LDQLLHY	0.000	0.000	0.000	0.001	4.580	2.195	0.000	3.902	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
LENPQPQAPARRLLPP	0.499	0.559	0.423	0.498	0.898	-0.156	0.498	0.303	tr|I3NGD6|I3NGD6_ICTTR	ProSAAS	PCSK1N
LGEAYHHVPESQRD	0.000	0.002	0.000	0.001	5.558	2.475	0.001	3.066	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
LGVLFNPYYDPLQWKSSRFE	0.070	0.119	0.059	0.098	0.536	-0.901	0.098	1.010	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
MYEENSRDNPF	0.010	0.026	0.021	0.047	1.880	0.911	0.026	1.581	tr|I3N6V4|I3N6V4_ICTTR	Secretogranin-3	SCG2
PEFRDDGAEPGPRE(-.98)	0.002	0.009	0.002	0.009	2.148	1.103	0.009	2.066	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
Q(-17.03)EPGAPAAGM	0.294	0.363	0.167	0.238	0.804	-0.316	0.238	0.623	tr|I3MIV6|I3MIV6_ICTTR	Carboxypeptidase E	CPE
Q(-17.03)HWSYGLRPG(-.98)	0.001	0.005	0.006	0.015	2.119	1.084	0.005	2.268	tr|I3MCA0|I3MCA0_ICTTR	Progonadoliberin-1	GNRH1
Q(-17.03)KIAEKFSQ	0.003	0.010	0.001	0.003	3.025	1.597	0.003	2.545	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
Q(-17.03)KIAEKFSQR(-.98)	0.001	0.004	0.000	0.001	5.746	2.523	0.001	2.890	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
Q(-17.03)LEDEAKELQ	0.095	0.147	0.167	0.238	1.922	0.942	0.147	0.833	tr|I3MGR7|I3MGR7_ICTTR	Proenkephalin-A	PENK
Q(-17.03)YDRVAELDQLLHY	0.045	0.087	0.046	0.079	0.418	-1.258	0.079	1.101	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
QYDRVAELDQLLHY	0.214	0.283	0.167	0.238	0.688	-0.540	0.238	0.623	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
RPVKVYPNGAEDES(+79.97)AESFPLEF	0.006	0.019	0.000	0.001	0.091	-3.466	0.001	3.066	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
RPVKVYPNGAEDESAESFPLEF	0.001	0.006	0.002	0.009	0.232	-2.109	0.006	2.213	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
RTLYGFGG	0.003	0.010	0.004	0.012	2.269	1.182	0.010	1.987	tr|I3N7M0|I3N7M0_ICTTR	Histone H4	LOC106145262
RVLTEEEEKELENL	0.936	0.966	0.888	0.898	0.933	-0.100	0.898	0.047	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
S(+42.01)YSMEHFRWGKPV(-.98)	0.023	0.052	0.021	0.047	0.633	-0.661	0.047	1.325	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
SANSNPAM(+15.99)APRE	0.095	0.147	0.074	0.119	3.343	1.741	0.119	0.926	tr|I3MP38|I3MP38_ICTTR	Somatostatin	SST
SANSNPAMAPRE	0.000	0.002	0.001	0.003	3.701	1.888	0.002	2.796	tr|I3MP38|I3MP38_ICTTR	Somatostatin	SST
SLDSPAGPA	0.968	0.978	0.743	0.785	1.050	0.071	0.785	0.105	tr|I3M5J1|I3M5J1_ICTTR	Corticoliberin	CRH
SLDSPAGPAE	0.003	0.010	0.004	0.012	3.484	1.801	0.010	1.987	tr|I3M5J1|I3M5J1_ICTTR	Corticoliberin	CRH
SPQLEDE	0.073	0.120	0.036	0.070	3.874	1.954	0.070	1.152	tr|I3MGR7|I3MGR7_ICTTR	Proenkephalin-A	PENK
SPQLEDEAKEL	0.007	0.019	0.006	0.015	1.794	0.843	0.015	1.817	tr|I3MGR7|I3MGR7_ICTTR	Proenkephalin-A	PENK
SPQLEDEAKELQ	0.010	0.026	0.011	0.029	1.542	0.625	0.026	1.581	tr|I3MGR7|I3MGR7_ICTTR	Proenkephalin-A	PENK
SQGGEPGAYLTPDT	0.056	0.102	0.074	0.119	1.755	0.812	0.102	0.990	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
SQGGEPGAYLTPDTREE	0.010	0.026	0.002	0.009	1.626	0.701	0.009	2.066	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
SQTPLVTLF	0.335	0.394	0.481	0.544	0.728	-0.458	0.394	0.404	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
SQTPLVTLFKNAIIKNAHKKGQ	0.155	0.221	0.200	0.272	3.215	1.685	0.221	0.656	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
SSQEGNPPSEERGHTLQEPEEAKVGE	0.113	0.168	0.236	0.308	1.409	0.494	0.168	0.774	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
SVNPYLQGQRL	0.000	0.001	0.001	0.003	6.020	2.590	0.001	3.031	tr|I3MEQ9|I3MEQ9_ICTTR	Neuroendocrine protein 7B2	SCG5
SVNPYLQGQRLDNVV	0.032	0.068	0.021	0.047	0.593	-0.755	0.047	1.325	tr|I3MEQ9|I3MEQ9_ICTTR	Neuroendocrine protein 7B3	SCG5
SVNPYLQGQRLDNVVA	0.893	0.932	0.888	0.898	1.062	0.087	0.898	0.047	tr|I3MEQ9|I3MEQ9_ICTTR	Neuroendocrine protein 7B4	SCG5
SVSTVLTSKYR	0.067	0.117	0.093	0.141	1.829	0.871	0.117	0.932	tr|A0A287D1R9|A0A287D1R9_ICTTR	Hemoglobin subunit alpha	LOC101956162
SYSMEHF	0.035	0.071	0.027	0.060	1.578	0.658	0.060	1.223	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
SYSMEHFRWGKPV(-.98)	0.188	0.252	0.167	0.238	1.477	0.562	0.238	0.623	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
TLGGSQDKSLHN	0.001	0.004	0.000	0.001	4.597	2.201	0.001	3.066	tr|I3MPL6|I3MPL6_ICTTR	Secretogranin-4	SCG3
TLYGFGG	0.984	0.984	1.000	1.000	1.062	0.087	0.984	0.007	tr|I3N7M0|I3N7M0_ICTTR	Histone H4	LOC106145262
TNEIVEEQYTPQSL	0.057	0.104	0.059	0.098	1.868	0.901	0.098	1.010	tr|I3N6V4|I3N6V4_ICTTR	Secretogranin-4	SCG2
TPLVTLFKNAIIKNAHKKGQ	0.171	0.240	0.200	0.272	2.089	1.063	0.240	0.620	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
VAELDQLLHY	0.731	0.781	0.888	0.898	1.103	0.142	0.781	0.107	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
VANALSHKYH	0.017	0.040	0.036	0.070	2.994	1.582	0.040	1.398	tr|A0A287CUK4|A0A287CUK4_ICTTR	Hemoglobin subunit beta-S/F-like	LOC101967841
VGRPEWWMDYQ	0.332	0.394	0.321	0.403	0.732	-0.450	0.394	0.404	tr|I3MGR7|I3MGR7_ICTTR	Proenkephalin-A	PENK
VKMALQQEGFD	0.001	0.004	0.002	0.006	0.397	-1.335	0.004	2.423	tr|I3MEG8|I3MEG8_ICTTR	Proprotein convertase subtilisin/kexin type 2	PCSK2
VKVGVNGF	0.027	0.059	0.046	0.079	1.819	0.863	0.059	1.227	tr|A0A287D887|A0A287D887_ICTTR	Glyceraldehyde-3-phosphate dehydrogenase	GAPDH
Y(+42.01)GGFM(+15.99)TSEKSQTPLVTLFKNAIIKNAH	0.000	0.001	0.000	0.001	0.224	-2.159	0.001	3.141	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
Y(+42.01)GGFMTSEKSQTPLVTLFKNAIIKNA	0.000	0.000	0.000	0.001	0.071	-3.819	0.000	4.971	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
Y(+42.01)GGFMTSEKSQTPLVTLFKNAIIKNAH	0.000	0.001	0.000	0.002	0.255	-1.974	0.001	2.848	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
Y(+42.01)GGFMTSEKSQTPLVTLFKNAIIKNAHKKGQ	0.304	0.371	0.200	0.272	0.561	-0.833	0.272	0.566	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
YDRVAELDQLLHY	0.230	0.296	0.321	0.403	0.783	-0.353	0.296	0.529	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
YGGFM(+15.99)RF	0.277	0.348	0.481	0.544	1.136	0.184	0.348	0.459	tr|I3MGR7|I3MGR7_ICTTR	Proenkephalin-A	PENK
YGGFM(+15.99)RGL	0.681	0.744	0.370	0.446	0.815	-0.295	0.446	0.350	tr|I3MGR7|I3MGR7_ICTTR	Proenkephalin-A	PENK
YGGFM(+15.99)TSEKSQTPLVTLFKNAIIKNAHKKGQ	0.946	0.966	0.673	0.727	1.365	0.449	0.727	0.138	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
YGGFMRF	0.038	0.074	0.021	0.047	1.673	0.743	0.047	1.325	tr|I3MGR7|I3MGR7_ICTTR	Proenkephalin-A	PENK
YGGFMRGL	0.011	0.026	0.027	0.060	1.972	0.979	0.026	1.581	tr|I3MGR7|I3MGR7_ICTTR	Proenkephalin-A	PENK
YGGFMTSE	0.348	0.403	0.481	0.544	1.249	0.321	0.403	0.394	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
YGGFMTSEKSQTPLVTLFKNAIIKNA	0.004	0.013	0.006	0.015	0.490	-1.030	0.013	1.885	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
YGGFMTSEKSQTPLVTLFKNAIIKNAH	0.008	0.023	0.006	0.015	0.386	-1.374	0.015	1.817	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
YGGFMTSEKSQTPLVTLFKNAIIKNAHKKGQ	0.034	0.071	0.046	0.079	1.735	0.795	0.071	1.152	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC

Significant stat

Code

# get DE features only tibble
uni_res_filt <- uni_res %>%
        dplyr::filter(BHT < FDR & BHW < FDR) %>%
        mutate(status = case_when(`FC(log2)` < - 0.6 ~ "Down",
                                  `FC(log2)` > 0.6 ~ "Up") )

if (nrow(uni_res_filt) == 0){
        print(paste("There is no significant differentiated features between",
                    temp_d$Label[2], "and", temp_d$Label[1], sep = " "))
} else {
        uni_res_filt = uni_res_filt %>%
                left_join(., feature_meta, by = c("variable" = "Peptide"))
        
        uni_res_filt %>% 
                knitr::kable(., digits = 3, "html") %>%
                kableExtra::kable_styling(bootstrap_options = c("striped", "hover", "condensed", "responsive")) %>%
                kableExtra::scroll_box(width = "100%", height = "500px")
}

variable	pT	BHT	pW	BHW	FC(lin)	FC(log2)	padj	-log10padj	status	Accession	Name	GN
AGAPEPAEPMELAKPSA	0.001	0.006	0.002	0.006	0.248	-2.011	0.006	2.213	Down	tr|I3MXH3|I3MXH3_ICTTR	Vasopressin-neurophysin 2-copeptin-like	AVP
EDREEASLQDRQYASHHT	0.007	0.021	0.021	0.047	3.147	1.654	0.021	1.677	Up	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
EDREEASLQDRQYASHHTTE	0.001	0.004	0.001	0.003	2.874	1.523	0.003	2.545	Up	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
EEASLQDRQYASHHT	0.000	0.002	0.001	0.003	3.938	1.978	0.002	2.680	Up	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
EEASLQDRQYASHHTTE	0.000	0.003	0.001	0.003	2.415	1.272	0.003	2.545	Up	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
FLGEAYHHVPE	0.011	0.027	0.011	0.029	2.349	1.232	0.027	1.563	Up	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
FLGEAYHHVPESQ	0.006	0.018	0.001	0.004	2.518	1.332	0.004	2.355	Up	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
FLGEAYHHVPESQRD	0.000	0.000	0.000	0.001	4.087	2.031	0.000	4.650	Up	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
FLGEAYHHVPESQRDKA	0.000	0.001	0.000	0.001	4.343	2.119	0.001	3.141	Up	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
FPTLGGSQDKSLHN	0.000	0.000	0.000	0.001	2.250	1.170	0.000	5.140	Up	tr|I3MPL6|I3MPL6_ICTTR	Secretogranin-3	SCG3
GEAYHHVPESQRD	0.001	0.004	0.001	0.004	3.056	1.611	0.004	2.423	Up	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
GEAYHHVPESQRDKA	0.001	0.004	0.000	0.001	4.852	2.279	0.001	3.066	Up	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
GPEFRDDGAEPGPRE(-.98)	0.002	0.006	0.002	0.009	1.895	0.923	0.006	2.209	Up	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
KVVAGVANALSHKYH	0.006	0.018	0.004	0.012	1.886	0.916	0.012	1.936	Up	tr|A0A287CUK4|A0A287CUK4_ICTTR	Hemoglobin subunit beta-S/F-like	LOC101967841
LAGAPEPA	0.000	0.001	0.000	0.001	0.043	-4.541	0.001	3.141	Down	tr|I3MXH3|I3MXH3_ICTTR	Vasopressin-neurophysin 2-copeptin-like	AVP
LDQLLHY	0.000	0.000	0.000	0.001	4.580	2.195	0.000	3.902	Up	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
LGEAYHHVPESQRD	0.000	0.002	0.000	0.001	5.558	2.475	0.001	3.066	Up	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
MYEENSRDNPF	0.010	0.026	0.021	0.047	1.880	0.911	0.026	1.581	Up	tr|I3N6V4|I3N6V4_ICTTR	Secretogranin-3	SCG2
PEFRDDGAEPGPRE(-.98)	0.002	0.009	0.002	0.009	2.148	1.103	0.009	2.066	Up	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
Q(-17.03)HWSYGLRPG(-.98)	0.001	0.005	0.006	0.015	2.119	1.084	0.005	2.268	Up	tr|I3MCA0|I3MCA0_ICTTR	Progonadoliberin-1	GNRH1
Q(-17.03)KIAEKFSQ	0.003	0.010	0.001	0.003	3.025	1.597	0.003	2.545	Up	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
Q(-17.03)KIAEKFSQR(-.98)	0.001	0.004	0.000	0.001	5.746	2.523	0.001	2.890	Up	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
RPVKVYPNGAEDES(+79.97)AESFPLEF	0.006	0.019	0.000	0.001	0.091	-3.466	0.001	3.066	Down	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
RPVKVYPNGAEDESAESFPLEF	0.001	0.006	0.002	0.009	0.232	-2.109	0.006	2.213	Down	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
RTLYGFGG	0.003	0.010	0.004	0.012	2.269	1.182	0.010	1.987	Up	tr|I3N7M0|I3N7M0_ICTTR	Histone H4	LOC106145262
SANSNPAMAPRE	0.000	0.002	0.001	0.003	3.701	1.888	0.002	2.796	Up	tr|I3MP38|I3MP38_ICTTR	Somatostatin	SST
SLDSPAGPAE	0.003	0.010	0.004	0.012	3.484	1.801	0.010	1.987	Up	tr|I3M5J1|I3M5J1_ICTTR	Corticoliberin	CRH
SPQLEDEAKEL	0.007	0.019	0.006	0.015	1.794	0.843	0.015	1.817	Up	tr|I3MGR7|I3MGR7_ICTTR	Proenkephalin-A	PENK
SPQLEDEAKELQ	0.010	0.026	0.011	0.029	1.542	0.625	0.026	1.581	Up	tr|I3MGR7|I3MGR7_ICTTR	Proenkephalin-A	PENK
SQGGEPGAYLTPDTREE	0.010	0.026	0.002	0.009	1.626	0.701	0.009	2.066	Up	tr|I3NBA1|I3NBA1_ICTTR	Secretogranin-1	CHGB
SVNPYLQGQRL	0.000	0.001	0.001	0.003	6.020	2.590	0.001	3.031	Up	tr|I3MEQ9|I3MEQ9_ICTTR	Neuroendocrine protein 7B2	SCG5
TLGGSQDKSLHN	0.001	0.004	0.000	0.001	4.597	2.201	0.001	3.066	Up	tr|I3MPL6|I3MPL6_ICTTR	Secretogranin-4	SCG3
VKMALQQEGFD	0.001	0.004	0.002	0.006	0.397	-1.335	0.004	2.423	Down	tr|I3MEG8|I3MEG8_ICTTR	Proprotein convertase subtilisin/kexin type 2	PCSK2
Y(+42.01)GGFM(+15.99)TSEKSQTPLVTLFKNAIIKNAH	0.000	0.001	0.000	0.001	0.224	-2.159	0.001	3.141	Down	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
Y(+42.01)GGFMTSEKSQTPLVTLFKNAIIKNA	0.000	0.000	0.000	0.001	0.071	-3.819	0.000	4.971	Down	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
Y(+42.01)GGFMTSEKSQTPLVTLFKNAIIKNAH	0.000	0.001	0.000	0.002	0.255	-1.974	0.001	2.848	Down	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
YGGFMTSEKSQTPLVTLFKNAIIKNA	0.004	0.013	0.006	0.015	0.490	-1.030	0.013	1.885	Down	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC
YGGFMTSEKSQTPLVTLFKNAIIKNAH	0.008	0.023	0.006	0.015	0.386	-1.374	0.015	1.817	Down	tr|I3M7B6|I3M7B6_ICTTR	Pro-opiomelanocortin	POMC

Volcano plot

Code

# volcano plot
volcano_plot(uni_res_annotation, feature_col = "variable", fdr = FDR, log2fc = LOG2FC, save = FALSE)

Heatmap

Code

# heatmap
anno <- data.frame(Label = as.factor(t(d5)[, "class"])) %>%
        filter(Label %in% contrast ) %>%
        droplevels()

hm <- d5[-1, d5[1,] %in% contrast ] %>%
        rownames_to_column("variable") %>%
        mutate(across(-variable, as.numeric)) #%>%
        #left_join(., feature_meta, by = c("variable" = "Feature"))

# select color palette of current pair
selected_colors = anno_colors[levels(anno[[1]])]

heatmap(hm, feature_col = "variable", sample_anno = anno, sample_colors = selected_colors, feature_anno = feature_meta, feature_label = "Peptide", rowname_switch = TRUE, save = FALSE)