Geographic variation in the polygenic score of height in Japan.

A geographical gradient of height has existed in Japan for approximately 100 years. People in northern Japan tend to be taller than those in southern Japan. The differences in annual temperature and day length between the northern and southern prefectures of Japan have been suggested as possible causes of the height gradient. Although height is well known to be a polygenic trait with high heritability, the genetic contributions to the gradient have not yet been explored. Polygenic score (PS) is calculated by aggregating the effects of genetic variants identified by genome-wide association studies (GWASs) to predict the traits of individual subjects. Here, we calculated the PS of height for 10,840 Japanese individuals from all 47 prefectures in Japan. The median height PS for each prefecture was significantly correlated with the mean height of females and males obtained from another independent Japanese nation-wide height dataset, suggesting genetic contribution to the observed height gradient. We also found that individuals and prefectures genetically closer to continental East Asian ancestry tended to have a higher PS; modern Japanese people are considered to have originated as result of admixture between indigenous Jomon people and immigrants from continental East Asia. Another PS analysis based on the GWAS using only the mainland Japanese was conducted to evaluate the effect of population stratification on PS. The result also supported genetic contribution to height, and indicated that the PS might be affected by a bias due to population stratification even in a relatively homogenous population like Japanese.

Prefecture-level population structure of the Japanese based on SNP genotypes of 11,069 individuals.

We analyzed genome-wide single-nucleotide polymorphism data of 11,069 Japanese individuals recruited from all 47 prefectures of Japan to clarify their genetic structure. The principal component analysis at the prefectural level enabled us to study the relationship between geographical location and genetic differentiation. The results revealed that the mainland Japanese were not genetically homogeneous, and the genetic structure could be explained mainly by the degree of Jomon ancestry and the geographical location. One of the interesting findings was that individuals in the Shikoku region (i.e., Tokushima Prefecture, Kagawa Prefecture, Ehime Prefecture, and Kochi Prefecture) were genetically close to Han Chinese. Therefore, the genetic components of immigrants from continental East Asia in the Yayoi period may have been well maintained in Shikoku. The present results will be useful for understanding the peopling of Japan, and also provide suggestions for recruiting subjects in genetic association studies.

Molecular evolution of the semenogelin 1 and 2 and mating system in gibbons.

OBJECTIVES: Semenogelin 1 and 2 (SEMG1 and SEMG2) are known as semen coagulating proteins in primates with a repetitive structure of 60-amino acids. The number of repeats varies among species and is hypothesized to be related to the level of primate sperm competition. Gibbons until recently were thought to be monogamous primates, but it is now known that gibbon social structure is flexible. Thus, hypotheses of the relationship between the SEMGs evolution and mating systems were tested. MATERIALS AND METHODS: The sequences of the exon 2 of the SEMG1 and SEMG2 were obtained from 50 captive gibbons comprising six species belonging to three genera (Hylobates, Symphalangus, and Nomascus). Then we quantified the levels of polymorphism and estimated rates of protein evolution by calculating d N /d S ratio. RESULTS: Several mutations that create a premature stop codon in the SEMG1 and a reduction in the repeats of the SEMG2 in the genus Hylobates were observed and may alter the coding properties for these proteins. We also found different level of nucleotide diversity in each gene and between genera. Strikingly, in Nomascus leucogenys we discovered a high d N /d S ratio in the SEMG1 and SEMG2. The Nomascus SEMG2 also showed significantly lower nucleotide diversity than the other two genera. DISCUSSION: These results are consistent with the presence of a strong positive selection in the Nomascus lineage even if the exact selective forces acting on these genes are not yet conclusively known. We were not able to demonstrate, among gibbons, unambiguous relationships between the SEMGs evolution and mating systems.

Admixture with indigenous people helps local adaptation: admixture-enabled selection in Polynesians.

BACKGROUND: Homo sapiens have experienced admixture many times in the last few thousand years. To examine how admixture affects local adaptation, we investigated genomes of modern Polynesians, who are shaped through admixture between Austronesian-speaking people from Southeast Asia (Asian-related ancestors) and indigenous people in Near Oceania (Papuan-related ancestors). METHODS: In this study local ancestry was estimated across the genome in Polynesians (23 Tongan subjects) to find the candidate regions of admixture-enabled selection contributed by Papuan-related ancestors. RESULTS: The mean proportion of Papuan-related ancestry across the Polynesian genome was estimated as 24.6% (SD = 8.63%), and two genomic regions, the extended major histocompatibility complex (xMHC) region on chromosome 6 and the ATP-binding cassette transporter sub-family C member 11 (ABCC11) gene on chromosome 16, showed proportions of Papuan-related ancestry more than 5 SD greater than the mean (> 67.8%). The coalescent simulation under the assumption of selective neutrality suggested that such signals of Papuan-related ancestry enrichment were caused by positive selection after admixture (false discovery rate = 0.045). The ABCC11 harbors a nonsynonymous SNP, rs17822931, which affects apocrine secretory cell function. The approximate Bayesian computation indicated that, in Polynesian ancestors, a strong positive selection (s = 0.0217) acted on the ancestral allele of rs17822931 derived from Papuan-related ancestors. CONCLUSIONS: Our results suggest that admixture with Papuan-related ancestors contributed to the rapid local adaptation of Polynesian ancestors. Considering frequent admixture events in human evolution history, the acceleration of local adaptation through admixture should be a common event in humans.

Association of HLA-DRB1*09:01 with severe COVID-19.

HLA-A, -C, -B, and -DRB1 genotypes were analyzed in 178 Japanese COVID-19 patients to investigate the association of HLA with severe COVID-19. Analysis of 32 common HLA alleles at four loci revealed a significant association between HLA-DRB1*09:01 and severe COVID-19 (odds ratio [OR], 3.62; 95% CI, 1.57-8.35; p = 0.00251 [permutation p value = 0.0418]) when age, sex, and other common HLA alleles at the DRB1 locus were adjusted. The DRB1*09:01 allele was more significantly associated with risk for severe COVID-19 compared to preexisting medical conditions such as hypertension, diabetes, and cardiovascular diseases. These results indicate a potential role for HLA in predisposition to severe COVID-19.

Admixture and natural selection shaped genomes of an Austronesian-speaking population in the Solomon Islands.

People in the Solomon Islands today are considered to have derived from Asian- and Papuan-related ancestors. Papuan-related ancestors colonized Near Oceania about 47,000 years ago, and Asian-related ancestors were Austronesian (AN)-speaking population, called Lapita, who migrated from Southeast Asia about 3,500 years ago. These two ancestral populations admixed in Near Oceania before the expansion of Lapita people into Remote Oceania. To understand the impact of the admixture on the adaptation of AN-speaking Melanesians in Near Oceania, we performed the genome-wide single nucleotide polymorphism (SNP) analysis of 21 individuals from Munda, the main town of the New Georgia Islands in the western Solomon Islands. Population samples from Munda were genetically similar to other Solomon Island population samples. The analysis of genetic contribution from the two different ancestries to the Munda genome revealed significantly higher proportions of Asian- and Papuan-related ancestries in the region containing the annexin A1 (ANXA1) gene (Asian component > 82.6%) and in the human leukocyte antigen (HLA) class II region (Papuan component > 85.4%), respectively. These regions were suspected to have undergone natural selection since the time of admixture. Our results suggest that admixture had affected adaptation of AN-speaking Melanesians in the Solomon Islands.

A novel variant of torque teno virus 7 identified in patients with Kawasaki disease.

Kawasaki disease (KD), first identified in 1967, is a pediatric vasculitis of unknown etiology that has an increasing incidence in Japan and many other countries. KD can cause coronary artery aneurysms. Its epidemiological characteristics, such as seasonality and clinical picture of acute systemic inflammation with prodromal intestinal/respiratory symptoms, suggest an infectious etiology for KD. Interestingly, multiple host genotypes have been identified as predisposing factors for KD. To explore experimental methodology for identifying etiological agent(s) for KD and to optimize epidemiological study design (particularly the sample size) for future studies, we conducted a pilot study. For a 1-year period, we prospectively enrolled 11 patients with KD. To each KD patient, we assigned two control individuals (one with diarrhea and the other with respiratory infections), matched for age, sex, and season of diagnosis. During the acute phase of disease, we collected peripheral blood, nasopharyngeal aspirate, and feces. We also determined genotypes, to identify those that confer susceptibility to KD. There was no statistically significant difference in the frequency of the risk genotypes between KD patients and control subjects. We also used unbiased metagenomic sequencing to analyze these samples. Metagenomic sequencing and PCR detected torque teno virus 7 (TTV7) in two patients with KD (18%), but not in control subjects (P = 0.111). Sanger sequencing revealed that the TTV7 found in the two KD patients contained almost identical variants in nucleotide and identical changes in resulting amino acid, relative to the reference sequence. Additionally, we estimated the sample size that would be required to demonstrate a statistical correlation between TTV7 and KD. Future larger scale studies with carefully optimized metagenomic sequencing experiments and adequate sample size are warranted to further examine the association between KD and potential pathogens, including TTV7.