Time series study of a 17-year record of (7)Be and (210)Pb fluxes in northern Taiwan using ensemble empirical mode decomposition.

Using the ensemble empirical mode decomposition (EEMD) and a significance test method, we have analyzed time series data on the fluxes of (7)Be and (210)Pb collected over a span of 17 y in Northern Taiwan. Among nine intrinsic mode functions (IMFs) extracted from the method five (IMF4-8) are non-trivial for (210)Pb and have adequate S/N with significant power in localized windows around the periodicities of 0.5 y, 1 y, 2 y, 5 y, and 11 y, respectively. For (7)Be, IMF5 and IMF8 with periods around 1 y and 11 y, respectively, have adequate S/N. The semi-annual and annual cycles represented by IMF4 and IMF5, respectively, are dominated by East Asian monsoon. The sum of IMF6 and IMF7 reveals an inter-annual cycle where both (7)Be and (210)Pb fluxes are well-correlated with the East Asian winter monsoon index (EAWMI). The close tracking of the (210)Pb and (7)Be in IMF8 cases may reflect an 11 y cycle; implying that it is caused by common climatologic factors, likely related to solar cycle, rather than their distinct production modes.

Application of chromosome 4q35-qter marker (pFR-1) for DNA rearrangement of facioscapulohumeral muscular dystrophy patients in Taiwan.

Facioscapulohumeral muscular dystrophy (FSHD) has been found to be linked to chromosome 4qter. A chromosome 4q35-ter marker, pFR-1 (subclone of the cosmid c51), has been recently isolated and used as a probe for mapping near, or within, the FSHD gene. To examine FSHD-associated DNA rearrangements in the Taiwan population, we used the pFR-1 probe to perform Southern blot analysis on 142 individuals, including 32 FSHD patients within 9 autosomal dominant families, five sporadic FSHD patients from 4 families (include one pair of twins), three sporadic scapuloperoneal syndrome (SPS) patients and two sporadic polymyositis patients with their unaffected parents, and 29 healthy controls. In 29 healthy individuals, 3 SPS and 2 polymyositis patients with their families, probe pFR-1 analysis revealed that all had polymorphic restriction fragments that were larger than 28 kb in length. All but 1 FSHD-affected individual had specific smaller EcoRI fragments (ranging in size from 10.5 to 27 kb). Two point linkage analysis between pFR-1 and the FSHD locus provided significant evidence for FSHD linkage (Z(max)=6.84). A similar smaller fragment was also present in 5 sporadic patients, while this smaller fragment could not be found in one of their parents. Identical EcoRI restriction fragment length polymorphism (RFLP) patterns linked to FSHD were shown in the monozygotic twins, even though they showed extreme variability in the expression of FSHD. We conclude that the pFR-1 probe is a tightly linked marker of FSHD and can be used to detect most DNA rearrangements associated with this disease in the Taiwan population. However, the same RFLP patterns may represent extreme variability in the expression of the FSHD gene.

Chinese population data on three tetrameric short tandem repeat loci--HUMTHO1, TPOX, and CSF1PO--derived using multiplex PCR and manual typing.

Allele and genotype frequencies for three tetrameric short tandem repeat (STR) loci were determined in a Chinese sample population using multiplex polymerase chain reaction (PCR), electrophoresis of the PCR products in DNA sequencing gels and subsequent detection by silver staining. The loci are HUMTHO1, TPOX, and CSF1PO. All loci meet Hardy-Weinberg expectations. In addition, there is no evidence for association of alleles among the three loci. The allelic frequency data can be used in human identity testing to estimate the frequency of a multiple STR locus DNA profile in the Chinese population.

Cytochrome b gene for species identification of the conservation animals.

A partial DNA sequence of cytochrome b gene was used to identify the remains of endangered animals and species endemic to Taiwan. The conservation of animals species included in this study were: the formosan gem-faced civets, leopard cats, tigers, clouded leopards, lion, formosan muntjacs, formosan sika deers, formosan sambars, formosan serows, water buffalo, formosan pangolins and formosan macaques. The control species used included domestic cats, domestic dogs, domestic sheeps, domestic cattles, domestic pigs and humans. Heteroplasmy was detected in the formosan macaque, domestic pig and domestic cats. The frequencies of heteroplasmy in these animals were about 0.25% (1 in 402bp). Sequences were aligned by Pileup program of GCG computer package, and the phylogenetic tree was constructed by the neighbor-joining method. The results of sequence comparison showed that the percentage range of sequence diversity in the same species was from 0.25 to 2.74%, and that between the different species was from 5.97 to 34.83%. The results of phylogenetic analysis showed that the genetic distance between the different species was from 6.33 to 40.59. Animals of the same species, both the endangered animal species and domestic animals, were clustered together in the neighbor-joining tree. Three unknown samples of animal remains were identified by this system. The partial sequence of cytochrome b gene adopted in this study proved to be usable for animal identification.

Fixed bin population data for the VNTR loci D1S7, D2S44, D4S139, D5S110, and D17S79 in Chinese from Taiwan.

Fixed bin frequencies for the VNTR loci D1S7, D2S44, D4S139, D5S110, and D17S79 were determined in a Chinese sample population. The data were generated by RFLP analysis of Hae III-digested genomic DNA and chemiluminescence detection. The five VNTR loci meet Hardy-Weinberg expectations in the Chinese sample population, and there is little evidence for association of alleles between the VNTR loci. The frequency data can be used in forensic analyses and paternity tests to estimate the frequency of a DNA profile in Chinese.

Noise and poise: Enhancement of postural complexity in the elderly with a stochastic-resonance-based therapy.

Pathologic states are associated with a loss of dynamical complexity. Therefore, therapeutic interventions that increase physiologic complexity may enhance health status. Using multiscale entropy analysis, we show that the postural sway dynamics of healthy young and healthy elderly subjects are more complex than that of elderly subjects with a history of falls. Application of subsensory noise to the feet has been demonstrated to improve postural stability in the elderly. We next show that this therapy significantly increases the multiscale complexity of sway fluctuations in healthy elderly subjects. Quantification of changes in dynamical complexity of biologic variability may be the basis of a new approach to assessing risk and to predicting the efficacy of clinical interventions, including noise-based therapies.

Chinese population data on the PCR-based loci HLA-DQ alpha, low-density-lipoprotein receptor, glycophorin A, hemoglobin gamma G, D7S8, and group-specific component.

Allele and genotype frequencies for six polymerase chain reaction (PCR)-based DNA genetic markers were determined in a Chinese sample population. The loci are HLA-DQ alpha, low-density-lipoprotein receptor, glycophorin A, hemoglobin gamma G, D7S8, and group-specific component. These results were obtained using the AmpliType HLA-DQ alpha forensic DNA amplification and typing kit and the AmpliType PM PCR amplification and typing kit. All loci met Hardy-Weinberg expectations and there was no evidence for association of alleles between the loci. The frequency data can be used in forensic analyses and paternity tests to estimate the frequency of a DNA profile in the Chinese population.

D1S80 allele frequencies in a Chinese population.

Allele frequencies for the VNTR locus D1S80 were determined in a Chinese population sample using the polymerase chain reaction and subsequent analysis of the amplified products by polyacrylamide gel electrophoresis and silver staining. A total of 18 nominal D1S80 alleles were observed in 105 unrelated Chinese. The data demonstrate that D1S80 is highly polymorphic in Chinese with a heterozygosity of 90.5%. The D1S80 frequency distribution meets Hardy-Weinberg expectations. This D1S80 data can be used in forensic analyses and paternity tests to estimate the frequency of a DNA profile in a Chinese population.

Nonlinear indicial response of complex nonstationary oscillations as pulmonary hypertension responding to step hypoxia.

This paper is devoted to the quantization of the degree of nonlinearity of the relationship between two biological variables when one of the variables is a complex nonstationary oscillatory signal. An example of the situation is the indicial responses of pulmonary blood pressure (P) to step changes of oxygen tension (DeltapO2) in the breathing gas. For a step change of DeltapO2 beginning at time t1, the pulmonary blood pressure is a nonlinear function of time and DeltapO2, which can be written as P(t-t1 | DeltapO2). An effective method does not exist to examine the nonlinear function P(t-t1 | DeltapO2). A systematic approach is proposed here. The definitions of mean trends and oscillations about the means are the keys. With these keys a practical method of calculation is devised. We fit the mean trends of blood pressure with analytic functions of time, whose nonlinearity with respect to the oxygen level is clarified here. The associated oscillations about the mean can be transformed into Hilbert spectrum. An integration of the square of the Hilbert spectrum over frequency yields a measure of oscillatory energy, which is also a function of time, whose mean trends can be expressed by analytic functions. The degree of nonlinearity of the oscillatory energy with respect to the oxygen level also is clarified here. Theoretical extension of the experimental nonlinear indicial functions to arbitrary history of hypoxia is proposed. Application of the results to tissue remodeling and tissue engineering of blood vessels is discussed.

Use of intrinsic modes in biology: examples of indicial response of pulmonary blood pressure to +/- step hypoxia.

Recently, a new method to analyze biological nonstationary stochastic variables has been presented. The method is especially suitable to analyze the variation of one biological variable with respect to changes of another variable. Here, it is illustrated by the change of the pulmonary blood pressure in response to a step change of oxygen concentration in the gas that an animal breathes. The pressure signal is resolved into the sum of a set of oscillatory intrinsic mode functions, which have zero "local mean," and a final nonoscillatory mode. With this device, we obtain a set of "mean trends," each of which represents a "mean" in a definitive sense, and together they represent the mean trend systematically with different degrees of oscillatory content. Correspondingly, the oscillatory content of the signal about any mean trend can be represented by a set of partial sums of intrinsic mode functions. When the concept of "indicial response function" is used to describe the change of one variable in response to a step change of another variable, we now have a set of indicial response functions of the mean trends and another set of indicial response functions to describe the energy or intensity of oscillations about each mean trend. Each of these can be represented by an analytic function whose coefficients can be determined by a least-squares curve-fitting procedure. In this way, experimental results are stated sharply by analytic functions.

Engineering analysis of biological variables: an example of blood pressure over 1 day.

Almost all variables in biology are nonstationarily stochastic. For these variables, the conventional tools leave us a feeling that some valuable information is thrown away and that a complex phenomenon is presented imprecisely. Here, we apply recent advances initially made in the study of ocean waves to study the blood pressure waves in the lung. We note first that, in a long wave train, the handling of the local mean is of predominant importance. It is shown that a signal can be described by a sum of a series of intrinsic mode functions, each of which has zero local mean at all times. The process of deriving this series is called the "empirical mode decomposition method." Conventionally, Fourier analysis represents the data by sine and cosine functions, but no instantaneous frequency can be defined. In the new way, the data are represented by intrinsic mode functions, to which Hilbert transform can be used. Titchmarsh [Titchmarsh, E. C. (1948) Introduction to the Theory of Fourier Integrals (Oxford Univ. Press, Oxford)] has shown that a signal and i times its Hilbert transform together define a complex variable. From that complex variable, the instantaneous frequency, instantaneous amplitude, Hilbert spectrum, and marginal Hilbert spectrum have been defined. In addition, the Gumbel extreme-value statistics are applied. We present all of these features of the blood pressure records here for the reader to see how they look. In the future, we have to learn how these features change with disease or interventions.

Identification of victims of the 1998 Taoyuan Airbus crash accident using DNA analysis.

In February 1998 a civilian aeroplane carrying 196 individuals crashed in Taiwan and killed another 6 people on the ground. Although there were dental and medical records, fingerprints, photographic evidence and personal effects to identify some of the victims, DNA analysis was required to further identify severely damaged remains. From the 202 people known to have perished in the plane crash, a total of 685 fragments of human remains were subjected to DNA analysis. The analysis was carried out using nine microsatellite loci, plus amelogenin to cluster the 685 fragments into 202 groups, accounting for all the victims. To establish genetic relatedness of the victims to other victims and living relatives, additional DNA loci were used. In this case the paternity index was increased by using HLA DQA1 plus Polymarker. The same 16 DNA loci were used to test blood samples from 201 relatives to establish parent/child and sibling relationships. With the exception of 19 victims identified by non-genetic evidence, 183 victims were successfully identified by DNA typing with relatively high values of paternity index by the direct or indirect comparison of relatives. The 202 victims were from 37 different families, ranging in size from 2 to 13 members and 74 individuals known to be unrelated to any other victim. The DNA from living relatives was used to identify one member of a family group, from which other victims of the family could be identified. ABO blood group information was further used to confirm genetic relatedness within families. A comparison of the DNA profiling results to the ABO blood group of the victims showed no discrepancies with the exception of two mutations in the FGA locus. In cases of severely damaged victims from a plane crash, DNA analysis proved to be the best choice to identify victims.

Haplotype frequencies of nine Y-chromosome STR loci in the Taiwanese Han population.

Haplotype frequencies of nine Y-chromosome STR loci (DYS19, DYS385, DYS388, DYS389I, DYS389II, DYS390, DYS391, DYS392 and DYS393) in the Taiwanese Han population were established. A total of 183 unrelated individuals produced 162 haplotypes, of which 146 were unique, 1 was found in 5 individuals, 2 were found in 3 individuals and 13 were found in 2 individuals. The haplotype diversity (99.99%) and discrimination capacity (88.5%) were calculated. A family study of 109 father/son pairs in 100 families showed 2 mutational events in the DYS389II locus and 1 in the DYS392 locus.