Spatially dense 3D facial heritability and modules of co-heritability in a father-offspring design

Introduction: The human face is a complex trait displaying a strong genetic component as illustrated by various studies on facial heritability. Most of these start from sparse descriptions of facial shape using a limited set of landmarks. Subsequently, facial features are preselected as univariate measurements or principal components and the heritability is estimated for each of these features separately. However, none of these studies investigated multivariate facial features, nor the co-heritability between different facial features. Here we report a spatially dense multivariate analysis of facial heritability and co-heritability starting from data from fathers and their children available within ALSPAC. Additionally, we provide an elaborate overview of related craniofacial heritability studies. Methods: In total, 3D facial images of 762 father-offspring pairs were retained after quality control. An anthropometric mask was applied to these images to establish spatially dense quasi-landmark configurations. Partial least squares regression was performed and the (co-)heritability for all quasi-landmarks (∼7160) was computed as twice the regression coefficient. Subsequently, these were used as input to a hierarchical facial segmentation, resulting in the definition of facial modules that are internally integrated through the biological mechanisms of inheritance. Finally, multivariate heritability estimates were obtained for each of the resulting modules. Results: Nearly all modular estimates reached statistical significance under 1,000,000 permutations and after multiple testing correction (p ≤ 1.3889 × 10-3), displaying low to high heritability scores. Particular facial areas showing the greatest heritability were similar for both sons and daughters. However, higher estimates were obtained in the former. These areas included the global face, upper facial part (encompassing the nasion, zygomas and forehead) and nose, with values reaching 82% in boys and 72% in girls. The lower parts of the face only showed low to moderate levels of heritability. Conclusion: In this work, we refrain from reducing facial variation to a series of individual measurements and analyze the heritability and co-heritability from spatially dense landmark configurations at multiple levels of organization. Finally, a multivariate estimation of heritability for global-to-local facial segments is reported. Knowledge of the genetic determination of facial shape is useful in the identification of genetic variants that underlie normal-range facial variation

Spatially Dense 3D Facial Heritability and Modules of Co-heritability in a Father-Offspring Design

Introduction: The human face is a complex trait displaying a strong genetic component as illustrated by various studies on facial heritability. Most of these start from sparse descriptions of facial shape using a limited set of landmarks. Subsequently, facial features are preselected as univariate measurements or principal components and the heritability is estimated for each of these features separately. However, none of these studies investigated multivariate facial features, nor the co-heritability between different facial features. Here we report a spatially dense multivariate analysis of facial heritability and co-heritability starting from data from fathers and their children available within ALSPAC. Additionally, we provide an elaborate overview of related craniofacial heritability studies.Methods: In total, 3D facial images of 762 father-offspring pairs were retained after quality control. An anthropometric mask was applied to these images to establish spatially dense quasi-landmark configurations. Partial least squares regression was performed and the (co-)heritability for all quasi-landmarks (∼7160) was computed as twice the regression coefficient. Subsequently, these were used as input to a hierarchical facial segmentation, resulting in the definition of facial modules that are internally integrated through the biological mechanisms of inheritance. Finally, multivariate heritability estimates were obtained for each of the resulting modules.Results: Nearly all modular estimates reached statistical significance under 1,000,000 permutations and after multiple testing correction (p ≤ 1.3889 × 10-3), displaying low to high heritability scores. Particular facial areas showing the greatest heritability were similar for both sons and daughters. However, higher estimates were obtained in the former. These areas included the global face, upper facial part (encompassing the nasion, zygomas and forehead) and nose, with values reaching 82% in boys and 72% in girls. The lower parts of the face only showed low to moderate levels of heritability.Conclusion: In this work, we refrain from reducing facial variation to a series of individual measurements and analyze the heritability and co-heritability from spatially dense landmark configurations at multiple levels of organization. Finally, a multivariate estimation of heritability for global-to-local facial segments is reported. Knowledge of the genetic determination of facial shape is useful in the identification of genetic variants that underlie normal-range facial variation

3D facial phenotyping by biometric sibling matching used in contemporary genomic methodologies

The analysis of contemporary genomic data typically operates on one-dimensional phenotypic measurements (e.g. standing height). Here we report on a data-driven, family-informed strategy to facial phenotyping that searches for biologically relevant traits and reduces multivariate 3D facial shape variability into amendable univariate measurements, while preserving its structurally complex nature. We performed a biometric identification of siblings in a sample of 424 children, defining 1,048 sib-shared facial traits. Subsequent quantification and analyses in an independent European cohort (n = 8,246) demonstrated significant heritability for a subset of traits (0.17–0.53) and highlighted 218 genome-wide significant loci (38 also study-wide) associated with facial variation shared by siblings. These loci showed preferential enrichment for active chromatin marks in cranial neural crest cells and embryonic craniofacial tissues and several regions harbor putative craniofacial genes, thereby enhancing our knowledge on the genetic architecture of normal-range facial variation

Population genetics and population ecology in management of endangered species

Abstract Knowledge of the determinants of the viability of populations is essential in order to undertake effective conservation and management of endangered species. In this study, long-term demographic data was combined with genetic data to study the viability of an endangered orchid species, Epipactis atrorubens. The genetic analyses revealed low levels of genetic variation and the presence of population genetic differentiation independent of the spatial scale. Low levels of seed-mediated gene flow, possibly linked to low seedling recruitment, is the likely cause of the low levels of gene flow. Indications of slow post-glacial colonisation rates were found, which together with the low gene flow predict a limited capacity of the species to shift its range to more suitable habitats after environmental change. Low genetic variation as a proxy for low evolutionary potential also suggests that the species has limited capacity to adapt to new environmental conditions. Furthermore, poor seedling recruitment lowers population viability in small populations, as highlighted by the low population growth rates. In addition, we found a strong effect of stochasticity that limits the viability of populations. Both the genetic and demographic analyses indicated low viability of the studied species and that seedling recruitment could be the main determinant for the viability.Tiivistelmä Luonnonsuojelun perusta on populaatioiden elinkykyyn vaikuttavien tekijöiden tuntemus. Tässä väitöskirjatyössä tutkittiin uhanalaisen orkidean, tummaneidonvaipan (Epipactis atrorubens), elinkykyyn vaikuttavia tekijöitä yhdistämällä pitkäaikaisseurannoilla kerätyt demografiset aineistot geneettisin menetelmin kerättyihin aineistoihin. Lajin populaatioiden geneettisen muuntelun määrän havaittiin olevan pieni ja populaatioiden todettiin olevan geneettisesti erilaistuneita maantieteellisestä skaalasta riippumatta. Geneettisen erilaistumisen syy voi olla alhainen geenivirta, joka on seurausta vähäisestä siemendispersaalista ja huonosta taimettumisesta. Populaatioiden evolutiivista historiaa tutkittaessa havaittiin merkkejä hitaasta jääkauden jälkeisestä kolonisaatiosta, mikä yhdessä alhaisen geenivirran kanssa ennustaa, että lajilla on huono kyky siirtyä sille sopivammille alueille, jos ympäristö muuttuu. Huonoa evolutiivista potentiaalia kuvastava vähäinen geneettinen muuntelu ennustaa, että lajilla on huono kyky sopeutua uusiin ympäristöoloihin. Tämän lisäksi huono taimettuminen laskee elinkykyä etenkin pienissä populaatioissa, mikä näkyy muun muassa pienten populaatioiden matalina kasvukertoimina. Stokastinen vaihtelu vaikutti elinkykyä alentavasti, mikä pitäisikin huomioida nykyistä paremmin elinkykyanalyyseissä. Sekä geneettiset että demografiset analyysit osoittivat taimettumisen mahdollisesti olevan määräävä tekijä tummaneidonvaipan populaatioiden elinkyvylle

Low population viability in small endangered orchid populations:genetic variation, seedling recruitment and stochasticity

Abstract There are only few studies that use both demographic and genetic data to assess population viability of plant species. We combined genetic and demographic data from 11 endangered perennial orchid populations of varying size in order to reveal determinants of viability. Small populations had substantially lower viability compared to large populations. Seedling recruitment rates were remarkably lower in small populations; this was not due to pollination limitation or inbreeding depression because the fruit set and heterozygosity were not correlated with population size, suggesting that there may be differences in successful germination. Low recruitment resulted in significantly lower predicted population growth rates in small populations. The impact of stochasticity on viability varied among populations and stochastic simulations indicated that only one large population was viable, whereas all the other large populations were predicted to go extinct within decades. While there was a positive correlation between the deterministic population growth rate and allelic richness, we did not find any other correlations between genetic variation and fitness or population size. The study populations are likely remnant populations of a once large meta-population that decreased in size due to unfavourable environmental conditions. Management should focus on the maintenance of large population size, which is needed to avoid negative consequences of stochasticity and to enhance seedling recruitment rates

Genetic structure of a regionally endangered orchid, the dark red helleborine (Epipactis atrorubens) at the edge of its distribution

Abstract The genetic structure and diversity of species is determined by both current population dynamics and historical processes. Population genetic structure at the edge of the distribution is often expected to differ substantially from populations at the centre, as these edge populations are often small and fragmented. In addition, populations located in regions that have experienced repeated glaciations throughout the Pleistocene, may still carry imprints from the genetic consequences of frequent distribution shifts. Using chloroplast DNA sequences and nuclear microsatellite markers we studied the genetic structure of Epipactis atrorubens at the northern edge of its distribution. Contrary to populations in the centre of the distribution, populations at the northern range are regionally endangered as they are small and disjunct. Sequence data of 2 chloroplast loci and allelic data from 6 nuclear microsatellite markers were obtained from 297 samples from Finland, Estonia and Russia. We sought for genetic indicators of past population processes, such as post-glacial colonisation history of E. atrorubens. As expected, we observed low genetic variation, in terms of numbers of substitutions, haplotypes and alleles, and significant levels of differentiation, especially pronounced in the chloroplast DNA. These features suggest that the edge populations could be prone to extinction