RESUMO

For decades, inbreeding in cattle has been evaluated using pedigree information. Nowadays, inbreeding coefficients can be obtained using genomic information such as runs of homozygosity (ROH). The aims of this study were to quantify ROH and heterozygosity-rich regions (HRR) in a subpopulation of Guzerá dual-purpose cattle, to examine ROH and HRR islands, and to compare inbreeding coefficients obtained by ROH with alternative genomic inbreeding coefficients. A subpopulation of 1733 Guzerá animals genotyped for 50k SNPs was used to obtain the ROH and HRR segments. Inbreeding coefficients by ROH (FROH ), by genomic relationship matrix based on VanRaden's method 1 using reference allele frequency in the population (FGRM ), by genomic relationship matrix based on VanRaden's method 1 using allele frequency fixed in 0.5 (FGRM_0.5 ), and by the proportion of homozygous loci (FHOM ) were calculated. A total of 15,660 ROH were identified, and the chromosome with the highest number of ROH was BTA6. A total of 4843 HRRs were identified, and the chromosome with the highest number of HRRs was BTA23. No ROH and HRR islands were identified according to established criteria, but the regions closest to the definition of an island were examined from 64 to 67 Mb of BTA6, from 36 to 37 Mb of BTA2 and from 0.50 to 1.25 Mb of BTA23. The genes identified in ROH islands have previously been associated with dairy and beef traits, while genes identified on HRR islands have previously been associated with reproductive traits and disease resistance. FROH was equal to 0.095 ± 0.084, and its Spearman correlation with FGRM was low (0.44) and moderate-high with FHOM (0.79) and with FGRM_0.5 (0.80). The inbreeding coefficients determined by ROH were higher than other cattle breeds' and higher than pedigree-based inbreeding in the Guzerá breed obtained in previous studies. It is recommended that future studies investigate the effects of inbreeding determined by ROH on the traits under selection in the subpopulation studied.

Assuntos

RESUMO

The analysis of the genomic landscape of inbreeding using runs of homozygosity (ROH) patterns is becoming an interesting tool to partially understand phenotypic differences among individuals. In this study, we analysed genome-wide ROH patterns in two groups of Florida goats. We first determined the inbreeding levels of each individual by calculating ROH-based inbreeding coefficients (FROH ). Then, the individuals were divided into two groups based on FROH : high inbreeding (HI, FROH >0.1) and low inbreeding (LI, FROH <0.03). Finally, we performed an extensive in-depth analysis of ROH distribution in each group separately. We found a higher abundance of short ROH in LI, whereas long ROH was more frequent in HI. Furthermore, ROH abundance was not evenly distributed among chromosomes within groups, with some chromosomes showing larger numbers of ROH, like CHI6, CHI7 and CHI27. A different landscape was observed in recent inbreeding (ROH >8 Mb), with significant increases in CHI6, CHI11 and CHI28. Determination of genomic regions with significantly increased ROH (ROH islands-ROHi) showed 13 ROHi related to whole inbreeding and five ROHi associated with recent inbreeding analysis. Within these genomic regions, 123 and 101 genes were identified in HI and LI, respectively, including 10 and seven candidate genes previously related to production, fertility and heat resistance in goats and livestock species.

Assuntos

RESUMO

One of the strategies to genetically improve reproductive traits, despite their low inheritability, has been the identification of candidate genes. Therefore, the objective of this study was to detect candidate genes associated with fecundity through the fixation index (FST) and runs of homozygosity (ROH) of selection signatures in Katahdin ewes. Productive and reproductive records from three years were used and the genotypes (OvineSNP50K) of 48 Katahdin ewes. Two groups of ewes were identified to carry out the genetic comparison: with high fecundity (1.3 ± 0.03) and with low fecundity (1.1 ± 0.06). This study shows for the first time evidence of the influence of the CNOT11, GLUD1, GRID1, MAPK8, and CCL28 genes in the fecundity of Katahdin ewes; in addition, new candidate genes were detected for fecundity that were not reported previously in ewes but that were detected for other species: ANK2 (sow), ARHGAP22 (cow and buffalo cow), GHITM (cow), HERC6 (cow), DPF2 (cow), and TRNAC-GCA (buffalo cow, bull). These new candidate genes in ewes seem to have a high expression in reproduction. Therefore, future studies are needed focused on describing the physiological basis of changes in the reproductive behavior influenced by these genes.

RESUMO

Admixture is known to greatly impact the genetic landscape of a population and, while genetic variation underlying human phenotypes has been shown to differ among populations, studies on admixed subjects are still scarce. Latin American populations are the result of complex demographic history, such as 2 or 3-way admixing events, bottlenecks and/or expansions, and adaptive events unique to the American continent. To explore the impact of these events on the genetic structure of Latino populations, we evaluated the following haplotype features: linkage disequilibrium, shared identity by descent segments, runs of homozygosity, and extended haplotype homozygosity (integrated haplotype score) in Latinos represented in the 1000 Genome Project along with array data from 171 Brazilians sampled in the South and Southeast regions of Brazil. We found that linkage disequilibrium decay relates to the amount of American and African ancestry. The extent of identity by descent sharing positively correlates with historical effective population sizes, which we found to be steady or growing, except for Puerto Ricans and Colombians. Long runs of homozygosity, a particular instance of autozygosity, was only enriched in Peruvians and Native Americans. We used simulations to account for random sampling and linkage disequilibrium to filter positive selection indexes and found 244 unique markers under selection, 26 of which are common to 2 or more populations. Some markers exhibiting positive selection signals had estimated time to the most recent common ancestor consistent with human adaptation to the American continent. In conclusion, Latino populations present highly divergent haplotype characteristics that impact genetic architecture and underlie complex phenotypes.

Assuntos

RESUMO

BACKGROUND: The detection of signatures of selection in genomic regions provides insights into the evolutionary process, enabling discoveries regarding complex phenotypic traits. In this research, we focused on identifying genomic regions affected by different selection pressures, mainly highlighting the recent positive selection, as well as understanding the candidate genes and functional pathways associated with the signatures of selection in the Mangalarga Marchador genome. Besides, we seek to direct the discussion about genes and traits of importance in this breed, especially traits related to the type and quality of gait, temperament, conformation, and locomotor system. RESULTS: Three different methods were used to search for signals of selection: Tajima's D (TD), the integrated haplotype score (iHS), and runs of homozygosity (ROH). The samples were composed of males (n = 62) and females (n = 130) that were initially chosen considering well-defined phenotypes for gait: picada (n = 86) and batida (n = 106). All horses were genotyped using a 670 k Axiom® Equine Genotyping Array​ (Axiom MNEC670). In total, 27, 104 (chosen), and 38 candidate genes were observed within the signatures of selection identified in TD, iHS, and ROH analyses, respectively. The genes are acting in essential biological processes. The enrichment analysis highlighted the following functions: anterior/posterior pattern for the set of genes (GLI3, HOXC9, HOXC6, HOXC5, HOXC4, HOXC13, HOXC11, and HOXC10); limb morphogenesis, skeletal system, proximal/distal pattern formation, JUN kinase activity (CCL19 and MAP3K6); and muscle stretch response (MAPK14). Other candidate genes were associated with energy metabolism, bronchodilator response, NADH regeneration, reproduction, keratinization, and the immunological system. CONCLUSIONS: Our findings revealed evidence of signatures of selection in the MM breed that encompass genes acting on athletic performance, limb development, and energy to muscle activity, with the particular involvement of the HOX family genes. The genome of MM is marked by recent positive selection. However, Tajima's D and iHS results point also to the presence of balancing selection in specific regions of the genome.

Assuntos

RESUMO

Blanco Orejinegro (BON) cattle have 500 years of adaptation to the Colombian tropic, but little is known about their genetic history. Our aim was to estimate levels of linkage disequilibrium (LD), effective population size (Ne), genomic inbreeding for runs of homozygosity (FROH ), genomic relation matrix (FGRM ), excess of homozygotes (FHOM ) and pedigree information (FPEDCOMP ) and to characterize the runs of homozygosity (ROH), searching for selection signatures. A total of 419 BON animals were genotyped, 70 with a 150K chip and 349 with a 50K chip. Next, an imputation to 50K was performed, and, after editing, databases of 40K were obtained. The PLINK v1.90 and R programs were used to estimate LD, ROH, FROH and FHOM . The SNeP v1.1 program was used to obtain Ne, and PreGSf90 was used to elaborate the scaled G matrix. The MTDFNRM program was used to estimate FPEDCOMP . The LD mean as r2 at 1 Mb was 0.21 (r2  > 0.30 at a distance of 96.72kb), and Ne was 123 ± 1. A total of 7,652 homozygous segments were obtained, with a mean of 18.35 ± 0.55 ROH/animal. Most of the genome was covered by long ROHs (ROH>8 Mb  = 4.86%), indicating significant recent inbreeding. The average inbreeding coefficient for FPEDCOM , FGRM , FHOM and FROH was 4.41%, 4.18%, 5.58% and 6.78%, respectively. The highest correlation was observed between FHOM and FROH (0.95). ROH hotspots/islands were defined using the extreme values of a box plot that was generated, and correspond to QTLs related to milk yield (55.11%), external appearance (13.47%), production (13.30%), reproduction (8.15%), health (5.24%) and meat carcass (4.74%).

Assuntos

RESUMO

With the availability of high-density SNP panels and the establishment of approaches for characterizing homozygosity and heterozygosity sites, it is possible to access fine-scale information regarding genomes, providing more than just comparisons of different inbreeding coefficients. This is the first study that seeks to access such information for the Mangalarga Marchador (MM) horse breed on a genomic scale. To this end, we aimed to assess inbreeding levels using different coefficients, as well as to characterize homozygous and heterozygous runs in the population. Using Axiom ® Equine Genotyping Array-670k SNP (Thermo Fisher), 192 horses were genotyped. Our results showed different estimates: inbreeding from genomic coefficients (FROH ) = 0.16; pedigree-based (FPED ) = 0.008; and a method based on excess homozygosity (FHOM ) = 0.010. The correlations between the inbreeding coefficients were low to moderate, and some comparisons showed negative correlations, being practically null. In total, 85,295 runs of homozygosity (ROH) and 10,016 runs of heterozygosity (ROHet) were characterized for the 31 horse autosomal chromosomes. The class with the highest percentage of ROH was 0-2 Mbps, with 92.78% of the observations. In the ROHet results, only the 0-2 class presented observations, with chromosome 11 highlighted in a region with high genetic variability. Three regions from the ROHet analyses showed genes with known functions: tripartite motif-containing 37 (TRIM37), protein phosphatase, Mg2+ /Mn2+ dependent 1E (PPM1E) and carbonic anhydrase 10 (CA10). Therefore, our findings suggest moderate inbreeding, possibly attributed to breed formation, annulling possible recent inbreeding. Furthermore, regions with high variability in the MM genome were identified (ROHet), associated with the recent selection and important events in the development and performance of MM horses over generations.

Assuntos

RESUMO

Characterization of autozygosity is relevant to monitor genetic diversity and manage inbreeding levels in breeding programs. Identification of autozygosity hotspots can unravel genomic regions targeted by selection for economically important traits and can help identify candidate genes for selection. In this study, we estimated the inbreeding levels of a Brazilian population of Murrah buffalo undergoing selection for milk production traits, particularly milk yield. We also studied the distribution of runs of homozygosity (ROH) islands and identified putative genes and quantitative trait loci (QTL) under selection. We genotyped 422 Murrah buffalo for 51,611 SNP; 350 of these had ROH longer than 10 Mb, indicating the occurrence of inbreeding in the last 5 generations. The mean length of the ROH per animal was 4.28 ± 1.85 Mb. Inbreeding coefficients were calculated from the genomic relationship matrix, the pedigree, and the ROH, with estimates varying between 0.242 and 0.035. Inbreeding estimates from the pedigree had a low correlation with the genomic estimates, and estimates from the genomic relationship matrix were much higher than those from the pedigree or the ROH. Signatures of selection were identified in 6 genomic regions, located on chromosomes 1, 2, 3, 5, 16, and 18, encompassing a total of 190 genes and 174 QTL. Many of the genes (e.g., APRT and ACSF3) and QTL identified are related to milk production traits, such as milk yield, milk fat yield and percentage, and milk protein yield and percentage. Other genes are associated with reproduction and immune response traits as well as morphological aspects of the buffalo species. Inbreeding levels in this population are still low but are increasing due to selection and should be managed to avoid future losses due to inbreeding depression. The proximity of genes linked to milk production traits with genes associated with reproduction and immune system traits suggests the need to include these latter genes in the breeding program to avoid negatively affecting them due to selection for production traits.

Assuntos

RESUMO

BACKGROUND: Currently, chromosomal microarrays (CMA) are recommended as first-tier test in the investigation of developmental disorders to examine copy number variations. The modern platforms also include probes for single nucleotide polymorphisms (SNPs) that detect homozygous regions in the genome, such as long contiguous stretches of homozygosity (LCSH) also named runs of homozygosity (ROH). LCHS are chromosomal segments resulting from complete or segmental chromosomal homozygosity, which may be indicative of uniparental disomy (UPD), consanguinity, as well as replicative DNA repair events, however also are common findings in normal populations. Knowing common LCSH of a population, which probably represent ancestral haplotypes of low-recombination regions in the genome, facilitates the interpretation of LCSH found in patients, allowing to prioritize those with possible clinical significance. However, population records of ancestral haplotype derived LCSH by SNP arrays are still scarce, particularly for countries such as Brazil where even for the clinic, microarrays that include SNPs are difficult to request due to their high cost. METHODS: In this study, we evaluate the frequencies and implications of LCSH detected by Affymetrix CytoScan® HD or 750 K platforms in 430 patients with neurodevelopmental disorders in southern Brazil. LCSH were analyzed in the context of pathogenic significance and also explored to identify ancestral haplotype derived LCSH. The criteria for considering a region as LCSH was homozygosis ≥3 Mbp on an autosome. RESULTS: In 95% of the patients, at least one LCSH was detected, a total of 1478 LCSH in 407 patients. In 2.6%, the findings were suggestive of UPD. For about 8.5% LCSH suggest offspring from first to fifth grade, more likely to have a clinical impact. Considering recurrent LCSH found at a frequency of 5% or more, we outline 11 regions as potentially representing ancestral haplotypes in our population. The region most involved with homozygosity was 16p11.2p11.1 (49%), followed by 1q21.2q21.3 (21%), 11p11.2p11.12 (19%), 3p21.31p21.2 (16%), 15q15 1q33p32.3 (12%), 2q11.1q12.1 (9%), 1p33p32.3 (6%), 20q11.21q11.23 (6%), 10q22.1q23.31 (5%), 6p22.2p22 (5%), and 7q11.22q11.23 (5%). CONCLUSIONS: In this work, we show the importance and usefulness of interpreting LCSH in the results of CMA wich incorporate SNPs.

Assuntos

RESUMO

Splicing-related gene mutations might affect the expression of a single gene or multiple genes and cause clinically heterogeneous diseases. With the advent of next-generation sequencing, several splicing gene mutations have been exposed, yet most major spliceosome genes have no reports of germline mutations and therefore, their effects are largely unknown. We describe the previously unreported concurrence of intellectual disability, short stature, poor speech, and minor craniofacial and hand anomalies in 2 female siblings with 3 homozygous missense variants in SNRPA (a component of the U1 small nuclear ribonucleoprotein complex) characterized by homozygosity mapping and whole exome sequencing. Combined, c.97A>G, c.98T>C, and c.100T>A, in exon 2 of SNRPA lead to p.Ile33Ala and p.Phe34Ile exchanges, which were predicted in silico to be deleterious. Although both patients exhibited some clinical features seen in other spliceosomal disorders, their complete clinical phenotype appears to be rather uncommon, a finding that may further support the notion that mutations in components of the major spliceosome do not strictly lead to the same syndromes/phenotypes.

Assuntos