Scientists without borders: lessons from Ukraine

Conflicts and natural disasters affect entire populations of the countries involved and, in addition to the thousands of lives destroyed, have a substantial negative impact on the scientific advances these countries provide. The unprovoked invasion of Ukraine by Russia, the devastating earthquake in Turkey and Syria, and the ongoing conflicts in the Middle East are just a few examples. Millions of people have been killed or displaced, their futures uncertain. These events have resulted in extensive infrastructure collapse, with loss of electricity, transportation, and access to services. Schools, universities, and research centers have been destroyed along with decades' worth of data, samples, and findings. Scholars in disaster areas face short- and long-term problems in terms of what they can accomplish now for obtaining grants and for employment in the long run. In our interconnected world, conflicts and disasters are no longer a local problem but have wide-ranging impacts on the entire world, both now and in the future. Here, we focus on the current and ongoing impact of war on the scientific community within Ukraine and from this draw lessons that can be applied to all affected countries where scientists at risk are facing hardship. We present and classify examples of effective and feasible mechanisms used to support researchers in countries facing hardship and discuss how these can be implemented with help from the international scientific community and what more is desperately needed. Reaching out, providing accessible training opportunities, and developing collaborations should increase inclusion and connectivity, support scientific advancements within affected communities, and expedite postwar and disaster recovery.

The Pioneer Advantage: Filling the blank spots on the map of genome diversity in Europe

Documenting genome diversity is important for the local biomedical communities and instrumental in developing precision and personalized medicine. Currently, tens of thousands of whole-genome sequences from Europe are publicly available, but most of these represent populations of developed countries of Europe. The uneven distribution of the available data is further impaired by the lack of data sharing. Recent whole-genome studies in Eastern Europe, one in Ukraine and one in Russia, demonstrated that local genome diversity and population structure from Eastern Europe historically had not been fully represented. An unexpected wealth of genomic variation uncovered in these studies was not so much a consequence of high variation within their population, but rather due to the "pioneer advantage." We discovered more variants because we were the first to prospect in the Eastern European genome pool. This simple comparison underscores the importance of removing the remaining geographic genome deserts from the rest of the world map of the human genome diversity.

Genetic diversity and selection in Puerto Rican horses

Since the first Spanish settlers brought horses to America centuries ago, several local varieties and breeds have been established in the New World. These were generally a consequence of the admixture of the different breeds arriving from Europe. In some instances, local horses have been selectively bred for specific traits, such as appearance, endurance, strength, and gait. We looked at the genetics of two breeds, the Puerto Rican Non-Purebred (PRNPB) (also known as the "Criollo") horses and the Puerto Rican Paso Fino (PRPF), from the Caribbean Island of Puerto Rico. While it is reasonable to assume that there was a historic connection between the two, the genetic link between them has never been established. In our study, we started by looking at the genetic ancestry and diversity of current Puerto Rican horse populations using a 668 bp fragment of the mitochondrial DNA D-loop (HVR1) in 200 horses from 27 locations on the island. We then genotyped all 200 horses in our sample for the "gait-keeper" DMRT3 mutant allele previously associated with the paso gait especially cherished in this island breed. We also genotyped a subset of 24 samples with the Illumina Neogen Equine Community genome-wide array (65,000 SNPs). This data was further combined with the publicly available PRPF genomes from other studies. Our analysis show an undeniable genetic connection between the two varieties in Puerto Rico, consistent with the hypothesis that PRNPB horses represent the descendants of the original genetic pool, a mix of horses imported from the Iberian Peninsula and elsewhere in Europe. Some of the original founders of PRNRB population must have carried the "gait-keeper" DMRT3 allele upon arrival to the island. From this admixture, the desired traits were selected by the local people over the span of centuries. We propose that the frequency of the mutant "gait-keeper" allele originally increased in the local horses due to the selection for the smooth ride and other characters, long before the PRPF breed was established. To support this hypothesis, we demonstrate that PRNPB horses, and not the purebred PRPF, carry a signature of selection in the genomic region containing the DMRT3 locus to this day. The lack of the detectable signature of selection associated with the DMRT3 in the PRPF would be expected if this native breed was originally derived from the genetic pool of PRNPB horses established earlier and most of the founders already had the mutant allele. Consequently, selection specific to PRPF later focused on allels in other genes (including CHRM5, CYP2E1, MYH7, SRSF1, PAM, PRN and others) that have not been previously associated with the prized paso gait phenotype in Puerto Rico or anywhere else.

Molecular Phylogeny and Evolution of Amazon Parrots in the Greater Antilles

Amazon parrots (Amazona spp.) colonized the islands of the Greater Antilles from the Central American mainland, but there has not been a consensus as to how and when this happened. Today, most of the five remaining island species are listed as endangered, threatened, or vulnerable as a consequence of human activity. We sequenced and annotated full mitochondrial genomes of all the extant Amazon parrot species from the Greater Antillean (A. leucocephala (Cuba), A. agilis, A. collaria (both from Jamaica), A. ventralis (Hispaniola), and A. vittata (Puerto Rico)), A. albifrons from mainland Central America, and A. rhodocorytha from the Atlantic Forest in Brazil. The assembled and annotated mitogenome maps provide information on sequence organization, variation, population diversity, and evolutionary history for the Caribbean species including the critically endangered A. vittata. Despite the larger number of available samples from the Puerto Rican Parrot Recovery Program, the sequence diversity of the A. vittata population in Puerto Rico was the lowest among all parrot species analyzed. Our data support the stepping-stone dispersal and speciation hypothesis that has started approximately 3.47 MYA when the ancestral population arrived from mainland Central America and led to diversification across the Greater Antilles, ultimately reaching the island of Puerto Rico 0.67 MYA. The results are presented and discussed in light of the geological history of the Caribbean and in the context of recent parrot evolution, island biogeography, and conservation. This analysis contributes to understating evolutionary history and empowers subsequent assessments of sequence variation and helps design future conservation efforts in the Caribbean.

Genome diversity in Ukraine

Background

The main goal of this collaborative effort is to provide genome-wide data for the previously underrepresented population in Eastern Europe, and to provide cross-validation of the data from genome sequences and genotypes of the same individuals acquired by different technologies. We collected 97 genome-grade DNA samples from consented individuals representing major regions of Ukraine that were consented for public data release. BGISEQ-500 sequence data and genotypes by an Illumina GWAS chip were cross-validated on multiple samples and additionally referenced to 1 sample that has been resequenced by Illumina NovaSeq6000 S4 at high coverage.

Results

The genome data have been searched for genomic variation represented in this population, and a number of variants have been reported: large structural variants, indels, copy number variations, single-nucletide polymorphisms, and microsatellites. To our knowledge, this study provides the largest to-date survey of genetic variation in Ukraine, creating a public reference resource aiming to provide data for medical research in a large understudied population.

Conclusions

Our results indicate that the genetic diversity of the Ukrainian population is uniquely shaped by evolutionary and demographic forces and cannot be ignored in future genetic and biomedical studies. These data will contribute a wealth of new information bringing forth a wealth of novel, endemic and medically related alleles.

Genome-wide sequence analyses of ethnic populations across Russia

The Russian Federation is the largest and one of the most ethnically diverse countries in the world, however no centralized reference database of genetic variation exists to date. Such data are crucial for medical genetics and essential for studying population history. The Genome Russia Project aims at filling this gap by performing whole genome sequencing and analysis of peoples of the Russian Federation. Here we report the characterization of genome-wide variation of 264 healthy adults, including 60 newly sequenced samples. People of Russia carry known and novel genetic variants of adaptive, clinical and functional consequence that in many cases show allele frequency divergence from neighboring populations. Population genetics analyses revealed six phylogeographic partitions among indigenous ethnicities corresponding to their geographic locales. This study presents a characterization of population-specific genomic variation in Russia with results important for medical genetics and for understanding the dynamic population history of the world's largest country.

Genomes of Three Closely Related Caribbean Amazons Provide Insight for Species History and Conservation

Islands have been used as model systems for studies of speciation and extinction since Darwin published his observations about finches found on the Galapagos. Amazon parrots inhabiting the Greater Antillean Islands represent a fascinating model of species diversification. Unfortunately, many of these birds are threatened as a result of human activity and some, like the Puerto Rican parrot, are now critically endangered. In this study we used a combination of de novo and reference-assisted assembly methods, integrating it with information obtained from related genomes to perform genome reconstruction of three amazon species. First, we used whole genome sequencing data to generate a new de novo genome assembly for the Puerto Rican parrot (Amazona vittata). We then improved the obtained assembly using transcriptome data from Amazona ventralis and used the resulting sequences as a reference to assemble the genomes Hispaniolan (A. ventralis) and Cuban (Amazona leucocephala) parrots. Finally, we, annotated genes and repetitive elements, estimated genome sizes and current levels of heterozygosity, built models of demographic history and provided interpretation of our findings in the context of parrot evolution in the Caribbean.

A Recurrent BRCA2 Mutation Explains the Majority of Hereditary Breast and Ovarian Cancer Syndrome Cases in Puerto Rico

Breast cancer is the most common cause of cancer diagnosis in women and is responsible for considerable mortality among the women of Puerto Rico. However, there are few studies in Puerto Rico on the genetic factors influencing risk. To determine the contribution of pathogenic mutations in BRCA1 and BRCA2, we sequenced these genes in 302 cases from two separate medical centers, who were not selected for age of onset or family history. We identified nine cases that are carriers of pathogenic germline mutation. This represents 2.9% of unselected cases and 5.6% of women meeting National Comprehensive Cancer Network (NCCN) criteria for BRCA testing. All of the identified pathogenic mutations were in the BRCA2 gene and the most common mutation is the p.Glu1308Ter (E1308X) mutation in BRCA2 found in eight out of nine cases, representing 89% of the pathogenic carriers. The E1308X mutation has been identified in breast and ovarian cancer families in Spain, and analysis of flanking DNA polymorphisms shows that all E1308X carriers occur on the same haplotype. This is consistent with BRCA2 E1308X being a founder mutation for the Puerto Rican population. These results will contribute to better inform genetic screening and counseling of breast and ovarian cancer cases in Puerto Rico and Puerto Rican populations in mainland United States.

Innovative assembly strategy contributes to understanding the evolution and conservation genetics of the endangered Solenodon paradoxus from the island of Hispaniola

Solenodons are insectivores that live in Hispaniola and Cuba. They form an isolated branch in the tree of placental mammals that are highly divergent from other eulipothyplan insectivores The history, unique biology, and adaptations of these enigmatic venomous species could be illuminated by the availability of genome data. However, a whole genome assembly for solenodons has not been previously performed, partially due to the difficulty in obtaining samples from the field. Island isolation and reduced numbers have likely resulted in high homozygosity within the Hispaniolan solenodon (Solenodon paradoxus). Thus, we tested the performance of several assembly strategies on the genome of this genetically impoverished species. The string graph-based assembly strategy seemed a better choice compared to the conventional de Bruijn graph approach due to the high levels of homozygosity, which is often a hallmark of endemic or endangered species. A consensus reference genome was assembled from sequences of 5 individuals from the southern subspecies (S. p. woodi). In addition, we obtained an additional sequence from 1 sample of the northern subspecies (S. p. paradoxus). The resulting genome assemblies were compared to each other and annotated for genes, with an emphasis on venom genes, repeats, variable microsatellite loci, and other genomic variants. Phylogenetic positioning and selection signatures were inferred based on 4,416 single-copy orthologs from 10 other mammals. We estimated that solenodons diverged from other extant mammals 73.6 million years ago. Patterns of single-nucleotide polymorphism variation allowed us to infer population demography, which supported a subspecies split within the Hispaniolan solenodon at least 300 thousand years ago.

Mitogenomic sequences support a north-south subspecies subdivision within Solenodon paradoxus

Solenodons are insectivores found only in Hispaniola and Cuba, with a Mesozoic divergence date versus extant mainland mammals. Solenodons are the oldest lineage of living eutherian mammal for which a mitogenome sequence has not been reported. We determined complete mitogenome sequences for six Hispaniolan solenodons (Solenodon paradoxus) using next-generation sequencing. The solenodon mitogenomes were 16,454-16,457 bp long and carried the expected repertoire of genes. A mitogenomic phylogeny confirmed the basal position of solenodons relative to shrews and moles, with solenodon mitogenomes estimated to have diverged from those of other mammals ca. 78 Mya. Control region sequences of solenodons from the northern (n = 3) and southern (n = 5) Dominican Republic grouped separately in a network, with F_ST = 0.72 (p = 0.036) between north and south. This regional genetic divergence supports previous morphological and genetic reports recognizing northern (S. p. paradoxus) and southern (S. p. woodi) subspecies in need of separate conservation plans.

Genomic legacy of the African cheetah, Acinonyx jubatus

BACKGROUND

Patterns of genetic and genomic variance are informative in inferring population history for human, model species and endangered populations.

RESULTS

Here the genome sequence of wild-born African cheetahs reveals extreme genomic depletion in SNV incidence, SNV density, SNVs of coding genes, MHC class I and II genes, and mitochondrial DNA SNVs. Cheetah genomes are on average 95 % homozygous compared to the genomes of the outbred domestic cat (24.08 % homozygous), Virunga Mountain Gorilla (78.12 %), inbred Abyssinian cat (62.63 %), Tasmanian devil, domestic dog and other mammalian species. Demographic estimators impute two ancestral population bottlenecks: one >100,000 years ago coincident with cheetah migrations out of the Americas and into Eurasia and Africa, and a second 11,084-12,589 years ago in Africa coincident with late Pleistocene large mammal extinctions. MHC class I gene loss and dramatic reduction in functional diversity of MHC genes would explain why cheetahs ablate skin graft rejection among unrelated individuals. Significant excess of non-synonymous mutations in AKAP4 (p<0.02), a gene mediating spermatozoon development, indicates cheetah fixation of five function-damaging amino acid variants distinct from AKAP4 homologues of other Felidae or mammals; AKAP4 dysfunction may cause the cheetah's extremely high (>80 %) pleiomorphic sperm.

CONCLUSIONS

The study provides an unprecedented genomic perspective for the rare cheetah, with potential relevance to the species' natural history, physiological adaptations and unique reproductive disposition.

Sequencing rare and common APOL1 coding variants to determine kidney disease risk

A third of African Americans with sporadic focal segmental glomerulosclerosis (FSGS) or HIV-associated nephropathy (HIVAN) do not carry APOL1 renal risk genotypes. This raises the possibility that other APOL1 variants may contribute to kidney disease. To address this question, we sequenced all APOL1 exons in 1, 437 Americans of African and European decent, including 464 patients with biopsy-proven FSGS/HIVAN. Testing for association with 33 common and rare variants with FSGS/HIVAN revealed no association independent of strong recessive G1 and G2 effects. Seeking additional variants that might have been under selection by pathogens and could represent candidates for kidney disease risk, we also sequenced an additional 1, 112 individuals representing 53 global populations. Except for G1 and G2, none of the 7 common codon-altering variants showed evidence of selection or could restore lysis against trypanosomes causing human African trypanosomiasis. Thus, only APOL1 G1 and G2 confer renal risk and other common and rare APOL1 missense variants, including the archaic G3 haplotype, do not contribute to sporadic FSGS and HIVAN in the United States population. Hence, in most potential clinical or screening applications, our study suggests that sequencing APOL1 exons is unlikely to bring additional information compared to genotyping only APOL1 G1 and G2 risk alleles.

SmileFinder: a resampling-based approach to evaluate signatures of selection from genome-wide sets of matching allele frequency data in two or more diploid populations

BACKGROUND

Adaptive alleles may rise in frequency as a consequence of positive selection, creating a pattern of decreased variation in the neighboring loci, known as a selective sweep. When the region containing this pattern is compared to another population with no history of selection, a rise in variance of allele frequencies between populations is observed. One challenge presented by large genome-wide datasets is the ability to differentiate between patterns that are remnants of natural selection from those expected to arise at random and/or as a consequence of selectively neutral demographic forces acting in the population.

FINDINGS

SmileFinder is a simple program that looks for diversity and divergence patterns consistent with selection sweeps by evaluating allele frequencies in windows, including neighboring loci from two or more populations of a diploid species against the genome-wide neutral expectation. The program calculates the mean of heterozygosity and FST in a set of sliding windows of incrementally increasing sizes, and then builds a resampled distribution (the baseline) of random multi-locus sets matched to the sizes of sliding windows, using an unrestricted sampling. Percentiles of the values in the sliding windows are derived from the superimposed resampled distribution. The resampling can easily be scaled from 1 K to 100 M; the higher the number, the more precise the percentiles ascribed to the extreme observed values.

CONCLUSIONS

The output from SmileFinder can be used to plot percentile values to look for population diversity and divergence patterns that may suggest past actions of positive selection along chromosome maps, and to compare lists of suspected candidate genes under random gene sets to test for the overrepresentation of these patterns among gene categories. Both applications of the algorithm have already been used in published studies. Here we present a publicly available, open source program that will serve as a useful tool for preliminary scans of selection using worldwide databases of human genetic variation, as well as population datasets for many non-human species, from which such data is rapidly emerging with the advent of new genotyping and sequencing technologies.

Reconstructing Native American migrations from whole-genome and whole-exome data

There is great scientific and popular interest in understanding the genetic history of populations in the Americas. We wish to understand when different regions of the continent were inhabited, where settlers came from, and how current inhabitants relate genetically to earlier populations. Recent studies unraveled parts of the genetic history of the continent using genotyping arrays and uniparental markers. The 1000 Genomes Project provides a unique opportunity for improving our understanding of population genetic history by providing over a hundred sequenced low coverage genomes and exomes from Colombian (CLM), Mexican-American (MXL), and Puerto Rican (PUR) populations. Here, we explore the genomic contributions of African, European, and especially Native American ancestry to these populations. Estimated Native American ancestry is in MXL, in CLM, and in PUR. Native American ancestry in PUR is most closely related to populations surrounding the Orinoco River basin, confirming the Southern America ancestry of the Taíno people of the Caribbean. We present new methods to estimate the allele frequencies in the Native American fraction of the populations, and model their distribution using a demographic model for three ancestral Native American populations. These ancestral populations likely split in close succession: the most likely scenario, based on a peopling of the Americas thousand years ago (kya), supports that the MXL Ancestors split kya, with a subsequent split of the ancestors to CLM and PUR kya. The model also features effective populations of in Mexico, in Colombia, and in Puerto Rico. Modeling Identity-by-descent (IBD) and ancestry tract length, we show that post-contact populations also differ markedly in their effective sizes and migration patterns, with Puerto Rico showing the smallest effective size and the earlier migration from Europe. Finally, we compare IBD and ancestry assignments to find evidence for relatedness among European founders to the three populations.

Populations of the Americas have a rich and heterogeneous genetic and cultural heritage that draws from a diversity of pre-Columbian Native American, European, and African populations. Characterizing this diversity facilitates the development of medical genetics research in diverse populations and the transfer of medical knowledge across populations. It also represents an opportunity to better understand the peopling of the Americas, from the crossing of Beringia to the post-Columbian era. Here, we take advantage sequencing of individuals of Colombian (CLM), Mexican (MXL), and Puerto Rican (PUR) origin by the 1000 Genomes project to improve our demographic models for the peopling of the Americas. The divergence among African, European, and Native American ancestors to these populations enables us to infer the continent of origin at each locus in the sampled genomes. The resulting patterns of ancestry suggest complex post-Columbian migration histories, starting later in CLM than in MXL and PUR. Whereas European ancestral segments show evidence of relatedness, a demographic model of synonymous variation suggests that the Native American Ancestors to MXL, PUR, and CLM panels split within a few hundred years over 12 thousand years ago. Together with early archeological sites in South America, these results support rapid divergence during the initial peopling of the Americas.

Colonization of islands in the Mona Passage by endemic dwarf geckoes (genus Sphaerodactylus) reconstructed with mitochondrial phylogeny

Little is known about the natural history of the Sphaerodactylus species endemic to the three islands located in the Mona Passage separating the Greater Antillean islands of Hispaniola and Puerto Rico. In this study, parts of two mitochondrial genes, 16S rRNA and 12S rRNA, were sequenced to determine the relationships between the sphaerodactylids that live in the Mona Passage and other Caribbean species from the same genus. While the main goal was to identify the biogeographical origin of these species, we also identified a genetically distinct type of dwarf gecko that warrants future evaluation as a possible new species. According to the reconstructed phylogenies, we propose a stepwise model of colonization wherein S. nicholsi from southwestern Puerto Rico or a very close ancestor gave rise through a founder event to Sphaerodactylus monensis on Mona Island. In a similar fashion, S. monensis or a very close ancestor on Mona Island gave rise to S. levinsi on Desecheo Island. This study also suggests that the most recent common ancestor between the species from the islands in the Mona Passage and Puerto Rico existed approximately 3 MYA.

The Sp185/333 immune response genes and proteins are expressed in cells dispersed within all major organs of the adult purple sea urchin

Purple sea urchins (Strongylocentrotus purpuratus) express a highly variable set of immune genes called Sp185/333 by two subtypes of coelomocytes: the polygonal and small phagocytes. We report that the Sp185/333 genes and their encoded proteins are also expressed in all of the major organs in the adult sea urchin, including the axial organ, pharynx, esophagus, intestine and gonads. After immune challenge, there is an increase in the level of Sp185/333 mRNA in cells associated with the intestine and axial organ. The Sp185/333 proteins increase in the axial organ, pharynx, esophagus and intestine after challenge. However, the proportion of Sp185/333-positive cells only increases in the axial organ, while there is no change in that proportion in the other organs after challenge. The size range of the major Sp185/333 proteins expressed by organs is broader (5 kDa to > 250 kDa) compared with those in coelomocytes (∼40 kDa to < 250 kDa). Images of the different organs do not clarify whether coelomocytes or parenchymal cells express the Sp185/333 proteins. The increase in levels of Sp185/333 transcripts, protein expression and Sp185/333-positive cells in the axial organ in response to challenge suggests that this organ may have an important role in immunity for this species.

Evidence for selection at HIV host susceptibility genes in a West Central African human population

Background

HIV-1 derives from multiple independent transfers of simian immunodeficiency virus (SIV) strains from chimpanzees to human populations. We hypothesized that human populations in west central Africa may have been exposed to SIV prior to the pandemic, and that previous outbreaks may have selected for genetic resistance to immunodeficiency viruses. To test this hypothesis, we examined the genomes of Biaka Western Pygmies, who historically resided in communities within the geographic range of the central African chimpanzee subspecies (Pan troglodytes troglodytes) that carries strains of SIV ancestral to HIV-1.

Results

SNP genotypes of the Biaka were compared to those of African human populations who historically resided outside the range of P. t. troglodytes, including the Mbuti Eastern Pygmies. Genomic regions showing signatures of selection were compared to the genomic locations of genes reported to be associated with HIV infection or pathogenesis. In the Biaka, a strong signal of selection was detected at CUL5, which codes for a component of the vif-mediated APOBEC3 degradation pathway. A CUL5 allele protective against AIDS progression was fixed in the Biaka. A signal of selection was detected at TRIM5, which codes for an HIV post-entry restriction factor. A protective mis-sense mutation in TRIM5 had the highest frequency in Biaka compared to other African populations, as did a protective allele for APOBEC3G, which codes for an anti-HIV-1 restriction factor. Alleles protective against HIV-1 for APOBEC3H, CXCR6 and HLA-C were at higher frequencies in the Biaka than in the Mbuti. Biaka genomes showed a strong signal of selection at TSG101, an inhibitor of HIV-1 viral budding.

Conclusions

We found protective alleles or evidence for selection in the Biaka at a number of genes associated with HIV-1 infection or progression. Pygmies have also been reported to carry genotypes protective against HIV-1 for the genes CCR5 and CCL3L1. Our hypothesis that HIV-1 may have shaped the genomes of some human populations in West Central Africa appears to merit further investigation.

A locally funded Puerto Rican parrot (Amazona vittata) genome sequencing project increases avian data and advances young researcher education

BACKGROUND

Amazona vittata is a critically endangered Puerto Rican endemic bird, the only surviving native parrot species in the United States territory, and the first parrot in the large Neotropical genus Amazona, to be studied on a genomic scale.

FINDINGS

In a unique community-based funded project, DNA from an A. vittata female was sequenced using a HiSeq Illumina platform, resulting in a total of ~42.5 billion nucleotide bases. This provided approximately 26.89x average coverage depth at the completion of this funding phase. Filtering followed by assembly resulted in 259,423 contigs (N50 = 6,983 bp, longest = 75,003 bp), which was further scaffolded into 148,255 fragments (N50 = 19,470, longest = 206,462 bp). This provided ~76% coverage of the genome based on an estimated size of 1.58 Gb. The assembled scaffolds allowed basic genomic annotation and comparative analyses with other available avian whole-genome sequences.

CONCLUSIONS

The current data represents the first genomic information from and work carried out with a unique source of funding. This analysis further provides a means for directed training of young researchers in genetic and bioinformatics analyses and will facilitate progress towards a full assembly and annotation of the Puerto Rican parrot genome. It also adds extensive genomic data to a new branch of the avian tree, making it useful for comparative analyses with other avian species. Ultimately, the knowledge acquired from these data will contribute to an improved understanding of the overall population health of this species and aid in ongoing and future conservation efforts.

Reconciling apparent conflicts between mitochondrial and nuclear phylogenies in African elephants

Conservation strategies for African elephants would be advanced by resolution of conflicting claims that they comprise one, two, three or four taxonomic groups, and by development of genetic markers that establish more incisively the provenance of confiscated ivory. We addressed these related issues by genotyping 555 elephants from across Africa with microsatellite markers, developing a method to identify those loci most effective at geographic assignment of elephants (or their ivory), and conducting novel analyses of continent-wide datasets of mitochondrial DNA. Results showed that nuclear genetic diversity was partitioned into two clusters, corresponding to African forest elephants (99.5% Cluster-1) and African savanna elephants (99.4% Cluster-2). Hybrid individuals were rare. In a comparison of basal forest "F" and savanna "S" mtDNA clade distributions to nuclear DNA partitions, forest elephant nuclear genotypes occurred only in populations in which S clade mtDNA was absent, suggesting that nuclear partitioning corresponds to the presence or absence of S clade mtDNA. We reanalyzed African elephant mtDNA sequences from 81 locales spanning the continent and discovered that S clade mtDNA was completely absent among elephants at all 30 sampled tropical forest locales. The distribution of savanna nuclear DNA and S clade mtDNA corresponded closely to range boundaries traditionally ascribed to the savanna elephant species based on habitat and morphology. Further, a reanalysis of nuclear genetic assignment results suggested that West African elephants do not comprise a distinct third species. Finally, we show that some DNA markers will be more useful than others for determining the geographic origins of illegal ivory. These findings resolve the apparent incongruence between mtDNA and nuclear genetic patterns that has confounded the taxonomy of African elephants, affirm the limitations of using mtDNA patterns to infer elephant systematics or population structure, and strongly support the existence of two elephant species in Africa.

History shaped the geographic distribution of genomic admixture on the island of Puerto Rico

Contemporary genetic variation among Latin Americans human groups reflects population migrations shaped by complex historical, social and economic factors. Consequently, admixture patterns may vary by geographic regions ranging from countries to neighborhoods. We examined the geographic variation of admixture across the island of Puerto Rico and the degree to which it could be explained by historic and social events. We analyzed a census-based sample of 642 Puerto Rican individuals that were genotyped for 93 ancestry informative markers (AIMs) to estimate African, European and Native American ancestry. Socioeconomic status (SES) data and geographic location were obtained for each individual. There was significant geographic variation of ancestry across the island. In particular, African ancestry demonstrated a decreasing East to West gradient that was partially explained by historical factors linked to the colonial sugar plantation system. SES also demonstrated a parallel decreasing cline from East to West. However, at a local level, SES and African ancestry were negatively correlated. European ancestry was strongly negatively correlated with African ancestry and therefore showed patterns complementary to African ancestry. By contrast, Native American ancestry showed little variation across the island and across individuals and appears to have played little social role historically. The observed geographic distributions of SES and genetic variation relate to historical social events and mating patterns, and have substantial implications for the design of studies in the recently admixed Puerto Rican population. More generally, our results demonstrate the importance of incorporating social and geographic data with genetics when studying contemporary admixed populations.

Association of trypanolytic ApoL1 variants with kidney disease in African Americans

African Americans have higher rates of kidney disease than European Americans. Here, we show that, in African Americans, focal segmental glomerulosclerosis (FSGS) and hypertension-attributed end-stage kidney disease (H-ESKD) are associated with two independent sequence variants in the APOL1 gene on chromosome 22 {FSGS odds ratio = 10.5 [95% confidence interval (CI) 6.0 to 18.4]; H-ESKD odds ratio = 7.3 (95% CI 5.6 to 9.5)}. The two APOL1 variants are common in African chromosomes but absent from European chromosomes, and both reside within haplotypes that harbor signatures of positive selection. ApoL1 (apolipoprotein L-1) is a serum factor that lyses trypanosomes. In vitro assays revealed that only the kidney disease-associated ApoL1 variants lysed Trypanosoma brucei rhodesiense. We speculate that evolution of a critical survival factor in Africa may have contributed to the high rates of renal disease in African Americans.

Worldwide distribution of the MYH9 kidney disease susceptibility alleles and haplotypes: evidence of historical selection in Africa

MYH9 was recently identified as renal susceptibility gene (OR 3-8, p < 10(-8)) for major forms of kidney disease disproportionately affecting individuals of African descent. The risk haplotype (E-1) occurs at much higher frequencies in African Americans (> or = 60%) than in European Americans (< 4%), revealing a genetic basis for a major health disparity. The population distributions of MYH9 risk alleles and the E-1 risk haplotype and the demographic and selective forces acting on the MYH9 region are not well explored. We reconstructed MYH9 haplotypes from 4 tagging single nucleotide polymorphisms (SNPs) spanning introns 12-23 using available data from HapMap Phase II, and by genotyping 938 DNAs from the Human Genome Diversity Panel (HGDP). The E-1 risk haplotype followed a cline, being most frequent within sub-Saharan African populations (range 50-80%), less frequent in populations from the Middle East (9-27%) and Europe (0-9%), and rare or absent in Asia, the Americas, and Oceania. The fixation indexes (F(ST)) for pairwise comparisons between the risk haplotypes for continental populations were calculated for MYH9 haplotypes; F(ST) ranged from 0.27-0.40 for Africa compared to other continental populations, possibly due to selection. Uniquely in Africa, the Yoruba population showed high frequency extended haplotype length around the core risk allele (C) compared to the alternative allele (T) at the same locus (rs4821481, iHs = 2.67), as well as high population differentiation (F(ST(CEU vs. YRI)) = 0.51) in HapMap Phase II data, also observable only in the Yoruba population from HGDP (F(ST) = 0.49), pointing to an instance of recent selection in the genomic region. The population-specific divergence in MYH9 risk allele frequencies among the world's populations may prove important in risk assessment and public health policies to mitigate the burden of kidney disease in vulnerable populations.

Genome-wide scans for footprints of natural selection

Detecting recent selected ‘genomic footprints’ applies directly to the discovery of disease genes and in the imputation of the formative events that molded modern population genetic structure. The imprints of historic selection/adaptation episodes left in human and animal genomes allow one to interpret modern and ancestral gene origins and modifications. Current approaches to reveal selected regions applied in genome-wide selection scans (GWSSs) fall into eight principal categories: (I) phylogenetic footprinting, (II) detecting increased rates of functional mutations, (III) evaluating divergence versus polymorphism, (IV) detecting extended segments of linkage disequilibrium, (V) evaluating local reduction in genetic variation, (VI) detecting changes in the shape of the frequency distribution (spectrum) of genetic variation, (VII) assessing differentiating between populations (F_ST), and (VIII) detecting excess or decrease in admixture contribution from one population. Here, we review and compare these approaches using available human genome-wide datasets to provide independent verification (or not) of regions found by different methods and using different populations. The lessons learned from GWSSs will be applied to identify genome signatures of historic selective pressures on genes and gene regions in other species with emerging genome sequences. This would offer considerable potential for genome annotation in functional, developmental and evolutionary contexts.

Extended IL10 haplotypes and their association with HIV progression to AIDS

Interleukin-10 (IL-10) is a pleiotropic cytokine with both immunosuppressive and immunostimulatory functions. Its roles in infections and autoimmunity may have resulted in selective pressures on polymorphisms within the gene, leading to genomic coexistence of several semi-conserved haplotypes involved with diverse pathogen interactions during genomic evolution. Previous studies focused either exclusively on promoter haplotypes or on individual SNPs. We genotyped 21 single nucleotide polymorphisms in the human IL10 gene and examined this variation compared to other mammalian species sequences. Haplotype heterogeneity in human populations is centered around 'classic' 'proximal' promoter polymorphisms: -592, -819 and -1082. High-producing GCC haplotypes are by far the most numerous and diverse group, the intermediate IL-10 producing ACC-inclusive haplotypes seem to be related most closely to the ancestral haplotype, and the ATA-inclusive haplotypes cluster a separate branch with strong bootstrap support. We looked at associations of corresponding haplotypes with HIV progression. A haplotype trend regression confirmed that individuals carrying the low-producing ATA-inclusive haplotypes in European Americans progress to AIDS faster, and most likely explain the role of IL10. Our findings are consistent with the hypothesis that existing polymorphisms in this gene may reflect a balance of historic adaptive responses to autoimmune, infectious and other disease agents.

Genome and gene alterations by insertions and deletions in the evolution of human and chimpanzee chromosome 22

BACKGROUND

Understanding structure and function of human genome requires knowledge of genomes of our closest living relatives, the primates. Nucleotide insertions and deletions (indels) play a significant role in differentiation that underlies phenotypic differences between humans and chimpanzees. In this study, we evaluated distribution, evolutionary history, and function of indels found by comparing syntenic regions of the human and chimpanzee genomes.

RESULTS

Specifically, we identified 6,279 indels of 10 bp or greater in a ~33 Mb alignment between human and chimpanzee chromosome 22. After the exclusion of those in repetitive DNA, 1,429 or 23% of indels still remained. This group was characterized according to the local or genome-wide repetitive nature, size, location relative to genes, and other genomic features. We defined three major classes of these indels, using local structure analysis: (i) those indels found uniquely without additional copies of indel sequence in the surrounding (10 Kb) region, (ii) those with at least one exact copy found nearby, and (iii) those with similar but not identical copies found locally. Among these classes, we encountered a high number of exactly repeated indel sequences, most likely due to recent duplications. Many of these indels (683 of 1,429) were in proximity of known human genes. Coding sequences and splice sites contained significantly fewer of these indels than expected from random expectations, suggesting that selection is a factor in limiting their persistence. A subset of indels from coding regions was experimentally validated and their impacts were predicted based on direct sequencing in several human populations as well as chimpanzees, bonobos, gorillas, and two subspecies of orangutans.

CONCLUSION

Our analysis demonstrates that while indels are distributed essentially randomly in intergenic and intronic genomic regions, they are significantly under-represented in coding sequences. There are substantial differences in representation of indel classes among genomic elements, most likely caused by differences in their evolutionary histories. Using local sequence context, we predicted origins and phylogenetic relationships of gene-impacting indels in primate species. These results suggest that genome plasticity is a major force behind speciation events separating the great ape lineages.

Evaluation of IL10, IL19 and IL20 gene polymorphisms and chronic hepatitis B infection outcome

Hepatitis B virus (HBV) infection remains a serious global health problem despite the availability of a highly effective vaccine. Approximately 5% of HBV-infected adults develop chronic hepatitis B, which may result in liver cirrhosis or hepatocellular carcinoma. Variants of interleukin-10 (IL10) have been previously associated with chronic hepatitis B infection and progression to hepatocellular carcinoma. Single nucleotide polymorphisms (SNP; n = 42) from the IL10, IL19 and IL20 gene regions were examined for an association with HBV infection outcome, either chronic or recovered, in a nested case-control study of African Americans and European Americans. Among African Americans, three nominally statistically significant SNP associations in IL10, two in IL20, and one haplotype association were observed with different HBV infection outcomes (P = 0.005-0.04). A SNP (rs1518108) in IL20 deviated significantly from Hardy-Weinberg equilibrium in African Americans, with a large excess of heterozygotes in chronic HBV-infected cases (P = 0.0006), which suggests a strong genetic effect. Among European Americans, a nominally statistically significant SNP association in IL20 and an IL20 haplotype were associated with HBV recovery (P = 0.01-0.04). These results suggest that IL10 and IL20 gene variants influence HBV infection outcome and encourage the pursuit of further studies of these cytokines in HBV pathogenesis.

Identifying selected regions from heterozygosity and divergence using a light-coverage genomic dataset from two human populations

When a selective sweep occurs in the chromosomal region around a target gene in two populations that have recently separated, it produces three dramatic genomic consequences: 1) decreased multi-locus heterozygosity in the region; 2) elevated or diminished genetic divergence (F_ST) of multiple polymorphic variants adjacent to the selected locus between the divergent populations, due to the alternative fixation of alleles; and 3) a consequent regional increase in the variance of F_ST (S²F_ST) for the same clustered variants, due to the increased alternative fixation of alleles in the loci surrounding the selection target. In the first part of our study, to search for potential targets of directional selection, we developed and validated a resampling-based computational approach; we then scanned an array of 31 different-sized moving windows of SNP variants (5–65 SNPs) across the human genome in a set of European and African American population samples with 183,997 SNP loci after correcting for the recombination rate variation. The analysis revealed 180 regions of recent selection with very strong evidence in either population or both. In the second part of our study, we compared the newly discovered putative regions to those sites previously postulated in the literature, using methods based on inspecting patterns of linkage disequilibrium, population divergence and other methodologies. The newly found regions were cross-validated with those found in nine other studies that have searched for selection signals. Our study was replicated especially well in those regions confirmed by three or more studies. These validated regions were independently verified, using a combination of different methods and different databases in other studies, and should include fewer false positives. The main strength of our analysis method compared to others is that it does not require dense genotyping and therefore can be used with data from population-based genome SNP scans from smaller studies of humans or other species.

Genetics of cattails in radioactively contaminated areas around Chornobyl

Research on populations from radioactively contaminated areas around Chornobyl has produced ambiguous results for the presence of radiation effects. More studies are needed to provide information on whether radiation exposure at Chornobyl significantly affected genetic diversity in natural populations of various taxa. Eleven and nine variable microsatellite loci were used to test for differences in genetic diversity between reference and Chornobyl populations of two cattail species (Typha angustifolia and Typha latifolia, respectively) from Ukraine. Our purpose was to determine whether radiation had a significant impact on genetic diversities of the Chornobyl Typha populations, or if their genetic composition might be better explained by species demography and/or changes in population dynamics, mainly in sexual and asexual reproduction. Populations closest to the reactor had increased genetic diversities and high number of genets, which likely were due to factors other than radiation including increased gene flow among Chornobyl populations, enhanced sexual reproduction within populations, and/or origin of the genets from seed bank. Both Typha species also demonstrated small but significant effects associated with latitude, geographical regions, and watersheds. Typha's demography in Ukraine possibly varies with these three factors, and the small difference between Chornobyl and reference populations of T. latifolia detected after partitioning the total genetic variance between them is probably due primarily to these factors. However, the positive correlations of several genetic characteristics with radionuclide concentrations suggest that radiation may have also affected genetics of Chornobyl Typha populations but much less than was expected considering massive contamination of the Chornobyl area.

Independent occurrences of multiple repeats in the control region of mitochondrial DNA of white-tailed deer

Deer in the genera Mazama and Odocoileus generally have two copies of a 75-base-pair (bp) repeat in the left domain of the control region of the mitochondrial DNA (mtDNA). Phylogenetic analyses further suggest an ancient origin for the duplication supporting a previously stated contention that this event occurred before the separation of Mazama and Odocoileus. However, white-tailed deer (Odocoileus virginianus) had three or four copies of a 75-bp repeat in the control region of their mtDNA in 7.8% of the individuals analyzed, and all of these animals were from the coastal plain of the southeastern United States. When copy 3 is present, it is very similar in sequence to copy 2, but variation suggests that copy 3 probably evolved multiple times from copy 2. The pattern of phylogenetic clustering of the haplotypes from across the coastal plain also suggests that phenotypes with three or four copies of the repeat have originated multiple times. The 44 observed haplotypes showed strong spatial subdivision across the area with subpopulations frequently showing complete shifts in haplotype frequencies from others taken from nearby areas. Many of the subpopulations right along the coast or on adjacent barrier islands have a limited number of haplotypes as would occur in populations undergoing drift because of small numbers of breeding females and limited female dispersal.

Single nucleotide polymorphisms and haplotypes in the IL10 region associated with HCV clearance

Hepatitis C virus (HCV) is an infectious blood-borne pathogen that usually persists as a chronic infection. However, approximately 15% of the time, patients can clear the virus, indicating that host differences could be critical in determining the course of HCV infection. The inflammatory response is crucial to resolving or failing to resolve an acute HCV infection. Some previous reports have implicated interleukin 10 (IL10) polymorphisms with successful anti-HCV therapy and natural viral clearance. We tested 54 single nucleotide polymorphisms (SNPs) in the IL10 region (+/-300 kb and 24 within the IL10 gene itself), which contains 13 genes including the IL10 immunomodulatory paralogs IL19, IL20, and IL24, for association with HCV clearance vs persistence. SNPs from two haplotype block regions, one at IL10 and the other from IL19/IL20, were associated with HCV clearance in African Americans (91 clearance cases and 183 chronically infected matched controls; P=0.05-0.002) while with expectation-maximization algorithm-reconstructed haplotypes, these associations remained (P=0.05-0.002). However, no significant associations were detected in European Americans (108 clearance and 245 chronic). Our results indicate that variants of the immunomodulatory IL10 and IL19/IL20 genes may be involved in natural clearance of HCV in the African-American population.

Evaluating association and transmission of eight inflammatory genes with Viliuisk encephalomyelitis susceptibility

Since the discovery of Viliuisk encephalomyelitis (VE) in 1887, scientists have tried to understand the natural history and aetiology of this endemic neurological disorder among the native Sakha population of Central Siberia. Familial aggregation and segregation analysis suggested a genetic influence on VE incidence. However, recent studies have implicated an unknown virus, possibly from the alpha herpesvirus family, as a possible cause for this disease. As VE is a neurological disease characterized by the inflammatory reactions systematically observed in the spinocerebellar fluid and in the brain tissue of deceased patients, we examined 17 single nucleotide polymorphisms (SNPs) across seven inflammation-related candidate gene regions, including chemokine receptors type 2 and 5 (CCR2/CCR5), interferon-gamma (IFN-gamma), interleukin-4 (IL-4), IL-6, IL-10, stromal cell-derived factor (SDF) and chemokine regulated upon activation, normal T-cell expressed and presumably secreted (RANTES). Our main objective was to analyse the degree of genetic association between VE and candidate genes that have been previously implicated in other inflammatory diseases. Samples were collected from 83 affected families comprising 88 verified VE cases, 156 family members, and an additional 69 unrelated, unaffected inhabitants of the same geographical area. This collection included substantially all of the cases that are currently on the VE Registry. The experimental design included both case-control and transmission/disequilibrium test (TDT)-based familial association analyses. None of 17 SNPs analysed was significantly associated with VE occurrence. Exclusion of these eight genes based on the lack of association has important implications for identifying the disease agent, as well as prescribing therapy and understanding Viliuisk encephalomyelitis.

High levels of fluctuating asymmetry in populations of Apodemus flavicollis from the most contaminated areas in Chornobyl

Random deviations from the perfect symmetry of normally bilaterally symmetrical characters for an individual with a given genotype occur during individual development due to the influence of multiple environmental factors. Fluctuating asymmetry (FA) is often used as a measure of developmental instability, and can be estimated as the variance of the distribution of differences between the left and right sides. We addressed the question of whether levels of FA were elevated in radioactively contaminated populations living around Chornobyl compared to those in reference populations of the yellow-necked mouse (Apodemus flavicollis). In addition, we studied amounts of directional asymmetry (DA) when one side is larger than the other on average. There was a significant difference among populations, including reference populations, in the amount of both FA and DA. A higher level of FA was documented for the contaminated populations in close proximity to the failed Chornobyl reactor for both the asymmetry of size and shape. The FAs of size and shape were highest in populations from the most contaminated locations in the Chornobyl exclusion zone. Although the directional asymmetry of shape was also highest in the contaminated populations, it was not significantly different from those in most of the reference populations. Populations from less contaminated areas inside the Chornobyl exclusion zone did not express FA values different from those of the reference populations outside the affected area. FA of skulls of A. flavicollis may indicate the degree to which the level of radioactive contamination affects the development of animals at Chornobyl. However, the mechanisms leading to these effects are not clear and probably vary from population to population. There were significant correlations between the overall right to left differences for the Procrustes aligned shape configurations, centroid sizes, and intramuscular (137)Cs. Detectable effects of radiation on developmental stability probably start to occur between 0.132 and 0.297 microGy/h.

Frequency distributions of 137Cs in fish and mammal populations

We collected fish and mammals in several radioactively contaminated locations in the Chornobyl Exclusion Zone and analyzed them for 137Cs content. Frequency distributions were built for populations of channel catfish, yellow-necked mice and bank voles. We combined our data with similar data from several other studies to demonstrate the relationship between the standard deviations and means of 137Cs of fish and mammal populations. The frequency distributions of 137Cs in populations of fish and mammals are not normal, as indicated by the strong relationship between standard deviation and mean. Distributions for mammals are more skewed than those for fish. Fish and mammals probably use their environments in fundamentally different ways. The highest concentrations and thus greatest risks are therefore confined to relatively few individuals in each population.