Hong-Xiang Zheng

A Mitochondrial Revelation of Early Human Migrations to the Tibetan Plateau Before and After the Last Glacial Maximum

As the highest plateau surrounded by towering mountain ranges, the Tibetan Plateau was once considered to be one of the last populated areas of modern humans. However, this view has been tremendously changed by archeological, linguistic, and genetic findings in the past 60 years. Nevertheless, the timing and routes of entry of modern humans into the Tibetan Plateau is still unclear. To make these problems clear, we carried out high-resolution mitochondrial-DNA (mtDNA) analyses on 562 Tibeto-Burman inhabitants from nine different regions across the plateau. By examining the mtDNA haplogroup distributions and their principal components, we demonstrated that maternal diversity on the plateau reflects mostly a northern East Asian ancestry. Furthermore, phylogeographic analysis of plateau-specific sublineages based on 31 complete mtDNA sequences revealed two primary components: pre-last glacial maximum (LGM) inhabitants and post-LGM immigrants. Also, the analysis of one major pre-LGM sublineage A10 showed a strong signal of post-LGM population expansion (about 15,000 years ago) and greater diversity in the southern part of the Tibetan Plateau, indicating the southern plateau as a refuge place when climate dramatically changed during LGM.

Y Chromosomes of 40% Chinese Descend from Three Neolithic Super-Grandfathers

Demographic change of human populations is one of the central questions for delving into the past of human beings. To identify major population expansions related to male lineages, we sequenced 78 East Asian Y chromosomes at 3.9 Mbp of the non-recombining region, discovered >4,000 new SNPs, and identified many new clades. The relative divergence dates can be estimated much more precisely using a molecular clock. We found that all the Paleolithic divergences were binary; however, three strong star-like Neolithic expansions at ∼6 kya (thousand years ago) (assuming a constant substitution rate of 1×10−9/bp/year) indicates that ∼40% of modern Chinese are patrilineal descendants of only three super-grandfathers at that time. This observation suggests that the main patrilineal expansion in China occurred in the Neolithic Era and might be related to the development of agriculture.

mtDNA Lineage Expansions in Sherpa Population Suggest Adaptive Evolution in Tibetan Highlands

Sherpa population is an ethnic group living in south mountainside of Himalayas for hundreds of years. They are famous as extraordinary mountaineers and guides, considered as a good example for successful adaptation to low oxygen environment in Tibetan highlands. Mitochondrial DNA (mtDNA) variations might be important in the highland adaption given its role in coding core subunits of oxidative phosphorylation in mitochondria. In this study, we sequenced the complete mtDNA genomes of 76 unrelated Sherpa individuals. Generally, Sherpa mtDNA haplogroup constitution was close to Tibetan populations. However, we found three lineage expansions in Sherpas, two of which (C4a3b1 and A4e3a) were Sherpa-specific. Both lineage expansions might begin within the past hundreds of years. Especially, nine individuals carry identical Haplogroup C4a3b1. According to the history of Sherpas and Bayesian skyline plot, we constructed various demographic models and found out that it is unlikely for these lineage expansions to occur in neutral models especially for C4a3b1. Nonsynonymous mutations harbored in C4a3b1 (G3745A) and A4e3a (T4216C) are both ND1 mutants (A147T and Y304H, respectively). Secondary structure predictions showed that G3745A were structurally closing to other pathogenic mutants, whereas T4216C itself was reported as the primary mutation for Lebers hereditary optic neuropathy. Thus, we propose that these mutations had certain effect on Complex I function and might be important in the high altitude adaptation for Sherpa people.

Common variants of four bilirubin metabolism genes and their association with serum bilirubin and coronary artery disease in Chinese Han population.

Objectives Studies have revealed an inverse relationship between serum total bilirubin (TBIL) levels and coronary artery disease (CAD). This study investigated the genetic variants of four bilirubin metabolism genes — heme oxygenase-1 (HMOX1), biliverdin reductase A (BLVRA), solute carrier organic anion transporter family member 1B1 (SLCO1B1), and uridine diphosphate glycosyltransferase 1A1 (UGT1A1) — in relation to TBIL levels and CAD. Methods and results Thirty-five common single nucleotide polymorphisms (SNPs) were genotyped in 2380 unrelated Han participants who underwent angiocardiography at hospitals in Shanghai, China. Only three genetic variants — rs4399719 (UGT1A1 T-2473G), rs887829 (UGT1A1 G-364A), and rs4148323 (UGT1A1 G211A) — were associated with TBIL levels (each P<0.001). Four significant associations with CAD were detected after controlling age and the false discovery rate at 15%: the recessive effect of SNP rs887829 (UGT1A1 G-364A) [age-adjusted odds ratio (OR): 0.24; 95% confidence interval (CI): 0.10–0.60; P=0.0014] and dominant effect of rs4149013 (SLCO1B1 A-12099G) (age-adjusted OR: 0.70; 95% CI: 0.55–0.91; P=0.0069) on male CAD, and the additive effects of rs2877262 (BLVRA G+1238/in6C) (age-adjusted OR: 0.73; 95% CI: 0.59–0.89; P=0.0021) and rs2690381 (BLVRA G+2613/in6A) (age-adjusted OR: 0.70; 95% CI: 0.56–0.86; P=0.0008) on female CAD. SNPs rs2877262 and rs2690381 were both in a linkage disequilibrium block within BLVRA with r2 greater than 0.750. Correspondingly, this block was identified to be associated with female CAD. Conclusion Our study provides genetic evidences for the difference in the impact of these four bilirubin metabolism genes on TBIL levels and CAD.

Mitochondrial genomes and exceptional longevity in a Chinese population: the Rugao longevity study

Genetic variants of whole mitochondrial DNA (mtDNA) that predispose to exceptional longevity need to be systematically identified and appraised. Here, we conducted a case-control study with 237 exceptional longevity subjects (aged 95–107) and 444 control subjects (aged 40–69) randomly recruited from a “longevity town”—the city of Rugao in China—to investigate the effects of mtDNA variants on exceptional longevity. We sequenced the entire mtDNA genomes of the 681 subjects using a next-generation platform and employed a complete mtDNA phylogenetic analytical strategy. We identified T3394C as a candidate that counteracts longevity, and we observed a higher load of private nonsynonymous mutations in the COX1 gene predisposing to female longevity. Additionally, for the first time, we identified several variants and new subhaplogroups related to exceptional longevity. Our results provide new clues for genetic mechanisms of longevity and shed light on strategies for evaluating rare mitochondrial variants that underlie complex traits.

HIF2A Variants Were Associated with Different Levels of High-Altitude Hypoxia among Native Tibetans.

Hypoxia inducible factors, including HIF1A and HIF2A, play central roles in response to high-altitude hypoxia and genetic variants of HIF1A or HIF2A were associated with high-altitude sickness or adaptation. However, it remains to determine whether they are associated with tolerance to different levels of high-altitude selection pressure among native Tibetans. We recruited 189 Tibetan subjects living at 2,700 meters (Low level of high altitude, LHA), 197 at 3,200 meters (Middle level of high altitude of high altitude, MHA), 249 at 3,700 meters (High level of high altitude, HHA) and 269 at 4,700 meters (Very high level of high altitude, VHA) and performed association analysis of twelve tSNPs (tagging SNPs) in HIF1A and HIF2A with high-altitude. We found (1) a increasing trend of HIF2A rs5621780-C(18.4%, 15.9%, 32.8% and 31.1%, respectively, in LHA, MHA, HHA and VHA)(P = 3.56E-9); (2) increasing trends of HIF2A rs6756667-A(68.7%, 73.4%, 79.9% and 89.6%), rs7589621- G(74.6%, 77.9%, 83.7%, and 92.1%) and rs1868092-A(64.1%, 67.3%, 75.1% and 84.4%) (P = 3.56E-9, 4.68E-16, 1.17E-13 and 7.09E-14, respectively); (3) a increasing trend of haplotype AG (68.7%, 73.1%, 79.9% and 89.6%) (P = 2.22E-7) which was constructed by rs6756667 and rs7589621; (4) a strong linear correlation between major alleles of rs6756667-A (R 2 = 0.997, P = 0.002), rs7589621-G (R 2 = 0.994, P = 0.003), rs1868092-A (R 2 = 0.985, P = 0.008) and altitude by linear correlation test. The associations between HIF2A variants and different level of high altitude support that extremely high-altitude hypoxia challenge imposes selective effects on HIF2A variants among native Tibetans.

Positive selection on mitochondrial M7 lineages among the Gelong people in Hainan.

Selections on human mitochondrial variations are difficult to examine. In this study, we found possible signs of selection on mitochondrial M7 lineages among the Gelong people who migrated from Guizhou to Hainan (the hottest province in China) throughout the last 1000 years. The genetic structure of the Gelong people shows an obvious sex-biased population admixture pattern with only 4.9% paternal contribution but 30.7% maternal contribution from indigenous Hlai people. According to frequency spectrum tests for deviation from neutrality and mismatch tests of demographic expansion, part of the maternal mitochondrial M7 lineages among the Gelong came from the Hlai had spread quickly and therefore might have undergone positive selection. In the future, whole mitochondrial genome sequencing might reveal the functional advantage of the M7 lineages.

A novel mutation (C1425Y) in the FBN2 gene in a father and son with congenital contractural arachnodactyly.

Congenital contractural arachnodactyly (Beals syndrome) is a rare autosomal dominantly inherited connective tissue disorder characterized by flexion contractures, arachnodactyly, crumpled ears, and mild muscular hypoplasia. Here, a father and son with congenital contractural arachnodactyly features were identified. After sequencing 15 exons (22 to 36) of the FBN2 gene, a novel mutation (C1425Y) was found in exon 33. This de novo mutation presented first in the father and was transmitted to his son, but not in the other 14 unaffected family members and 365 normal people. The C1425Y mutation occurs at the 19th cbEGF domain. Cysteines in this cbEGF domain are rather conserved in species, from human down to ascidian. The cbEGF12-13 in human FBN1 was employed as the template to perform homology modeling of cbEGF18-19 of human FBN2 protein. The mutation has also been evaluated by further prediction tools, for example, SIFT, Blosum62, biochemical Yus matrice, and UMD-Predictor tool. In all analysis, the mutation is predicted to be pathogenic. Thus, the structure destabilization by C1425Y might be the cause of the disorder.

Resolving ambiguity in the phylogenetic relationship of genotypes A, B, and C of hepatitis B virus

BackgroundHepatitis B virus (HBV) is an important infectious agent that causes widespread concern because billions of people are infected by at least 8 different HBV genotypes worldwide. However, reconstruction of the phylogenetic relationship between HBV genotypes is difficult. Specifically, the phylogenetic relationships among genotypes A, B, and C are not clear from previous studies because of the confounding effects of genotype recombination. In order to clarify the evolutionary relationships, a rigorous approach is required that can effectively explore genetic sequences with recombination.ResultIn the present study, phylogenetic relationship of the HBV genotypes was reconstructed using a consensus phylogeny of phylogenetic trees of HBV genome segments. Reliability of the reconstructed phylogeny was extensively evaluated in agreements of local phylogenies of genome segments.The reconstructed phylogenetic tree revealed that HBV genotypes B and C had a closer phylogenetic relationship than genotypes A and B or A and C. Evaluations showed the consensus method was capable to reconstruct reliable phylogenetic relationship in the presence of recombinants.ConclusionThe consensus method implemented in this study provides an alternative approach for reconstructing reliable phylogenetic relationships for viruses with possible genetic recombination. Our approach revealed the phylogenetic relationships of genotypes A, B, and C of HBV.

Human mitochondrial DNA haplogroup M8a influences the penetrance of m.8684C>T in Han Chinese men with non-obstructive azoospermia

RESEARCH QUESTION What is the role of mitochondrial DNA (mtDNA) in the pathogenesis of non-obstructive azoospermia (NOA)? DESIGN mtDNA genome sequencing followed by an independent population validation were performed in 628 NOA cases and 584 healthy controls. Antioxidant capacity of serum was evaluated in 54 randomly selected cases out of 536 and 49 out of 489 controls. RESULTS In the screening stage, 13 mtDNA haplogroups (hg) were ascertained, and 10 susceptible variants were observed. In the validation stage, hg M8* in individuals was found to be associated with increased risk of NOA [odds ratio (OR) 2.61, 95% confidence interval (CI) 1.47-4.61] (P=0.001). Unexpectedly, the frequency of m.8684C>T, the defining marker for hg M8a, was also higher in NOA (OR 4.14, 95% CI 1.56-11.03) (P=0.002). Subsequently, the frequency distributions were compared among the sub-hg of hg M8* (including hg M8a, C and Z) and, intriguingly, no significance was found in hg C and Z. Additionally, the level of total antioxidant capacity was significantly decreased (P<0.05) compared with the control group. CONCLUSIONS hg M8a background in general played an active role in the penetrance of 8684C>T in NOA, and mtDNA genetic variants (causing low antioxidant levels) might increase mtDNA damage and impair normal spermatogenesis.