The need to diversify genomic studies: Insights from Andean highlanders and Amazonians.

More globally diverse perspectives are needed in genomic studies and precision medicine practices on non-Europeans. Here, we illustrate this by discussing the distribution of clinically actionable genetic variants involved in drug response in Andean highlanders and Amazonians, considering their environment, history, genetic structure, and historical biases in the perception of biological diversity of Native Americans.

Population genetics of PDE4B (phosphodiesterase-4B) in neglected Native Americans: Implications for cancer pharmacogenetics.

PDE4B (phosphodiesterase-4B) has an important role in cancer and in pharmacology of some disorders, such as inflammatory diseases. Remarkably in Native Americans, PDE4B variants are associated with acute lymphoblastic leukemia (ALL) relapse, as this gene modulates sensitivity of glucocorticoids used in ALL chemotherapy. PDE4B allele rs6683977.G, associated with genomic regions of Native American origin in US-Hispanics (admixed among Native Americans, Europeans, and Africans), increases ALL relapse risk, contributing to an association between Native American ancestry and ALL relapse that disappeared with an extra-phase of chemotherapy. This result insinuates that indigenous populations along the Americas may have high frequencies of rs6683977.G, but this has never been corroborated. We studied ancestry and PDE4B diversity in 951 healthy individuals from nine Latin American populations. In non-admixed Native American populations rs6683977.G has frequencies greater than 90%, is in linkage disequilibrium with other ALL relapse associated and regulatory variants in PDE4B-intron-7, conforming haplotypes showing their highest worldwide frequencies in Native Americans (>0.82). Our findings inform the discussion on the pertinence of an extra-phase of chemotherapy in Native American populations, and exemplifies how knowledge generated in US-Hispanics is relevant for their even more neglected and vulnerable Native American ancestors along the American continent.

Tracing the Distribution of European Lactase Persistence Genotypes Along the Americas.

In adulthood, the ability to digest lactose, the main sugar present in milk of mammals, is a phenotype (lactase persistence) observed in historically herder populations, mainly Northern Europeans, Eastern Africans, and Middle Eastern nomads. As the -13910∗T allele in the MCM6 gene is the most well-characterized allele responsible for the lactase persistence phenotype, the -13910C > T (rs4988235) polymorphism is commonly evaluated in lactase persistence studies. Lactase non-persistent adults may develop symptoms of lactose intolerance when consuming dairy products. In the Americas, there is no evidence of the consumption of these products until the arrival of Europeans. However, several American countries' dietary guidelines recommend consuming dairy for adequate human nutrition and health promotion. Considering the extensive use of dairy and the complex ancestry of Pan-American admixed populations, we studied the distribution of -13910C > T lactase persistence genotypes and its flanking haplotypes of European origin in 7,428 individuals from several Pan-American admixed populations. We found that the -13910∗T allele frequency in Pan-American admixed populations is directly correlated with allele frequency of the European sources. Moreover, we did not observe any overrepresentation of European haplotypes in the -13910C > T flanking region, suggesting no selective pressure after admixture in the Americas. Finally, considering the dominant effect of the -13910∗T allele, our results indicate that Pan-American admixed populations are likely to have higher frequency of lactose intolerance, suggesting that general dietary guidelines deserve further evaluation across the continent.

Human-SARS-CoV-2 interactome and human genetic diversity: TMPRSS2-rs2070788, associated with severe influenza, and its population genetics caveats in Native Americans.

For human/SARS-CoV-2 interactome genes ACE2, TMPRSS2 and BSG, there is a convincing evidence of association in Asians with influenza-induced SARS for TMPRSS2-rs2070788, tag-SNP of the eQTL rs383510. This case illustrates the importance of population genetics and of sequencing data in the design of genetic association studies in different human populations: the high linkage disequilibrium (LD) between rs2070788 and rs383510 is Asian-specific. Leveraging on a combination of genotyping and sequencing data for Native Americans (neglected in genetic studies), we show that while their frequencies of the Asian tag-SNP rs2070788 is, surprisingly, the highest worldwide, it is not in LD with the eQTL rs383510, that therefore, should be directly genotyped in genetic association studies of SARS in populations with Native American ancestry.

Characterizing the Genetic Architecture of Parkinson's Disease in Latinos.

OBJECTIVE: This work was undertaken in order to identify Parkinson's disease (PD) risk variants in a Latino cohort, to describe the overlap in the genetic architecture of PD in Latinos compared to European-ancestry subjects, and to increase the diversity in PD genome-wide association (GWAS) data. METHODS: We genotyped and imputed 1,497 PD cases and controls recruited from nine clinical sites across South America. We performed a GWAS using logistic mixed models; variants with a p-value <1 × 10-5 were tested in a replication cohort of 1,234 self-reported Latino PD cases and 439,522 Latino controls from 23andMe, Inc. We also performed an admixture mapping analysis where local ancestry blocks were tested for association with PD status. RESULTS: One locus, SNCA, achieved genome-wide significance (p-value <5 × 10-8 ); rs356182 achieved genome-wide significance in both the discovery and the replication cohorts (discovery, G allele: 1.58 OR, 95% CI 1.35-1.86, p-value 2.48 × 10-8 ; 23andMe, G allele: 1.26 OR, 95% CI 1.16-1.37, p-value 4.55 × 10-8 ). In our admixture mapping analysis, a locus on chromosome 14, containing the gene STXBP6, achieved significance in a joint test of ancestries and in the Native American single-ancestry test (p-value <5 × 10-5 ). A second locus on chromosome 6, containing the gene RPS6KA2, achieved significance in the African single-ancestry test (p-value <5 × 10-5 ). INTERPRETATION: This study demonstrated the importance of the SNCA locus for the etiology of PD in Latinos. By leveraging the demographic history of our cohort via admixture mapping, we identified two potential PD risk loci that merit further study. ANN NEUROL 2021;90:353-365.

Admixture/fine-mapping in Brazilians reveals a West African associated potential regulatory variant (rs114066381) with a strong female-specific effect on body mass and fat mass indexes.

BACKGROUND/OBJECTIVES: Admixed populations are a resource to study the global genetic architecture of complex phenotypes, which is critical, considering that non-European populations are severely underrepresented in genomic studies. Here, we study the genetic architecture of BMI in children, young adults, and elderly individuals from the admixed population of Brazil. SUBJECTS/METHODS: Leveraging admixture in Brazilians, whose chromosomes are mosaics of fragments of Native American, European, and African origins, we used genome-wide data to perform admixture mapping/fine-mapping of body mass index (BMI) in three Brazilian population-based cohorts from Northeast (Salvador), Southeast (Bambuí), and South (Pelotas). RESULTS: We found significant associations with African-associated alleles in children from Salvador (PALD1 and ZMIZ1 genes), and in young adults from Pelotas (NOD2 and MTUS2 genes). More importantly, in Pelotas, rs114066381, mapped in a potential regulatory region, is significantly associated only in females (p = 2.76e-06). This variant is rare in Europeans but with frequencies of ~3% in West Africa and has a strong female-specific effect (95% CI: 2.32-5.65 kg/m2 per each A allele). We confirmed this sex-specific association and replicated its strong effect for an adjusted fat mass index in the same Pelotas cohort, and for BMI in another Brazilian cohort from São Paulo (Southeast Brazil). A meta-analysis confirmed the significant association. Remarkably, we observed that while the frequency of rs114066381-A allele ranges from 0.8 to 2.1% in the studied populations, it attains ~9% among women with morbid obesity from Pelotas, São Paulo, and Bambuí. The effect size of rs114066381 is at least five times higher than the FTO SNPs rs9939609 and rs1558902, already emblematic for their high effects. CONCLUSIONS: We identified six candidate SNPs associated with BMI. rs114066381 stands out for its high effect that was replicated and its high frequency in women with morbid obesity. We demonstrate how admixed populations are a source of new relevant phenotype-associated genetic variants.

The genetic structure and adaptation of Andean highlanders and Amazonians are influenced by the interplay between geography and culture.

Western South America was one of the worldwide cradles of civilization. The well-known Inca Empire was the tip of the iceberg of an evolutionary process that started 11,000 to 14,000 years ago. Genetic data from 18 Peruvian populations reveal the following: 1) The between-population homogenization of the central southern Andes and its differentiation with respect to Amazonian populations of similar latitudes do not extend northward. Instead, longitudinal gene flow between the northern coast of Peru, Andes, and Amazonia accompanied cultural and socioeconomic interactions revealed by archeology. This pattern recapitulates the environmental and cultural differentiation between the fertile north, where altitudes are lower, and the arid south, where the Andes are higher, acting as a genetic barrier between the sharply different environments of the Andes and Amazonia. 2) The genetic homogenization between the populations of the arid Andes is not only due to migrations during the Inca Empire or the subsequent colonial period. It started at least during the earlier expansion of the Wari Empire (600 to 1,000 years before present). 3) This demographic history allowed for cases of positive natural selection in the high and arid Andes vs. the low Amazon tropical forest: in the Andes, a putative enhancer in HAND2-AS1 (heart and neural crest derivatives expressed 2 antisense RNA1, a noncoding gene related to cardiovascular function) and rs269868-C/Ser1067 in DUOX2 (dual oxidase 2, related to thyroid function and innate immunity) genes and, in the Amazon, the gene encoding for the CD45 protein, essential for antigen recognition by T and B lymphocytes in viral-host interaction.

A positively selected FBN1 missense variant reduces height in Peruvian individuals.

On average, Peruvian individuals are among the shortest in the world1. Here we show that Native American ancestry is associated with reduced height in an ethnically diverse group of Peruvian individuals, and identify a population-specific, missense variant in the FBN1 gene (E1297G) that is significantly associated with lower height. Each copy of the minor allele (frequency of 4.7%) reduces height by 2.2 cm (4.4 cm in homozygous individuals). To our knowledge, this is the largest effect size known for a common height-associated variant. FBN1 encodes the extracellular matrix protein fibrillin 1, which is a major structural component of microfibrils. We observed less densely packed fibrillin-1-rich microfibrils with irregular edges in the skin of individuals who were homozygous for G1297 compared with individuals who were homozygous for E1297. Moreover, we show that the E1297G locus is under positive selection in non-African populations, and that the E1297 variant shows subtle evidence of positive selection specifically within the Peruvian population. This variant is also significantly more frequent in coastal Peruvian populations than in populations from the Andes or the Amazon, which suggests that short stature might be the result of adaptation to factors that are associated with the coastal environment in Peru.

Lung Ultrasound Volume Sweep Imaging for Pneumonia Detection in Rural Areas: Piloting Training in Rural Peru.

OBJECTIVE: Pneumonia is the leading cause of pediatric mortality worldwide among children 0-5 years old. Lung ultrasound can be used to diagnose pneumonia in rural areas as it is a portable and relatively economic imaging modality with ~95% sensitivity and specificity for pneumonia in children. Lack of trained sonographers is the current limiting factor to its deployment in rural areas. In this study, we piloted training of a volume sweep imaging (VSI) ultrasound protocol for pneumonia detection in Peru with rural health workers. VSI may be taught to individuals with limited medical/ultrasound experience as it requires minimal anatomical knowledge and technical skill. In VSI, the target organ is imaged with a series of sweeps and arcs of the ultrasound probe in relation to external body landmarks. METHODS: Rural health workers in Peru were trained on a VSI ultrasound protocol for pneumonia detection. Subjects were given a brief didactic session followed by hands-on practice with the protocol. Each attempt was timed and mistakes were recorded. Participants performed the protocol until they demonstrated two mistake-free attempts. RESULTS: It took participants a median number of three attempts (range 1-6) to perform the VSI protocol correctly. Time to mastery took 51.4 ± 17.7 min. There were no significant differences among doctors, nurses, and technicians in total training time (P = 0.43) or number of attempts to success (P = 0.72). Trainee age was not found to be significantly correlated with training time (P = 0.50) or number of attempts to success (P = 0.40). CONCLUSION: Rural health workers learned a VSI protocol for pneumonia detection with relative ease in a short amount of time. Future studies should investigate the clinical efficacy of this VSI protocol for pneumonia detection. KEY MESSAGE: A volume sweep imaging (VSI) protocol for pneumonia detection can be taught with minimal difficulty to rural health workers without prior ultrasound experience. No difference was found in training performance related to education level or age. VSI involves no significant knowledge of anatomy or technical skill.

Native American Genomic Diversity through Ancient DNA.

Ancient DNA is a powerful tool to understand the evolutionary dynamics of both current and ancestral populations. Posth et al. use ancient DNA to elucidate important questions surrounding the peopling of Central and South America, giving us greater insights into the ancestry of genetically understudied populations.

Evolutionary genomic dynamics of Peruvians before, during, and after the Inca Empire.

Native Americans from the Amazon, Andes, and coastal geographic regions of South America have a rich cultural heritage but are genetically understudied, therefore leading to gaps in our knowledge of their genomic architecture and demographic history. In this study, we sequence 150 genomes to high coverage combined with an additional 130 genotype array samples from Native American and mestizo populations in Peru. The majority of our samples possess greater than 90% Native American ancestry, which makes this the most extensive Native American sequencing project to date. Demographic modeling reveals that the peopling of Peru began ∼12,000 y ago, consistent with the hypothesis of the rapid peopling of the Americas and Peruvian archeological data. We find that the Native American populations possess distinct ancestral divisions, whereas the mestizo groups were admixtures of multiple Native American communities that occurred before and during the Inca Empire and Spanish rule. In addition, the mestizo communities also show Spanish introgression largely following Peruvian Independence, nearly 300 y after Spain conquered Peru. Further, we estimate migration events between Peruvian populations from all three geographic regions with the majority of between-region migration moving from the high Andes to the low-altitude Amazon and coast. As such, we present a detailed model of the evolutionary dynamics which impacted the genomes of modern-day Peruvians and a Native American ancestry dataset that will serve as a beneficial resource to addressing the underrepresentation of Native American ancestry in sequencing studies.

Relationships between phyllosphere bacterial communities and plant functional traits in a neotropical forest.

The phyllosphere--the aerial surfaces of plants, including leaves--is a ubiquitous global habitat that harbors diverse bacterial communities. Phyllosphere bacterial communities have the potential to influence plant biogeography and ecosystem function through their influence on the fitness and function of their hosts, but the host attributes that drive community assembly in the phyllosphere are poorly understood. In this study we used high-throughput sequencing to quantify bacterial community structure on the leaves of 57 tree species in a neotropical forest in Panama. We tested for relationships between bacterial communities on tree leaves and the functional traits, taxonomy, and phylogeny of their plant hosts. Bacterial communities on tropical tree leaves were diverse; leaves from individual trees were host to more than 400 bacterial taxa. Bacterial communities in the phyllosphere were dominated by a core microbiome of taxa including Actinobacteria, Alpha-, Beta-, and Gammaproteobacteria, and Sphingobacteria. Host attributes including plant taxonomic identity, phylogeny, growth and mortality rates, wood density, leaf mass per area, and leaf nitrogen and phosphorous concentrations were correlated with bacterial community structure on leaves. The relative abundances of several bacterial taxa were correlated with suites of host plant traits related to major axes of plant trait variation, including the leaf economics spectrum and the wood density-growth/mortality tradeoff. These correlations between phyllosphere bacterial diversity and host growth, mortality, and function suggest that incorporating information on plant-microbe associations will improve our ability to understand plant functional biogeography and the drivers of variation in plant and ecosystem function.

Serum alpha-fetoprotein screening for hepatoblastoma in children with Beckwith-Wiedemann syndrome or isolated hemihyperplasia.

An elevated risk of hepatoblastoma for children with Beckwith-Wiedemann syndrome or isolated hemihyperplasia is well established. We describe five children with Beckwith-Wiedemann syndrome or isolated hemihyperplasia for whom serial serum alpha-fetoprotein screening, usually in combination with abdominal ultrasound, led to early detection of hepatoblastoma.