Regions of Homozygosity in the Porcine Genome: Consequence of Demography and the Recombination Landscape

Inbreeding has long been recognized as a primary cause of fitness reduction in both wild and domesticated populations. Consanguineous matings cause inheritance of haplotypes that are identical by descent (IBD) and result in homozygous stretches along the genome of the offspring. Size and position of regions of homozygosity (ROHs) are expected to correlate with genomic features such as GC content and recombination rate, but also direction of selection. Thus, ROHs should be non-randomly distributed across the genome. Therefore, demographic history may not fully predict the effects of inbreeding. The porcine genome has a relatively heterogeneous distribution of recombination rate, making Sus scrofa an excellent model to study the influence of both recombination landscape and demography on genomic variation. This study utilizes next-generation sequencing data for the analysis of genomic ROH patterns, using a comparative sliding window approach. We present an in-depth study of genomic variation based on three different parameters: nucleotide diversity outside ROHs, the number of ROHs in the genome, and the average ROH size. We identified an abundance of ROHs in all genomes of multiple pigs from commercial breeds and wild populations from Eurasia. Size and number of ROHs are in agreement with known demography of the populations, with population bottlenecks highly increasing ROH occurrence. Nucleotide diversity outside ROHs is high in populations derived from a large ancient population, regardless of current population size. In addition, we show an unequal genomic ROH distribution, with strong correlations of ROH size and abundance with recombination rate and GC content. Global gene content does not correlate with ROH frequency, but some ROH hotspots do contain positive selected genes in commercial lines and wild populations. This study highlights the importance of the influence of demography and recombination on homozygosity in the genome to understand the effects of inbreeding.     

Evidence of Inbreeding in Hodgkin Lymphoma

Genome-wide association studies (GWASs) have identified several, mainly co-dominantly acting, single-nucleotide polymorphisms (SNPs) associated with Hodgkin lymphoma (HL). We searched for recessively acting disease loci by performing an analysis of runs of homozygosity (ROH) based on windows of homozygous SNP-blocks and by calculating genomic inbreeding coefficients on a SNP-wise basis. We used data from a previous GWAS with 906 cases and 1217 controls from a population with a long history of no matings between relatives. Ten recurrent ROHs were identified among 25 055 ROHs across all individuals but their association with HL was not genome-wide significant. All recurrent ROHs showed significant evidence for natural selection. As a novel finding genomic inbreeding among cases was significantly higher than among controls (P = 2.11*10⁻¹⁴) even after correcting for covariates. Higher inbreeding among the cases was mainly based on a group of individuals with a higher average length of ROHs per person. This result suggests a correlation of higher levels of inbreeding with higher cancer incidence and might reflect the existence of recessive alleles causing HL. Genomic inbreeding may result in a higher expression of deleterious recessive genes within a population.     

Conservation genetics of the annual hemiparasitic plant <Emphasis Type="Italic">Melampyrum sylvaticum</Emphasis> (Orobanchaceae) in the UK and Scandinavia

Melampyrum sylvaticum is an endangered annual hemiparasitic plant that is found in only 19 small and isolated populations in the United Kingdom (UK). To evaluate the genetic consequences of this patchy distribution we compared levels of diversity, inbreeding and differentiation from ten populations from the UK with eight relatively large populations from Sweden and Norway where the species is more continuously distributed. We demonstrate that in both the UK and Scandinavia, the species is highly inbreeding (global F _IS = 0.899). Levels of population differentiation were high (F’_ST = 0.892) and significantly higher amongst UK populations (F’_ST = 0.949) than Scandinavian populations (F’_ST = 0.762; P < 0.01). The isolated populations in the UK have, on average, lower genetic diversity (allelic richness, proportion of loci that are polymorphic, gene diversity) than Scandinavian populations, and this diversity difference is associated with the smaller census size and population area of UK populations. From a conservation perspective, the naturally inbreeding nature of the species may buffer the species against immediate effects of inbreeding depression, but the markedly lower levels of genetic diversity in UK populations may represent a genetic constraint to evolutionary change. In addition, the high levels of population differentiation suggest that gene flow among populations will not be effective at replenishing lost variation. We thus recommend supporting in situ conservation management with ex situ populations and human-mediated seed dispersal among selected populations in the UK.     

Inbreeding depression in intra-provenance crosses driven by founder relatedness in white spruce

Genetic means for height growth differed between intra- and inter-provenance crosses, which we hypothesized was due partly to unidentified relatedness among intra-provenance base parents resulting in mild inbreeding and inbreeding depression among their offspring. A dense array of 5,844 single nucleotide polymorphisms was used to directly construct a genomic relationship matrix (G) that had four elements ranging from 0.17 to 0.24, between five intra-provenance base parents. Adjusting the numerator relationship matrix for this relatedness among base parents produced inbreeding coefficients of F _i ≈ 0.1 in their offspring, which displayed depressed height growth. Accounting for inbreeding level as a covariate in a mixed model decreased grossly overestimated (up to 2 ×) dominance variance in models without the covariate adjustment. Height growth decreased 39 cm (～ 6 %) for every 0.1 increase in F _i .     

The pattern of runs of homozygosity and genomic inbreeding in world-wide sheep populations

Genome-wide pattern of runs of homozygosity (ROH) across ovine genome can provide a useful resource for studying diversity and demography history in sheep. We analyzed 50 k SNPs chip data of 2536 animals to identify pattern, distribution and level of ROHs in 68 global sheep populations. A total of 60,301 ROHs were detected in all breeds. The majority of the detected ROHs were <16 Mb and the average total number of ROHs per individual was 23.8 ± 13.8. The ROHs greater than 1 Mb covered on average 8.2% of the sheep autosomes, 1% of which was related to the ROHs with 1–4 Mb of length. The mean sum of ROH length in two-thirds of the populations was less than 250 Mb ranging from 21.7 to near 570 Mb. The level of genomic inbreeding was relatively low. The average of the inbreeding coefficients based on ROH (F_ROH) was 0.09 ± 0.05. It was rising in a stepwise manner with distance from Southwest Asia and maximum values were detected in North European breeds. A total of 465 ROH hotspots were detected in 25 different autosomes which partially surrounding 257 Refseq genes across the genome. Most of the detected genes were related to growth, body weight, meat production and quality, wool production and pigmentation. In conclusion, our analysis showed that the sheep genome, compared with other livestock species such as cattle and pig, displays low levels of homozygosity and appropriate genetic diversity for selection response and genetic merit gain.     

Surviving in isolation: genetic variation,bottlenecks and reproductive strategies in the Canarian endemic <Emphasis Type="Italic">Limonium macrophyllum</Emphasis> (Plumbaginaceae)

Oceanic archipelagos are typically rich in endemic taxa, because they offer ideal conditions for diversification and speciation in isolation. One of the most remarkable evolutionary radiations on the Canary Islands comprises the 16 species included in Limonium subsection Nobiles, all of which are subject to diverse threats, and legally protected. Since many of them are single-island endemics limited to one or a few populations, there exists a risk that a loss of genetic variation might limit their long-term survival. In this study, we used eight newly developed microsatellite markers to characterize the levels of genetic variation and inbreeding in L. macrophyllum, a species endemic to the North-east of Tenerife that belongs to Limonium subsection Nobiles. We detected generally low levels of genetic variation over all populations (H _T = 0.363), and substantial differentiation among populations (F _ST = 0.188; R _ST = 0.186) coupled with a negligible degree of inbreeding (F?=?0.042). Obligate outcrossing may have maintained L. macrophyllum relatively unaffected by inbreeding despite the species’ limited dispersal ability and the genetic bottlenecks likely caused by a prolonged history of grazing. Although several factors still constitute a risk for the conservation of L. macrophyllum, the lack of inbreeding and the recent positive demographic trends observed in the populations of this species are factors that favour its future persistence.     

利用高密度SNP标记分析中国荷斯坦牛基因组近交

畜禽育种中传统上利用系谱信息评估群体近交程度?近年来随着高通量单核苷酸多态(single nucleotide polymorphism, SNP)检测成本降低,使利用基因组信息分析真实的基因组近交程度成为可能?本研究利用牛54 K SNP 芯片数据统计了北京地区2107头荷斯坦牛基因组上的长纯合片段(runs of homozygosity, ROH)的频率和分布,计算了2种基因组近交系数,即染色体上ROH的长度占基因组总长度的比例(Froh)及个体所有标记基因型中纯合子所占比例,即基因组纯合度(Fhom),进而分析了两种基因组近交系数之间的相关性以及基因组近交与系谱近交系数之间的相关性?结果表明,共检测到44 676个ROH片段,其长度主要分布在1~10 Mb之间?不同长度的ROH散布于个体基因组内,短ROH较长ROH更为常见?ROH在染色体上并非均匀分布,ROH频率最高的区域为10号染色体中部?两种基因组近交系数之间相关性很高(91%以上),但基因组近交与系谱近交之间的相关性较低(低于50%)?系谱完整性是影响基因组近交与系谱近交结果一致的重要因素,基因组近交系数能够反映个体真实的近交,本研究为评估群体近交水平提供了有力工具?     

Assessing genetic structure in a rare clonal eucalypt as a basis for augmentation and introduction translocations

Genetics can provide information that is critical for planning translocations for conservation, such as levels of diversity and divergence of target and source populations. For clonal plants, assessing population characteristics (size, diversity, mortality, gene flow) that influence conservation values also requires identification of different genetic individuals. We used 12 microsatellite markers to guide germplasm source recommendations for augmentation and introduction translocations to conserve the critically endangered Eucalyptus cuprea that occurs in fragmented populations in the semi-arid shrublands of Western Australia. Ramet clumps with identical multilocus genotypes were identified in all populations but clonal richness (R = 0 ? 0.86) and heterogeneity (D = 0 ? 0.98) varied among populations. Genetic diversity was low to moderate in all populations (mean H _o = 0.61, mean A = 3.78) and did not differ significantly between localities. There was evidence of inbreeding in some populations but outcrossing (t _m = 0.495) in the small number of families available for study (N = 4) and genotypic diversity of the larger extant populations suggest the generation of novel genotypes is a component of the reproductive strategy. Most diversity was within populations and differentiation among populations was moderate (F _ST = 0.100) suggesting mixing of source population for translocation is unlikely to lead to outbreeding depression. Principal Co-ordinate and Bayesian analyses indicated the Northern population is distinct from Central/Southern populations. We recommend use of mixed germplasm to conserve the moderate diversity characterising larger remnant populations and to enable the production of recombinants through sexual reproduction. But given seed availability and the distinction of the Northern population, an initial precautionary approach to a translocation proposed for south of the geographical range may be to source germplasm from the Central/Southern locality.     

Influence of ethnic traditional cultures on genetic diversity of rice landraces under on-farm conservation in southwest China

Background

Crop genetic resources are important components of biodiversity. However, with the large-scale promotion of mono-cropping, genetic diversity has largely been lost. Ex-situ conservation approaches were widely used to protect traditional crop varieties worldwide. However, this method fails to maintain the dynamic evolutionary processes of crop genetic resources in their original habitats, leading to genetic diversity reduction and even loss of the capacity of resistance to new diseases and pests. Therefore, on-farm conservation has been considered a crucial complement to ex-situ conservation. This study aimed at clarifying the genetic diversity differences between ex-situ conservation and on-farm conservation and to exploring the influence of traditional cultures on genetic diversity of rice landraces under on-farm conservation.

Methods

The conservation status of rice landrace varieties, including Indica and Japonica, non-glutinous rice (Oryza sativa) and glutinous rice (Oryza sativa var. glutinosa Matsum), was obtained through ethno-biology investigation method in 12 villages of ethnic groups from Guizhou, Yunnan and Guangxi provinces of China. The genetic diversity between 24 pairs of the same rice landraces from different times were compared using simple sequence repeat (SSR) molecular markers technology. The landrace paris studied were collected in 1980 and maintained ex-situ, while 2014 samples were collected on-farm in southwest of China.

Results

The results showed that many varieties of rice landraces have been preserved on-farm by local farmers for hundreds or thousands of years. The number of alleles (Na), effective number of alleles (Ne), Nei genetic diversity index (He) and Shannon information index (I) of rice landraces were significantly higher by 12.3–30.4 % under on-farm conservation than under ex-situ conservation. Compared with the ex-situ conservation approach, rice landraces under on-farm conservation programs had more alleles and higher genetic diversity. In every site we investigated, ethnic traditional cultures play a positive influence on rice landrace variety diversity and genetic diversity.

Conclusion

Most China’s rice landraces were conserved in the ethnic areas of southwest China. On-farm conservation can effectively promote the allelic variation and increase the genetic diversity of rice landraces over the past 35 years. Moreover, ethnic traditional culture practices are a crucial foundation to increase genetic diversity of rice landraces and implement on-farm conservation.

     

A genetic assessment of seed production areas (SPAs) for restoration

The conservation and restoration of degraded landscapes continues to be a major activity with the demand for seed increasing to meet growing global targets. Seed collected from natural plant populations is often irregular and unpredictable, limiting the diversity of species that can be conserved through restoration programs and the area that can be restored. Seed production areas (SPAs; seed orchards, seed increase) are one option for improving seed supply and quality to meet the demand from conservation and restoration programs. We evaluated genetic diversity, inbreeding and mating system parameters in adults and seed from remnants and SPAs of Acacia montana and Dodonaea viscosa subsp. cuneata, two key restoration species in south-eastern Australia. Overall, we found no significant differences in genetic diversity between remnants and SPAs of both species or between adult and seed cohorts. In contrast, we found significant inbreeding in many remnants and their seed crops suggesting that some inbred seed were used to establish these SPAs. We also found significant inbreeding in some SPAs and their seed crops and that SPAs were often biased towards one to two source remnants. Additional germplasm is now required to broaden the representativeness of remnants in the SPAs and the genetic base of seed being produced for restoration programs.     

Diversity of fungi associated with roots of <Emphasis Type="Italic">Calanthe</Emphasis> orchid species in Korea

While symbiotic fungi play a key role in the growth of endangered Calanthe orchid species, the relationship between fungal diversity and Calanthe species remains unclear. Here, we surveyed root associated fungal diversity of six Calanthe orchid species by sequencing the internal transcribed spacer (ITS) region using 454 pyrosequencing. Our results revealed that Paraboeremia and Coprinopsis are dominant fungal genera among Calanthe species. In terms of overall relative abundance, Paraboeremia was the most common fungal genus associated with Calanthe roots, followed by Coprinopsis. Overall fungal diversity showed a significant degree of variation depending on both location and Calanthe species. In terms of number of different fungal genera detected within Calanthe species, C. discolor had the most diverse fungal community, with 10 fungal genera detected. This study will contribute toward a better understanding of those fungi that are required for successful cultivation and conservation of Korean Calanthe species.     

Limited diversity associated with duplicated class II MHC-<Emphasis Type="Italic">DRB</Emphasis> genes in the red squirrel population in the United Kingdom compared with continental Europe

The red squirrel (Sciurus vulgaris) population in the United Kingdom has declined over the last century and is now on the UK endangered species list. This is the result of competition from the eastern grey squirrel (S. carolinensis) which was introduced in the 19th century. However, recent evidence suggests that the rate of population decline is enhanced by squirrelpox disease, caused by a viral infection carried asymptomatically by grey squirrels but to which red squirrels are highly susceptible. Population genetic diversity provides some resilience to rapidly evolving or exotic pathogens. There is currently no data on genetic diversity of extant UK squirrel populations with respect to genes involved in disease resistance. Diversity is highest at loci involved in the immune response including genes clustered within the major histocompatibility complex (MHC). Using the class II DRB locus as a marker for diversity across the MHC region we genotyped 110 red squirrels from locations in the UK and continental Europe. Twenty-four Scvu-DRB alleles at two functional loci; Scvu-DRB1 and Scvu-DRB2, were identified. High levels of diversity were identified at both loci in the continental populations. In contrast, no diversity was observed at the Scvu-DRB2 locus in the mainland UK population while a high level of homozygosity was observed at the Scvu-DRB1 locus. The red squirrel population in the UK appears to lack the extensive MHC diversity associated with continental populations, a feature which may have contributed to their rapid decline.     

Evaluation of artificial selection in Standard Poodles using whole-genome sequencing

Identifying regions of artificial selection within dog breeds may provide insights into genetic variation that underlies breed-specific traits or diseases—particularly if these traits or disease predispositions are fixed within a breed. In this study, we searched for runs of homozygosity (ROH) and calculated the d _i statistic (which is based upon F _ST) to identify regions of artificial selection in Standard Poodles using high-coverage, whole-genome sequencing data of 15 Standard Poodles and 49 dogs across seven other breeds. We identified consensus ROH regions ≥1 Mb in length and common to at least ten Standard Poodles covering 0.6 % of the genome, and d _i regions that most distinguish Standard Poodles from other breeds covering 3.7 % of the genome. Within these regions, we identified enriched gene pathways related to olfaction, digestion, and taste, as well as pathways related to adrenal hormone biosynthesis, T cell function, and protein ubiquitination that could contribute to the pathogenesis of some Poodle-prevalent autoimmune diseases. We also validated variants related to hair coat and skull morphology that have previously been identified as being under selective pressure in Poodles, and flagged additional polymorphisms in genes such as ITGA2B, CBX4, and TNXB that may represent strong candidates for other common Poodle disorders.     

Genetic structure of island and mainland populations of a Neotropical bumble bee species

As a consequence of founder effects, small population size and demographic constraints, island populations are often characterized by low genetic diversity and high inbreeding. The effects of inbreeding are more pronounced in haplo-diploid insects like bees than in similar diploid species, because their method of sex determination requires heterozygosity at a sex locus. Inbreeding leads to homozygosity at the sex locus and the production of non-viable diploid males. This means that island populations of bees are particularly prone to extinction. Here we determine the levels of diversity and isolation between islands and mainland populations of the bumble bee Bombus morio in southeast Brazil. We analyzed 659 individuals from 24 populations, sequencing two mitochondrial genes (COI and Cytb) and genotyping all individuals at 14 microsatellite loci. Surprisingly, genetic diversity was high and genetic isolation was low in all populations except Teodoro Sampaio (mainland) and Ilha da Vitória (island). Genetic diversity is not significantly correlated with island area, but is lower in populations that are more distant from the mainland. Except perhaps for Ilha da Vitória, we suggest that the island populations are unlikely to go extinct due to genetic factors. Finally, based on its genetic distance from all other populations, we identify a putative new subspecies in the Teodoro Sampaio region.     

Genome-wide single nucleotide polymorphisms (SNPs) for a model invasive ascidian <Emphasis Type="Italic">Botryllus schlosseri</Emphasis>

Invasive species cause huge damages to ecology, environment and economy globally. The comprehensive understanding of invasion mechanisms, particularly genetic bases of micro-evolutionary processes responsible for invasion success, is essential for reducing potential damages caused by invasive species. The golden star tunicate, Botryllus schlosseri, has become a model species in invasion biology, mainly owing to its high invasiveness nature and small well-sequenced genome. However, the genome-wide genetic markers have not been well developed in this highly invasive species, thus limiting the comprehensive understanding of genetic mechanisms of invasion success. Using restriction site-associated DNA (RAD) tag sequencing, here we developed a high-quality resource of 14,119 out of 158,821 SNPs for B. schlosseri. These SNPs were relatively evenly distributed at each chromosome. SNP annotations showed that the majority of SNPs (63.20%) were located at intergenic regions, and 21.51% and 14.58% were located at introns and exons, respectively. In addition, the potential use of the developed SNPs for population genomics studies was primarily assessed, such as the estimate of observed heterozygosity (H _O ), expected heterozygosity (H _E ), nucleotide diversity (π), Wright’s inbreeding coefficient (F _IS ) and effective population size (Ne). Our developed SNP resource would provide future studies the genome-wide genetic markers for genetic and genomic investigations, such as genetic bases of micro-evolutionary processes responsible for invasion success.     

Genetic structure of <Emphasis Type="Italic">Pseudococcus microcirculus</Emphasis> (Hemiptera: Pseudococcidae) populations on epiphytic orchids in south Florida

In 2012, the orchid mealy bug Pseudococcus microcirculuswas first detected in situ in North America’s more diverse orchid region, the Big Cypress Basin (Collier Co FL). A follow-up survey showed that the mealy bug is more widespread and found on epiphytic orchids in two locations, in both the Fakahatchee Strand State Preserve (sites B and F) and the Florida Panther National Wildlife Refuge (sites M and C). There, we collected mealy bugs (n=54) from 35 orchid individuals and screened allelic variation at seven microsatellite loci. We estimated genetic diversity and differentiation among all sites and compared the variation among individuals collected on the same plant. Genetic differentiation between sites M and C (F _ST=0.03, P<0.01) and, M and B (F _ST=0.04, P<0.01) was detected. We also detected significantly lower mean pairwise relatedness among individuals from site B compared to all the other locations, and this population had the lowest inbreeding coefficient. Genetic diversity and mean pairwise relatedness were highly variable among plants with multiple individuals; however, plants from sites F and M tend to have collections of individuals with higher mean pairwise relatedness compared to sites B and C. Our results indicate that there is genetic diversity and differentiation among mealy bugs in these locations, and that collections of individuals on the same plant are genetically diverse. As such, the mealy bugs throughout these areas are likely to be genetically diverse and exist in multiple distinct populations.     

Persistence of the gypsophile <Emphasis Type="Italic">Lepidospartum burgessii</Emphasis> (Asteraceae) through clonal growth and limited gene flow

Lepidospartum burgessii is a rare gypsophilic shrub with limited distribution in New Mexico and Texas. Most of the known plants are restricted to two large populations, with a few small, isolated populations scattered in the surrounding area. The low recruitment observed in the two largest populations may be due to low seed set resulting from high inbreeding and/or self-incompatibility. We used eight microsatellite loci to quantify diversity, relatedness, inbreeding, population structure, and frequency of clonal reproduction. Seed collections were made to quantify seed set and germination rates. Overall, there was a moderate level of clonal diversity within patches of L. burgessii indicating asexual growth is important for population persistence. Inbreeding coefficients were high both between and within populations. Most patches showed a significant level of relatedness between individuals. At a fine scale, patches within populations were significantly different from each other, however when all patches were combined, the two populations of L. burgessii were not genetically distinct. Compared to a population of its common congener, Lepidospartum latisquamum, L. burgessii populations had similar measures of diversity, more clonal reproduction, and lower germination rates. High relatedness and inbreeding may explain the low seed set and recruitment in L. burgessii, however factors such as insect herbivory and precipitation changes may further depress recruitment.     

Ecomorph Evolution in <Emphasis Type="Italic">Myotis</Emphasis> (Vespertilionidae,Chiroptera)

The genus Myotis is unique among mammals in its high taxonomic diversity and global distribution. Their phylogenetic relationships reflect biogeographic affinities rather than phenotypes. Myotis diverged from other bats in the early Miocene, with a subsequent split between Old and New World lineages about 19 million years ago. Similar ecomorphs (‘Leuconoe’ [near-water hunters], ‘Myotis’ [gleaners], ‘Selysius’ [aerial hawkers]) emerged independently in different lineages of Myotis. We retrieved the probable ancestral ecomorph for each lineage. Phenetic diversity was estimated from the analysis of body and skull traits. It seems that evolution of Myotis fluctuated between ‘Leuconoe’, ‘Selysius’, and larger ‘Myotis’.     

Population viability analysis predicts decreasing genetic diversity in <Emphasis Type="Italic">ex situ</Emphasis> populations of the Itasenpara bitterling <Emphasis Type="Italic">Acheilognathus longipinnis</Emphasis> from the Kiso River,Japan

To conserve endangered species, the maintenance of ex situ captive populations with sustainable genetic diversity is often required, in combination with population viability analysis (PVA). Since 2010, the threatened Itasenpara bitterling Acheilognathus longipinnis lineages in the Kiso region, Japan, have been maintained in ex situ rearing facilities to allow for conservation efforts. In this study, we obtained microsatellite data from DNA extracted from these captive populations to elucidate their genetic diversity and effective population size. The populations of several initial generations indicated a deviation from Hardy–Weinberg equilibrium, probably due to the limited number of extracted founder individuals analyzed. The effective population size of the captive population tended to increase over the course of generations, although the degree of genetic diversity tended to decrease highlighting the concern for the progression of inbreeding. Our prediction based on the PVA suggests that the maintenance of the captive population under the current conditions could lead to extinction of the Itasenpara bitterling in 50 years. In contrast, simultaneously increasing the carrying capacity and individual exchange among populations appears to enhance the effective management of captive Itasenpara bitterling populations.