Refining the evolutionary tree of the horse Y chromosome.

The Y chromosome carries information about the demography of paternal lineages, and thus, can prove invaluable for retracing both the evolutionary trajectory of wild animals and the breeding history of domesticates. In horses, the Y chromosome shows a limited, but highly informative, sequence diversity, supporting the increasing breeding influence of Oriental lineages during the last 1500 years. Here, we augment the primary horse Y-phylogeny, which is currently mainly based on modern horse breeds of economic interest, with haplotypes (HT) segregating in remote horse populations around the world. We analyze target enriched sequencing data of 5 Mb of the Y chromosome from 76 domestic males, together with 89 whole genome sequenced domestic males and five Przewalski's horses from previous studies. The resulting phylogeny comprises 153 HTs defined by 2966 variants and offers unprecedented resolution into the history of horse paternal lineages. It reveals the presence of a remarkable number of previously unknown haplogroups in Mongolian horses and insular populations. Phylogenetic placement of HTs retrieved from 163 archaeological specimens further indicates that most of the present-day Y-chromosomal variation evolved after the domestication process that started around 4200 years ago in the Western Eurasian steppes. Our comprehensive phylogeny significantly reduces ascertainment bias and constitutes a robust evolutionary framework for analyzing horse population dynamics and diversity.

First report of cucumber mosaic virus infecting Azorina vidalii in Azores Islands.

Virus diseases that occur in crops pose a major threat not only to global food security but also to wild plant communities growing in natural ecosystems (Jones, 2020, and references therein). In Azores (Portugal) little is known about viruses present on native flora and therefore, they have not been taken into consideration in conservation programs. Considering this, we selected Azorina vidalii (Campanulaceae), an endangered (IUCN) plant, endemic to Azores (Bilz, 2011), to survey for plant viruses. A. vidalii, the sole species of its genus, is often found in crevices with no soil accumulation on coastal cliffs, exposed to storms and sea spray, and is used as an ornamental. Leaves from 53 plants of A. vidalii from three populations from Terceira Island and three populations from Flores Island were randomly collected without obvious symptoms of virus infection, between the summer of 2021 and fall of 2022. RNA extraction was performed using the Plant/Fungi Total RNA Purification Kit (Norgen Biotek, Canada). RNA extracts from each population were pooled into six distinct composite samples (AvT1, AvT2, AvT3, AvF1, AvF4 and AvF5) and sent to Lexogen (Austria) for small RNA library preparation and High-Throughput Sequencing. Single-end RNA sequencing using Illumina NextSeq2000 system yielded between 10.1 M and 33.8 M raw reads. Adaptors and low-quality reads were removed with Trim Galore! and PRINSEQ. Trimmed reads were mapped to the genome phylogenetically nearest to A. vidalii available at the NCBI database (Adenophora triphylla). The resulting 2.5 M - 13.5 M unmapped reads were analysed with VirusDetect online version (database v248) (Zheng et al., 2017) for virus detection and identification. Sequences of cucumber mosaic virus (CMV) (contigs of up to 3045 nt for RNA1, 2917 nt for RNA2 and 2086 nt for RNA3) were identified in five (AvT1, AvT2, AvT3, AvF1 and AvF5) of the six composite samples and CMV satellite sequences (two contigs with 145 and 197 nt) were identified in only one (AvT1) of the composite samples. To confirm the presence of CMV, all samples were tested by two-step RT-PCR using primers targeting CMV-specific RdRp gene (513 bp) (Grieco et al., 2000), which retrieved 18 positive samples (34%). Nine samples were selected for Sanger sequencing (six out of 13 from Terceira and three out of five from Flores) based on digestion profile obtained with AluI and MboI. The resulting sequences (OQ176229-OQ176233, OQ732757-OQ732760) share an identity of 97.2-100% and BLASTn showed them to have 98.3-99.6% identity to CMV strain TN (AB176848). A Neighbour-Joining tree (Supplementary material) inferred in MEGA11 (Tamura et al., 2021) with 237 additional CMV-RdRp sequences, showed that A. vidalii CMV-derived isolates clustered together with reference strains of subgroup II, as those used by Roossinck (2002) for phylogenetic analysis of the 2a ORF. Besides CMV, tomato spotted wilt virus and polerovirus-associated RNAs sequences were found in one of the A. vidalii populations, but with lower coverage, and need to be further investigated. To the best of our knowledge, this is the first report of CMV infecting A. vidalli. CMV, genus Cucumovirus, is an agriculturally important virus and one of the most successful viruses known, infecting over 1,200 species of plants (Palukaitis & García-Arenal, 2003). In addition to A. vidalii being a CMV reservoir, which may have implications on adjacent crop fields, further research is needed to investigate the impact of CMV on A. vidalli fitness.

Developing a 670k genotyping array to tag ~2M SNPs across 24 horse breeds.

BACKGROUND: To date, genome-scale analyses in the domestic horse have been limited by suboptimal single nucleotide polymorphism (SNP) density and uneven genomic coverage of the current SNP genotyping arrays. The recent availability of whole genome sequences has created the opportunity to develop a next generation, high-density equine SNP array. RESULTS: Using whole genome sequence from 153 individuals representing 24 distinct breeds collated by the equine genomics community, we cataloged over 23 million de novo discovered genetic variants. Leveraging genotype data from individuals with both whole genome sequence, and genotypes from lower-density, legacy SNP arrays, a subset of ~5 million high-quality, high-density array candidate SNPs were selected based on breed representation and uniform spacing across the genome. Considering probe design recommendations from a commercial vendor (Affymetrix, now Thermo Fisher Scientific) a set of ~2 million SNPs were selected for a next-generation high-density SNP chip (MNEc2M). Genotype data were generated using the MNEc2M array from a cohort of 332 horses from 20 breeds and a lower-density array, consisting of ~670 thousand SNPs (MNEc670k), was designed for genotype imputation. CONCLUSIONS: Here, we document the steps taken to design both the MNEc2M and MNEc670k arrays, report genomic and technical properties of these genotyping platforms, and demonstrate the imputation capabilities of these tools for the domestic horse.

Analysis of copy number variants by three detection algorithms and their association with body size in horses.

BACKGROUND: Copy number variants (CNVs) have been shown to play an important role in genetic diversity of mammals and in the development of many complex phenotypic traits. The aim of this study was to perform a standard comparative evaluation of CNVs in horses using three different CNV detection programs and to identify genomic regions associated with body size in horses. RESULTS: Analysis was performed using the Illumina Equine SNP50 genotyping beadchip for 854 horses. CNVs were detected by three different algorithms, CNVPartition, PennCNV and QuantiSNP. Comparative analysis revealed 50 CNVs that affected 153 different genes mainly involved in sensory perception, signal transduction and cellular components. Genome-wide association analysis for body size showed highly significant deleted regions on ECA1, ECA8 and ECA9. Homologous regions to the detected CNVs on ECA1 and ECA9 have also been shown to be correlated with human height. CONCLUSIONS: Comparative analysis of CNV detection algorithms was useful to increase the specificity of CNV detection but had certain limitations dependent on the detection tool. GWAS revealed genome-wide associated CNVs for body size in horses.

Genetic diversity in the modern horse illustrated from genome-wide SNP data.

Horses were domesticated from the Eurasian steppes 5,000-6,000 years ago. Since then, the use of horses for transportation, warfare, and agriculture, as well as selection for desired traits and fitness, has resulted in diverse populations distributed across the world, many of which have become or are in the process of becoming formally organized into closed, breeding populations (breeds). This report describes the use of a genome-wide set of autosomal SNPs and 814 horses from 36 breeds to provide the first detailed description of equine breed diversity. F(ST) calculations, parsimony, and distance analysis demonstrated relationships among the breeds that largely reflect geographic origins and known breed histories. Low levels of population divergence were observed between breeds that are relatively early on in the process of breed development, and between those with high levels of within-breed diversity, whether due to large population size, ongoing outcrossing, or large within-breed phenotypic diversity. Populations with low within-breed diversity included those which have experienced population bottlenecks, have been under intense selective pressure, or are closed populations with long breed histories. These results provide new insights into the relationships among and the diversity within breeds of horses. In addition these results will facilitate future genome-wide association studies and investigations into genomic targets of selection.

Genome-wide analysis reveals selection for important traits in domestic horse breeds.

Intense selective pressures applied over short evolutionary time have resulted in homogeneity within, but substantial variation among, horse breeds. Utilizing this population structure, 744 individuals from 33 breeds, and a 54,000 SNP genotyping array, breed-specific targets of selection were identified using an F(ST)-based statistic calculated in 500-kb windows across the genome. A 5.5-Mb region of ECA18, in which the myostatin (MSTN) gene was centered, contained the highest signature of selection in both the Paint and Quarter Horse. Gene sequencing and histological analysis of gluteal muscle biopsies showed a promoter variant and intronic SNP of MSTN were each significantly associated with higher Type 2B and lower Type 1 muscle fiber proportions in the Quarter Horse, demonstrating a functional consequence of selection at this locus. Signatures of selection on ECA23 in all gaited breeds in the sample led to the identification of a shared, 186-kb haplotype including two doublesex related mab transcription factor genes (DMRT2 and 3). The recent identification of a DMRT3 mutation within this haplotype, which appears necessary for the ability to perform alternative gaits, provides further evidence for selection at this locus. Finally, putative loci for the determination of size were identified in the draft breeds and the Miniature horse on ECA11, as well as when signatures of selection surrounding candidate genes at other loci were examined. This work provides further evidence of the importance of MSTN in racing breeds, provides strong evidence for selection upon gait and size, and illustrates the potential for population-based techniques to find genomic regions driving important phenotypes in the modern horse.