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.

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.

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.

A course-based research experience: how benefits change with increased investment in instructional time.

There is widespread agreement that science, technology, engineering, and mathematics programs should provide undergraduates with research experience. Practical issues and limited resources, however, make this a challenge. We have developed a bioinformatics project that provides a course-based research experience for students at a diverse group of schools and offers the opportunity to tailor this experience to local curriculum and institution-specific student needs. We assessed both attitude and knowledge gains, looking for insights into how students respond given this wide range of curricular and institutional variables. While different approaches all appear to result in learning gains, we find that a significant investment of course time is required to enable students to show gains commensurate to a summer research experience. An alumni survey revealed that time spent on a research project is also a significant factor in the value former students assign to the experience one or more years later. We conclude: 1) implementation of a bioinformatics project within the biology curriculum provides a mechanism for successfully engaging large numbers of students in undergraduate research; 2) benefits to students are achievable at a wide variety of academic institutions; and 3) successful implementation of course-based research experiences requires significant investment of instructional time for students to gain full benefit.

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.

Reconstructing the population genetic history of the Caribbean.

The Caribbean basin is home to some of the most complex interactions in recent history among previously diverged human populations. Here, we investigate the population genetic history of this region by characterizing patterns of genome-wide variation among 330 individuals from three of the Greater Antilles (Cuba, Puerto Rico, Hispaniola), two mainland (Honduras, Colombia), and three Native South American (Yukpa, Bari, and Warao) populations. We combine these data with a unique database of genomic variation in over 3,000 individuals from diverse European, African, and Native American populations. We use local ancestry inference and tract length distributions to test different demographic scenarios for the pre- and post-colonial history of the region. We develop a novel ancestry-specific PCA (ASPCA) method to reconstruct the sub-continental origin of Native American, European, and African haplotypes from admixed genomes. We find that the most likely source of the indigenous ancestry in Caribbean islanders is a Native South American component shared among inland Amazonian tribes, Central America, and the Yucatan peninsula, suggesting extensive gene flow across the Caribbean in pre-Columbian times. We find evidence of two pulses of African migration. The first pulse--which today is reflected by shorter, older ancestry tracts--consists of a genetic component more similar to coastal West African regions involved in early stages of the trans-Atlantic slave trade. The second pulse--reflected by longer, younger tracts--is more similar to present-day West-Central African populations, supporting historical records of later transatlantic deportation. Surprisingly, we also identify a Latino-specific European component that has significantly diverged from its parental Iberian source populations, presumably as a result of small European founder population size. We demonstrate that the ancestral components in admixed genomes can be traced back to distinct sub-continental source populations with far greater resolution than previously thought, even when limited pre-Columbian Caribbean haplotypes have survived.

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 48% in MXL, 25% in CLM, and 13% in PUR. Native American ancestry in PUR is most closely related to populations surrounding the Orinoco River basin, confirming the Southern American 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 16 thousand years ago (kya), supports that the MXL Ancestors split 12.2kya, with a subsequent split of the ancestors to CLM and PUR 11.7kya. The model also features effective populations of 62,000 in Mexico, 8,700 in Colombia, and 1,900 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.

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.

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.

Estudios sobre ADN mitocondrial sugieren un linaje predominante en la cordillera Oriental de Colombia y un vínculo suramericano para los arcaicos de Puerto Rico / Mitochondrial DNA studies suggest a predominant lineage in the Eastern Cordillera of Colombia and South American link for the archaics of Puerto Rico

Objetivo: Este trabajo integra información de la secuencia del ADNmt del norte de Suramérica con Puerto Rico, con el fin de comprender el poblamiento del Caribe, especialmente de los taínos. De paso, arroja información sobre hechos demográficos en la Colombia precolombina. Metodología: Se obtuvieron 59 muestras de Colombia y Venezuela, las cuales fueron analizadas junto a otras dos pertenecientes a los indios warao y disponibles en el Genbank. Se alinearon secuencias HVR-I y II (Hypervariable Region) y se compararon con el rCRS. El 93,4% de las muestras resultaron ser de origen amerindio. Resultados: Un venezolano exhibió mutaciones relacionadas con el linaje antiguo C-II de Puerto Rico, el cual se estima que arribó a Puerto Rico en la era prearahuaca. Mediante secuenciación completa del ADNmt se demostró que esta muestra, VE6, pertenece al cladoamericano nativo C1b. Dos personas de Colombia y Venezuela presentaban la transición 16129 que define el linaje A-VIII de Puerto Rico. Dicha transición dentro del haplogrupo A también se haencontrado en los ciboneyes de Cuba y en otras tribus americanas. La deleción de un par de bases –498d– define el linaje B-I de Colombia (Bogotá y Villa de Leyva, Boyacá), un polimorfismo encontrado en los departamentos correspondientes a la cordillera Oriental y que se extiende al Valle del Cauca y a Panamá. Conclusión: Este linaje experimentó una expansión demográfica en la cordillera Oriental que lo llevó a expandirse geográficamente hasta Panamá. Sería recomendable ampliar el muestreo de la costa norte de Colombia y Venezuela, para encontrar más conexiones precolombinas con Puerto Rico. Además, sería conveniente verificar la distribución geográfica de 498d con un muestreo más numeroso y que cubra una zona más amplia de Colombia.

Objetive: This work integrates sequencing information of mtDNA from Northern South America with Puerto Rico, to reach an understanding of the peopling of the Caribbean, especially the Tainos. At the same time, it sheds light on demographic events in the Pre-Columbian Colombia. Methodology: Fifty nine samples from Colombia and Venezuela were obtained, and then analyzed along with two others from Warao Indians available in Genbank. HVR (Hypervariable Region) I andII sequences were aligned and compared to the rCRS. Fully 93,4% of the mtDNA samples were shown to be of Amerindian origin. Results: A Venezuelan exhibited ancient mutations related to lineage C-II of Puerto Rico, which has been estimated to have arrived to this island in pre-Arawak times. Through complete mtDNA sequencing, it was shown that this sample, VE6, belongs to the Native American C1b clade. Two individuals from Colombia and Venezuela showed the 16129 transition that defines lineage A-VIII of Puerto Rico. This transition has also been found in the Cuban Ciboneys and in various American tribes. A one base pair deletion –498d–defines lineage B-I from Colombia (Bogotá and “Villa de Leyva”, Boyacá), a polymorphism found in the departments belonging to the Eastern cordillera and extending to the Cauca Valley and Panamá. Conclusions: This lineage went through a demographic expansion in the Eastern Cordillera that may have triggered its geographic expansion to Panamá. It would be recommendable to expand the sampling of the Northern Coast of Colombia and Venezuela to find more pre-Columbian connections with Puerto Rico. Furthermore, it would be convenient to verify the geographic distribution of 498d with a bigger sample covering a wider region of Colombia.

Ancestry-environment interactions and asthma risk among Puerto Ricans.

BACKGROUND: Puerto Ricans, an admixed population of African, European, and Native American ancestries, have the highest asthma prevalence, morbidity, and mortality rates of any United States' population. Although socioeconomic status (SES) is negatively correlated with asthma incidence in most populations, no such relationship has been identified among Puerto Ricans. We hypothesized that, in this admixed population, the association between SES and asthma may interact with genetic ancestry. METHODS: We analyzed 135 Puerto Rican subjects with asthma and 156 control subjects recruited from six different recruitment centers in Puerto Rico. Individual ancestry for each subject was estimated using 44 ancestry informative markers. SES was assigned using the census tracts' median family income. Analyses of SES were based on the SES of the clinic site from which the subjects were recruited and on a subset of individuals on whom home address-based SES was available. RESULTS: In the two (independent) analyses, we found a significant interaction between SES, ancestry, and asthma disease status. At lower SES, European ancestry was associated with increased risk of asthma, whereas African ancestry was associated with decreased risk. The opposite was true for their higher SES counterparts. CONCLUSIONS: The observed interaction may help to explain the unique pattern of risk for asthma in Puerto Ricans and the lack of association with SES observed in previous studies when not accounting for varying proportions of ancestry.