CN109273051B - Identity information coding method based on telomere length - Google Patents

Identity information coding method based on telomere length Download PDF

Info

Publication number
CN109273051B
CN109273051B CN201811001443.7A CN201811001443A CN109273051B CN 109273051 B CN109273051 B CN 109273051B CN 201811001443 A CN201811001443 A CN 201811001443A CN 109273051 B CN109273051 B CN 109273051B
Authority
CN
China
Prior art keywords
code
telomere
subject
age
telomere length
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201811001443.7A
Other languages
Chinese (zh)
Other versions
CN109273051A (en
Inventor
夏茂
徐鑫
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Individual
Original Assignee
Individual
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Individual filed Critical Individual
Priority to CN201811001443.7A priority Critical patent/CN109273051B/en
Publication of CN109273051A publication Critical patent/CN109273051A/en
Application granted granted Critical
Publication of CN109273051B publication Critical patent/CN109273051B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30024Cell structures in vitro; Tissue sections in vitro

Landscapes

  • Engineering & Computer Science (AREA)
  • Quality & Reliability (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Health & Medical Sciences (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Image Processing (AREA)
  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)

Abstract

The invention belongs to the field of information coding, and particularly relates to an identity information coding method based on telomere length. The method at least comprises a static code and a dynamic code, wherein: the static code comprises a country code, a language code and a face image acquisition sequence code of a subject; the dynamic codes include a telomere information code, a credit code, and a health code of the subject. According to the method, besides the national, language, face identification and other information necessary for the traditional identity information coding method, the cell age and telomere length annual average change degree information of the subject is embedded, so that automatic information comparison can be rapidly carried out on the people in a specific sample area and even specific information people through the coding, and finally the overall health level of the society is favorably improved.

Description

Identity information coding method based on telomere length
Technical Field
The invention belongs to the field of information coding, and particularly relates to an identity information coding method based on telomere length.
Background
With the development of social science and technology, people pay more and more attention to self health, and meanwhile, the health information is gradually combined with modern science and technology, so that the self health information is timely collected, evaluated and exchanged, the defects in the health information are pertinently improved in the later period through the datamation of the information, and finally the health level of the society is integrally improved. Telomeres are located at the ends of chromosomes and serve to protect the ends of chromosomes from fusion and degeneration and to prevent DNA shedding during cell division. After each cell division, telomeres shorten because they are lost in part during each replication cycle; when telomeres can not be shortened any more, the cells reach an aging state; at this point, the cell irreversibly stops dividing and immediately activates the apoptotic mechanism, i.e., the cell goes towards apoptosis. By measuring the length of telomeres, our physiological age is clearly known. Whether telomere information can be coded so as to be combined with common codes or not, and then characteristic information such as the physiological age of a current information person can be obtained by the codes during information exchange so as to realize the aim of informationalized comparison of a group in a specific sample area and even a specific information person by utilizing the efficient transmission and exchange of information streams and finally realize the targeted health defect improvement effect, which is a technical problem to be solved urgently in recent years in the field.
Disclosure of Invention
The invention aims to overcome the defects of the prior art and provide an identity information coding method based on the telomere length, which is characterized in that the cell age and the average change degree of the telomere length per year of a subject are embedded in addition to the national, language, face identification and other information which are necessary in the traditional identity information coding method, so that automatic information comparison can be rapidly carried out on a group in a specific sample area and even a specific information person through the coding, and the whole health level of the society is improved.
In order to achieve the purpose, the invention adopts the following technical scheme:
an identity information coding method based on telomere length is characterized in that: the method at least comprises a static code and a dynamic code, wherein:
the static code comprises a country code, a language code and a face image acquisition sequence code of a subject;
the dynamic codes comprise a telomere information code, a credit code and a health code of the subject; the telomere information code is two bits, and the acquisition process comprises the following steps:
1) calculating the length of the telomere of the sample of the subject;
2) establishing a sample database according to the existing telomere lengths of the subjects with different ages and different sexes;
dividing the sample database into two groups according to gender, wherein the first group is male, and the second group is female; and the age is recorded as x0Wherein
Figure GDA0003241410480000021
Figure GDA0003241410480000022
The maximum age at which humans can theoretically survive; the hour notes
Figure GDA0003241410480000023
For male group data, applying statistical software to perform regression analysis on the telomere length of the samples of the subjects and the corresponding age of the subjects to obtain the following formula:
TRF1=a1x0+b1
for the female group data, applying statistical software to perform regression analysis on the telomere length of the subject samples and the corresponding age of the subject to obtain the following formula:
TRF0=a0x0+b0
3) and respectively carrying out telomere length detection once when the subject is at the age of t and t + k, wherein the telomere lengths of corresponding samples are TRF(t)、TRF(t+k)And sequentially executing the following substeps:
A. calculating cell age of the sample
Cellular age CA at t + k years of age of subject(t+k)Comprises the following steps:
Figure GDA0003241410480000024
B. calculating the average change degree of telomere length of the sample
Mean degree of change of telomere length in the subject from t to t + k years of aget,kComprises the following steps:
Figure GDA0003241410480000031
C. method for establishing telomere information code
Recording the Gaussian function as the maximum integer of Gauss (x) not more than x;
Figure GDA0003241410480000032
then m is5Is an integer of not less than 0 and not more than 35;
recording min (x, y) as the smaller number of x and y;
Figure GDA0003241410480000033
then m is4Is an integer of not less than 0 and not more than 35;
calculate m5And m4Then m is added5And m4Is converted into a thirty-six system symbol according to a thirty-six system, and is converted into a thirty-six system symbol according to m5 m4The two-bit character string is obtained when the two-bit character string is spliced together, and the two-bit character string is the telomere information code.
Preferably, in the step 1), calculating the length of the telomere of the sample of the subject comprises the following sub-steps:
a. 3-5 ml of peripheral blood is extracted and put into 0.5mol/L EDTAK2Freezing and storing in an anticoagulant centrifuge tube at-20 ℃;
b. extracting peripheral blood leukocyte genome DNA by adopting a saturated phenol-chloroform method, selecting a sample with an OD260/OD280 value of 1.8-2.0, and observing the integrity of the DNA under ultraviolet rays after 0.8% agarose gel electrophoresis;
c. performing mixed enzyme digestion on genomic DNA for 2h by using a Telo TAGGG Telomere Length Assay kit of Roche, performing electrophoresis for 4h in 0.8% agarose gel at 5V/cm for 5 h, denaturing and neutralizing the gel, transferring the DNA to a positively charged nylon membrane by a siphon blotting method, and baking for 20min at 120 ℃; pre-hybridizing at 42 ℃ for 1h, and adding a digoxin-labeled telomere probe for hybridization for 3 h; and (3) carrying out chemiluminescence development after washing the membrane, carrying out optical density scanning imaging on the obtained X-ray film, analyzing by using telomere length analysis software Telomeric1.2, and calculating the telomere length of each sample according to the following formula:
Figure GDA0003241410480000041
wherein L isiIs the i-point marker length;
ODithe optical density value of point i.
Preferably, the country code is two digits and refers to ISO 3166-1; the language code is two digits, and refers to ISO 693-1.
Preferably, the face image acquisition sequence code acquisition process comprises the following sub-steps:
(1) acquiring time of a subject when the face image is acquired, wherein the time comprises four parts of year, month, date, clock and minute; the ones of each part are complemented into two digits on the left side by 0;
(2) subtracting 2000 from year to obtain four digits, and if the four digits are not enough, the four digits are compensated by 0 on the left side; arranging the quartic, month, date, clock and minute in the original sequence;
(3) supplementing a two-digit sequence number so as to form 14-digit numbers together with the numbers sequentially arranged in the substep (2); the sequence number is determined according to the sequence of the images collected by the testee in the same minute and the sex of the images; specifically, if the subject is male, the last digit of the order number is odd, and the left side of the two-digit deficiency is supplemented with 0 to form two digits; if the subject is female, the last digit is an even number, and the left side of the two digits is supplemented with 0 to form two digits;
(4) and carrying out thirty-six system conversion on the 14-bit numbers in the substep (3) to form nine-bit face image acquisition sequence codes, wherein the left side of less than nine bits is complemented with 0 to form nine bits.
Preferably, the credit code is one bit, and the health code is one bit, and the forming steps are as follows:
when the subject credits are Cre and the credit full is Cremax(ii) a Subjects were classified as Hea healthy and Hea full healthymaxWhen the current is over;
calculating m by calculating Gauss function as the maximum integer not exceeding x3
Figure GDA0003241410480000042
Then m is put3Converting into symbols in thirty-six system according to thirty-six system to obtain corresponding credit codes;
and m is calculated according to the following formula2
Figure GDA0003241410480000051
Then m is put2And converting into symbols in thirty-six system according to the thirty-six system to obtain the corresponding health codes.
Preferably, the method further comprises a check code for detecting whether the static code and the dynamic code are mistakenly transmitted or not after the static code and the dynamic code are sequentially arranged from left to right; the check code is obtained as follows:
when the dynamic code and the static code are all known, setting the code as the dynamic code
M18M17M16M15M14M13M12M11M10M9M8M7M6M5M4M3M2
Wherein M isjRepresents a character, j is more than or equal to 2 and less than or equal to 18;
step1 character according to thirty-six systemConversion table will MjReverse conversion to a number mj
Step2 for each mjCalculating the weight ω according to the following formulaj
ωj=mod(2j-1,37)
Wherein mod (x,37) represents the remainder of the integer x divided by 37;
step3 calculating m according to the following formula1
Figure GDA0003241410480000052
Step4 converting the number m according to the thirty-six system character conversion table1Conversion to character M1The M being1Namely the check code.
The invention has the beneficial effects that:
1) by the method, the disordered original information of the subjects, namely the specific information persons or the subjects, namely the sample area crowd can be combed in an informationized manner, and the comparison function of the health data of the specific information persons and the people in the sample area crowd in a specified age range is realized. Even through the establishment of the database, the average value of the average change degree of the cell age and the telomere length year in the current society can be obtained, and when a new sample enters the database, the information can be coded, and then the two-bit characters of the telomere information code are compared with the thirty-six-digit characters converted from the average value of the average change degree of the cell age and the telomere length year, so that the aging degree and the vitality degree of the person who belongs to the current sample, namely the subject, can be judged, and whether to perform post-recuperation and exercise in a targeted manner can be determined. The coding method is purely information character coding, so that the coding method can be even applied to a block chain, thereby becoming the basis of a block chain BaaS platform, and finally realizing the functions of convenience and high efficiency of information communication of different people and even different crowds by depending on the BaaS platform.
Drawings
FIG. 1 is a layout comparison table of the present invention;
FIG. 2 is a table of thirty-six digits and characters;
fig. 3 is a 14-digit sequence table for calculating a face image capture sequence code.
Detailed Description
For ease of understanding, the specific structure and operation of the present invention is further described herein with reference to FIGS. 1-3:
an identity information encoding method based on telomere length, as shown in fig. 1, includes a static code, a dynamic code and a check code, which are arranged in sequence from left to right, wherein:
the static code comprises a country code, a language code and a face image acquisition sequence code of the subject; wherein:
the country code is two digits, and refers to ISO 3166-1; the language code is two digits, and refers to ISO 693-1. In the specific operation, the face image acquisition sequence code can be formed by adding the last two-bit allocation code to the time (accurate to minutes) when the face image is acquired by the subject. The acquisition process of the face image acquisition sequence code comprises the following steps:
1) firstly, acquiring the time of a subject for acquiring a face image, wherein the time is usually Beijing time; including year, month, date, clock and minute, the single digit must be complemented to the left by 0 to two digits.
2) Subtracting 2000 from year to obtain four digits, supplementing 0 to four digits on left side, and arranging the four digits with month, date, clock and minute in original sequence.
3) And a two digit sequence number is added to form a 14 digit number, as shown in figure 3. The order number depends on the sequence of images acquired by the subject in the same minute and the sex, so that the number of the order number can be increased as appropriate in subsequent use. Taking two-bit sequence numbers as an example: if the testee is a male, the last digit of the sequence number is an odd number, and the left side with less than two digits is complemented into two digits by 0, namely, the male with the 1 st collected image at a certain time point 01 is a male with the 2 nd collected image 03, and the like. If the testee is a woman, the last digit of the sequence number is an even number, and the left side with less than two digits is complemented into two digits with 0, namely, the woman whose image is acquired at the 1 st digit at a certain time point 00, the woman whose image is acquired at the 2 nd digit 02, and so on. For example, a subject was a 22 nd female subject who was 45 points 23/7/1/2018, and the corresponding 14 digits were 00180701234542. This is similar to the ID card number, with five to twelve digits from the last being the birth date and one to four digits from the last being the sex number.
4) And carrying out trihexa-hexadecimal conversion on the 14-bit digital in the step 3) to form nine-bit face image acquisition sequence codes, wherein the left side of less than nine bits is complemented into nine bits by 0.
The dynamic codes include a telomere information code, a credit code, and a health code of the subject.
The credit code is one bit, the credit score is prior art or available using prior third party technology, and the corresponding credit score is Cre (Cre)maxDivision, i.e. credit full into Cremax);
Calculating m by calculating Gauss function as the maximum integer not exceeding x3
Figure GDA0003241410480000071
Then m is put3And converting into symbols in thirty-six system according to the thirty-six system to obtain the corresponding credit codes.
The health code is one bit, the health score is also prior art, or obtained by applying the prior third party technology, and the corresponding health score is Hea (Hea)maxDivision, i.e. full health as Heamax);
M is calculated according to the following formula2
Figure GDA0003241410480000072
Then m is put2And converting into symbols in thirty-six system according to the thirty-six system to obtain the corresponding health codes.
The telomere information code is two bits, and the acquisition process comprises the following steps:
1) calculating the length of a sample telomere of the subject, comprising the following sub-steps:
a. 3-5 ml of peripheral blood is extracted and put into 0.5mol/L EDTAK2Freezing and storing in an anticoagulant centrifuge tube at-20 ℃;
b. extracting peripheral blood leukocyte genome DNA by adopting a saturated phenol-chloroform method, selecting a sample with an OD260/OD280 value of 1.8-2.0, and observing the integrity of the DNA under ultraviolet rays after 0.8% agarose gel electrophoresis;
c. performing mixed enzyme digestion on genomic DNA for 2h by using a Telo TAGGG Telomere Length Assay kit of Roche, performing electrophoresis for 4h in 0.8% agarose gel at 5V/cm for 5 h, denaturing and neutralizing the gel, transferring the DNA to a positively charged nylon membrane by a siphon blotting method, and baking for 20min at 120 ℃; pre-hybridizing at 42 ℃ for 1h, and adding a digoxin-labeled telomere probe for hybridization for 3 h; and (3) carrying out chemiluminescence development after washing the membrane, carrying out optical density scanning imaging on the obtained X-ray film, analyzing by using telomere length analysis software Telomeric1.2, and calculating the telomere length of each sample according to the following formula:
Figure GDA0003241410480000081
wherein L isiIs the i-point marker length;
ODithe optical density value of point i.
2) Establishing a sample database according to the existing telomere lengths of the subjects with different ages and different sexes;
dividing the sample database into two groups according to gender, wherein the first group is male, and the second group is female; and the age is recorded as x0Wherein
Figure GDA0003241410480000082
Figure GDA0003241410480000083
The maximum age at which humans can theoretically survive; the hour notes
Figure GDA0003241410480000084
For male group data, applying statistical software to perform regression analysis on the telomere length of the samples of the subjects and the corresponding age of the subjects to obtain the following formula:
TRF1=a1x0+b1
for the female group data, applying statistical software to perform regression analysis on the telomere length of the subject samples and the corresponding age of the subject to obtain the following formula:
TRF0=a0x0+b0
3) and respectively carrying out telomere length detection once when the subject is at the age of t and t + k, wherein the telomere lengths of corresponding samples are TRF(t)、TRF(t+k)And sequentially executing the following substeps:
A. calculating cell age of the sample
Cellular age CA at t + k years of age of subject(t+k)Comprises the following steps:
Figure GDA0003241410480000091
when actually calculating, if CA(t+k)T + k, i.e., the cell age is less than the actual age, means that the subject is younger at the cellular level than the peer, i.e., the subject is more viable than the peer. If CA(t+k)> t + k, i.e. cell age greater than actual age, means that the subject is older than the peer at the "cell level", i.e. the subject is not viable, when targeted conditioning and post-exercise needs to be considered, in order to expect the cell age CA at the next test(t+k)The value can be as small as or equal to the current age value.
B. Calculating the average change degree of telomere length of the sample
Mean degree of change of telomere length in the subject from t to t + k years of aget,kComprises the following steps:
Figure GDA0003241410480000092
in actual calculation, if | Δt,kA larger value of | indicates a faster rate of telomere shortening, indicating a faster subject aging. And if | Δt,k|≤|Δt,0If yes, the subject is shown to be senescent slower than the peer, otherwise, the subject is shown to be senescent faster than the peer. Mean degree of telomere length change in age t to t + k by subjectt,kIt can also be determined whether the subject needs to be considered for targeted conditioning and late exercise to reduce the rate of telomere shortening.
C. Method for establishing telomere information code
Note the book
Figure GDA0003241410480000101
Then m is5Is an integer of not less than 0 and not more than 35;
recording min (x, y) as the smaller number of x and y;
Figure GDA0003241410480000102
then m is4Is an integer of not less than 0 and not more than 35.
Calculate m5And m4Then m is added5And m4Is converted into a thirty-six system symbol according to a thirty-six system, and is converted into a thirty-six system symbol according to m5 m4The two-bit character string is obtained when the two-bit character string is spliced together, and the two-bit character string is the telomere information code.
During decoding of the session information code, if the first bit character is located earlier (or m) in the thirty-six system character conversion table5Smaller values) indicating better subject cell viability (or "younger" subject cells); if the second digit character is positioned more forward in the thirty-six system character conversion table (or in other wordsm4Smaller values) indicating that the subject is aging more slowly.
The check code is obtained as follows:
since the dynamic code and the static code are all known at this time, that is, the 2 nd to 17 th bits from right to left as shown in fig. 1 are known, the 2 th to 17 th bits are set as:
M18M17M16M15M14M13M12M11M10M9M8M7M6M5M4M3M2
wherein M isjRepresents a character, j is more than or equal to 2 and less than or equal to 18;
step1 converting M according to thirty-six system character conversion tablejReverse conversion to a number mj
Step2 for each mjCalculating the weight ω according to the following formulaj
ωj=mod(2j-1,37)
Wherein mod (x,37) represents the remainder of the integer x divided by 37;
step3 calculating m according to the following formula1
Figure GDA0003241410480000103
Step4 converting the number m according to the thirty-six system character conversion table1Conversion to character M1The M being1Namely the check code.
To facilitate a further understanding of the invention, the following specific examples are given herein:
example 1:
the subject A is a male, the nationality is China, the native language is Chinese, the face image is collected in Beijing at 2017, 5, 12, 9 and 8 minutes, and the 14 th male is collected at the time point. The telomere length measured at age 45 was 11.59kb, and at age 47 it was 11.45 kb. At the age of 47, the first was assigned a credit rating of 85 percent and a health rating of 76.3 percent.
According to the steps of the invention, the following are obtained:
the country code is CN;
the language code is ZH;
the 14-digit number in the image acquisition sequence code calculation step is 00170512090827, and the corresponding image acquisition sequence code is 026 BYHVBF;
m in the telomere code calculation step5=10,m 416, the corresponding telomere code is AG;
m in the step of calculating a credit code 330, the corresponding credit code is U;
m in the health code calculation step 227, the corresponding credit code is R;
m in the check code calculation step 126, the corresponding check code is Q;
the full code of the nail at the age of 47 years is CNZH026 BYHVBFAGURQ.
If the first full code is input into the system to be CNZH2LSBKYOZFAGURQ, the system calculates the first 17 bits of the full code according to the check code calculation rule to obtain the check code bit which is U, and the system prompts that the input is wrong different from Q, thereby realizing the check function.
Example 2:
subject B, Canada, English, was collected in 23 o' clock 23/7/1/2018 Beijing as the face image, which was the 22 nd female subject at this time point. Telomeres of 12.03kb were measured at age 44, and 11.66kb at age 46. At the age of 47, the first was rated 79 (percent) for credit and 90 (percent) for health.
According to the steps of the invention, the following are obtained:
the country code is CA;
the language code is EN;
the 14-bit number in the image acquisition sequence code calculation step is 00180701234542, and the corresponding image acquisition sequence code is 02B0 GUKXA;
m in the telomere code calculation step5=9,m 435, the corresponding telomere code is 9Z;
credit codeM in the calculation step 328, the corresponding credit code is S;
m in the health code calculation step 232, the corresponding credit code is W;
m in the check code calculation step 112, the corresponding check code is C;
the full code obtained for nail at age 47 is CAEN02B0GUKXA9 ZSHC.
If the first full code is CAEM02B0GUKXA9 ZSVC, the system calculates the first 17 bits of the full code according to the check code calculation rule to obtain that the check code bit is G, and the system prompts that the input is wrong, which is different from C, thereby realizing the check function.
Example 3 (when only the telomere information code and the credit code are used within the dynamic code, and no health code is included):
the subject A is a male, the nationality is China, the native language is Chinese, the face image is collected in Beijing at 2017, 5, 12, 9 and 8 minutes, and the 14 th male is collected at the time point. The telomere length measured at age 45 was 11.59kb, and at age 47 it was 11.45 kb. The credits for the nail were 85 (percent) at the age of 47.
According to the steps of the invention, the following are obtained:
the country code is CN;
the language code is ZH;
the 14-digit number in the image acquisition sequence code calculation step is 00170512090827, and the corresponding image acquisition sequence code is 026 BYHVBF;
m in the telomere code calculation step5=10,m 416, the corresponding telomere code is AG;
m in the step of calculating a credit code 330, the corresponding credit code is U;
m in the check code calculation step 14, the corresponding check code is 4;
the full code of the nail at the age of 47 years is obtained as CNZH026BYHVBFAGU 4.
If the first full code is input in the system as CNZH27 TZLL 37AGU4, the system calculates the first 16 bits of the full code according to the check code calculation rule to obtain that the check code bit is O, and the system prompts that the input is wrong different from 4, thereby realizing the check function.
Example 4 (when only the telomere information code and the health code are used within the dynamic code, and no credit code is included):
the first subject is male, the nationality is China, the native language is Chinese, and the human face image is collected at 5, 12 and 9 in 2017 of Beijing, which is the 14 th male subject at the time point. The telomere length measured at age 45 was 11.59kb, and at age 47 it was 11.45 kb. The nail health score was 76.3 (percent) at 47 years of age.
According to the steps of the invention, the following are obtained:
the country code is CN;
the language code is ZH;
the 11-digit number in the image acquisition sequence code calculation step is 00170512090827, and the corresponding image acquisition sequence code is 026 BYHVBF;
m in the telomere code calculation step5=10,m 416, the corresponding telomere code is AG;
m in the health code calculation step 227, the corresponding credit code is R;
m in the check code calculation step 134, the corresponding check code is Y;
the full code of the nail at the age of 47 years is obtained as CNZH026 BYHVBFAGRY.
If the input full code of the nail is CNZH27 TZLL 37AGRY in the system, the system calculates the first 16 bits of the full code according to the calculation rule of the check code to obtain that the check code bit is I, and the system prompts that the input is wrong when the check code bit is different from Y, thereby realizing the check function.

Claims (6)

1. An identity information coding method based on telomere length is characterized in that: the method at least comprises a static code and a dynamic code, wherein:
the static code comprises a country code, a language code and a face image acquisition sequence code of a subject;
the dynamic codes comprise a telomere information code, a credit code and a health code of the subject; the telomere information code is two bits, and the acquisition process comprises the following steps:
1) calculating the length of the telomere of the sample of the subject;
2) establishing a sample database according to the existing telomere lengths of the subjects with different ages and different sexes;
dividing the sample database into two groups according to gender, wherein the first group is male, and the second group is female; and the age is recorded as x0Wherein
Figure FDA0003241410470000011
Figure FDA0003241410470000012
The maximum age at which humans can theoretically survive; the hour notes
Figure FDA0003241410470000013
For male group data, applying statistical software to perform regression analysis on the telomere length of the samples of the subjects and the corresponding age of the subjects to obtain the following formula:
TRF1=a1x0+b1
for the female group data, applying statistical software to perform regression analysis on the telomere length of the subject samples and the corresponding age of the subject to obtain the following formula:
TRF0=a0x0+b0
3) and respectively carrying out telomere length detection once when the subject is at the age of t and t + k, wherein the telomere lengths of corresponding samples are TRF(t)、TRF(t+k)And sequentially executing the following substeps:
A. calculating cell age of the sample
Cellular age CA at t + k years of age of subject(t+k)Comprises the following steps:
Figure FDA0003241410470000014
B. calculating the average change degree of telomere length of the sample
Mean degree of change of telomere length in the subject from t to t + k years of aget,kComprises the following steps:
Figure FDA0003241410470000021
C. method for establishing telomere information code
Recording the Gaussian function as the maximum integer of Gauss (x) not more than x;
Figure FDA0003241410470000022
then m is5Is an integer of not less than 0 and not more than 35;
recording min (x, y) as the smaller number of x and y;
Figure FDA0003241410470000023
then m is4Is an integer of not less than 0 and not more than 35;
calculate m5And m4Then m is added5And m4Is converted into a thirty-six system symbol according to a thirty-six system, and is converted into a thirty-six system symbol according to m5 m4The two-bit character string is obtained when the two-bit character string is spliced together, and the two-bit character string is the telomere information code.
2. The identity information encoding method based on telomere length as claimed in claim 1, wherein: in the step 1), calculating the length of the telomere of the sample of the subject comprises the following sub-steps:
a. 3-5 ml of peripheral blood is extracted and put into 0.5mol/L EDTAK2Freezing and storing in an anticoagulant centrifuge tube at-20 ℃;
b. extracting peripheral blood leukocyte genome DNA by adopting a saturated phenol-chloroform method, selecting a sample with an OD260/OD280 value of 1.8-2.0, and observing the integrity of the DNA under ultraviolet rays after 0.8% agarose gel electrophoresis;
c. performing mixed enzyme digestion on genomic DNA for 2h by using a Telo TAGGG Telomere Length Assay kit of Roche, performing electrophoresis for 4h in 0.8% agarose gel at 5V/cm for 5 h, denaturing and neutralizing the gel, transferring the DNA to a positively charged nylon membrane by a siphon blotting method, and baking for 20min at 120 ℃; pre-hybridizing at 42 ℃ for 1h, and adding a digoxin-labeled telomere probe for hybridization for 3 h; and (3) carrying out chemiluminescence development after washing the membrane, carrying out optical density scanning imaging on the obtained X-ray film, analyzing by using telomere length analysis software Telomeric1.2, and calculating the telomere length of each sample according to the following formula:
Figure FDA0003241410470000031
wherein L isiIs the i-point marker length;
ODithe optical density value of point i.
3. The identity information encoding method based on telomere length as claimed in claim 1, wherein: the country code is two digits and refers to ISO 3166-1; the language code is two digits, and refers to ISO 693-1.
4. The identity information coding method based on the telomere length as claimed in claim 1, 2 or 3, wherein: the acquisition process of the face image acquisition sequence code comprises the following substeps:
(1) acquiring time including five parts of year, month, date, clock and minute when the human face image is acquired by the subject; the ones of each part are complemented into two digits on the left side by 0;
(2) subtracting 2000 from year to obtain four digits, and if the four digits are not enough, the four digits are compensated by 0 on the left side; arranging the quartic, month, date, clock and minute in the original sequence;
(3) supplementing a two-digit sequence number so as to form 14-digit numbers together with the numbers sequentially arranged in the substep (2); the sequence number is determined according to the sequence of the images collected by the testee in the same minute and the sex of the images; specifically, if the subject is male, the last digit of the order number is odd, and the left side of the two-digit deficiency is supplemented with 0 to form two digits; if the subject is female, the last digit is an even number, and the left side of the two digits is supplemented with 0 to form two digits;
(4) and carrying out thirty-six system conversion on the 14-bit numbers in the substep (3) to form nine-bit face image acquisition sequence codes, wherein the left side of less than nine bits is complemented with 0 to form nine bits.
5. The identity information coding method based on the telomere length as claimed in claim 4, wherein: the credit code is one bit, the health code is one bit, and the formation steps are as follows:
when the subject credits are Cre and the credit full is Cremax(ii) a Subjects were classified as Hea healthy and Hea full healthymaxWhen the current is over;
calculating m by calculating Gauss function as the maximum integer not exceeding x3
Figure FDA0003241410470000041
Then m is put3Converting into symbols in thirty-six system according to thirty-six system to obtain corresponding credit codes;
and m is calculated according to the following formula2
Figure FDA0003241410470000042
Then m is put2And converting into symbols in thirty-six system according to the thirty-six system to obtain the corresponding health codes.
6. The identity information coding method based on the telomere length as claimed in claim 4, wherein: the method also comprises a check code for detecting whether the static code and the dynamic code are mistakenly transmitted or not after the static code and the dynamic code are sequentially arranged from left to right; the check code is obtained as follows:
when the dynamic code and the static code are all known, setting the code as the dynamic code
M18M17M16M15M14M13M12M11M10M9M8M7M6M5M4M3M2
Wherein M isjRepresents a character, j is more than or equal to 2 and less than or equal to 18;
step1 converting M according to thirty-six system character conversion tablejReverse conversion to a number mj
Step2 for each mjCalculating the weight ω according to the following formulaj
ωj=mod(2j-1,37)
Wherein mod (x,37) represents the remainder of the integer x divided by 37;
step3 calculating m according to the following formula1
Figure FDA0003241410470000043
Step4 converting the number m according to the thirty-six system character conversion table1Conversion to character M1The M being1Namely the check code.
CN201811001443.7A 2018-08-30 2018-08-30 Identity information coding method based on telomere length Active CN109273051B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811001443.7A CN109273051B (en) 2018-08-30 2018-08-30 Identity information coding method based on telomere length

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811001443.7A CN109273051B (en) 2018-08-30 2018-08-30 Identity information coding method based on telomere length

Publications (2)

Publication Number Publication Date
CN109273051A CN109273051A (en) 2019-01-25
CN109273051B true CN109273051B (en) 2022-01-18

Family

ID=65154809

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811001443.7A Active CN109273051B (en) 2018-08-30 2018-08-30 Identity information coding method based on telomere length

Country Status (1)

Country Link
CN (1) CN109273051B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110334260B (en) * 2019-06-03 2024-05-07 平安科技(深圳)有限公司 Data analysis method, device, computer equipment and storage medium
CN110408684A (en) * 2019-08-22 2019-11-05 陕西九州医学检验有限公司 A kind of telomere length detection kit and method and biological age evaluation method
CN110974177B (en) * 2019-12-23 2022-09-09 南昌五元生物科技有限公司 Biological age calculation method and system based on telomere length

Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1306266A (en) * 2001-02-26 2001-08-01 华西医科大学法医学院 Gene identity card and its preparing process
CN102592158A (en) * 2011-01-11 2012-07-18 郝明学 Identity card coding method for global pets and other artificial fed animals
CN102618633A (en) * 2012-02-03 2012-08-01 常州亚当生物技术有限公司 Method for detecting length of the shortest telomere in cells by using improved STELA method
CN103233071A (en) * 2013-04-27 2013-08-07 南京优而生物科技发展有限公司 Method for measuring telomere absolute length
CN103412854A (en) * 2013-08-01 2013-11-27 汪风珍 Chinese identity card number coding system
CN103413221A (en) * 2012-07-01 2013-11-27 宾劲松 Member coding method with identity feature identification and application thereof
CN104694641A (en) * 2015-02-13 2015-06-10 冯文峰 Method for predicating gene age and disease susceptibility and kit
CN104769134A (en) * 2012-09-11 2015-07-08 赛拉诺斯股份有限公司 Information management systems and methods using a biological signature
CN104899521A (en) * 2015-06-08 2015-09-09 深圳市华傲数据技术有限公司 Methods for bleaching and reverse bleaching of 18-bits identity card number
CN105184083A (en) * 2015-09-14 2015-12-23 南京延长科技有限公司 Digital identity code used for medical instrument whole-process quality tracing and coding method thereof
CN106548435A (en) * 2015-09-22 2017-03-29 联邦应用基因股份有限公司 Cell Age analytical integration system
CN106755429A (en) * 2016-12-27 2017-05-31 上海三誉生物科技有限公司 A kind of method of accurate measurement crowd telomere length
CN107058591A (en) * 2017-06-16 2017-08-18 张晓� Telomere length detection method and kit

Patent Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1306266A (en) * 2001-02-26 2001-08-01 华西医科大学法医学院 Gene identity card and its preparing process
CN102592158A (en) * 2011-01-11 2012-07-18 郝明学 Identity card coding method for global pets and other artificial fed animals
CN102618633A (en) * 2012-02-03 2012-08-01 常州亚当生物技术有限公司 Method for detecting length of the shortest telomere in cells by using improved STELA method
CN103413221A (en) * 2012-07-01 2013-11-27 宾劲松 Member coding method with identity feature identification and application thereof
CN104769134A (en) * 2012-09-11 2015-07-08 赛拉诺斯股份有限公司 Information management systems and methods using a biological signature
CN103233071A (en) * 2013-04-27 2013-08-07 南京优而生物科技发展有限公司 Method for measuring telomere absolute length
CN103412854A (en) * 2013-08-01 2013-11-27 汪风珍 Chinese identity card number coding system
CN104694641A (en) * 2015-02-13 2015-06-10 冯文峰 Method for predicating gene age and disease susceptibility and kit
CN104899521A (en) * 2015-06-08 2015-09-09 深圳市华傲数据技术有限公司 Methods for bleaching and reverse bleaching of 18-bits identity card number
CN105184083A (en) * 2015-09-14 2015-12-23 南京延长科技有限公司 Digital identity code used for medical instrument whole-process quality tracing and coding method thereof
CN106548435A (en) * 2015-09-22 2017-03-29 联邦应用基因股份有限公司 Cell Age analytical integration system
CN106755429A (en) * 2016-12-27 2017-05-31 上海三誉生物科技有限公司 A kind of method of accurate measurement crowd telomere length
CN107058591A (en) * 2017-06-16 2017-08-18 张晓� Telomere length detection method and kit

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
"Estimation of Human Age According to Telomere Shortening in Peripheral Blood Leukocytes of Tibetan";Fu Ren 等;《ResearchGate》;20090930;第252-255页 *
"Telo TAGGG Telomere Length Assay";lele-kele;《豆丁网》;20140813;第1-22页 *
"The relationship between telomere length and beekeeping among Malaysians";Nurul Fatihah Mohamad Nasir 等;《Springer》;20151231;第1-6页 *
"人外周血白细胞端粒DNA长度与年龄的定量关系";任甫 等;《中国法医学杂志》;20071231;第22卷(第3期);第160-162页 *
"人外周血白细胞端粒DNA长度与年龄相关性的法医学应用研究";薄云峰;《万方数据知识服务平台》;20070921;第1-31页 *

Also Published As

Publication number Publication date
CN109273051A (en) 2019-01-25

Similar Documents

Publication Publication Date Title
CN109273051B (en) Identity information coding method based on telomere length
Berry et al. The politics of tax increases in the states
KURTE SEX DIMORPHISM AND SIZE TRENDS IN THE CA VE BEAR, URSUS SPELAEUS ROSENMÜLLER AND HEINROTH
Gao et al. The all-data-based evolutionary hypothesis of ciliated protists with a revised classification of the phylum Ciliophora (Eukaryota, Alveolata)
CN110246577B (en) Method for assisting gestational diabetes genetic risk prediction based on artificial intelligence
Orysiak et al. The association between ace gene variation and aerobic capacity in winter endurance disciplines
CN105349645A (en) Roe and fry quick classification and determination methods based on CoI DNA bar codes
CN109234358B (en) Telomere information coding method
Balakirev et al. Taxonomic revision of Niviventer (Rodentia, Muridae) from Vietnam: a morphological and molecular approach
Marshall et al. Social environment and breast cancer. A cohort analysis of patient survival
CN110874409A (en) Disease grading prediction system, method, electronic device and readable storage medium
CN115115620B (en) Pneumonia lesion simulation method and system based on deep learning
CN113593714A (en) Method, system, equipment and medium for detecting multi-classification new coronary pneumonia cases
CN112288749A (en) Skull image segmentation method based on depth iterative fusion depth learning model
Grismer et al. A new species of crocodile newt Tylototriton (Caudata: Salamandridae) from northern Myanmar (Burma)
CN102063568A (en) Individual diabetes mellitus prediction model
Salzano Demographic and genetic interrelationships among the Cayapo Indians of Brazil
Weiss et al. Temperature and barometric pressure are related to running speed and pacing of the fastest runners in the'Berlin Marathon'.
BR112020008566A2 (en) method to estimate inflammation area of periodontal pockets
CN112802605A (en) Prediction model for survival benefit of metastatic renal cancer patient after receiving system treatment and establishment method and application thereof
Poulain et al. Anthropometric traits at military medical examinations associated with demographic family characteristics in Sardinia at the turn of twentieth century
CN106446601A (en) Method for labeling lncRNA functions in large scale
CN114093427B (en) Antiviral peptide prediction method based on deep learning and machine learning
Wang Convergence and Prospects
Sagayama et al. Total energy expenditure in elite open-water swimmers

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant